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