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

Diff of /MITgcm/model/src/initialise_varia.F

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

revision 1.11 by heimbach, Tue Apr 10 22:35:25 2001 UTC revision 1.18 by heimbach, Mon Sep 10 16:35:27 2001 UTC
# Line 22  C     == Global variables == Line 22  C     == Global variables ==
22  #include "SIZE.h"  #include "SIZE.h"
23  #include "EEPARAMS.h"  #include "EEPARAMS.h"
24  #include "PARAMS.h"  #include "PARAMS.h"
 #ifdef ALLOW_OBCS  
25  #include "DYNVARS.h"  #include "DYNVARS.h"
 #endif  
26    
27  C     == Routine arguments ==  C     == Routine arguments ==
28        INTEGER myThid        INTEGER myThid
# Line 33  CEndOfInterface Line 31  CEndOfInterface
31  C     == Local variables ==  C     == Local variables ==
32        INTEGER bi,bj,K,iMin,iMax,jMin,jMax        INTEGER bi,bj,K,iMin,iMax,jMin,jMax
33    
34    #ifdef ALLOW_TAMC_CHECKPOINTING
35    
36          nIter0 = INT( startTime/deltaTClock )
37    
38    C--   Set Bo_surf => define the Linear Relation: Phi_surf(eta)
39          CALL INI_LINEAR_PHISURF( myThid )
40    
41    C--   Set coriolis operators
42          CALL INI_CORI( myThid )
43    
44    C--   Set laplace operators for use in 2D conjugate gradient solver.
45          CALL INI_CG2D( myThid )
46    
47    #ifdef ALLOW_NONHYDROSTATIC
48    C--   Set laplace operators for use in 3D conjugate gradient solver.
49          CALL INI_CG3D( myThid )
50    #endif
51    
52    #endif /* ALLOW_TAMC_CHECKPOINTING */
53          _BARRIER
54    
55    C--   Initialise 3-dim. diffusivities
56          CALL INI_MIXING( myThid )
57        _BARRIER        _BARRIER
58    
59    C--   Initialize DYNVARS arrays (state fields + G terms: Gu,Gv,...) to zero [always]
60          CALL INI_DYNVARS( myThid )
61    
62  C--   Initialise model fields.  C--   Initialise model fields.
63  C     Starting values of U, V, W, temp., salt. and tendency terms  C     Starting values of U, V, W, temp., salt. and tendency terms
64  C     are set here. Fields are either set to default or read from  C     are set here. Fields are either set to default or read from
# Line 41  C     stored files. Line 66  C     stored files.
66        CALL INI_FIELDS( myThid )        CALL INI_FIELDS( myThid )
67        _BARRIER        _BARRIER
68    
69    #ifdef ALLOW_PASSIVE_TRACER
70    C--   Initialise passive tracer(s)
71          CALL INI_TR1( myThid )
72          _BARRIER
73    #endif
74    
75    #ifndef ALLOW_AUTODIFF_TAMC
76        IF ( usePickupBeforeC35 ) THEN        IF ( usePickupBeforeC35 ) THEN
77  C-- IMPORTANT : Need to activate the following call to restart from  C-- IMPORTANT : Need to activate the following call to restart from
78  C     a pickup file written by MITgcmUV_checkpoint34 or earlier.  C     a pickup file written by MITgcmUV_checkpoint34 or earlier.
# Line 48  C     a pickup file written by MITgcmUV_ Line 80  C     a pickup file written by MITgcmUV_
80         CALL THE_CORRECTION_STEP(startTime, nIter0, myThid)         CALL THE_CORRECTION_STEP(startTime, nIter0, myThid)
81        ENDIF        ENDIF
82        ENDIF        ENDIF
83    #endif
84    
85  C--   Initial conditions are convectively adjusted (for historical reasons)  C--   Initial conditions are convectively adjusted (for historical reasons)
86        IF ( startTime .EQ. 0. ) THEN        IF ( startTime .EQ. 0. ) THEN
87    CADJ loop = parallel
88          DO bj = myByLo(myThid), myByHi(myThid)          DO bj = myByLo(myThid), myByHi(myThid)
89    CADJ loop = parallel
90           DO bi = myBxLo(myThid), myBxHi(myThid)           DO bi = myBxLo(myThid), myBxHi(myThid)
91             iMin=1-Olx             iMin=1-Olx
92             iMax=sNx+Olx             iMax=sNx+Olx
# Line 68  C--   Initial conditions are convectivel Line 103  C--   Initial conditions are convectivel
103  C--   Initialize variable data for packages  C--   Initialize variable data for packages
104        CALL PACKAGES_INIT_VARIABLES( myThid )        CALL PACKAGES_INIT_VARIABLES( myThid )
105    
106    #ifdef NONLIN_FRSURF
107    C--   Compute the surface level thickness, function of eta_n-1
108    C     and modify hFac(C,W,S) accordingly :
109          IF (nonlinFreeSurf.GT.0) THEN
110            CALL CALC_SURF_DR(etaNm1, startTime, nIter0, myThid )
111            CALL UPDATE_SURF_DR( startTime, nIter0, myThid )
112          ENDIF
113    C-    update also CG2D matrix (and preconditioner)
114          IF ( nonlinFreeSurf.GT.2) THEN
115            CALL UPDATE_CG2D( startTime, nIter0, myThid )
116          ENDIF
117    #endif
118    
119  C--   Finally summarise the model state  C--   Finally summarise the model state
120        CALL STATE_SUMMARY( myThid )        CALL STATE_SUMMARY( myThid )
121    
# Line 82  C--   initialise time-average arrays wit Line 130  C--   initialise time-average arrays wit
130        ENDIF        ENDIF
131  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
132    
133    #ifdef EXACT_CONSERV
134          IF (exactConserv) THEN
135    C--   Compute the initial Barotropic Flow Divergence :
136            DO bj=myByLo(myThid),myByHi(myThid)
137             DO bi=myBxLo(myThid),myBxHi(myThid)
138               CALL CALC_EXACT_ETA( bi,bj, uVel,vVel,
139         I                          startTime, nIter0, myThid )
140             ENDDO
141            ENDDO
142          ENDIF
143    #endif /* EXACT_CONSERV */
144          
145        RETURN        RETURN
146        END        END
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22