C $Header: /u/gcmpack/MITgcm/model/src/ini_fields.F,v 1.27 2017/03/25 01:28:26 jmc Exp $ C $Name: $ #include "PACKAGES_CONFIG.h" #include "CPP_OPTIONS.h" C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| CBOP C !ROUTINE: INI_FIELDS C !INTERFACE: SUBROUTINE INI_FIELDS( myThid ) C !DESCRIPTION: C Initial conditions of the model state variables are set here. For C continuation runs a restart file will be used. For original runs C custom rules can be included here. C !USES: IMPLICIT NONE #include "SIZE.h" #include "EEPARAMS.h" #include "PARAMS.h" C !INPUT PARAMETERS: INTEGER myThid CEOP C Set model variables to initial/restart conditions IF ( startTime .EQ. baseTime .AND. nIter0 .EQ. 0 & .AND. pickupSuff .EQ. ' ' ) THEN CALL INI_VEL( myThid ) CALL INI_THETA( myThid ) CALL INI_SALT( myThid ) CALL INI_PSURF( myThid ) CALL INI_PRESSURE( myThid ) #ifdef INCLUDE_EP_FORCING_CODE CALL INI_EP( myThid ) #endif ELSEIF ( .NOT.useOffLine .OR. nonlinFreeSurf.GT.0 ) THEN CALL READ_PICKUP( nIter0, myThid ) ENDIF #ifdef ALLOW_NONHYDROSTATIC IF ( nonHydrostatic ) THEN CALL INI_NH_FIELDS( myThid ) ENDIF #endif /* ALLOW_NONHYDROSTATIC */ #ifdef NONLIN_FRSURF # ifndef DISABLE_SIGMA_CODE IF ( selectSigmaCoord.NE.0 ) THEN CALL UPDATE_ETAWS( startTime, nIter0, myThid ) ENDIF # endif /* DISABLE_SIGMA_CODE */ #endif /* NONLIN_FRSURF */ RETURN END