22c -------------------------------------------------------------
33c
44 subroutine stepgrid (q ,fm ,fp ,gm ,gp ,mitot ,mjtot ,mbc ,dt ,dtnew ,dx ,dy ,
5- & nvar ,xlow ,ylow ,time ,mptr ,maux ,aux , actualstep )
5+ & nvar ,xlow ,ylow ,time ,mptr ,maux ,aux )
66c
77c
88c ::::::::::::::::::: STEPGRID ::::::::::::::::::::::::::::::::::::
@@ -28,7 +28,7 @@ subroutine stepgrid(q,fm,fp,gm,gp,mitot,mjtot,mbc,dt,dtnew,dx,dy,
2828
2929 use geoclaw_module
3030 use amr_module
31- use fgout_module, only: FGOUT_num_grids, FGOUT_fgrids,
31+ use fgout_module, only: FGOUT_num_grids, FGOUT_fgrids,
3232 & FGOUT_tcfmax, fgout_interp, fgout_grid,
3333 & FGOUT_ttol
3434 implicit double precision (a- h,o- z)
@@ -47,14 +47,13 @@ subroutine stepgrid(q,fm,fp,gm,gp,mitot,mjtot,mbc,dt,dtnew,dx,dy,
4747 logical :: debug = .false.
4848 logical :: dump = .false.
4949 !logical :: dump = .true.
50- logical , intent (in ) :: actualstep
5150 type(fgout_grid), pointer :: fgout
5251 !logical , allocatable :: fgout_interp_needed(:)
5352 logical :: fgout_interp_needed(FGOUT_num_grids)
5453 real (kind= 8 ) :: fgout_tnext
5554c
5655#ifdef WHERE_AM_I
57- write (* ,* ) " starting stepgrid grid " ,mptr
56+ write (* ,* ) " starting stepgrid grid " ,mptr
5857#endif
5958 FGOUT_tcfmax = - rinfinity
6059 level = node(nestlevel,mptr)
@@ -135,7 +134,7 @@ subroutine stepgrid(q,fm,fp,gm,gp,mitot,mjtot,mbc,dt,dtnew,dx,dy,
135134 fgout_interp_needed(ng) = .false.
136135 else
137136 fgout_tnext = fgout%output_times(fgout%next_output_index)
138- fgout_interp_needed(ng) =
137+ fgout_interp_needed(ng) =
139138 & ((fgout%x_low < xlowmbc + mx * dx) .and.
140139 & (fgout%x_hi > xlowmbc) .and.
141140 & (fgout%y_low < ylowmbc + my * dy) .and.
@@ -153,14 +152,14 @@ subroutine stepgrid(q,fm,fp,gm,gp,mitot,mjtot,mbc,dt,dtnew,dx,dy,
153152
154153
155154c ::::::::::::::::::::::::fgout Output:::::::::::::::::::::::::::::::::
156- c This has been moved to tick.f, after advancing all patches on
155+ c This has been moved to tick.f, after advancing all patches on
157156c finest level. No need to check on each patch separately.
158157c ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
159158
160159c This call has been moved out to advanc
161160c call b4step2(mbc,mx,my,nvar,q,
162161c & xlowmbc,ylowmbc,dx,dy,time,dt,maux,aux,actualstep)
163-
162+
164163c ::::::::::::::::::::::::FGOUT DATA before step:::::::::::::::::::::::
165164c # fill in values at fgout points affected at time tc0
166165 do ng= 1 ,FGOUT_num_grids
@@ -191,7 +190,7 @@ subroutine stepgrid(q,fm,fp,gm,gp,mitot,mjtot,mbc,dt,dtnew,dx,dy,
191190 & q,aux,dx,dy,dt,cflgrid,
192191 & fm,fp,gm,gp,rpn2,rpt2)
193192c
194- c
193+ c
195194 mptr_level = node(nestlevel,mptr)
196195
197196c write(outunit,811) mptr, mptr_level, cflgrid
@@ -324,9 +323,7 @@ subroutine stepgrid(q,fm,fp,gm,gp,mitot,mjtot,mbc,dt,dtnew,dx,dy,
324323 endif
325324
326325#ifdef WHERE_AM_I
327- write (* ,* ) " ending stepgrid grid " ,mptr
326+ write (* ,* ) " ending stepgrid grid " ,mptr
328327#endif
329328 return
330329 end
331-
332-
0 commit comments