model/src/ini_fields.F

C $Header: /u/gcmpack/MITgcm/model/src/ini_fields.F,v 1.24 2009/12/11 13:56:28 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 ) THEN
        CALL READ_PICKUP( nIter0, myThid )

c       IF (pickup_write_immed) THEN
c         CALL WRITE_PICKUP(
c    I                       .TRUE., startTime, nIter0, myThid )
c       ENDIF

      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
1	C $Header: /u/gcmpack/MITgcm/model/src/ini_fields.F,v 1.24 2009/12/11 13:56:28 jmc Exp $
2	C $Name: $
3
4	#include "PACKAGES_CONFIG.h"
5	#include "CPP_OPTIONS.h"
6
7	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
8	CBOP
9	C !ROUTINE: INI_FIELDS
10
11	C !INTERFACE:
12	SUBROUTINE INI_FIELDS( myThid )
13
14	C !DESCRIPTION:
15	C Initial conditions of the model state variables are set here. For
16	C continuation runs a restart file will be used. For original runs
17	C custom rules can be included here.
18
19	C !USES:
20	IMPLICIT NONE
21	#include "SIZE.h"
22	#include "EEPARAMS.h"
23	#include "PARAMS.h"
24
25	C !INPUT PARAMETERS:
26	INTEGER myThid
27	CEOP
28
29	C Set model variables to initial/restart conditions
30	IF ( startTime .EQ. baseTime .AND. nIter0 .EQ. 0
31	& .AND. pickupSuff .EQ. ' ' ) THEN
32	CALL INI_VEL( myThid )
33	CALL INI_THETA( myThid )
34	CALL INI_SALT( myThid )
35	CALL INI_PSURF( myThid )
36	CALL INI_PRESSURE( myThid )
37	#ifdef INCLUDE_EP_FORCING_CODE
38	CALL INI_EP( myThid )
39	#endif
40	ELSEIF ( .NOT.useOffLine ) THEN
41	CALL READ_PICKUP( nIter0, myThid )
42
43	c IF (pickup_write_immed) THEN
44	c CALL WRITE_PICKUP(
45	c I .TRUE., startTime, nIter0, myThid )
46	c ENDIF
47
48	ENDIF
49
50	#ifdef ALLOW_NONHYDROSTATIC
51	IF ( nonHydrostatic ) THEN
52	CALL INI_NH_FIELDS( myThid )
53	ENDIF
54	#endif /* ALLOW_NONHYDROSTATIC */
55
56	#ifdef NONLIN_FRSURF
57	# ifndef DISABLE_SIGMA_CODE
58	IF ( selectSigmaCoord.NE.0 ) THEN
59	CALL UPDATE_ETAWS( startTime, nIter0, myThid )
60	ENDIF
61	# endif /* DISABLE_SIGMA_CODE */
62	#endif /* NONLIN_FRSURF */
63
64	RETURN
65	END