[cig-commits] r16822 - in seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN: . DATA
yangl at geodynamics.org
yangl at geodynamics.org
Fri May 28 08:11:05 PDT 2010
Author: yangl
Date: 2010-05-28 08:11:05 -0700 (Fri, 28 May 2010)
New Revision: 16822
Added:
seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/NOISE_clean_create
seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/Z_submit
Modified:
seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/DATA/CMTSOLUTION
seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/DATA/Par_file
seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/DATA/STATIONS
seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/Makefile.in
seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/constants.h.in
seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/create_header_file.f90
seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/create_movie_shakemap_AVS_DX_GMT.f90
seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/flags.guess
seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/read_parameter_file.f90
seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/specfem3D.f90
Log:
noise cross-correlation capability
Modified: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/DATA/CMTSOLUTION
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/DATA/CMTSOLUTION 2010-05-28 09:53:29 UTC (rev 16821)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/DATA/CMTSOLUTION 2010-05-28 15:11:05 UTC (rev 16822)
@@ -1,13 +1,13 @@
-PDE 2001 9 9 23 59 17.78 34.0745 -118.3792 6.4 4.2 4.2 HOLLYWOOD POLET SEM
-event name: 9703873
+XXXX 2005 8 6 5 40 33.30 36.1488 -118.0652 3.7 3.9 3.9 14169456_SEMm12_Lin
+event name: 14169456
time shift: 0.0000
half duration: 0.0000
-latitude: 34.0745
-longitude: -118.3792
-depth: 5.4000
-Mrr: -0.002000e+23
-Mtt: -0.064000e+23
-Mpp: 0.066000e+23
-Mrt: -0.090000e+23
-Mrp: -0.002000e+23
-Mtp: 0.188000e+23
+latitude: 36.1488
+longitude: -118.0652
+depth: 3.6694
+Mrr: 0.0e+21
+Mtt: 0.0e+21
+Mpp: 0.0e+21
+Mrt: 0.0e+21
+Mrp: 0.0e+21
+Mtp: 0.0e+21
Modified: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/DATA/Par_file 2010-05-28 09:53:29 UTC (rev 16821)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/DATA/Par_file 2010-05-28 15:11:05 UTC (rev 16822)
@@ -1,13 +1,13 @@
# forward or adjoint simulation
SIMULATION_TYPE = 1 # 1 = forward, 2 = adjoint, 3 = both simultaneously
-SAVE_FORWARD = .false.
+SAVE_FORWARD = .true.
# coordinates of mesh block in latitude/longitude and depth in km
-LATITUDE_MIN = 33.8d0
-LATITUDE_MAX = 34.1d0
-LONGITUDE_MIN = -118.4d0
-LONGITUDE_MAX = -118.1d0
+LATITUDE_MIN = 32.2d0
+LATITUDE_MAX = 36.8d0
+LONGITUDE_MIN = -121.6d0
+LONGITUDE_MAX = -114.7d0
DEPTH_BLOCK_KM = 60.d0
UTM_PROJECTION_ZONE = 11
SUPPRESS_UTM_PROJECTION = .false.
@@ -15,54 +15,54 @@
# number of elements at the surface along edges of the mesh at the surface
# (must be 8 * multiple of NPROC below if mesh is not regular and contains mesh doublings)
# (must be multiple of NPROC below if mesh is regular)
-NEX_XI = 16
-NEX_ETA = 16
+NEX_XI = 336
+NEX_ETA = 288
# number of MPI processors along xi and eta (can be different)
-NPROC_XI = 2
-NPROC_ETA = 2
+NPROC_XI = 14
+NPROC_ETA = 12
# model (SoCal, Harvard_LA, Min_Chen_anisotropy)
-MODEL = Harvard_LA
+MODEL = SoCal
# parameters describing the model
-OCEANS = .true.
-TOPOGRAPHY = .true.
-ATTENUATION = .true.
+OCEANS = .false.
+TOPOGRAPHY = .false.
+ATTENUATION = .false.
USE_OLSEN_ATTENUATION = .false.
# absorbing boundary conditions for a regional simulation
ABSORBING_CONDITIONS = .true.
# record length in seconds
-RECORD_LENGTH_IN_SECONDS = 30.0
+RECORD_LENGTH_IN_SECONDS = 200.0
# save AVS or OpenDX movies
-MOVIE_SURFACE = .false.
+MOVIE_SURFACE = .true.
MOVIE_VOLUME = .false.
-NTSTEP_BETWEEN_FRAMES = 200
+NTSTEP_BETWEEN_FRAMES = 2000
CREATE_SHAKEMAP = .false.
SAVE_DISPLACEMENT = .false.
-USE_HIGHRES_FOR_MOVIES = .false.
+USE_HIGHRES_FOR_MOVIES = .true.
HDUR_MOVIE = 0.0
# save AVS or OpenDX mesh files to check the mesh
SAVE_MESH_FILES = .false.
# path to store the local database file on each node
-LOCAL_PATH = /scratch/komatits/DATABASES_MPI
+LOCAL_PATH = /scratch/yangl/DATABASES_MPI
# machine file for MPI
# this is not needed for new cluster and is in go_mesher/go_solver
# MACHINE_FILE = mymachines
# interval at which we output time step info and max of norm of displacement
-NTSTEP_BETWEEN_OUTPUT_INFO = 100
+NTSTEP_BETWEEN_OUTPUT_INFO = 2000
# interval in time steps for writing of seismograms
NTSTEP_BETWEEN_OUTPUT_SEISMOS = 200000
# print source time function
-PRINT_SOURCE_TIME_FUNCTION = .false.
+PRINT_SOURCE_TIME_FUNCTION = .true.
Modified: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/DATA/STATIONS
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/DATA/STATIONS 2010-05-28 09:53:29 UTC (rev 16821)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/DATA/STATIONS 2010-05-28 15:11:05 UTC (rev 16822)
@@ -1,1112 +1,13 @@
- BVA1 AZ 33.326500 -116.365500 0.00 0.00
- BVA2 AZ 33.326500 -116.365500 0.00 0.00
- BVDA1 AZ 33.326500 -116.365500 0.00 0.00
- BVDA2 AZ 33.326500 -116.365500 0.00 0.00
- BZN AZ 33.491500 -116.667000 0.00 0.00
- CPE AZ 32.888900 -117.105100 0.00 0.00
- CRY AZ 33.565400 -116.737300 0.00 0.00
- CSLB AZ 33.779300 -118.112600 0.00 0.00
- DHL AZ 33.388900 -115.787500 0.00 0.00
- FLV AZ 32.911400 -117.061900 0.00 0.00
- FLV2 AZ 32.911400 -117.061900 0.00 0.00
- FRD AZ 33.494700 -116.602200 0.00 0.00
- FSE AZ 33.005700 -116.592600 0.00 0.00
- GLA AZ 33.051200 -114.827000 0.00 0.00
- HWB AZ 33.026200 -116.959600 0.00 0.00
- KNW AZ 33.714100 -116.711900 0.00 0.00
- LVA2 AZ 33.351600 -116.561500 0.00 0.00
- MONP AZ 32.892700 -116.422500 0.00 0.00
- MTRP AZ 32.820300 -117.056700 0.00 0.00
- PFO AZ 33.611700 -116.459400 0.00 0.00
- RDM AZ 33.630000 -116.847800 0.00 0.00
- SCI2 AZ 32.915000 -118.487900 0.00 0.00
- SCV1 AZ 33.611800 -116.459700 0.00 0.00
- SCV2 AZ 33.611800 -116.459700 0.00 0.00
- SMER AZ 33.457700 -117.170800 0.00 0.00
- SMER2 AZ 33.457700 -117.170800 0.00 0.00
- SMTC AZ 32.944900 -115.799900 0.00 0.00
- SND AZ 33.551900 -116.612900 0.00 0.00
- SOL AZ 32.841000 -117.248000 0.00 0.00
- THSB AZ 32.878800 -117.226900 0.00 0.00
- TRO AZ 33.523400 -116.425700 0.00 0.00
- WMC AZ 33.573600 -116.674700 0.00 0.00
- YAQ AZ 33.166600 -116.353900 0.00 0.00
- ARC BK 40.877700 -124.077400 0.00 0.00
- BBEB BK 37.821600 -122.329800 0.00 0.00
- BDM BK 37.953970 -121.865540 0.00 0.00
- BKS BK 37.876200 -122.235600 0.00 0.00
- BRIB BK 37.918900 -122.151800 0.00 0.00
- BRK BK 37.873500 -122.261000 0.00 0.00
- CMB BK 38.034550 -120.386510 0.00 0.00
- CMSB BK 37.871900 -122.251700 0.00 0.00
- CRQB BK 38.055800 -122.224900 0.00 0.00
- CVS BK 38.345260 -122.458400 0.00 0.00
- ELFS BK 40.618300 -120.727900 0.00 0.00
- FARB BK 37.697820 -123.001100 0.00 0.00
- GASB BK 39.654710 -122.715950 0.00 0.00
- HAST BK 36.388700 -121.551400 0.00 0.00
- HATC BK 40.817300 -121.470500 0.00 0.00
- HELL BK 36.680100 -119.022800 0.00 0.00
- HERB BK 38.012400 -122.262200 0.00 0.00
- HOPS BK 38.993490 -123.072340 0.00 0.00
- HUMO BK 42.607100 -122.956690 0.00 0.00
- JCC BK 40.817450 -124.029550 0.00 0.00
- JRSC BK 37.403730 -122.238680 0.00 0.00
- KCC BK 37.323630 -119.318700 0.00 0.00
- MCCM BK 38.144780 -122.880180 0.00 0.00
- MHC BK 37.341600 -121.642600 0.00 0.00
- MHDL BK 37.842300 -122.494300 0.00 0.00
- MIN BK 40.346000 -121.606600 0.00 0.00
- MNRC BK 38.878740 -122.442770 0.00 0.00
- MOBB BK 36.690600 -122.166000 0.00 0.00
- MOD BK 41.902460 -120.302950 0.00 0.00
- OHLN BK 38.006300 -122.273000 0.00 0.00
- ORV BK 39.554510 -121.500360 0.00 0.00
- OXMT BK 37.499400 -122.424300 0.00 0.00
- PACP BK 37.008000 -121.287000 0.00 0.00
- PKD BK 35.945170 -120.541600 0.00 0.00
- PKD1 BK 35.889400 -120.426100 0.00 0.00
- PKD2 BK 35.888200 -120.430500 0.00 0.00
- POTR BK 38.202630 -121.935350 0.00 0.00
- RAMR BK 35.636000 -120.869800 0.00 0.00
- RFSB BK 37.916100 -122.336100 0.00 0.00
- RSRB BK 37.935700 -122.446500 0.00 0.00
- SAO BK 36.764000 -121.447200 0.00 0.00
- SBRN BK 37.685600 -122.411300 0.00 0.00
- SCCB BK 37.287400 -121.864200 0.00 0.00
- SM2B BK 37.838700 -122.110200 0.00 0.00
- SMCB BK 37.838800 -122.111600 0.00 0.00
- SUTB BK 39.229100 -121.786100 0.00 0.00
- SVIN BK 38.033200 -122.526300 0.00 0.00
- VALB BK 38.121500 -122.275300 0.00 0.00
- W02B BK 37.791100 -122.386300 0.00 0.00
- WDC BK 40.579880 -122.541130 0.00 0.00
- WENL BK 37.622110 -121.756970 0.00 0.00
- YBH BK 41.732040 -122.710390 0.00 0.00
- YBIB BK 37.814200 -122.359200 0.00 0.00
- CCRB BP 35.957160 -120.551610 0.00 0.00
- EADB BP 35.895220 -120.422600 0.00 0.00
- FROB BP 35.910950 -120.486880 0.00 0.00
- GHIB BP 35.832250 -120.347280 0.00 0.00
- JCNB BP 35.938970 -120.431120 0.00 0.00
- JCSB BP 35.921200 -120.434000 0.00 0.00
- LCCB BP 35.980007 -120.514230 0.00 0.00
- MMNB BP 35.956500 -120.496000 0.00 0.00
- RMNB BP 36.000870 -120.477720 0.00 0.00
- SCYB BP 36.009420 -120.536610 0.00 0.00
- SMNB BP 35.972970 -120.579880 0.00 0.00
- VARB BP 35.926080 -120.447050 0.00 0.00
- VCAB BP 35.921620 -120.533920 0.00 0.00
- 14405 CE 33.787000 -118.356000 0.00 0.00
- 14901 CE 33.562000 -118.117000 0.00 0.00
- 23732 CE 34.190000 -117.330000 0.00 0.00
- 24400 CE 34.037000 -118.178000 0.00 0.00
- 24851 CE 34.070000 -118.346000 0.00 0.00
- 24853 CE 34.077000 -118.286000 0.00 0.00
- 25131 CE 34.151000 -119.198000 0.00 0.00
- 26903 CE 34.611000 -120.731000 0.00 0.00
- 57202 CE 37.638000 -121.497000 0.00 0.00
- G405 CE 33.787000 -118.356000 0.00 0.00
- J732 CE 34.190000 -117.330000 0.00 0.00
- K851 CE 34.070000 -118.346000 0.00 0.00
- K853 CE 34.077000 -118.286000 0.00 0.00
- ADO CI 34.550460 -117.433910 0.00 0.00
- AGA CI 33.638380 -116.401100 0.00 0.00
- AGO CI 34.146470 -118.766990 0.00 0.00
- ALP CI 34.687080 -118.299460 0.00 0.00
- ARG CI 36.083700 -117.480600 0.00 0.00
- ARV CI 35.126900 -118.830090 0.00 0.00
- BAK CI 35.344440 -119.104450 0.00 0.00
- BAR CI 32.680050 -116.672150 0.00 0.00
- BBA CI 34.195510 -118.353400 0.00 0.00
- BBB CI 33.352620 -115.733250 0.00 0.00
- BBR CI 34.262300 -116.920750 0.00 0.00
- BBS CI 33.921390 -116.980580 0.00 0.00
- BC3 CI 33.655150 -115.453660 0.00 0.00
- BCC CI 33.575080 -117.261190 0.00 0.00
- BEL CI 34.000600 -115.998200 0.00 0.00
- BFS CI 34.238830 -117.658530 0.00 0.00
- BKR CI 35.269300 -116.070300 0.00 0.00
- BLA CI 34.069480 -116.388960 0.00 0.00
- BLY CI 33.750340 -114.523750 0.00 0.00
- BOR CI 33.268200 -116.417160 0.00 0.00
- BR2 CI 32.978990 -115.549360 0.00 0.00
- BRE CI 33.807760 -117.981160 0.00 0.00
- BTC CI 33.012130 -115.219870 0.00 0.00
- BTP CI 34.682240 -118.573980 0.00 0.00
- BVH CI 34.076180 -118.395900 0.00 0.00
- CAB CI 34.155730 -118.640930 0.00 0.00
- CAC CI 34.136780 -118.121990 0.00 0.00
- CAP CI 33.388220 -117.194950 0.00 0.00
- CAR CI 35.308190 -119.845830 0.00 0.00
- CBC CI 34.140110 -118.127500 0.00 0.00
- CCC CI 35.524950 -117.364530 0.00 0.00
- CFS CI 34.105390 -117.282870 0.00 0.00
- CGO CI 36.550400 -117.802950 0.00 0.00
- CHF CI 34.333410 -118.025850 0.00 0.00
- CHN CI 33.998790 -117.680440 0.00 0.00
- CIA CI 33.401860 -118.413720 0.00 0.00
- CLC CI 35.815740 -117.597510 0.00 0.00
- CLT CI 34.092840 -117.316870 0.00 0.00
- COO CI 33.896040 -118.216390 0.00 0.00
- CPP CI 34.060200 -117.809000 0.00 0.00
- CRN CI 33.875740 -117.561060 0.00 0.00
- CRP CI 34.136240 -118.127050 0.00 0.00
- CTC CI 33.655160 -115.990060 0.00 0.00
- CTD CI 33.354860 -118.444360 0.00 0.00
- CWC CI 36.439880 -118.080160 0.00 0.00
- DAN CI 34.637450 -115.381150 0.00 0.00
- DEC CI 34.253530 -118.333830 0.00 0.00
- DEV CI 33.935970 -116.577940 0.00 0.00
- DGR CI 33.650010 -117.009470 0.00 0.00
- DJJ CI 34.106180 -118.455050 0.00 0.00
- DJJB CI 34.106180 -118.455050 0.00 0.00
- DLA CI 33.848220 -118.096240 0.00 0.00
- DNR CI 33.566670 -116.630560 0.00 0.00
- DPP CI 32.998620 -116.941530 0.00 0.00
- DRC CI 32.805400 -115.446540 0.00 0.00
- DRE CI 32.805320 -115.446790 0.00 0.00
- DSC CI 35.142550 -116.103950 0.00 0.00
- DVT CI 32.659150 -116.100610 0.00 0.00
- EDW CI 34.883030 -117.991060 0.00 0.00
- EDW2 CI 34.881100 -117.993880 0.00 0.00
- ELC CI 32.781400 -115.535600 0.00 0.00
- EML CI 32.890830 -116.845660 0.00 0.00
- ERR CI 33.116750 -115.822670 0.00 0.00
- EW2 CI 33.940590 -116.408390 0.00 0.00
- FIG CI 34.728320 -119.988030 0.00 0.00
- FLL CI 34.397270 -118.918070 0.00 0.00
- FLS CI 34.970360 -117.039090 0.00 0.00
- FMP CI 33.712640 -118.293810 0.00 0.00
- FON CI 34.099570 -117.438760 0.00 0.00
- FPC CI 35.082000 -117.582670 0.00 0.00
- FUL CI 33.871700 -117.922510 0.00 0.00
- FUR CI 36.467030 -116.863220 0.00 0.00
- GLA CI 33.051490 -114.827060 0.00 0.00
- GMR CI 34.784570 -115.659940 0.00 0.00
- GOR CI 33.153700 -117.229210 0.00 0.00
- GR2 CI 34.118300 -118.299400 0.00 0.00
- GRA CI 36.996080 -117.366210 0.00 0.00
- GRF CI 34.119200 -118.300400 0.00 0.00
- GRH CI 34.308030 -118.559540 0.00 0.00
- GSA CI 34.136770 -118.128300 0.00 0.00
- GSC CI 35.301770 -116.805740 0.00 0.00
- GVR CI 34.049700 -118.119970 0.00 0.00
- HEC CI 34.829400 -116.335000 0.00 0.00
- HLL CI 34.176430 -118.359670 0.00 0.00
- HLN CI 34.121280 -117.218970 0.00 0.00
- IRM CI 34.157380 -115.145130 0.00 0.00
- ISA CI 35.662780 -118.474030 0.00 0.00
- JCS CI 33.085900 -116.595900 0.00 0.00
- JRC CI 35.982300 -117.807600 0.00 0.00
- JRC2 CI 35.982490 -117.808850 0.00 0.00
- JVA CI 34.366220 -116.612660 0.00 0.00
- KIK CI 34.150370 -118.101560 0.00 0.00
- KML CI 33.597380 -117.879090 0.00 0.00
- LAF CI 33.868890 -118.331430 0.00 0.00
- LAX CI 33.943820 -118.413900 0.00 0.00
- LBW1 CI 33.798000 -118.088370 0.00 0.00
- LBW2 CI 33.798000 -118.088370 0.00 0.00
- LCG CI 34.000330 -118.377940 0.00 0.00
- LCP CI 34.735510 -120.279960 0.00 0.00
- LDF CI 35.130660 -115.184160 0.00 0.00
- LDR CI 34.990600 -118.341560 0.00 0.00
- LEV CI 34.614620 -118.291040 0.00 0.00
- LFP CI 34.305290 -118.488050 0.00 0.00
- LGB CI 33.975300 -118.149180 0.00 0.00
- LGU CI 34.108190 -119.065870 0.00 0.00
- LJR CI 34.807620 -118.867750 0.00 0.00
- LKL CI 34.615940 -117.824930 0.00 0.00
- LLS CI 33.684470 -117.943040 0.00 0.00
- LMR CI 34.934260 -117.696610 0.00 0.00
- LMR2 CI 34.934390 -117.696290 0.00 0.00
- LRL CI 35.479540 -117.682120 0.00 0.00
- LTP CI 33.881100 -118.175680 0.00 0.00
- LTR CI 34.521100 -117.990300 0.00 0.00
- LUG CI 34.365600 -117.366830 0.00 0.00
- MAG CI 35.361650 -118.922820 0.00 0.00
- MCT CI 34.226450 -116.040730 0.00 0.00
- MGE CI 33.818410 -116.368740 0.00 0.00
- MIK CI 34.137070 -118.126080 0.00 0.00
- MIS CI 33.740300 -118.335170 0.00 0.00
- MLAC CI 37.630190 -118.836050 0.00 0.00
- MLS CI 34.004600 -117.561620 0.00 0.00
- MOP CI 34.280850 -118.904900 0.00 0.00
- MPI CI 34.812580 -119.145300 0.00 0.00
- MPM CI 36.057990 -117.489010 0.00 0.00
- MPP CI 34.888480 -119.813620 0.00 0.00
- MSJ CI 33.808010 -116.967890 0.00 0.00
- MTL CI 34.268800 -118.237400 0.00 0.00
- MTP CI 35.484340 -115.553200 0.00 0.00
- MUR CI 33.599990 -117.195430 0.00 0.00
- MWC CI 34.223620 -118.058320 0.00 0.00
- NBS CI 34.780350 -116.557980 0.00 0.00
- NEE CI 34.824900 -114.599410 0.00 0.00
- NEE2 CI 34.767590 -114.618810 0.00 0.00
- NJQ CI 34.534120 -120.177370 0.00 0.00
- NOT CI 34.228690 -118.558290 0.00 0.00
- NSS CI 33.556100 -115.946500 0.00 0.00
- NSS2 CI 33.555530 -115.945860 0.00 0.00
- OGC CI 33.788160 -117.844000 0.00 0.00
- OKV CI 34.396970 -119.298740 0.00 0.00
- OLI CI 33.945390 -117.923720 0.00 0.00
- OLP CI 32.607830 -116.930360 0.00 0.00
- OSI CI 34.614500 -118.723500 0.00 0.00
- PAS CI 34.148440 -118.171170 0.00 0.00
- PASA CI 34.148440 -118.171130 0.00 0.00
- PASC CI 34.171410 -118.185230 0.00 0.00
- PDE CI 34.441990 -118.582150 0.00 0.00
- PDM CI 34.303360 -114.141520 0.00 0.00
- PDR CI 33.962730 -118.437020 0.00 0.00
- PDU CI 34.120700 -117.638080 0.00 0.00
- PER CI 33.861640 -117.205290 0.00 0.00
- PHL CI 35.407730 -120.545560 0.00 0.00
- PLC CI 33.824360 -116.511950 0.00 0.00
- PLM CI 33.353610 -116.862650 0.00 0.00
- PLS CI 33.795300 -117.609060 0.00 0.00
- PMD CI 33.647850 -116.377690 0.00 0.00
- PSD CI 33.823930 -116.550260 0.00 0.00
- PSR CI 34.091810 -117.807090 0.00 0.00
- QUG CI 34.395580 -118.498370 0.00 0.00
- RCT CI 36.305230 -119.243840 0.00 0.00
- RIN CI 34.281960 -118.479250 0.00 0.00
- RINB CI 34.281960 -118.479250 0.00 0.00
- RIO CI 34.104730 -117.979560 0.00 0.00
- RPV CI 33.743460 -118.404120 0.00 0.00
- RRS CI 33.882170 -117.366460 0.00 0.00
- RRX CI 34.875330 -116.996840 0.00 0.00
- RSB CI 33.973100 -117.327200 0.00 0.00
- RSS CI 33.973270 -117.326740 0.00 0.00
- RUS CI 34.050730 -118.080850 0.00 0.00
- RVR CI 33.993510 -117.375450 0.00 0.00
- RXH CI 33.183130 -115.622570 0.00 0.00
- SAL CI 33.280100 -115.985850 0.00 0.00
- SAN CI 33.704320 -117.885780 0.00 0.00
- SBB2 CI 34.688500 -117.824480 0.00 0.00
- SBC CI 34.440760 -119.714920 0.00 0.00
- SBI CI 33.480460 -119.029860 0.00 0.00
- SBP CI 34.232360 -117.234880 0.00 0.00
- SBPX CI 34.232400 -117.234840 0.00 0.00
- SCI CI 32.979900 -118.546970 0.00 0.00
- SCI2 CI 32.979900 -118.546970 0.00 0.00
- SCZ CI 33.995320 -119.634350 0.00 0.00
- SCZ2 CI 33.995430 -119.635100 0.00 0.00
- SDD CI 33.552590 -117.661710 0.00 0.00
- SDG CI 32.784000 -117.138050 0.00 0.00
- SDP CI 34.565470 -120.501370 0.00 0.00
- SDR CI 32.735610 -116.942410 0.00 0.00
- SES CI 34.436920 -119.137500 0.00 0.00
- SHO CI 35.899530 -116.275300 0.00 0.00
- SIO CI 34.293000 -119.164600 0.00 0.00
- SJU CI 33.487250 -117.681140 0.00 0.00
- SLA CI 35.890950 -117.283320 0.00 0.00
- SLR CI 33.833590 -116.797370 0.00 0.00
- SMB CI 34.902180 -120.447090 0.00 0.00
- SMF CI 34.021590 -118.446750 0.00 0.00
- SMM CI 35.314200 -119.995810 0.00 0.00
- SMS CI 34.014380 -118.456170 0.00 0.00
- SMV CI 34.270910 -118.744070 0.00 0.00
- SNCC CI 33.248000 -119.524000 0.00 0.00
- SOT CI 34.416000 -118.449000 0.00 0.00
- SPF CI 34.059330 -118.646140 0.00 0.00
- SPG CI 36.135500 -118.810990 0.00 0.00
- SPG2 CI 36.200570 -118.766240 0.00 0.00
- SRN CI 33.828540 -117.789380 0.00 0.00
- SSW CI 33.177470 -115.615640 0.00 0.00
- STC CI 34.303020 -119.186760 0.00 0.00
- STG CI 33.664000 -117.768560 0.00 0.00
- STS CI 33.790330 -118.198780 0.00 0.00
- SVD CI 34.106470 -117.098220 0.00 0.00
- SWS CI 32.940800 -115.795800 0.00 0.00
- SYP CI 34.527750 -119.978340 0.00 0.00
- TA2 CI 34.382030 -117.678220 0.00 0.00
- TCF CI 34.083770 -118.599000 0.00 0.00
- TEH CI 35.291300 -118.420790 0.00 0.00
- TFT CI 35.145920 -119.419460 0.00 0.00
- THX CI 33.634950 -116.164020 0.00 0.00
- TIN CI 37.054220 -118.230090 0.00 0.00
- TOR CI 33.575260 -116.225840 0.00 0.00
- TOV CI 34.156070 -118.820390 0.00 0.00
- TUQ CI 35.435840 -115.923890 0.00 0.00
- USB CI 34.413000 -119.842700 0.00 0.00
- USC CI 34.019190 -118.286310 0.00 0.00
- VCS CI 34.483640 -118.117830 0.00 0.00
- VES CI 35.840890 -119.084690 0.00 0.00
- VRD CI 34.214590 -118.279640 0.00 0.00
- VTV CI 34.560650 -117.329600 0.00 0.00
- WBS CI 35.536640 -118.140350 0.00 0.00
- WER CI 35.060530 -119.027110 0.00 0.00
- WES CI 32.759040 -115.731610 0.00 0.00
- WGR CI 34.510850 -119.274070 0.00 0.00
- WLT CI 34.009480 -117.950770 0.00 0.00
- WNS CI 34.115000 -118.380000 0.00 0.00
- WSS CI 34.171700 -118.649710 0.00 0.00
- WTT CI 33.948690 -118.255470 0.00 0.00
- FABS FA 34.066670 -118.571113 0.00 0.00
- TUC IU 32.309600 -110.784600 0.00 0.00
- BMN LB 40.431470 -117.221780 0.00 0.00
- DAC LB 36.277000 -117.593700 0.00 0.00
- MVU LB 38.503750 -112.212330 0.00 0.00
- TPH LB 38.075000 -117.222500 0.00 0.00
- LTH LI 46.461700 -119.417680 0.00 0.00
- LTL LI 30.537400 -90.766000 0.00 0.00
- BJOB NC 36.618083 -121.314683 0.00 0.00
- CADB NC 37.163460 -121.626450 0.00 0.00
- CAG NC 37.857270 -122.422660 0.00 0.00
- CAL NC 37.450348 -121.800367 0.00 0.00
- CBP NC 37.747665 -121.708427 0.00 0.00
- CCH1 NC 37.743318 -122.096570 0.00 0.00
- CDOB NC 37.729400 -121.833460 0.00 0.00
- CGP1 NC 37.645448 -122.011278 0.00 0.00
- CHR NC 37.372883 -121.768867 0.00 0.00
- CLCB NC 37.737867 -122.064867 0.00 0.00
- CMOB NC 37.810417 -121.802717 0.00 0.00
- CMW1 NC 37.540530 -121.887425 0.00 0.00
- CPI NC 37.988903 -122.215605 0.00 0.00
- CRH NC 37.858840 -121.992640 0.00 0.00
- CRPB NC 37.912330 -121.907160 0.00 0.00
- CSU1 NC 37.643032 -121.940202 0.00 0.00
- CTA NC 38.026908 -122.015992 0.00 0.00
- CYD1 NC 37.562893 -122.096695 0.00 0.00
- GAXB NC 38.710733 -122.757333 0.00 0.00
- GCVB NC 38.769267 -123.014867 0.00 0.00
- GDXB NC 38.807775 -122.794742 0.00 0.00
- HCAB NC 37.025367 -121.484533 0.00 0.00
- JBG NC 37.343048 -122.339760 0.00 0.00
- JBMB NC 37.318607 -122.153145 0.00 0.00
- JBR NC 37.172170 -122.011210 0.00 0.00
- JECB NC 37.049225 -121.807237 0.00 0.00
- JELB NC 36.927450 -121.827633 0.00 0.00
- JFP NC 37.360447 -122.189857 0.00 0.00
- JGR NC 37.516042 -122.458153 0.00 0.00
- JHU NC 37.214270 -122.171680 0.00 0.00
- JJO NC 37.206310 -121.977260 0.00 0.00
- JLAB NC 37.154320 -121.734670 0.00 0.00
- JMGB NC 37.637190 -122.474820 0.00 0.00
- JPC NC 37.270800 -122.212760 0.00 0.00
- JPR NC 37.795270 -122.474950 0.00 0.00
- JPSB NC 37.198550 -122.346710 0.00 0.00
- JSA NC 37.582248 -122.420058 0.00 0.00
- JSFB NC 37.403720 -122.175960 0.00 0.00
- JSGB NC 37.283972 -122.050168 0.00 0.00
- JSP NC 37.586442 -122.502168 0.00 0.00
- KBO NC 42.212358 -124.225983 0.00 0.00
- KCPB NC 39.686308 -123.582422 0.00 0.00
- KCT NC 40.475550 -124.337393 0.00 0.00
- KEB NC 42.872210 -124.334248 0.00 0.00
- KHBB NC 40.659902 -123.219660 0.00 0.00
- KHMB NC 40.874833 -123.732500 0.00 0.00
- KMPB NC 40.417195 -124.120755 0.00 0.00
- KMR NC 40.202177 -123.708632 0.00 0.00
- KRMB NC 41.522958 -123.907972 0.00 0.00
- KRP NC 41.157653 -124.023297 0.00 0.00
- KSXB NC 41.830378 -123.876883 0.00 0.00
- LDH NC 40.511900 -121.581600 0.00 0.00
- MBS1 NC 37.761023 -118.944542 0.00 0.00
- MCB NC 37.644395 -118.896835 0.00 0.00
- MDP1 NC 37.632460 -119.079720 0.00 0.00
- MMLB NC 37.659093 -118.979065 0.00 0.00
- MMX1 NC 37.617115 -118.959752 0.00 0.00
- NAPC NC 38.439507 -122.252717 0.00 0.00
- NBO NC 37.921440 -122.718055 0.00 0.00
- NBRB NC 38.260208 -122.551788 0.00 0.00
- NCC NC 38.001983 -122.474667 0.00 0.00
- NEA NC 38.309883 -122.943002 0.00 0.00
- NEH NC 38.336560 -122.810820 0.00 0.00
- NFV NC 38.467690 -122.896510 0.00 0.00
- NGVB NC 38.280430 -122.215820 0.00 0.00
- NHM NC 38.154863 -121.804385 0.00 0.00
- NHS NC 38.656517 -122.614417 0.00 0.00
- NHV NC 38.164880 -122.712960 0.00 0.00
- NLH NC 38.122815 -122.149213 0.00 0.00
- NMH NC 38.669370 -122.633212 0.00 0.00
- NMI NC 38.076013 -122.258502 0.00 0.00
- NOLB NC 38.042480 -122.796837 0.00 0.00
- NPRB NC 37.996422 -123.019493 0.00 0.00
- NSM NC 38.348787 -122.577952 0.00 0.00
- NSP NC 38.200252 -122.464403 0.00 0.00
- NTAB NC 37.923070 -122.598290 0.00 0.00
- NTAC NC 37.923717 -122.596483 0.00 0.00
- NTO NC 38.143608 -122.449692 0.00 0.00
- NTR NC 38.240390 -122.798420 0.00 0.00
- NTYB NC 38.389210 -122.662970 0.00 0.00
- PAGB NC 35.730650 -120.249860 0.00 0.00
- PHOB NC 35.866610 -120.479603 0.00 0.00
- PHSB NC 35.824040 -121.053930 0.00 0.00
- PMPB NC 36.215905 -120.801263 0.00 0.00
- PSM NC 36.068807 -120.596180 0.00 0.00
- ADH NN 37.968300 -118.715300 0.00 0.00
- ANT NN 37.917800 -118.564000 0.00 0.00
- BAB NN 39.601300 -120.104000 0.00 0.00
- BCK NN 37.697500 -118.372000 0.00 0.00
- BEK NN 39.866700 -120.358600 0.00 0.00
- BEN NN 37.715500 -118.573100 0.00 0.00
- BFC NN 38.894000 -119.607700 0.00 0.00
- BHP NN 37.299500 -118.487300 0.00 0.00
- BMN NN 40.431300 -117.221700 0.00 0.00
- BMR NN 40.108700 -120.291000 0.00 0.00
- BON NN 37.955200 -118.301700 0.00 0.00
- CAS NN 37.574800 -118.551500 0.00 0.00
- CC12 NN 36.133700 -115.170900 0.00 0.00
- CC16 NN 36.161600 -115.061800 0.00 0.00
- CC22 NN 36.113400 -115.238200 0.00 0.00
- CCAD NN 39.164200 -119.778600 0.00 0.00
- CCHI NN 39.168100 -119.747300 0.00 0.00
- CF01 NN 39.159300 -119.764700 0.00 0.00
- CF02 NN 39.190700 -119.741800 0.00 0.00
- CF03 NN 39.125700 -119.767600 0.00 0.00
- CHYN NN 36.222300 -115.105700 0.00 0.00
- CLK NN 37.590700 -118.824200 0.00 0.00
- CLRK NN 36.084600 -115.204700 0.00 0.00
- CWC NN 36.439900 -118.080200 0.00 0.00
- DIX NN 39.802200 -118.082000 0.00 0.00
- DNY NN 41.086300 -119.277800 0.00 0.00
- DOEL NN 36.209700 -115.139300 0.00 0.00
- DON NN 39.351700 -120.320500 0.00 0.00
- DSP NN 37.369500 -117.974200 0.00 0.00
- EBP NN 38.582800 -119.806300 0.00 0.00
- EGLV NN 39.158000 -119.719600 0.00 0.00
- EMB NN 38.974800 -120.101900 0.00 0.00
- EUR NN 39.754200 -120.710800 0.00 0.00
- FLV NN 37.814000 -117.909500 0.00 0.00
- FMT NN 36.637500 -116.782600 0.00 0.00
- FPK NN 39.224700 -118.151000 0.00 0.00
- GMN NN 37.300400 -117.259800 0.00 0.00
- GMR NN 34.784600 -115.659900 0.00 0.00
- GNO NN 38.929200 -119.852800 0.00 0.00
- GSLV NN 36.037900 -115.044200 0.00 0.00
- GVN NN 36.999900 -117.343500 0.00 0.00
- GWY NN 36.186000 -116.668900 0.00 0.00
- GZY NN 39.962000 -120.650200 0.00 0.00
- HCK NN 38.075500 -118.592200 0.00 0.00
- HF05 NN 36.030700 -115.089400 0.00 0.00
- HTC NN 37.529800 -118.771000 0.00 0.00
- HVGC NN 39.489700 -119.720900 0.00 0.00
- HYDE NN 36.160900 -115.192300 0.00 0.00
- HYX NN 39.772800 -117.763300 0.00 0.00
- IND NN 39.434300 -120.291700 0.00 0.00
- KBF NN 39.506800 -120.211800 0.00 0.00
- KVN NN 39.051000 -118.100000 0.00 0.00
- LCH NN 37.233000 -117.645000 0.00 0.00
- LHV NN 38.251200 -118.504500 0.00 0.00
- LOVE NN 39.519900 -119.806200 0.00 0.00
- LOY NN 39.658700 -120.235800 0.00 0.00
- LUL NN 38.052300 -119.180300 0.00 0.00
- LV04 NN 36.160200 -115.131700 0.00 0.00
- LV07 NN 36.181600 -115.314600 0.00 0.00
- LV09 NN 36.247500 -115.241900 0.00 0.00
- LVO NN 39.871700 -120.871800 0.00 0.00
- MCA NN 36.646600 -117.280500 0.00 0.00
- MCC NN 37.919500 -119.024300 0.00 0.00
- MGN NN 37.813300 -118.695500 0.00 0.00
- MIL NN 38.024800 -118.185800 0.00 0.00
- MLN NN 37.593800 -118.342700 0.00 0.00
- MMC NN 38.360800 -119.128300 0.00 0.00
- MNA NN 38.432800 -118.156700 0.00 0.00
- MND1 NN 38.906800 -119.730400 0.00 0.00
- MOHS NN 36.254700 -115.138400 0.00 0.00
- MPK NN 39.295700 -120.030200 0.00 0.00
- MPT NN 38.062500 -118.775300 0.00 0.00
- MZP NN 37.700700 -117.383400 0.00 0.00
- NELL NN 36.220300 -115.034700 0.00 0.00
- NEN NN 35.649500 -114.843300 0.00 0.00
- NMHS NN 39.530900 -119.775400 0.00 0.00
- NOAA NN 39.568100 -119.795800 0.00 0.00
- OMM NN 37.621400 -118.999400 0.00 0.00
- ORC NN 37.635300 -118.656000 0.00 0.00
- PAH NN 39.706500 -119.384200 0.00 0.00
- PAN NN 36.393200 -117.099000 0.00 0.00
- PEA NN 39.607600 -119.960300 0.00 0.00
- PICO NN 39.431200 -119.775600 0.00 0.00
- PNT NN 39.089200 -119.598700 0.00 0.00
- POC NN 37.362300 -118.279200 0.00 0.00
- POW NN 38.409500 -118.631800 0.00 0.00
- PRN NN 37.406500 -115.050300 0.00 0.00
- QSM NN 35.965000 -116.868200 0.00 0.00
- RCC NN 37.487700 -118.721700 0.00 0.00
- RF05 NN 39.509500 -119.836000 0.00 0.00
- RF08 NN 39.542700 -119.856100 0.00 0.00
- RFMA NN 39.518800 -119.901000 0.00 0.00
- RFNV NN 39.574200 -119.827500 0.00 0.00
- RRK NN 36.065700 -115.466000 0.00 0.00
- RUB NN 39.051800 -120.153900 0.00 0.00
- RYN NN 38.628100 -118.523000 0.00 0.00
- SAT NN 39.602500 -120.447000 0.00 0.00
- SBT NN 39.626800 -120.665800 0.00 0.00
- SCH NN 37.365800 -118.687000 0.00 0.00
- SF02 NN 39.555700 -119.733400 0.00 0.00
- SHP NN 36.505500 -115.160200 0.00 0.00
- SJC NN 38.349200 -119.439200 0.00 0.00
- SKYF NN 39.482500 -119.834000 0.00 0.00
- SLK NN 37.834000 -119.128700 0.00 0.00
- SMI NN 39.867300 -120.529500 0.00 0.00
- SMRN NN 39.537200 -119.727400 0.00 0.00
- SPHI NN 39.543600 -119.757300 0.00 0.00
- SSFS NN 39.637000 -119.708300 0.00 0.00
- SWTP NN 39.516400 -119.703500 0.00 0.00
- TAH NN 39.150000 -120.163000 0.00 0.00
- TMO NN 36.804800 -117.403500 0.00 0.00
- TNK NN 39.267500 -120.235800 0.00 0.00
- TNP NN 38.082000 -117.217700 0.00 0.00
- TRC NN 38.349700 -115.585200 0.00 0.00
- UNLV NN 36.108200 -115.141600 0.00 0.00
- UNRN NN 39.527000 -119.818000 0.00 0.00
- UNRX NN 39.514100 -119.749300 0.00 0.00
- VIP NN 39.754000 -119.460800 0.00 0.00
- VOTK NN 36.080900 -115.069900 0.00 0.00
- VPK NN 39.474700 -120.037300 0.00 0.00
- WAK NN 38.504400 -119.437200 0.00 0.00
- WCK NN 37.527800 -118.650700 0.00 0.00
- WCN NN 39.301700 -119.756300 0.00 0.00
- WCT NN 36.792800 -116.625700 0.00 0.00
- WGLF NN 39.498600 -119.821600 0.00 0.00
- WHR NN 40.036700 -118.361200 0.00 0.00
- WIL NN 38.810400 -119.908300 0.00 0.00
- WMD NN 37.443500 -118.637000 0.00 0.00
- WTNK NN 36.191600 -115.009400 0.00 0.00
- WVA NN 39.944500 -119.824000 0.00 0.00
- WYRD NN 39.492700 -119.762800 0.00 0.00
- YER NN 38.984500 -119.240700 0.00 0.00
- 0655 NP 34.313130 -118.499180 0.00 0.00
- 1662 NP 37.840000 -122.296667 0.00 0.00
- 1676 NP 37.791667 -122.415000 0.00 0.00
- 1745 NP 37.456910 -122.170010 0.00 0.00
- 1761 NP 38.290000 -122.456667 0.00 0.00
- 1767 NP 38.440670 -122.704558 0.00 0.00
- 1768 NP 38.233635 -122.636712 0.00 0.00
- 1770 NP 37.760207 -122.415235 0.00 0.00
- 1771 NP 37.777368 -122.396888 0.00 0.00
- 1772 NP 37.750942 -122.490712 0.00 0.00
- 1773 NP 37.763863 -122.473602 0.00 0.00
- 1774 NP 37.761463 -122.504473 0.00 0.00
- 1775 NP 37.407013 -122.094607 0.00 0.00
- 1776 NP 37.474828 -122.239235 0.00 0.00
- 1777 NP 37.541873 -122.304710 0.00 0.00
- 1778 NP 37.302493 -121.979782 0.00 0.00
- 1779 NP 37.267475 -121.866562 0.00 0.00
- 1780 NP 37.263518 -121.916168 0.00 0.00
- 1781 NP 37.301433 -121.903900 0.00 0.00
- 1782 NP 37.318157 -121.943535 0.00 0.00
- 1783 NP 37.306598 -121.849152 0.00 0.00
- 1785 NP 37.660808 -121.873663 0.00 0.00
- 1786 NP 37.680305 -121.739572 0.00 0.00
- 1788 NP 37.417018 -121.915405 0.00 0.00
- 1789 NP 37.392333 -121.984667 0.00 0.00
- 1790 NP 37.590942 -122.367405 0.00 0.00
- 1792 NP 37.747660 -122.425277 0.00 0.00
- 1793 NP 37.346535 -121.807777 0.00 0.00
- 1794 NP 37.407182 -122.001070 0.00 0.00
- 1795 NP 37.746650 -122.386748 0.00 0.00
- 1796 NP 37.394145 -122.060413 0.00 0.00
- 1798 NP 37.376015 -121.843742 0.00 0.00
- 1799 NP 37.271133 -121.986237 0.00 0.00
- 1801 NP 37.450958 -122.428703 0.00 0.00
- 1802 NP 37.592678 -122.074493 0.00 0.00
- 1803 NP 37.717368 -121.727005 0.00 0.00
- 1804 NP 37.679253 -122.484045 0.00 0.00
- 1805 NP 37.784333 -122.300833 0.00 0.00
- 1806 NP 37.910328 -122.536818 0.00 0.00
- 1807 NP 37.887667 -122.299167 0.00 0.00
- 1813 NP 38.034333 -122.270167 0.00 0.00
- 1815 NP 37.489080 -122.263690 0.00 0.00
- 1816 NP 37.744833 -122.401333 0.00 0.00
- 1817 NP 37.779333 -122.486000 0.00 0.00
- 1818 NP 37.716333 -122.431167 0.00 0.00
- 1819 NP 37.810760 -122.292522 0.00 0.00
- 1820 NP 37.825098 -122.276065 0.00 0.00
- 1821 NP 37.798600 -122.268400 0.00 0.00
- 1823 NP 37.994480 -122.114720 0.00 0.00
- 1824 NP 37.724163 -122.448410 0.00 0.00
- 1825 NP 37.413953 -121.870850 0.00 0.00
- 1826 NP 37.680600 -121.920977 0.00 0.00
- 1827 NP 37.860592 -122.293267 0.00 0.00
- 1828 NP 37.858345 -122.252518 0.00 0.00
- 1829 NP 38.290272 -122.461250 0.00 0.00
- 1830 NP 37.374177 -122.078102 0.00 0.00
- 1831 NP 37.364082 -122.018175 0.00 0.00
- 1834 NP 37.425487 -122.077717 0.00 0.00
- 1835 NP 38.442097 -122.606800 0.00 0.00
- 1836 NP 37.691692 -122.169825 0.00 0.00
- 1837 NP 37.232935 -121.771860 0.00 0.00
- 1838 NP 37.383317 -121.910883 0.00 0.00
- 1839 NP 37.285373 -121.936020 0.00 0.00
- 1841 NP 37.358340 -121.977750 0.00 0.00
- 1842 NP 37.688730 -121.783512 0.00 0.00
- 1843 NP 37.726883 -122.150917 0.00 0.00
- 1844 NP 37.885200 -122.032167 0.00 0.00
- 1845 NP 38.578050 -122.578383 0.00 0.00
- 1846 NP 37.942150 -121.930467 0.00 0.00
- 1847 NP 38.012983 -122.134167 0.00 0.00
- 1848 NP 38.367433 -122.524033 0.00 0.00
- 1863 NP 37.666717 -121.857133 0.00 0.00
- 1865 NP 36.882867 -121.390717 0.00 0.00
- 286 NP 32.954768 -115.824303 0.00 0.00
- 5030 NP 34.521333 -117.991000 0.00 0.00
- 5062 NP 33.177470 -115.615640 0.00 0.00
- 5081 NP 34.084000 -118.600000 0.00 0.00
- 5271 NP 33.352667 -115.733000 0.00 0.00
- 5398 NP 34.194833 -118.353667 0.00 0.00
- 5399 NP 33.943000 -118.411000 0.00 0.00
- 5401 NP 34.397000 -118.918000 0.00 0.00
- 5402 NP 34.076000 -118.395000 0.00 0.00
- 5403 NP 34.396833 -119.298667 0.00 0.00
- 5404 NP 34.155000 -118.641000 0.00 0.00
- 5442 NP 34.049048 -116.940323 0.00 0.00
- 5443 NP 33.844470 -116.310790 0.00 0.00
- 5444 NP 33.634610 -116.095340 0.00 0.00
- 5445 NP 35.061310 -119.027940 0.00 0.00
- 5470 NP 35.060810 -119.026440 0.00 0.00
- 5471 NP 34.199880 -118.362010 0.00 0.00
- BBA NP 34.195510 -118.353400 0.00 0.00
- BBB NP 33.352620 -115.733250 0.00 0.00
- BVH NP 34.076180 -118.395900 0.00 0.00
- FLL NP 34.397270 -118.918070 0.00 0.00
- GRF NP 34.119200 -118.300400 0.00 0.00
- JAB NP 34.310920 -118.496700 0.00 0.00
- JGB NP 34.313130 -118.499180 0.00 0.00
- LAX NP 33.943820 -118.413900 0.00 0.00
- LT2 NP 34.521490 -117.990990 0.00 0.00
- LTR NP 34.521100 -117.990300 0.00 0.00
- OKV NP 34.396970 -119.298740 0.00 0.00
- TCF NP 34.083770 -118.599000 0.00 0.00
- NE70 NR 32.420950 -115.260780 0.00 0.00
- NE71 NR 31.689730 -115.905260 0.00 0.00
- NE72 NR 30.848430 -116.058570 0.00 0.00
- NE73 NR 30.065100 -115.348470 0.00 0.00
- NE74 NR 28.007510 -114.013800 0.00 0.00
- NE75 NR 27.293340 -112.856490 0.00 0.00
- NE76 NR 26.888940 -111.999050 0.00 0.00
- NE77 NR 26.015780 -111.361330 0.00 0.00
- NE78 NR 24.398200 -111.106430 0.00 0.00
- NE79 NR 23.119370 -109.756110 0.00 0.00
- NE80 NR 30.500000 -112.319930 0.00 0.00
- NE81 NR 28.918340 -109.636260 0.00 0.00
- NE82 NR 26.915660 -109.230840 0.00 0.00
- NE83 NR 24.730880 -107.739330 0.00 0.00
- B081 PB 33.711170 -116.714170 0.00 0.00
- B082 PB 33.598180 -116.596010 0.00 0.00
- B084 PB 33.611570 -116.456370 0.00 0.00
- B086 PB 33.557500 -116.530000 0.00 0.00
- B087 PB 33.495500 -116.602670 0.00 0.00
- DCD PG 35.212167 -120.840833 0.00 0.00
- DPD PG 35.233000 -120.781667 0.00 0.00
- SHD PG 35.168167 -120.761333 0.00 0.00
- WRD PG 35.458833 -120.881667 0.00 0.00
- BAHB RB 28.944000 -113.562000 0.00 0.00
- CHXB RB 31.472000 -115.051000 0.00 0.00
- GUAY RB 27.898660 -100.872720 0.00 0.00
- PPXB RB 31.335000 -113.632000 0.00 0.00
- TOPB RB 25.609000 -109.535000 0.00 0.00
- GVDA SB 33.668900 -116.673798 0.00 0.00
- WLA SB 33.097598 -115.530602 0.00 0.00
- AL5 SN 36.859600 -116.454700 0.00 0.00
- AL5S SN 36.859600 -116.454700 0.00 0.00
- AMD SN 36.452600 -116.280900 0.00 0.00
- BTW SN 36.997800 -116.566500 0.00 0.00
- CAF SN 36.839100 -116.337700 0.00 0.00
- CRF SN 36.811800 -116.534000 0.00 0.00
- DOM SN 37.002100 -116.408600 0.00 0.00
- ECO SN 37.210800 -116.329600 0.00 0.00
- EXHS SN 36.849500 -116.429400 0.00 0.00
- FMW SN 36.902100 -116.368800 0.00 0.00
- FRG SN 36.816900 -116.419500 0.00 0.00
- HEL SN 36.724600 -116.975000 0.00 0.00
- LEC SN 36.562700 -116.689600 0.00 0.00
- LSC SN 36.730700 -116.325500 0.00 0.00
- LWLS SN 36.644500 -116.397600 0.00 0.00
- MDVS SN 36.851900 -116.421400 0.00 0.00
- NCF SN 36.889900 -116.568200 0.00 0.00
- NI5 SN 36.852300 -116.461900 0.00 0.00
- NI5S SN 36.852300 -116.461900 0.00 0.00
- PIT SN 36.679800 -116.493700 0.00 0.00
- PUV SN 36.949400 -115.963300 0.00 0.00
- RED SN 36.689500 -116.093000 0.00 0.00
- RF13D SN 36.852900 -116.424900 0.00 0.00
- RF13S SN 36.852900 -116.424900 0.00 0.00
- RF15D SN 36.853700 -116.426700 0.00 0.00
- RF15S SN 36.853700 -116.426700 0.00 0.00
- RF16D SN 36.851700 -116.425800 0.00 0.00
- RF16S SN 36.851700 -116.425800 0.00 0.00
- RPY SN 36.851500 -116.456300 0.00 0.00
- SCF SN 36.756800 -116.544000 0.00 0.00
- SGR SN 36.980500 -117.032700 0.00 0.00
- SME SN 36.828600 -116.446800 0.00 0.00
- SMES SN 36.828600 -116.446800 0.00 0.00
- SPC SN 36.674600 -116.203000 0.00 0.00
- STC SN 37.293900 -116.435800 0.00 0.00
- STH SN 36.645700 -116.337500 0.00 0.00
- STO SN 36.860300 -116.474200 0.00 0.00
- SYM SN 36.741600 -116.446000 0.00 0.00
- SYMS SN 36.837700 -116.472300 0.00 0.00
- TAR SN 36.868000 -116.632200 0.00 0.00
- TIM SN 37.066700 -116.469400 0.00 0.00
- TPW SN 36.901600 -116.251900 0.00 0.00
- TWP SN 37.204700 -116.123400 0.00 0.00
- TYM SN 37.144100 -116.720800 0.00 0.00
- TYMS SN 36.839400 -116.467500 0.00 0.00
- WLD SN 36.792700 -116.625700 0.00 0.00
- YCW SN 36.922400 -116.475600 0.00 0.00
- YFT SN 37.076200 -115.973500 0.00 0.00
- 109C TA 32.888900 -117.105100 0.00 0.00
- 112A TA 32.535600 -114.580400 0.00 0.00
- 113A TA 32.768300 -113.766700 0.00 0.00
- 115A TA 32.700600 -112.227900 0.00 0.00
- 116A TA 32.561800 -111.704200 0.00 0.00
- 117A TA 32.571600 -110.739300 0.00 0.00
- 118A TA 32.639900 -109.969700 0.00 0.00
- 119A TA 32.766300 -109.302900 0.00 0.00
- 214A TA 31.955900 -112.811500 0.00 0.00
- 216A TA 32.002200 -111.457400 0.00 0.00
- 217A TA 31.774800 -110.816200 0.00 0.00
- 218A TA 31.973700 -110.046400 0.00 0.00
- 219A TA 31.998900 -109.259200 0.00 0.00
- 318A TA 31.439000 -109.990700 0.00 0.00
- 319A TA 31.375700 -109.280900 0.00 0.00
- A04A TA 48.719700 -122.707000 0.00 0.00
- A05A TA 48.998300 -122.085200 0.00 0.00
- A06A TA 49.098300 -121.480400 0.00 0.00
- A07A TA 49.048000 -120.383700 0.00 0.00
- A08A TA 48.953400 -119.272500 0.00 0.00
- A09A TA 48.975200 -118.585100 0.00 0.00
- A10A TA 48.981300 -117.558600 0.00 0.00
- A11A TA 48.958000 -116.363100 0.00 0.00
- A12A TA 48.933900 -115.652600 0.00 0.00
- A13A TA 48.933000 -114.413500 0.00 0.00
- B04A TA 48.057500 -123.504000 0.00 0.00
- B05A TA 48.264100 -122.096000 0.00 0.00
- B06A TA 48.518300 -121.484500 0.00 0.00
- B07A TA 48.460700 -120.119700 0.00 0.00
- B08A TA 48.358400 -119.333200 0.00 0.00
- B09A TA 48.422600 -118.148800 0.00 0.00
- B10A TA 48.299500 -117.225200 0.00 0.00
- B11A TA 48.436600 -116.367400 0.00 0.00
- B12A TA 48.469400 -115.588900 0.00 0.00
- B13A TA 48.373900 -114.467600 0.00 0.00
- BNLO TA 37.131100 -122.172900 0.00 0.00
- C03A TA 47.947500 -124.566200 0.00 0.00
- C04A TA 47.717200 -122.971800 0.00 0.00
- C05A TA 47.694700 -121.689500 0.00 0.00
- C06A TA 47.923000 -120.894400 0.00 0.00
- C07A TA 47.690200 -120.061300 0.00 0.00
- C08A TA 47.784900 -119.045400 0.00 0.00
- C09A TA 47.800600 -118.274100 0.00 0.00
- C10A TA 47.819200 -117.307700 0.00 0.00
- C11A TA 47.839700 -116.255700 0.00 0.00
- C12A TA 47.703900 -115.481100 0.00 0.00
- C12B TA 47.709200 -115.475000 0.00 0.00
- C13A TA 47.679500 -114.573500 0.00 0.00
- C14A TA 47.768700 -113.746800 0.00 0.00
- D03A TA 47.116100 -123.771400 0.00 0.00
- D04A TA 47.110900 -122.803500 0.00 0.00
- D05A TA 47.188100 -121.988800 0.00 0.00
- D06A TA 47.193500 -120.844500 0.00 0.00
- D07A TA 47.193600 -119.973200 0.00 0.00
- D08A TA 47.057300 -118.921000 0.00 0.00
- D09A TA 47.061500 -118.309200 0.00 0.00
- D10A TA 47.054700 -117.275900 0.00 0.00
- D11A TA 47.047700 -116.335300 0.00 0.00
- D12A TA 47.052800 -115.352500 0.00 0.00
- D13A TA 47.094300 -114.458700 0.00 0.00
- D14A TA 47.081600 -113.506200 0.00 0.00
- D15A TA 47.041800 -112.520400 0.00 0.00
- E03A TA 46.545900 -123.563200 0.00 0.00
- E04A TA 46.593400 -122.720400 0.00 0.00
- E05A TA 46.560500 -121.760800 0.00 0.00
- E06A TA 46.543300 -120.979100 0.00 0.00
- E07A TA 46.558500 -119.854800 0.00 0.00
- E08A TA 46.491100 -119.059500 0.00 0.00
- E09A TA 46.513900 -118.145500 0.00 0.00
- E10A TA 46.485000 -117.110200 0.00 0.00
- E11A TA 46.355900 -116.209200 0.00 0.00
- E12A TA 46.415100 -115.570500 0.00 0.00
- E13A TA 46.442300 -114.188000 0.00 0.00
- E14A TA 46.416200 -113.492600 0.00 0.00
- E15A TA 46.424600 -112.641200 0.00 0.00
- ELFS TA 40.618300 -120.727900 0.00 0.00
- F03A TA 45.930600 -123.559100 0.00 0.00
- F04A TA 45.932500 -122.419300 0.00 0.00
- F05A TA 45.883900 -121.458800 0.00 0.00
- F06A TA 45.766400 -120.781700 0.00 0.00
- F07A TA 45.895100 -119.927800 0.00 0.00
- F08A TA 45.797100 -118.776800 0.00 0.00
- F09A TA 45.708700 -117.909400 0.00 0.00
- F10A TA 45.973400 -117.227700 0.00 0.00
- F11A TA 45.888300 -116.154900 0.00 0.00
- F12A TA 45.756500 -115.254700 0.00 0.00
- F13A TA 45.789300 -114.331700 0.00 0.00
- F14A TA 45.812200 -113.369900 0.00 0.00
- F15A TA 45.840900 -112.492900 0.00 0.00
- G03A TA 45.315300 -123.281100 0.00 0.00
- G04A TA 45.205900 -122.478400 0.00 0.00
- G05A TA 45.242200 -121.316700 0.00 0.00
- G06A TA 45.235700 -120.635100 0.00 0.00
- G07A TA 45.266400 -119.669200 0.00 0.00
- G08A TA 45.290400 -118.959500 0.00 0.00
- G09A TA 45.278100 -117.780200 0.00 0.00
- G10A TA 45.292500 -117.120400 0.00 0.00
- G11A TA 45.399700 -116.268000 0.00 0.00
- G13A TA 45.093100 -114.232900 0.00 0.00
- G14A TA 45.243200 -113.460400 0.00 0.00
- G15A TA 45.166000 -112.488700 0.00 0.00
- H02A TA 44.676400 -123.999700 0.00 0.00
- H03A TA 44.676500 -123.292300 0.00 0.00
- H04A TA 44.683700 -122.186200 0.00 0.00
- H05A TA 44.647300 -121.227100 0.00 0.00
- H06A TA 44.734400 -120.334600 0.00 0.00
- H07A TA 44.591300 -119.564600 0.00 0.00
- H08A TA 44.519100 -118.670400 0.00 0.00
- H09A TA 44.665300 -117.664300 0.00 0.00
- H10A TA 44.589000 -116.747400 0.00 0.00
- H11A TA 44.703500 -116.012700 0.00 0.00
- H12A TA 44.549400 -114.855400 0.00 0.00
- H13A TA 44.564200 -114.254500 0.00 0.00
- H14A TA 44.616500 -113.367400 0.00 0.00
- H15A TA 44.617300 -112.643900 0.00 0.00
- H16A TA 44.703800 -111.247800 0.00 0.00
- HAST TA 36.388700 -121.551400 0.00 0.00
- HATC TA 40.816100 -121.461200 0.00 0.00
- HELL TA 36.680100 -119.022800 0.00 0.00
- I02A TA 44.003500 -123.829900 0.00 0.00
- I03A TA 43.972600 -123.277700 0.00 0.00
- I04A TA 43.794100 -122.411300 0.00 0.00
- I05A TA 44.162600 -121.267500 0.00 0.00
- I06A TA 43.943700 -120.210700 0.00 0.00
- I07A TA 44.082000 -119.504200 0.00 0.00
- I08A TA 43.910400 -118.569200 0.00 0.00
- I09A TA 43.973300 -117.740700 0.00 0.00
- I10A TA 44.086000 -116.802900 0.00 0.00
- I11A TA 43.912100 -115.957800 0.00 0.00
- I13A TA 43.914600 -114.116900 0.00 0.00
- I14A TA 43.928600 -113.451800 0.00 0.00
- I15A TA 43.999700 -112.485000 0.00 0.00
- I16A TA 43.875600 -111.486800 0.00 0.00
- J02A TA 43.365400 -123.574700 0.00 0.00
- J03A TA 43.371700 -122.964600 0.00 0.00
- J04A TA 43.240600 -122.109300 0.00 0.00
- J05A TA 43.284300 -121.235700 0.00 0.00
- J06A TA 43.251500 -120.152800 0.00 0.00
- J07A TA 43.373600 -119.310500 0.00 0.00
- J08A TA 43.358000 -118.474300 0.00 0.00
- J09A TA 43.347100 -117.754200 0.00 0.00
- J10A TA 43.427500 -116.767000 0.00 0.00
- J11A TA 43.415100 -115.827800 0.00 0.00
- J12A TA 43.250000 -115.098000 0.00 0.00
- J13A TA 43.397900 -114.174200 0.00 0.00
- J14A TA 43.323400 -113.517800 0.00 0.00
- J15A TA 43.399800 -112.433400 0.00 0.00
- J16A TA 43.274100 -111.611900 0.00 0.00
- K01A TA 42.808800 -124.469200 0.00 0.00
- K02A TA 42.766700 -123.489800 0.00 0.00
- K04A TA 42.613300 -121.730700 0.00 0.00
- K05A TA 42.725800 -120.893400 0.00 0.00
- K06A TA 42.799100 -120.251000 0.00 0.00
- K07A TA 42.691200 -119.246900 0.00 0.00
- K08A TA 42.730800 -118.486400 0.00 0.00
- K09A TA 42.699600 -117.725100 0.00 0.00
- K10A TA 42.777900 -116.870500 0.00 0.00
- K11A TA 42.771300 -116.032300 0.00 0.00
- K12A TA 42.636000 -114.902900 0.00 0.00
- K13A TA 42.649300 -114.084000 0.00 0.00
- K14A TA 42.545200 -113.176000 0.00 0.00
- K15A TA 42.685200 -112.530500 0.00 0.00
- K16A TA 42.832100 -111.588400 0.00 0.00
- L02A TA 42.156000 -123.601800 0.00 0.00
- L04A TA 42.175000 -121.891000 0.00 0.00
- L05A TA 42.047200 -120.833600 0.00 0.00
- L07A TA 42.018700 -119.339800 0.00 0.00
- L08A TA 42.190100 -118.344500 0.00 0.00
- L09A TA 42.019200 -117.666700 0.00 0.00
- L10A TA 42.077300 -116.471100 0.00 0.00
- L11A TA 42.166900 -115.754100 0.00 0.00
- L12A TA 42.146000 -115.016200 0.00 0.00
- L13A TA 42.088600 -113.944400 0.00 0.00
- L14A TA 42.034300 -113.239800 0.00 0.00
- L15A TA 42.004100 -112.386000 0.00 0.00
- L16A TA 42.014900 -111.431900 0.00 0.00
- LAVA TA 38.755200 -120.739500 0.00 0.00
- M01C TA 41.847300 -124.122100 0.00 0.00
- M02C TA 41.392000 -122.853800 0.00 0.00
- M03C TA 41.274200 -122.122000 0.00 0.00
- M04C TA 41.782600 -121.839300 0.00 0.00
- M05C TA 41.359300 -121.145700 0.00 0.00
- M06C TA 41.204700 -120.477200 0.00 0.00
- M07A TA 41.388400 -119.171100 0.00 0.00
- M08A TA 41.448300 -118.379200 0.00 0.00
- M09A TA 41.423000 -117.449700 0.00 0.00
- M10A TA 41.522000 -116.539600 0.00 0.00
- M11A TA 41.431100 -115.791200 0.00 0.00
- M12A TA 41.415700 -114.915100 0.00 0.00
- M13A TA 41.360200 -114.165500 0.00 0.00
- M14A TA 41.503000 -113.347100 0.00 0.00
- M15A TA 41.463200 -112.447700 0.00 0.00
- M16A TA 41.314600 -111.629800 0.00 0.00
- N02C TA 40.822000 -123.305700 0.00 0.00
- N06A TA 40.748400 -119.834600 0.00 0.00
- N07A TA 40.772100 -118.971600 0.00 0.00
- N07B TA 40.779700 -118.971100 0.00 0.00
- N08A TA 40.781100 -118.133700 0.00 0.00
- N09A TA 40.852000 -117.524400 0.00 0.00
- N10A TA 40.718600 -116.508000 0.00 0.00
- N11A TA 40.816700 -115.736000 0.00 0.00
- N12A TA 40.852200 -115.038700 0.00 0.00
- N13A TA 40.855600 -114.204500 0.00 0.00
- N14A TA 40.851300 -113.186700 0.00 0.00
- N15A TA 40.890300 -112.520100 0.00 0.00
- N16A TA 40.886900 -111.437000 0.00 0.00
- N17A TA 40.942500 -110.833500 0.00 0.00
- N18A TA 40.976300 -109.673100 0.00 0.00
- N19A TA 40.893600 -109.177200 0.00 0.00
- O01C TA 40.139800 -123.819500 0.00 0.00
- O02C TA 40.176700 -122.788300 0.00 0.00
- O03C TA 39.997100 -122.031900 0.00 0.00
- O04C TA 40.320000 -121.086000 0.00 0.00
- O05C TA 39.962000 -120.917600 0.00 0.00
- O06A TA 40.165000 -119.827500 0.00 0.00
- O07A TA 40.161400 -118.877200 0.00 0.00
- O08A TA 40.290300 -118.155000 0.00 0.00
- O09A TA 40.169700 -117.189900 0.00 0.00
- O10A TA 40.291700 -116.499900 0.00 0.00
- O11A TA 40.131300 -115.657000 0.00 0.00
- O12A TA 40.267900 -114.745400 0.00 0.00
- O13A TA 40.131400 -113.981400 0.00 0.00
- O15A TA 40.280700 -112.469400 0.00 0.00
- O16A TA 40.206700 -111.502400 0.00 0.00
- O17A TA 40.194800 -110.737800 0.00 0.00
- O18A TA 40.265500 -110.008200 0.00 0.00
- O19A TA 40.298200 -109.123600 0.00 0.00
- P01C TA 39.469000 -123.337500 0.00 0.00
- P05C TA 39.302800 -120.607700 0.00 0.00
- P06A TA 39.678500 -119.898300 0.00 0.00
- P07A TA 39.539900 -118.889300 0.00 0.00
- P08A TA 39.694600 -118.080000 0.00 0.00
- P09A TA 39.551600 -117.139500 0.00 0.00
- P10A TA 39.620200 -116.463900 0.00 0.00
- P11A TA 39.553000 -115.753600 0.00 0.00
- P12A TA 39.473100 -114.907500 0.00 0.00
- P13A TA 39.455000 -114.015600 0.00 0.00
- P14A TA 39.590600 -113.068700 0.00 0.00
- P15A TA 39.570800 -112.278600 0.00 0.00
- P16A TA 39.609200 -111.659500 0.00 0.00
- P17A TA 39.473000 -110.740000 0.00 0.00
- P18A TA 39.627900 -110.245700 0.00 0.00
- P19A TA 39.633300 -108.981200 0.00 0.00
- Q03C TA 38.633300 -122.014500 0.00 0.00
- Q04C TA 38.839000 -121.376900 0.00 0.00
- Q07A TA 38.938300 -118.807800 0.00 0.00
- Q08A TA 38.860600 -117.931600 0.00 0.00
- Q09A TA 38.834000 -117.181600 0.00 0.00
- Q10A TA 38.824700 -116.399900 0.00 0.00
- Q11A TA 38.845500 -115.654100 0.00 0.00
- Q12A TA 39.040000 -114.829900 0.00 0.00
- Q13A TA 38.955100 -114.020200 0.00 0.00
- Q14A TA 38.988200 -113.276900 0.00 0.00
- Q15A TA 38.999500 -112.379300 0.00 0.00
- Q16A TA 38.917600 -111.171600 0.00 0.00
- Q18A TA 39.102500 -110.132800 0.00 0.00
- Q19A TA 38.955300 -109.263000 0.00 0.00
- R04C TA 38.257100 -120.936300 0.00 0.00
- R05C TA 38.702800 -120.075500 0.00 0.00
- R06C TA 38.522600 -119.450900 0.00 0.00
- R07C TA 38.089000 -119.046900 0.00 0.00
- R08A TA 38.348900 -118.106400 0.00 0.00
- R09A TA 38.239700 -117.071800 0.00 0.00
- R10A TA 38.288600 -116.302100 0.00 0.00
- R11A TA 38.348900 -115.585400 0.00 0.00
- R12A TA 38.328100 -114.607600 0.00 0.00
- R13A TA 38.180200 -113.969100 0.00 0.00
- R14A TA 38.298700 -113.021300 0.00 0.00
- R15A TA 38.210800 -112.276700 0.00 0.00
- R16A TA 38.284300 -111.482700 0.00 0.00
- R17A TA 38.418700 -110.710900 0.00 0.00
- R18A TA 38.386200 -109.894200 0.00 0.00
- R19A TA 38.291600 -109.260700 0.00 0.00
- S04C TA 37.504900 -121.327800 0.00 0.00
- S05C TA 37.346000 -120.329600 0.00 0.00
- S06C TA 37.881700 -119.848800 0.00 0.00
- S08C TA 37.499300 -118.171100 0.00 0.00
- S09A TA 37.724300 -117.224600 0.00 0.00
- S10A TA 37.923000 -116.594800 0.00 0.00
- S11A TA 37.644400 -115.747200 0.00 0.00
- S12A TA 37.608300 -114.848900 0.00 0.00
- S13A TA 37.580800 -113.860400 0.00 0.00
- S14A TA 37.760100 -113.168400 0.00 0.00
- S15A TA 37.676100 -112.362700 0.00 0.00
- S16A TA 37.722200 -111.596100 0.00 0.00
- S17A TA 37.635400 -110.801800 0.00 0.00
- S18A TA 37.690300 -109.994800 0.00 0.00
- S19A TA 37.745900 -109.136600 0.00 0.00
- SUTB TA 39.229100 -121.786100 0.00 0.00
- T05C TA 36.895300 -120.673700 0.00 0.00
- T06C TA 37.006900 -119.708900 0.00 0.00
- T11A TA 37.240800 -115.220200 0.00 0.00
- T12A TA 36.725600 -114.714700 0.00 0.00
- T13A TA 37.019500 -113.907300 0.00 0.00
- T14A TA 37.061500 -113.084000 0.00 0.00
- T15A TA 37.018300 -112.382400 0.00 0.00
- T16A TA 36.983900 -111.506000 0.00 0.00
- T17A TA 36.996900 -110.804100 0.00 0.00
- T18A TA 37.136400 -109.874100 0.00 0.00
- T19A TA 36.829900 -109.024900 0.00 0.00
- U04C TA 36.362700 -120.782800 0.00 0.00
- U05C TA 36.335600 -120.120500 0.00 0.00
- U10A TA 36.419300 -116.329700 0.00 0.00
- U11A TA 36.423000 -115.383500 0.00 0.00
- U12A TA 36.432100 -114.538800 0.00 0.00
- U13A TA 36.415100 -113.965300 0.00 0.00
- U14A TA 36.418200 -113.180500 0.00 0.00
- U15A TA 36.428000 -112.291500 0.00 0.00
- U16A TA 36.142800 -111.129700 0.00 0.00
- U17A TA 36.599600 -110.662400 0.00 0.00
- U18A TA 36.419900 -109.869600 0.00 0.00
- V03C TA 36.021400 -121.235600 0.00 0.00
- V04C TA 35.636000 -120.869900 0.00 0.00
- V05C TA 35.866700 -119.902800 0.00 0.00
- V11A TA 35.838400 -115.430500 0.00 0.00
- V12A TA 35.726600 -114.851100 0.00 0.00
- V13A TA 35.852200 -113.984000 0.00 0.00
- V14A TA 35.633900 -113.105300 0.00 0.00
- V15A TA 35.819100 -112.173100 0.00 0.00
- V18A TA 35.711300 -109.932700 0.00 0.00
- V19A TA 35.714700 -109.045600 0.00 0.00
- W12A TA 35.301000 -114.870100 0.00 0.00
- W13A TA 35.099000 -113.885400 0.00 0.00
- W14A TA 35.212700 -113.083400 0.00 0.00
- W15A TA 35.178700 -112.266600 0.00 0.00
- W16A TA 35.095100 -111.531800 0.00 0.00
- W17A TA 35.078900 -110.712700 0.00 0.00
- W18A TA 35.118100 -109.735700 0.00 0.00
- W19A TA 35.111800 -109.387900 0.00 0.00
- X13A TA 34.593500 -113.830200 0.00 0.00
- X14A TA 34.469200 -112.890600 0.00 0.00
- X15A TA 34.486600 -112.236800 0.00 0.00
- X16A TA 34.417800 -111.441100 0.00 0.00
- X17A TA 34.337100 -110.805800 0.00 0.00
- X18A TA 34.529300 -109.950100 0.00 0.00
- X19A TA 34.427500 -109.290100 0.00 0.00
- Y12C TA 33.750300 -114.523800 0.00 0.00
- Y13A TA 33.814200 -113.828700 0.00 0.00
- Y14A TA 33.938300 -113.004800 0.00 0.00
- Y15A TA 33.953500 -112.333100 0.00 0.00
- Y16A TA 33.879800 -111.478300 0.00 0.00
- Y17A TA 33.695300 -110.844400 0.00 0.00
- Y18A TA 33.777600 -110.034100 0.00 0.00
- Y19A TA 33.957100 -109.254100 0.00 0.00
- Y22C TA 34.074100 -106.921100 0.00 0.00
- Z14A TA 33.362700 -112.945800 0.00 0.00
- Z15A TA 33.289300 -112.158100 0.00 0.00
- Z16A TA 33.341000 -111.427300 0.00 0.00
- Z17A TA 33.296900 -110.472300 0.00 0.00
- Z18A TA 33.085100 -110.036200 0.00 0.00
- Z19A TA 33.292100 -109.265700 0.00 0.00
- BOZ US 45.647000 -111.629600 0.00 0.00
- DUG US 40.195000 -112.813300 0.00 0.00
- EGMT US 48.024000 -109.754700 0.00 0.00
- ELK US 40.744840 -115.238780 0.00 0.00
- HLID US 43.562500 -114.413800 0.00 0.00
- MNV US 38.432800 -118.153110 0.00 0.00
- MVCO US 37.210300 -108.498600 0.00 0.00
- NLWA US 47.390000 -123.870000 0.00 0.00
- TPNV US 36.948800 -116.249500 0.00 0.00
- WUAZ US 35.516900 -111.373900 0.00 0.00
- WVOR US 42.433950 -118.636717 0.00 0.00
- SRU UU 39.110830 -110.523830 0.00 0.00
+A1 II 34.5 -118.00 0.0 0.0
+A2 II 34.5 -117.75 0.0 0.0
+A3 II 34.5 -117.50 0.0 0.0
+A4 II 34.5 -117.25 0.0 0.0
+A5 II 34.5 -117.00 0.0 0.0
+A6 II 34.5 -116.75 0.0 0.0
+A7 II 34.5 -116.50 0.0 0.0
+A8 II 34.5 -116.25 0.0 0.0
+A9 II 34.5 -116.00 0.0 0.0
+A10 II 34.5 -115.75 0.0 0.0
+A11 II 34.5 -115.50 0.0 0.0
+A12 II 34.5 -115.25 0.0 0.0
+A13 II 34.5 -115.00 0.0 0.0
Modified: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/Makefile.in
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/Makefile.in 2010-05-28 09:53:29 UTC (rev 16821)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/Makefile.in 2010-05-28 15:11:05 UTC (rev 16822)
@@ -138,6 +138,7 @@
$O/assemble_MPI_vector.o \
$O/read_arrays_solver.o \
$O/specfem3D.o \
+ $O/noise_tomography.o \
$(EMPTY_MACRO)
# objects toggled between the parallel and serial version
@@ -283,6 +284,9 @@
$O/parallel.o: constants.h parallel.f90
${MPIFCCOMPILE_CHECK} -c -o $O/parallel.o parallel.f90
+$O/noise_tomography.o: constants.h OUTPUT_FILES/values_from_mesher.h noise_tomography.f90
+ ${MPIFCCOMPILE_NO_CHECK} -c -o $O/noise_tomography.o ${FCFLAGS_f90} noise_tomography.f90
+
###
### serial compilation without optimization
###
Added: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/NOISE_clean_create
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/NOISE_clean_create (rev 0)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/NOISE_clean_create 2010-05-28 15:11:05 UTC (rev 16822)
@@ -0,0 +1,7 @@
+#!/bin/bash
+ulimit -s unlimited
+
+LOCAL_PATH=$1
+
+rm -rf $LOCAL_PATH
+mkdir -p $LOCAL_PATH
Property changes on: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/NOISE_clean_create
___________________________________________________________________
Name: svn:executable
+ *
Added: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/Z_submit
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/Z_submit (rev 0)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/Z_submit 2010-05-28 15:11:05 UTC (rev 16822)
@@ -0,0 +1,91 @@
+#!/bin/bash
+#PBS -S /bin/bash
+#PBS -N NOISE_BASIN_12
+#PBS -q tromp
+#PBS -o zzz_job_info/job.o_all
+#PBS -e zzz_job_info/job.e_all
+
+#PBS -l nodes=21:ppn=8,walltime=100:00:00
+
+
+cd $PBS_O_WORKDIR
+
+######################################## pre-simulation ###########################################
+## gather simulation infomation
+NPROC_XI=`grep NPROC_XI DATA/Par_file | cut -d = -f 2 `
+NPROC_ETA=`grep NPROC_ETA DATA/Par_file | cut -d = -f 2`
+numnodes=$(( $NPROC_XI * $NPROC_ETA ))
+cat $PBS_NODEFILE > zzz_job_info/compute_nodes_all
+echo "$PBS_JOBID" > zzz_job_info/jobid_all
+
+### clean local nodes
+d=`date`
+LOCAL_PATH_ALL=`grep LOCAL_PATH DATA/Par_file | cut -d = -f 2 | sed 's/ //g'`
+echo "Start cleaning local disk on nodes, $d"
+#for (( ievent=1;ievent<=$numnodes;ievent=ievent+1 ))
+#do
+# let line=$ievent
+# machine_all=$(awk 'NR=='$line'{print}' zzz_job_info/compute_nodes_all)
+# mpiexec -np 1 -host $machine_all ./NOISE_clean_create $LOCAL_PATH_ALL < /dev/null &
+#done
+#wait
+for (( ievent=1;ievent<=$numnodes;ievent=ievent+8 ))
+do
+ let line=$ievent
+ machine_all=$(awk 'NR=='$line'{print}' zzz_job_info/compute_nodes_all)
+ mpiexec -np 1 -host $machine_all ./NOISE_clean_create $LOCAL_PATH_ALL < /dev/null &
+done
+wait
+d=`date`
+echo "Finish cleaning local disk on nodes, $d"
+
+## clean OUTPUT_FILES
+rm -rf OUTPUT_FILES
+mkdir OUTPUT_FILES
+
+#### backup Fortran source files
+##rm -rf ZZZ_src
+##mkdir -p ZZZ_src
+##cp -p *.f90 ZZZ_src/
+##cp -p *.c ZZZ_src/
+##cp -p *.in ZZZ_src/
+##cp -p *.h ZZZ_src/
+######################################## simulation ###############################################
+################################################################### step 1: ensemble forward source
+sleep 2
+mpiexec -np $numnodes -hostfile zzz_job_info/compute_nodes_all $PWD/xmeshfem3D_NOISE1
+sleep 2
+mpiexec -np $numnodes -hostfile zzz_job_info/compute_nodes_all $PWD/xspecfem3D_NOISE1
+# backup step 1 results
+rm -rf ZZZ_1
+mkdir -p ZZZ_1
+mv OUTPUT_FILES/* ZZZ_1/
+cp DATA/Par_file ZZZ_1/
+cp DATA/CMTSOLUTION ZZZ_1/
+cp DATA/STATIONS ZZZ_1/
+################################################################ step 2: ensemble forward wavefield
+cp DATA/Par_file_NOISE_2_attenuation DATA/Par_file
+sleep 2
+mpiexec -np $numnodes -hostfile zzz_job_info/compute_nodes_all $PWD/xmeshfem3D_NOISE2
+sleep 2
+mpiexec -np $numnodes -hostfile zzz_job_info/compute_nodes_all $PWD/xspecfem3D_NOISE2
+# backup step 2 results
+rm -rf ZZZ_2
+mkdir -p ZZZ_2
+mv OUTPUT_FILES/* ZZZ_2/
+cp DATA/Par_file ZZZ_2/
+cp DATA/CMTSOLUTION ZZZ_2/
+cp DATA/STATIONS ZZZ_2/
+
+
+d=`date`
+echo "Start cleaning local disk on nodes, $d"
+for (( ievent=1;ievent<=$numnodes;ievent=ievent+8 ))
+do
+ let line=$ievent
+ machine_all=$(awk 'NR=='$line'{print}' zzz_job_info/compute_nodes_all)
+ mpiexec -np 1 -host $machine_all ./NOISE_clean_create $LOCAL_PATH_ALL < /dev/null &
+done
+wait
+d=`date`
+echo "Finish cleaning local disk on nodes, $d"
Property changes on: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/Z_submit
___________________________________________________________________
Name: svn:executable
+ *
Modified: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/constants.h.in
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/constants.h.in 2010-05-28 09:53:29 UTC (rev 16821)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/constants.h.in 2010-05-28 15:11:05 UTC (rev 16822)
@@ -61,6 +61,8 @@
! input, output and main MPI I/O files
integer, parameter :: ISTANDARD_OUTPUT = 6
integer, parameter :: IIN = 40,IOUT = 41
+ integer, parameter :: IIN_NOISE = 43,IOUT_NOISE = 44
+ integer, parameter :: NOISE_TOMOGRAPHY = 1
! uncomment this to write messages to a text file
integer, parameter :: IMAIN = 42
! uncomment this to write messages to the screen
Modified: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/create_header_file.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/create_header_file.f90 2010-05-28 09:53:29 UTC (rev 16821)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/create_header_file.f90 2010-05-28 15:11:05 UTC (rev 16822)
@@ -98,8 +98,12 @@
print *
print *,'on NEC SX, make sure "loopcnt=" parameter'
print *,'in Makefile is greater than max vector length = ',NGLOB_AB
-
print *
+ print *, '********************************'
+ print *,'number of time steps = ',NSTEP
+ print *,'time-stepping of the solver will be: ',DT
+ print *, '********************************'
+ print *
print *,'done'
print *
Modified: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/create_movie_shakemap_AVS_DX_GMT.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/create_movie_shakemap_AVS_DX_GMT.f90 2010-05-28 09:53:29 UTC (rev 16821)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/create_movie_shakemap_AVS_DX_GMT.f90 2010-05-28 15:11:05 UTC (rev 16822)
@@ -40,7 +40,7 @@
real(kind=CUSTOM_REAL), parameter :: THRESHOLD = 1._CUSTOM_REAL / 100._CUSTOM_REAL
! coefficient of power law used for non linear scaling
- logical, parameter :: NONLINEAR_SCALING = .true.
+ logical, parameter :: NONLINEAR_SCALING = .false.
real(kind=CUSTOM_REAL), parameter :: POWER_SCALING = 0.25_CUSTOM_REAL
integer it,it1,it2,ivalue,nspectot_AVS_max,ispec
Modified: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/flags.guess
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/flags.guess 2010-05-28 09:53:29 UTC (rev 16821)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/flags.guess 2010-05-28 15:11:05 UTC (rev 16822)
@@ -26,7 +26,7 @@
# Intel ifort Fortran90 for Linux
#
if test x"$FLAGS_CHECK" = x; then
- FLAGS_CHECK="-O3 -vec-report0 -e03 -std03 -implicitnone -warn truncated_source -warn argument_checking -warn unused -warn declarations -warn alignments -warn ignore_loc -warn usage -check nobounds -align sequence -assume byterecl -fpe0 -ftz -traceback -ftrapuv" # -mcmodel=medium
+ FLAGS_CHECK="-O3 -vec-report0 -implicitnone -warn truncated_source -warn argument_checking -warn unused -warn declarations -warn alignments -warn ignore_loc -warn usage -check nobounds -align sequence -assume byterecl -fpe0 -ftz -traceback -ftrapuv" # -mcmodel=medium
fi
if test x"$FLAGS_NO_CHECK" = x; then
# standard options (leave option -ftz, which is *critical* for performance)
@@ -35,7 +35,7 @@
#
# use -xP or -xS with Intel version 10, -xSSE4.2 with Intel version 11 (for Intel Nehalem processors; otherwise -xS is fine)
#
- FLAGS_NO_CHECK="-O3 -xS -vec-report0 -e03 -std03 -implicitnone -warn truncated_source -warn argument_checking -warn unused -warn declarations -warn alignments -warn ignore_loc -warn usage -check nobounds -align sequence -assume byterecl -fpe3 -ftz" # -mcmodel=medium
+ FLAGS_NO_CHECK="-O3 -xS -vec-report0 -implicitnone -warn truncated_source -warn argument_checking -warn unused -warn declarations -warn alignments -warn ignore_loc -warn usage -check nobounds -align sequence -assume byterecl -fpe3 -ftz" # -mcmodel=medium
fi
#MPI_LIBS = -Vaxlib
;;
Modified: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/read_parameter_file.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/read_parameter_file.f90 2010-05-28 09:53:29 UTC (rev 16821)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/read_parameter_file.f90 2010-05-28 15:11:05 UTC (rev 16822)
@@ -221,6 +221,7 @@
! compute total number of time steps, rounded to next multiple of 100
NSTEP = 100 * (int(RECORD_LENGTH_IN_SECONDS / (100.d0*DT)) + 1)
+ if ( NOISE_TOMOGRAPHY /= 0 ) NSTEP = 2*NSTEP-1 ! time steps needs to be doubled, due to +/- branches
! compute the total number of sources in the CMTSOLUTION file
! there are NLINES_PER_CMTSOLUTION_SOURCE lines per source in that file
Modified: seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/specfem3D.f90
===================================================================
--- seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/specfem3D.f90 2010-05-28 09:53:29 UTC (rev 16821)
+++ seismo/3D/SPECFEM3D_SESAME/tags/v1.4.4_last_BASIN/specfem3D.f90 2010-05-28 15:11:05 UTC (rev 16822)
@@ -345,7 +345,7 @@
real(kind=CUSTOM_REAL) b_additional_term,b_force_normal_comp, kappa_k, mu_k
real(kind=CUSTOM_REAL), dimension(NDIM,NGLOB_ADJOINT) :: b_displ, b_veloc, b_accel
real(kind=CUSTOM_REAL), dimension(NGLLX,NGLLY,NGLLZ,NSPEC_ADJOINT) :: rho_kl, mu_kl, kappa_kl, &
- rhop_kl, beta_kl, alpha_kl
+ rhop_kl, beta_kl, alpha_kl,Sigma_kl_crust_mantle
real(kind=CUSTOM_REAL) dsxx,dsxy,dsxz,dsyy,dsyz,dszz
real(kind=CUSTOM_REAL) b_duxdxl,b_duxdyl,b_duxdzl,b_duydxl,b_duydyl,b_duydzl,b_duzdxl,b_duzdyl,b_duzdzl
real(kind=CUSTOM_REAL) b_duxdxl_plus_duydyl,b_duxdxl_plus_duzdzl,b_duydyl_plus_duzdzl
@@ -553,6 +553,15 @@
real(kind=CUSTOM_REAL), dimension(NGLLX,NGLLY,NGLLZ) :: dvxdxl,dvxdyl,dvxdzl,dvydxl,dvydyl,dvydzl,dvzdxl,dvzdyl,dvzdzl
real(kind=CUSTOM_REAL), dimension(:,:,:,:),allocatable:: div, curl_x, curl_y, curl_z
+!<YANGL
+! NOISE_TOMOGRAPHY
+ real(kind=CUSTOM_REAL), dimension(:,:,:,:,:), allocatable :: noise_sourcearray
+ integer :: irec_master_noise
+ real(kind=CUSTOM_REAL), dimension(:), allocatable :: &
+ normal_x_noise,normal_y_noise,normal_z_noise, mask_noise
+!>YANGL
+
+
! ************** PROGRAM STARTS HERE **************
! sizeprocs returns number of processes started
@@ -1383,12 +1392,12 @@
rho_kl(:,:,:,:) = 0._CUSTOM_REAL
mu_kl(:,:,:,:) = 0._CUSTOM_REAL
kappa_kl(:,:,:,:) = 0._CUSTOM_REAL
-
+ if (NOISE_TOMOGRAPHY == 3) Sigma_kl_crust_mantle(:,:,:,:) = 0._CUSTOM_REAL
endif
! allocate files to save movies and shaking map
- if(MOVIE_SURFACE .or. CREATE_SHAKEMAP) then
- if (USE_HIGHRES_FOR_MOVIES) then
+ if(MOVIE_SURFACE .or. CREATE_SHAKEMAP .or. NOISE_TOMOGRAPHY /= 0) then ! for noise tomography, store_val_x/y/z/ux/uy/uz needed for 'surface movie'
+ if (USE_HIGHRES_FOR_MOVIES .or. NOISE_TOMOGRAPHY /= 0) then ! for noise tomography, must NOT be coarse
nmovie_points = NGLLX * NGLLY * NSPEC2D_TOP
else
nmovie_points = NGNOD2D * NSPEC2D_TOP
@@ -1410,14 +1419,15 @@
allocate(store_val_norm_displ(nmovie_points))
allocate(store_val_norm_veloc(nmovie_points))
allocate(store_val_norm_accel(nmovie_points))
+ if (MOVIE_SURFACE .or. CREATE_SHAKEMAP) then ! those arrays are not neccessary for noise tomography, so only allocate them in MOVIE_SURFACE case
+ allocate(store_val_x_all(nmovie_points,0:NPROC-1))
+ allocate(store_val_y_all(nmovie_points,0:NPROC-1))
+ allocate(store_val_z_all(nmovie_points,0:NPROC-1))
+ allocate(store_val_ux_all(nmovie_points,0:NPROC-1))
+ allocate(store_val_uy_all(nmovie_points,0:NPROC-1))
+ allocate(store_val_uz_all(nmovie_points,0:NPROC-1))
+ endif
- allocate(store_val_x_all(nmovie_points,0:NPROC-1))
- allocate(store_val_y_all(nmovie_points,0:NPROC-1))
- allocate(store_val_z_all(nmovie_points,0:NPROC-1))
- allocate(store_val_ux_all(nmovie_points,0:NPROC-1))
- allocate(store_val_uy_all(nmovie_points,0:NPROC-1))
- allocate(store_val_uz_all(nmovie_points,0:NPROC-1))
-
! to compute max of norm for shaking map
store_val_norm_displ(:) = -1.
store_val_norm_veloc(:) = -1.
@@ -1513,6 +1523,35 @@
endif
close(27)
+!<YANGL
+ ! NOISE TOMOGRAPHY
+ if ( NOISE_TOMOGRAPHY /= 0 ) then
+ allocate(noise_sourcearray(NDIM,NGLLX,NGLLY,NGLLZ,NSTEP))
+ allocate(normal_x_noise(nmovie_points))
+ allocate(normal_y_noise(nmovie_points))
+ allocate(normal_z_noise(nmovie_points))
+ allocate(mask_noise(nmovie_points))
+ noise_sourcearray(:,:,:,:,:) = 0._CUSTOM_REAL
+ normal_x_noise(:) = 0._CUSTOM_REAL
+ normal_y_noise(:) = 0._CUSTOM_REAL
+ normal_z_noise(:) = 0._CUSTOM_REAL
+ mask_noise(:) = 0._CUSTOM_REAL
+
+ call read_parameters_noise(myrank,nrec,NSTEP,nmovie_points, &
+ islice_selected_rec,xi_receiver,eta_receiver,gamma_receiver,nu, &
+ noise_sourcearray,xigll,yigll,zigll,NSPEC2D_TOP_VAL, &
+ 1, ibool, ibelm_top, &
+ xstore,ystore,zstore, &
+ irec_master_noise,normal_x_noise,normal_y_noise,normal_z_noise,mask_noise)
+
+ if (myrank == 0) &
+ call check_parameters_noise(myrank,NOISE_TOMOGRAPHY,SIMULATION_TYPE,SAVE_FORWARD, &
+ 1, 1, .false., &
+ .false., .false., USE_HIGHRES_FOR_MOVIES)
+ endif
+!>YANGL
+
+
!
! s t a r t t i m e i t e r a t i o n s
!
@@ -3030,6 +3069,42 @@
endif
+!<YANGL
+ ! NOISE_TOMOGRAPHY
+ if ( NOISE_TOMOGRAPHY == 1 ) then
+ ! the first step of noise tomography is to use |S(\omega)|^2 as a point force source at one of the receivers.
+ ! hence, instead of a moment tensor 'sourcearrays', a 'noise_sourcearray' for a point force is needed.
+ ! furthermore, the CMTSOLUTION needs to be zero, i.e., no earthquakes.
+ ! now this must be manually set in DATA/CMTSOLUTION, by USERS.
+ call add_source_master_rec_noise(myrank,nrec, &
+ NSTEP,accel,noise_sourcearray, &
+ ibool,islice_selected_rec,ispec_selected_rec, &
+ it,irec_master_noise)
+ elseif ( NOISE_TOMOGRAPHY == 2 ) then
+ ! second step of noise tomography, i.e., read the surface movie saved at every timestep
+ ! use the movie to drive the unsemble forward wavefield
+ call noise_read_add_surface_movie(myrank,nmovie_points,accel, &
+ normal_x_noise,normal_y_noise,normal_z_noise,mask_noise, &
+ store_val_ux,store_val_uy,store_val_uz, &
+ ibelm_top,ibool,NSPEC2D_TOP, &
+ 1,NSTEP-it+1,LOCAL_PATH,jacobian2D_top,wgllwgll_xy)
+ ! be careful, since ensemble forward sources are reversals of generating wavefield "eta"
+ ! hence the "NSTEP-it+1", i.e., start reading from the last timestep
+ ! note the ensemble forward sources are generally distributed on the surface of the earth
+ ! that's to say, the ensemble forward source is kind of a surface force density, not a body force density
+ ! therefore, we must add it here, before applying the inverse of mass matrix
+ elseif ( NOISE_TOMOGRAPHY == 3 ) then
+ ! third step of noise tomography, i.e., read the surface movie saved at every timestep
+ ! use the movie to reconstruct the ensemble forward wavefield
+ ! the ensemble adjoint wavefield is done as usual
+ ! note instead of "NSTEP-it+1", now we us "it", since reconstruction is a reversal of reversal
+ call noise_read_add_surface_movie(myrank,nmovie_points,b_accel, &
+ normal_x_noise,normal_y_noise,normal_z_noise,mask_noise, &
+ store_val_ux,store_val_uy,store_val_uz, &
+ ibelm_top,ibool,NSPEC2D_TOP, &
+ 1,it,LOCAL_PATH,jacobian2D_top,wgllwgll_xy)
+ endif
+!>YANGL
! assemble all the contributions between slices using MPI
call assemble_MPI_vector(accel,iproc_xi,iproc_eta,addressing, &
@@ -3333,6 +3408,15 @@
enddo
enddo
+!<YANGL
+ ! NOISE TOMOGRAPHY --- source strength kernel
+ if (NOISE_TOMOGRAPHY == 3) &
+ call compute_kernels_strength_noise(myrank,ibool, &
+ Sigma_kl_crust_mantle,displ,deltat,it, &
+ nmovie_points,normal_x_noise,normal_y_noise,normal_z_noise, &
+ NSPEC2D_TOP,ibelm_top,LOCAL_PATH)
+!>YANGL
+
if (SAVE_MOHO_MESH) then
do ispec2D = 1, nspec2D_moho
ispec_top = ibelm_moho_top(ispec2D)
@@ -3362,6 +3446,20 @@
endif
endif
+!<YANGL
+ ! first step of noise tomography, i.e., save a surface movie at every time step
+ ! modified from the subroutine 'write_movie_surface'
+ if ( NOISE_TOMOGRAPHY == 1 ) then
+ call noise_save_surface_movie(myrank,nmovie_points,displ, &
+ xstore,ystore,zstore, &
+ store_val_x,store_val_y,store_val_z, &
+ store_val_ux,store_val_uy,store_val_uz, &
+ ibelm_top,ibool, &
+ NSPEC2D_TOP, &
+ 1,it,LOCAL_PATH)
+ endif
+!>YANGL
+
! save MOVIE on the SURFACE
if(MOVIE_SURFACE .and. mod(it,NTSTEP_BETWEEN_FRAMES) == 0) then
@@ -3692,6 +3790,13 @@
write(27) moho_kl
close(27)
endif
+!<YANGL
+ ! noise strength kernel
+ if (NOISE_TOMOGRAPHY == 3) then
+ call save_kernels_strength_noise(myrank,LOCAL_PATH, &
+ Sigma_kl_crust_mantle,1.0,1.0)
+ endif
+!>YANGL
endif
More information about the CIG-COMMITS
mailing list