From 83d53b959614c89399e67099ca2ba3baef05b655 Mon Sep 17 00:00:00 2001 From: Fabio Silva Date: Wed, 5 Jun 2024 00:01:54 -0700 Subject: [PATCH] Increased mm array as per RWG's instructions to support larger ruptures. Also added RWG's updates to hb_high_v6.1.1_sub.f --- bbp/src/gp/StochSim/V6.0/hb_high_v6.1.1_sub.f | 14 +++++++++++--- bbp/src/gp/StochSim/V6.0/params.h | 2 +- 2 files changed, 12 insertions(+), 4 deletions(-) diff --git a/bbp/src/gp/StochSim/V6.0/hb_high_v6.1.1_sub.f b/bbp/src/gp/StochSim/V6.0/hb_high_v6.1.1_sub.f index 483333c3..a0cc01d5 100644 --- a/bbp/src/gp/StochSim/V6.0/hb_high_v6.1.1_sub.f +++ b/bbp/src/gp/StochSim/V6.0/hb_high_v6.1.1_sub.f @@ -169,7 +169,7 @@ subroutine hb_high(cap, stlon, stlat, velfile, outname, duration, dt, ds, nevnt, dimension astop(lv),dtop(lv),shyp(lv),dhyp(lv),nx(lv),nw(lv),dx(lv),dw(lv) - DIMENSION stdd(mmv,3),dfr(mm),fgrand(mmv),fgrands(mm) + DIMENSION stdd(mm,3),dfr(mm),fgrand(mmv),fgrands(mm) complex cs(mm,3),sgc dimension fn(nlaymax),an(nlaymax) @@ -296,13 +296,21 @@ subroutine hb_high(cap, stlon, stlat, velfile, outname, duration, dt, ds, nevnt, fcfac = fcfac_default fcfac = 0.0 +c RWG 20240220: initialization of zhyp_max, nstot, farea_in added + zhyp_max = 0.0 + nstot = 0 + farea_in = 0.0 + do 132 iv=1,nevnt nstot = nx(iv)*nw(iv) + nstot farea_in = nx(iv)*dx(iv)*nw(iv)*dw(iv) + farea_in astop(iv) = 0.5*nx(iv)*dx(iv) -c RWG 201241219: '*' should have been '/' in below, fixed 21041219 - zhyp = dtop(iv) + dhyp(iv)/sin(dipq(iv)*pu) +c##- c RWG 201241219: '*' should have been '/' in below, fixed 21041219 +c##- zhyp = dtop(iv) + dhyp(iv)/sin(dipq(iv)*pu) +c RWG 20240220: '/' above should have been '*', below is correct +c I have no idea why above was done + zhyp = dtop(iv) + dhyp(iv)*sin(dipq(iv)*pu) if(zhyp.gt.zhyp_max) zhyp_max = zhyp 132 continue diff --git a/bbp/src/gp/StochSim/V6.0/params.h b/bbp/src/gp/StochSim/V6.0/params.h index fdc0286a..276d94b6 100644 --- a/bbp/src/gp/StochSim/V6.0/params.h +++ b/bbp/src/gp/StochSim/V6.0/params.h @@ -1,5 +1,5 @@ C old way with large 'lv', this is nevent and needed this for ShakeOut & Hayward sims C parameter (mm=65392,mmv=60000,mst=6,lv=1000) C more general way with small 'lv' - parameter (mm=65392,mmv=60000,mst=6,lv=10) + parameter (mm=131072,mmv=60000,mst=6,lv=10) parameter (nq=2000,np=500,nlaymax=500)