/[MITgcm]/MITgcm/model/src/ini_fields.F
ViewVC logotype

Annotation of /MITgcm/model/src/ini_fields.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.9 - (hide annotations) (download)
Sun Mar 25 22:33:52 2001 UTC (23 years, 2 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint38, c37_adj, checkpoint39
Changes since 1.8: +3 -2 lines
Modifications and additions to enable automatic differentiation.
Detailed info's in doc/notes_c37_adj.txt

1 heimbach 1.9 C $Header: /u/gcmpack/models/MITgcmUV/model/src/ini_fields.F,v 1.8 2001/03/12 20:41:02 jmc Exp $
2     C $Name: checkpoint37 $
3 cnh 1.1
4 cnh 1.4 #include "CPP_OPTIONS.h"
5 cnh 1.1
6     CStartOfInterface
7     SUBROUTINE INI_FIELDS( myThid )
8     C /==========================================================\
9     C | SUBROUTINE INI_FIELDS |
10     C | o Initialise model fields |
11     C |==========================================================|
12     C | Initial conditions of the model state variables are |
13     C | set here. For continuation runs a restart file will be |
14     C | used. For original runs custom rules can be included |
15     C | here. |
16     C \==========================================================/
17 adcroft 1.5 IMPLICIT NONE
18 cnh 1.1
19     C === Global variables ===
20     #include "SIZE.h"
21     #include "EEPARAMS.h"
22     #include "PARAMS.h"
23    
24     C == Routine arguments ==
25     C myThid - Number of this instance of INI_FIELDS
26     INTEGER myThid
27     CEndOfInterface
28    
29     C-- Set model variables to initial/restart conditions
30 cnh 1.4 IF ( startTime .EQ. 0. .AND. nIter0 .EQ. 0 ) THEN
31 cnh 1.1 CALL INI_UVEL( myThid )
32     CALL INI_VVEL( myThid )
33 heimbach 1.9 CALL INI_WVEL( myThid )
34 cnh 1.1 CALL INI_THETA( myThid )
35     CALL INI_SALT( myThid )
36 adcroft 1.6 CALL INI_PSURF( myThid )
37 jmc 1.8 #ifdef ALLOW_NONHYDROSTATIC
38     IF ( nonHydrostatic ) THEN
39     CALL INI_GW( myThid )
40     ENDIF
41     #endif
42 cnh 1.1 ELSE
43 cnh 1.3 IF ( nIter0 .EQ. 0 ) THEN
44     nIter0 = INT(startTime/deltaTClock )
45     ENDIF
46     CALL READ_CHECKPOINT( nIter0, myThid )
47 cnh 1.1 ENDIF
48    
49     RETURN
50     END

  ViewVC Help
Powered by ViewVC 1.1.22