C- program hazFXnga13l.f;; Use with NGA relations, or others.c March 16, 2015: remove a repeated line of code in the GR-without-uncert loop.c This repeated line caused CEUS GMM calculations to be skipped for GR-distr.csource magnitudes. Discovery by Morgan Moschetti, 3/16/2015.c Feb 12, 2015: modify code to handle PGV in CEUS (for two Atkinson relations only)cPGV units are cm/s. There is a clamp at 400 cm/s for these relations.c Oct 30 2014: repair /gail3/ common to have freq(14) dimension.c Mar 14, 2014: A08', AB06', and Pez11 all have added 1.5-s table entries. Tablesc with additional 1.5s end in.rev instead of.datc jan 27, 2015: getSomer: add 1.5s coeffs. Repair 1.5s Jan 27 2015(CEUS relation used for finite src)c Jan 24, 2014: getAB06: Add coeffs for 1.5s jan 24 2014 SH. See AB06.1p5s.f for the details.c Jan 24, 2014: getCampCEUS: include interpo.
This is the software package (zip file) for SFEdit2 and supports R-WearStudio and MetazaStudio software only. 1.83 MB VersaWorks 6.5.2 Offline Updater: Use this updater if you are running VersaWorks version. Click the 'Download Software' button below for your version of the free PowerScan software. 'Save' the download to your computer. Exit all programs before installing the software. Locate the.zip folder that you downloaded. ScanPro 2000 latest software ver3.345a2019-07-22c Download Software: ScanPro 1100 latest software ver6.345a.
Coeffs for 1.5s spectral period see 'campCEUS.1p5s.f' for detailsc Jan 24, 2014: getToro, getSilva: include interpo. Coeffs for 1.5s spectral period see 'silva.1p5s.f' for detailsc 1/22/2014: correct aleatory sigma in Idriss NGA-W(2).c: Use same median clamps in amean11 as are used in other CEUS GMMs. PPowersc noticed some discrepencies in earlier versions of this code.c: Update ASK2013 a1 and vlin for some short periods. From Sanaz R emailc OCT 1, 2013: use Mcap of 7.5 on s.d.
Computation in Idriss2013.c 9/17/2013In ASK13, Ry0 is always -1. (previously some cases Ry0=dmin(5)=-1)c 8/28/2013CB13: 'Save' philny and always calculate philny(22) on first time thruc 8/27/2013CB13: Update c0 vector for PGA and several short-period SA Bozorgnia emailc 8/27/2013changed z1km to z1ref in call to ASKc 8/21/2013CB13: Improve Zhyp further for downdip ruptures.c 8/20/2013 Inline CB13 coeffs.
Begin standardizing Zhyp in CB13.c Correction Aug 15 2013: initialize PI in CB13 subroutine (previously wasn't).c Also, calculate A1100 first time through CB13cc 8/13/2013: Use the Rx0 version of ASK13. Set z1rock = -1 for the hardrock calls of ASKc 8/05/2013: Minor repair in vicinity of line 1095 (some compiler subrange errors in prior vers)c 8/02/2013: The correlation coeff vector rho was updated in CB13. It is nowc stored as a data vector rather than being read in.c: ASK13: Use estimated Vs30 setting if index is 36. Use measured Vs30 if index is -36.c: Incorporate CY2013 vers of July 29. Use estimated Vs30 if index is 35.c Use measured Vs30 if index is -35 (new 8/01/2013).c: Incorporate ASKNGA2013V11.f which is the July 24 version of A S &Kamai. Uses Ry0 distance.ccc July 22, 2013: Using latest BSSA coeffs. Clin and Vclin for soil site computations.c These postdate the Eq Spectra coeffs.c July 16, 2013: Increase the number of clustered-event groups from 5 to 8.
To handlec NMSZ SSC clustered-event models. The M-uncertainty density functionc read in in relwt has to be the same for all eqs in a group. This is ac fundamental coding limitation. It is not checked internally.c 6/27/2013: For the Pezeshk 2012 relation, force rkm to be max(1., rkm). Veryc close hypocenters will cause Pezeshk GMPE to issue NaN median value.c For the AB06' GMPE(index 26), force rkm to be max(1.8,rkm)c 6/11/2013: correct z1 units when calling bssa2013drv.
Use Chiou Z1cal.c 5/21/2013: correct c11 term in gksa13v2c 5/20/2013: Add May 18 BSSA.c 5/17/2013: correct some e5 coeffs in bssa2013drv. Make cDPP = 0 in CY2013. (no directivity)c 5/16/2013: Update AS13 to May version. Norm A recommends using this instead of AS08 as of May 15.c 5/15/2013: update CY NGA-W to May version. 5/16: correct small typo from Chiou emailc 5/14/2013: Add anelastic attenuation in CB13, for R80 km.
Bozorgnia email May 2013.c 4/16/2013: bssa corrections completed. All 107 periods workc 4/12/2013: bssa2013 updated.
4/15: BSSA corrected some.c 4/12/2013: CY2013 updated. Corrected a line about HWtaper3 in AS-2013 routine (their only update).c 4/11/2013: Update Idriss NGA-W GMPE coeffs to latest available (emailed 4/10)c 4/09/2013: revise CB coeff file. The GMPE has not changed from feb, only the coeffs.c 4/09/2013: GK13 routine modified to version 2. Index is still 38. Basin and Range Q reset to 205cc 4/05/2013: correct the ia limit in wus02 ceus02 and some other arrays.
IAMX is the limit, now set to 10cc 4/04/2013: increase iamax, number of atten models, to 10 (per spectral period)cc 4/04/2013: Reduce basin effect in GK13 if Vs30600 m/s. Basin bottom is depth where Vs1.5 km/s. Cant be 0 for BCcsite condition.
Update sigmaaleatory model in gk13.c 3/29/2013: include GK13 model which replaces GK12. Index 38.c 3/19/2013: increase nfltmx to 1000.c 3/14/2013: revise AS08 to use our Rx. 3/15:Use vs30class=1 (measured).
Index 16C 3/11/2013: working on BSSA13. Simplified the coding some. Perhaps more needed (SH).c 3/06/2013: working on CY2013. Dirctivity, cDPP, is turned off for initial exercise.c A long document discusses some suggested ways to implement.c 3/01/2013: update A&S GMPE. 3/29/2013: The A&S subr. Is expecting z1 to come in with units km/s not m/s.cc 2/25/2013: Update CB-NGA2 GMPE Index is 34c 2/19/2013: Add option to output inelastic spectral displacement, SDi.c This option is invoked if the second argument is 'dy' or 'dY' and the third is thec value of dy to use (units: cm).
Expected use, with WUS crustal-src GMPEs.c 1/07/2013: can use a fixed aleatory sigma for all NGA-W gmpes.c To invoke this option:c hazFXnga13l.test input.file fixsigma 0.6c That is make arg(2) 'fixsigma' and make arg(3) the value you want (natural log sigma)c: update the CB2012 input file to 'v2'. Update affects the tau part of sigmaaleatory.c: add and repair the AS2012 GMPE for NGA-W. Use index 36. Implementation may have the Ry0 term.c To get the Ry0 version, set icode(ip,ia)=1 in input line.
July 31, 2013: also use ASK2013 if index is -36c: add the latest Graizer Kalkan gmpe. Includes basin effect. See dgkbasin. Index 38cc: add the BSSA 2012 GMPE: index 33c: CB12 with v3 update. Previously CB12 had v2 update.
(removed 2/25/2013)c: Idriss 2012 relation added index 37. Updated to Idriss 2013c: the Grazier GMPE is included for use in tectonically active areas. Index 39.c: the CY2012 GMPE is included for use in tectonically active areas.
Use thec sigma model which assumes Vs30 is known. From Chiou email c 3/21/2011: add 3 2011 GMPES for CENA from G. Atkinson email. Indexes 25, 26, 27.c These 2011 CENA GMPEs are set up for hard rock and BC(firm rock) only. Don't use them for the soil site condition.c: use Pezeshk11 sigma from article.
Not Gail's 'stopgap'.c: repair indexc problem when accessing sigPez11.c: Add Silva single-corner, m-dependent stress, index 28. BC or Ac FOr variable Vs30 do not use this code. Instead use hazFXnga7vs.f.c from hazFXnga7x.f. This version has event-clustering option. If this option isc selected, all events in input file must be part of the cluster model. Initialc work on this was for NMSZ scenarios. These scenarios dictate several limitsc on the available model.
Do not mix unclustered sources (faults) with clusteredc in any given input file,c but you can use this code on either clustered or independent events. Clustered eventsc are set up for Characteristic without uncertainty. Models need more work.c- Solaris compile: cc -c iosubs128.cc f95 -o hazFXnga13l hazFXnga13l.f iosubs128.o -fast -e -ftrap=%nonec -fast has been known to produce segmentation errors whereas -O is not known to produce them all else equalc -fast definitely fails (seg error) when the clustered source is invoked.c for debug:c f95 -o hazFXnga13l hazFXnga13l.f iosubs128.o -C -e -ftrap=%nonec for typical runs:c f95 -o hazFXnga13l hazFXnga13l.f iosubs128.o -O -e -ftrap=%nonec -e to extend linelength. Several statements go past column 72. Flag -fast no longer works due to somec excess stuff in the BSSA initialization. Best might be -Oc The ftrap flag is required to prevent crashing when NaNs such as log(0) arec attempted. Code has been modified to avoid some of these where known.c for detailed source listing -Xlist is the argument.
Comments are closed.
|
AuthorWrite something about yourself. No need to be fancy, just an overview. Archives
March 2023
Categories |