model/src/ini_fields.F

C $Header: /u/gcmpack/MITgcm/model/src/ini_fields.F,v 1.23 2009/05/17 20:15:22 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 */

      RETURN
      END
1	jmc	1.24	C $Header: /u/gcmpack/MITgcm/model/src/ini_fields.F,v 1.23 2009/05/17 20:15:22 jmc Exp $
2	adcroft	1.10	C $Name: $
3	cnh	1.1
4	edhill	1.16	#include "PACKAGES_CONFIG.h"
5	cnh	1.4	#include "CPP_OPTIONS.h"
6	cnh	1.1
7	edhill	1.19	C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
8	cnh	1.12	CBOP
9			C !ROUTINE: INI_FIELDS
10	edhill	1.19
11	cnh	1.12	C !INTERFACE:
12	cnh	1.1	SUBROUTINE INI_FIELDS( myThid )
13	edhill	1.19
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	cnh	1.12
19			C !USES:
20	adcroft	1.5	IMPLICIT NONE
21	cnh	1.1	#include "SIZE.h"
22			#include "EEPARAMS.h"
23			#include "PARAMS.h"
24
25	edhill	1.19	C !INPUT PARAMETERS:
26	cnh	1.1	INTEGER myThid
27	cnh	1.12	CEOP
28	cnh	1.1
29	edhill	1.19	C Set model variables to initial/restart conditions
30	jmc	1.20	IF ( startTime .EQ. baseTime .AND. nIter0 .EQ. 0
31	adcroft	1.15	& .AND. pickupSuff .EQ. ' ' ) THEN
32	edhill	1.19	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	heimbach	1.13	#ifdef INCLUDE_EP_FORCING_CODE
38	edhill	1.19	CALL INI_EP( myThid )
39	heimbach	1.13	#endif
40	jmc	1.21	ELSEIF ( .NOT.useOffLine ) THEN
41	jmc	1.22	CALL READ_PICKUP( nIter0, myThid )
42	edhill	1.19
43	jmc	1.23	c IF (pickup_write_immed) THEN
44			c CALL WRITE_PICKUP(
45			c I .TRUE., startTime, nIter0, myThid )
46			c ENDIF
47	edhill	1.19
48	heimbach	1.17	ENDIF
49	cnh	1.1
50	jmc	1.24	#ifdef ALLOW_NONHYDROSTATIC
51			IF ( nonHydrostatic ) THEN
52			CALL INI_NH_FIELDS( myThid )
53			ENDIF
54			#endif /* ALLOW_NONHYDROSTATIC */
55
56	cnh	1.1	RETURN
57			END