31 |
INTEGER bi,bj,K |
INTEGER bi,bj,K |
32 |
#endif |
#endif |
33 |
|
|
|
#ifdef ALLOW_GMREDI |
|
|
C-- Initialize GM/Redi parameterization |
|
|
IF (use_GMRedi) CALL GMREDI_INIT( myThid ) |
|
|
#endif |
|
|
|
|
|
#ifdef ALLOW_KPP |
|
|
C-- Initialize KPP vertical mixing scheme. |
|
|
IF (use_KPPmixing) THEN |
|
|
CALL KPP_INIT( myThid ) |
|
|
CALL KPP_OPEN_DIAGS( myThid ) |
|
|
ENDIF |
|
|
#endif |
|
|
|
|
34 |
C-- Initialise model fields. |
C-- Initialise model fields. |
35 |
C Starting values of U, V, W, temp., salt. and tendency terms |
C Starting values of U, V, W, temp., salt. and tendency terms |
36 |
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 |
38 |
CALL INI_FIELDS( myThid ) |
CALL INI_FIELDS( myThid ) |
39 |
_BARRIER |
_BARRIER |
40 |
|
|
41 |
C-- Initialise time-stepping (predictor variables are carried in G) |
C-- Check dependances between packages |
42 |
IF ( startTime .EQ. 0. ) |
CALL PACKAGES_INIT( myThid ) |
|
& CALL INI_PREDICTOR( myThid ) |
|
|
_BARRIER |
|
43 |
|
|
44 |
C-- Custom routine to set forcing fields |
C-- Initialise time-stepping (predictor variables are carried in G) |
45 |
CALL INI_FORCING( myThid ) |
IF ( startTime .EQ. 0. ) THEN |
46 |
|
CALL INI_PREDICTOR( myThid ) |
47 |
|
_BARRIER |
48 |
|
END IF |
49 |
|
|
50 |
#ifdef ALLOW_OBCS |
#ifdef ALLOW_OBCS |
51 |
C-- Calculate values on open boundaries at iter=nIter0 |
C-- Calculate values on open boundaries at iter=nIter0 |