natl_12/code/ini_fields.F

C $Header: /u/gcmpack/MITgcm/model/src/ini_fields.F,v 1.15 2003/06/05 16:03:05 adcroft Exp $
C $Name: checkpoint51c_post $

#include "CPP_OPTIONS.h"

CBOP
C     !ROUTINE: INI_FIELDS
C     !INTERFACE:
      SUBROUTINE INI_FIELDS( myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE INI_FIELDS                                     
C     | o Initialise model fields                                 
C     *==========================================================*
C     | Initial conditions of the model state variables are       
C     | set here. For continuation runs a restart file will be    
C     | used. For original runs custom rules can be included      
C     | here.                                                     
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     myThid -  Number of this instance of INI_FIELDS
      INTEGER myThid
CEOP

C--   Set model variables to initial/restart conditions
      IF ( startTime .EQ. 0. .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
       IF (nonHydrostatic) CALL INI_PNH( myThid )
      ELSE

CcnhTemporaryFixStarts
C      Call ini_theta and salt so that sponge thetaRef and saltRef get set
C      for pickup. Values of theta and salt are overwritten from pickup 
C      'CALL READ_CHECKPOINT' immediately below
       CALL INI_THETA( myThid )
       CALL INI_SALT( myThid )
CcnhTemporaryFixEnds
       IF ( nIter0 .EQ. 0 ) THEN
        nIter0 = INT(startTime/deltaTClock )
       ENDIF
       CALL READ_CHECKPOINT( nIter0, myThid )
      ENDIF

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/model/src/ini_fields.F,v 1.15 2003/06/05 16:03:05 adcroft Exp $
2	C $Name: checkpoint51c_post $
3
4	#include "CPP_OPTIONS.h"
5
6	CBOP
7	C !ROUTINE: INI_FIELDS
8	C !INTERFACE:
9	SUBROUTINE INI_FIELDS( myThid )
10	C !DESCRIPTION: \bv
11	C ==========================================================
12	C \| SUBROUTINE INI_FIELDS
13	C \| o Initialise model fields
14	C ==========================================================
15	C \| Initial conditions of the model state variables are
16	C \| set here. For continuation runs a restart file will be
17	C \| used. For original runs custom rules can be included
18	C \| here.
19	C ==========================================================
20	C \ev
21
22	C !USES:
23	IMPLICIT NONE
24	C === Global variables ===
25	#include "SIZE.h"
26	#include "EEPARAMS.h"
27	#include "PARAMS.h"
28
29	C !INPUT/OUTPUT PARAMETERS:
30	C == Routine arguments ==
31	C myThid - Number of this instance of INI_FIELDS
32	INTEGER myThid
33	CEOP
34
35	C-- Set model variables to initial/restart conditions
36	IF ( startTime .EQ. 0. .AND. nIter0 .EQ. 0
37	& .AND. pickupSuff .EQ. ' ' ) THEN
38	CALL INI_VEL( myThid )
39	CALL INI_THETA( myThid )
40	CALL INI_SALT( myThid )
41	CALL INI_PSURF( myThid )
42	CALL INI_PRESSURE( myThid )
43	#ifdef INCLUDE_EP_FORCING_CODE
44	CALL INI_EP( myThid )
45	#endif
46	IF (nonHydrostatic) CALL INI_PNH( myThid )
47	ELSE
48
49	CcnhTemporaryFixStarts
50	C Call ini_theta and salt so that sponge thetaRef and saltRef get set
51	C for pickup. Values of theta and salt are overwritten from pickup
52	C 'CALL READ_CHECKPOINT' immediately below
53	CALL INI_THETA( myThid )
54	CALL INI_SALT( myThid )
55	CcnhTemporaryFixEnds
56	IF ( nIter0 .EQ. 0 ) THEN
57	nIter0 = INT(startTime/deltaTClock )
58	ENDIF
59	CALL READ_CHECKPOINT( nIter0, myThid )
60	ENDIF
61
62	RETURN
63	END