22 |
#include "SIZE.h" |
#include "SIZE.h" |
23 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
24 |
#include "PARAMS.h" |
#include "PARAMS.h" |
|
#ifdef ALLOW_OBCS |
|
25 |
#include "DYNVARS.h" |
#include "DYNVARS.h" |
|
#endif |
|
26 |
|
|
27 |
C == Routine arguments == |
C == Routine arguments == |
28 |
INTEGER myThid |
INTEGER myThid |
32 |
INTEGER bi,bj,K,iMin,iMax,jMin,jMax |
INTEGER bi,bj,K,iMin,iMax,jMin,jMax |
33 |
|
|
34 |
_BARRIER |
_BARRIER |
35 |
|
C-- Initialize DYNVARS arrays (state fields + G terms: Gu,Gv,...) to zero [always] |
36 |
|
CALL INI_DYNVARS( myThid ) |
37 |
|
|
38 |
C-- Initialise model fields. |
C-- Initialise model fields. |
39 |
C Starting values of U, V, W, temp., salt. and tendency terms |
C Starting values of U, V, W, temp., salt. and tendency terms |
40 |
C are set here. Fields are either set to default or read from |
C are set here. Fields are either set to default or read from |
83 |
ENDIF |
ENDIF |
84 |
#endif /* ALLOW_TIMEAVE */ |
#endif /* ALLOW_TIMEAVE */ |
85 |
|
|
86 |
|
#ifdef EXACT_CONSERV |
87 |
|
IF (exactConserv) THEN |
88 |
|
C-- Compute the initial Barotropic Flow Divergence : |
89 |
|
DO bj=myByLo(myThid),myByHi(myThid) |
90 |
|
DO bi=myBxLo(myThid),myBxHi(myThid) |
91 |
|
CALL CALC_EXACT_ETA( bi,bj, uVel,vVel, |
92 |
|
I startTime, nIter0, myThid ) |
93 |
|
ENDDO |
94 |
|
ENDDO |
95 |
|
ENDIF |
96 |
|
#endif /* EXACT_CONSERV */ |
97 |
|
|
98 |
RETURN |
RETURN |
99 |
END |
END |