pkg/fizhi/step_fizhi_fg.F

       subroutine step_fizhi_fg (myTime, myIter, myThid, dt)
c-----------------------------------------------------------------------
c  Subroutine step_fizhi_fg - 'Wrapper' routine to advance
c        the physics state and make a 'first guess' at the new
c        value. At this point, increment with the physics 
c        tendency only. 
c        Also: Set up "bi, bj loop" and some timers and clocks here.
c Call: step_physics
c-----------------------------------------------------------------------
       implicit none
#include "CPP_OPTIONS.h"
#include "SIZE.h"
#include "fizhi_SIZE.h"
#include "land_SIZE.h"
#include "fizhi_coms.h"
#include "EEPARAMS.h"

       integer myTime, myIter, myThid 

       integer i, j, L, bi, bj
       integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2
       _RL dt

       im1 = 1-OLx
       im2 = sNx+OLx
       jm1 = 1-OLy
       jm2 = sNy+OLy
       idim1 = 1
       idim2 = sNx
       jdim1 = 1
       jdim2 = sNy

       do bj = myByLo(myThid), myByHi(myThid)
       do bi = myBxLo(myThid), myBxHi(myThid)

c  Step forward the physics state using physics tendencies
        call step_physics(uphy,vphy,thphy,sphy,dt,im1,im2,jm1,jm2,
     .      Nrphys,Nsx,Nsy,1,sNx,1,sNy,bi,bj,duphy,dvphy,dthphy,dsphy)

       enddo
       enddo

       return
       end
1	subroutine step_fizhi_fg (myTime, myIter, myThid, dt)
2	c-----------------------------------------------------------------------
3	c Subroutine step_fizhi_fg - 'Wrapper' routine to advance
4	c the physics state and make a 'first guess' at the new
5	c value. At this point, increment with the physics
6	c tendency only.
7	c Also: Set up "bi, bj loop" and some timers and clocks here.
8	c Call: step_physics
9	c-----------------------------------------------------------------------
10	implicit none
11	#include "CPP_OPTIONS.h"
12	#include "SIZE.h"
13	#include "fizhi_SIZE.h"
14	#include "land_SIZE.h"
15	#include "fizhi_coms.h"
16	#include "EEPARAMS.h"
17
18	integer myTime, myIter, myThid
19
20	integer i, j, L, bi, bj
21	integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2
22	_RL dt
23
24	im1 = 1-OLx
25	im2 = sNx+OLx
26	jm1 = 1-OLy
27	jm2 = sNy+OLy
28	idim1 = 1
29	idim2 = sNx
30	jdim1 = 1
31	jdim2 = sNy
32
33	do bj = myByLo(myThid), myByHi(myThid)
34	do bi = myBxLo(myThid), myBxHi(myThid)
35
36	c Step forward the physics state using physics tendencies
37	call step_physics(uphy,vphy,thphy,sphy,dt,im1,im2,jm1,jm2,
38	. Nrphys,Nsx,Nsy,1,sNx,1,sNy,bi,bj,duphy,dvphy,dthphy,dsphy)
39
40	enddo
41	enddo
42
43	return
44	end