/[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.18 by heimbach, Mon Sep 10 16:35:27 2001 UTC revision 1.22 by jmc, Wed Jan 30 04:12:12 2002 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  CStartOfInterface  CBOP
7    C     !ROUTINE: INITIALISE_VARIA
8    C     !INTERFACE:
9        SUBROUTINE INITIALISE_VARIA(myThid)        SUBROUTINE INITIALISE_VARIA(myThid)
10  C     /==========================================================\  C     !DESCRIPTION: \bv
11  C     | SUBROUTINE INITIALISE_VARIA                              |  C     *==========================================================*
12  C     | o Set the initial conditions for dynamics variables      |  C     | SUBROUTINE INITIALISE_VARIA                              
13  C     |   and time dependent arrays                              |  C     | o Set the initial conditions for dynamics variables      
14  C     |==========================================================|  C     |   and time dependent arrays                              
15  C     | This routine reads/writes  data from an input file and   |  C     *==========================================================*
16  C     | from various binary files.                               |  C     | This routine reads/writes  data from an input file and    
17  C     | Each thread invokes an instance of this routine as does  |  C     | from various binary files.                                
18  C     | each process in a multi-process parallel environment like|  C     | Each thread invokes an instance of this routine as does  
19  C     | MPI.                                                     |  C     | each process in a multi-process parallel environment like
20  C     \==========================================================/  C     | MPI.                                                      
21        IMPLICIT NONE  C     *==========================================================*
22    C     \ev
23    
24    C     !CALLING SEQUENCE:
25    C     INITIALISE_VARIA
26    C       |
27    C       |-- INI_LINEAR_PHISURF
28    C       |
29    C       |-- INI_CORI
30    C       |
31    C       |-- INI_CG2D
32    C       |
33    C       |-- INI_CG3D
34    C       |
35    C       |-- INI_MIXING
36    C       |
37    C       |-- INI_DYNVARS
38    C       |
39    C       |-- INI_FIELDS
40    C       |
41    C       |-- INI_TR1
42    C       |
43    C       |-- THE_CORRECTION_STEP
44    C       |
45    C       |-- CALL CONVECTIVE_ADJUSTMENT_INI
46    C       |
47    C       |-- PACKAGES_INIT_VARIABLES
48    C       |
49    C       |-- CALC_SURF_DR
50    C       |
51    C       |-- UPDATE_SURF_DR
52    C       |
53    C       |-- UPDATE_CG2D
54    C       |
55    C       |-- STATE_SUMMARY
56    C       |
57    C       |-- TIMEAVE_STATVARS
58    
59    C     !USES:
60          IMPLICIT NONE
61  C     == Global variables ==  C     == Global variables ==
62  #include "SIZE.h"  #include "SIZE.h"
63  #include "EEPARAMS.h"  #include "EEPARAMS.h"
64  #include "PARAMS.h"  #include "PARAMS.h"
65  #include "DYNVARS.h"  #include "DYNVARS.h"
66    
67    C     !INPUT/OUTPUT PARAMETERS:
68  C     == Routine arguments ==  C     == Routine arguments ==
69        INTEGER myThid        INTEGER myThid
 CEndOfInterface  
70    
71    C     !LOCAL VARIABLES:
72  C     == Local variables ==  C     == Local variables ==
73        INTEGER bi,bj,K,iMin,iMax,jMin,jMax        INTEGER bi,bj,K,iMin,iMax,jMin,jMax
74    CEOP
75    
76  #ifdef ALLOW_TAMC_CHECKPOINTING  #ifdef ALLOW_TAMC_CHECKPOINTING
77    
# Line 72  C--   Initialise passive tracer(s) Line 114  C--   Initialise passive tracer(s)
114        _BARRIER        _BARRIER
115  #endif  #endif
116    
117    #ifdef ALLOW_AUTODIFF_TAMC
118    C--   Initialise active fields to help TAMC
119          CALL INI_AUTODIFF( myThid )
120          _BARRIER
121    #endif
122    
123  #ifndef ALLOW_AUTODIFF_TAMC  #ifndef ALLOW_AUTODIFF_TAMC
124        IF ( usePickupBeforeC35 ) THEN        IF ( usePickupBeforeC35 ) THEN
125  C-- IMPORTANT : Need to activate the following call to restart from  C-- IMPORTANT : Need to activate the following call to restart from
# Line 104  C--   Initialize variable data for packa Line 152  C--   Initialize variable data for packa
152        CALL PACKAGES_INIT_VARIABLES( myThid )        CALL PACKAGES_INIT_VARIABLES( myThid )
153    
154  #ifdef NONLIN_FRSURF  #ifdef NONLIN_FRSURF
155  C--   Compute the surface level thickness, function of eta_n-1  C--   Compute the surface level thickness <-- function of etaH(n)
156  C     and modify hFac(C,W,S) accordingly :  C     and modify hFac(C,W,S) accordingly :
157        IF (nonlinFreeSurf.GT.0) THEN        IF (nonlinFreeSurf.GT.0) THEN
158          CALL CALC_SURF_DR(etaNm1, startTime, nIter0, myThid )          CALL CALC_SURF_DR(etaH, startTime, -1 , myThid )
159          CALL UPDATE_SURF_DR( startTime, nIter0, myThid )          CALL UPDATE_SURF_DR( startTime, nIter0, myThid )
160        ENDIF        ENDIF
161  C-    update also CG2D matrix (and preconditioner)  C-    update also CG2D matrix (and preconditioner)
# Line 130  C--   initialise time-average arrays wit Line 178  C--   initialise time-average arrays wit
178        ENDIF        ENDIF
179  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
180    
 #ifdef EXACT_CONSERV  
       IF (exactConserv) THEN  
 C--   Compute the initial Barotropic Flow Divergence :  
         DO bj=myByLo(myThid),myByHi(myThid)  
          DO bi=myBxLo(myThid),myBxHi(myThid)  
            CALL CALC_EXACT_ETA( bi,bj, uVel,vVel,  
      I                          startTime, nIter0, myThid )  
          ENDDO  
         ENDDO  
       ENDIF  
 #endif /* EXACT_CONSERV */  
         
181        RETURN        RETURN
182        END        END
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22