model/src/ini_fields.F

a5bd9cbb61 Ed H*0001 #include "PACKAGES_CONFIG.h"
1dbaea09ee Chri*0002 #include "CPP_OPTIONS.h"
924557e60a Chri*0003
a30418b6b9 Ed H*0004 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
9366854e02 Chri*0005 CBOP
                0006 C     !ROUTINE: INI_FIELDS
a30418b6b9 Ed H*0007
9366854e02 Chri*0008 C     !INTERFACE:
924557e60a Chri*0009       SUBROUTINE INI_FIELDS( myThid )
a30418b6b9 Ed H*0010
                0011 C     !DESCRIPTION:
                0012 C     Initial conditions of the model state variables are set here. For
                0013 C     continuation runs a restart file will be used. For original runs
                0014 C     custom rules can be included here.
924557e60a Chri*0015
9366854e02 Chri*0016 C     !USES:
                0017       IMPLICIT NONE
924557e60a Chri*0018 #include "SIZE.h"
                0019 #include "EEPARAMS.h"
                0020 #include "PARAMS.h"
                0021
a30418b6b9 Ed H*0022 C     !INPUT PARAMETERS:
924557e60a Chri*0023       INTEGER myThid
9366854e02 Chri*0024 CEOP
924557e60a Chri*0025
a30418b6b9 Ed H*0026 C     Set model variables to initial/restart conditions
51d88069bf Jean*0027       IF ( startTime .EQ. baseTime .AND.  nIter0 .EQ. 0
b7bf24ec60 Alis*0028      &     .AND. pickupSuff .EQ. ' ' ) THEN
a30418b6b9 Ed H*0029         CALL INI_VEL( myThid )
                0030         CALL INI_THETA( myThid )
                0031         CALL INI_SALT( myThid )
                0032         CALL INI_PSURF( myThid )
                0033         CALL INI_PRESSURE( myThid )
b275747e24 Patr*0034 #ifdef INCLUDE_EP_FORCING_CODE
a30418b6b9 Ed H*0035         CALL INI_EP( myThid )
b275747e24 Patr*0036 #endif
b8c40edd31 Jean*0037       ELSEIF ( .NOT.useOffLine .OR. nonlinFreeSurf.GT.0 ) THEN
931983c6ab Jean*0038         CALL READ_PICKUP( nIter0, myThid )
c0d4d5aecc Patr*0039       ENDIF
a5bd9cbb61 Ed H*0040
5700f22e40 Jean*0041 #ifdef ALLOW_NONHYDROSTATIC
                0042       IF ( nonHydrostatic ) THEN
                0043         CALL INI_NH_FIELDS( myThid )
                0044       ENDIF
                0045 #endif /* ALLOW_NONHYDROSTATIC */
                0046
9313308f9b Jean*0047 #ifdef NONLIN_FRSURF
                0048 # ifndef DISABLE_SIGMA_CODE
                0049       IF ( selectSigmaCoord.NE.0 ) THEN
                0050         CALL UPDATE_ETAWS( startTime, nIter0, myThid )
                0051       ENDIF
                0052 # endif /* DISABLE_SIGMA_CODE */
                0053 #endif /* NONLIN_FRSURF */
                0054
924557e60a Chri*0055       RETURN
                0056       END