Increased mm array as per RWG's instructions to support larger ruptur…

…es. 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

@@ -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
 

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)