fake advance a step to calculate magnetic field
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer | :: | igrid | ||||
| integer | :: | idimmin | ||||
| integer | :: | idimmax | ||||
| type(state) | :: | s |
subroutine fake_advance(igrid,idimmin,idimmax,s) use mod_global_parameters use mod_fix_conserve integer :: igrid,idimmin,idimmax type(state) :: s double precision :: dx1,dx2,dx3 double precision :: fC(ixGlo1:ixGhi1,ixGlo2:ixGhi2,& ixGlo3:ixGhi3,1:nwflux,1:ndim) double precision :: fE(ixGlo1:ixGhi1,ixGlo2:ixGhi2,& ixGlo3:ixGhi3,sdim:3) dx1=rnode(rpdx1_,igrid);dx2=rnode(rpdx2_,igrid);dx3=rnode(rpdx3_,igrid); call fake_update(ixGlo1,ixGlo2,ixGlo3,ixGhi1,ixGhi2,ixGhi3,s,fC,fE,dx1,dx2,& dx3) call store_flux(igrid,fC,idimmin,idimmax,3) call store_edge(igrid,ixGlo1,ixGlo2,ixGlo3,ixGhi1,ixGhi2,ixGhi3,fE,idimmin,& idimmax) end subroutine fake_advance