21 |
#include "SIZE.h" |
#include "SIZE.h" |
22 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
23 |
#include "PARAMS.h" |
#include "PARAMS.h" |
24 |
|
#ifdef ALLOW_OBCS |
25 |
|
#include "DYNVARS.h" |
26 |
|
#endif |
27 |
|
|
28 |
C == Routine arguments == |
C == Routine arguments == |
29 |
INTEGER myThid |
INTEGER myThid |
30 |
CEndOfInterface |
CEndOfInterface |
31 |
|
|
32 |
C == Local variables == |
C == Local variables == |
33 |
#ifdef ALLOW_OBCS |
INTEGER bi,bj,K,iMin,iMax,jMin,jMax |
|
INTEGER bi,bj,K |
|
|
#endif |
|
34 |
|
|
35 |
|
_BARRIER |
36 |
C-- Initialise model fields. |
C-- Initialise model fields. |
37 |
C Starting values of U, V, W, temp., salt. and tendency terms |
C Starting values of U, V, W, temp., salt. and tendency terms |
38 |
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 |
40 |
CALL INI_FIELDS( myThid ) |
CALL INI_FIELDS( myThid ) |
41 |
_BARRIER |
_BARRIER |
42 |
|
|
43 |
C-- Initialise time-stepping (predictor variables are carried in G) |
C-- Initial conditions are convectively adjusted (for historical reasons) |
44 |
IF ( startTime .EQ. 0. ) THEN |
IF ( startTime .EQ. 0. ) THEN |
|
CALL INI_PREDICTOR( myThid ) |
|
|
_BARRIER |
|
|
END IF |
|
|
|
|
|
C-- Check dependances between packages |
|
|
CALL PACKAGES_INIT( myThid ) |
|
|
|
|
|
#ifdef ALLOW_OBCS |
|
|
C-- Calculate values on open boundaries at iter=nIter0 |
|
|
IF (openBoundaries) THEN |
|
45 |
DO bj = myByLo(myThid), myByHi(myThid) |
DO bj = myByLo(myThid), myByHi(myThid) |
46 |
DO bi = myBxLo(myThid), myBxHi(myThid) |
DO bi = myBxLo(myThid), myBxHi(myThid) |
47 |
DO K=1,Nr |
iMin=1-Olx |
48 |
CALL SET_OBCS( K, bi, bj, startTime, myThid ) |
iMax=sNx+Olx |
49 |
Caja CALL CYCLE_OBCS( K, bi, bj, myThid ) |
jMin=1-Oly |
50 |
ENDDO |
jMax=sNy+Oly |
51 |
|
CALL CONVECTIVE_ADJUSTMENT( |
52 |
|
I bi, bj, iMin, iMax, jMin, jMax, |
53 |
|
I startTime, 0, myThid ) |
54 |
ENDDO |
ENDDO |
55 |
ENDDO |
ENDDO |
56 |
ENDIF |
_BARRIER |
57 |
#endif |
END IF |
58 |
|
|
59 |
|
C-- Initialize variable data for packages |
60 |
|
CALL PACKAGES_INIT_VARIABLES( myThid ) |
61 |
|
|
62 |
C-- Finally summarise the model state |
C-- Finally summarise the model state |
63 |
CALL STATE_SUMMARY( myThid ) |
CALL STATE_SUMMARY( myThid ) |