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

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

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

Revision 1.27 - (hide annotations) (download)
Tue Jan 20 20:46:55 2015 UTC (9 years, 4 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65i, checkpoint65n, checkpoint65l, checkpoint65m, checkpoint66a, checkpoint65o
Changes since 1.26: +1 -5 lines 
- move ALLOW_EDDYPSI block out of DYNVARS.h and merge it into FFIELDS.h
- rename uMean,vMean --> uEulerMean,vEulerMean (+ change name in pickup file)
- add frictionHeating field to pickup-files (for synchronous time-stepping)
1 jmc 1.27 C $Header: /u/gcmpack/MITgcm/model/src/ini_dynvars.F,v 1.26 2014/08/15 19:20:57 jmc Exp $
2 mlosch 1.7 C $Name: $
3 adcroft 1.2
4     #include "CPP_OPTIONS.h"
5    
6 cnh 1.5 CBOP
7     C !ROUTINE: INI_DYNVARS
8     C !INTERFACE:
9 adcroft 1.2 SUBROUTINE INI_DYNVARS( myThid )
10 cnh 1.5 C !DESCRIPTION: \bv
11     C *==========================================================*
12 jmc 1.20 C | SUBROUTINE INI_DYNVARS
13     C | o Initialise to zero all DYNVARS.h arrays
14 cnh 1.5 C *==========================================================*
15     C | Sets all the State variables to zero.
16     C | Sets all the Gs (arrays used for multi-level time-stepping)
17     C | to zero.
18     C *==========================================================*
19     C \ev
20    
21     C !USES:
22 adcroft 1.2 IMPLICIT NONE
23     C === Global variables ===
24     #include "SIZE.h"
25     #include "EEPARAMS.h"
26     #include "PARAMS.h"
27     #include "DYNVARS.h"
28    
29 cnh 1.5 C !INPUT/OUTPUT PARAMETERS:
30 adcroft 1.2 C == Routine arguments ==
31 jmc 1.20 C myThid :: my Thread Id. number
32 adcroft 1.2 INTEGER myThid
33    
34 cnh 1.5 C !LOCAL VARIABLES:
35 adcroft 1.2 C == Local variables ==
36 jmc 1.20 C bi,bj :: tile indices
37 jmc 1.24 C i,j,k :: loop index
38 adcroft 1.2 INTEGER bi, bj
39 jmc 1.24 INTEGER i, j, k
40 cnh 1.5 CEOP
41 adcroft 1.2
42     C-- Over all tiles
43     DO bj = myByLo(myThid), myByHi(myThid)
44     DO bi = myBxLo(myThid), myBxHi(myThid)
45    
46     C- 3D arrays
47 jmc 1.24 DO k=1,Nr
48     DO j=1-OLy,sNy+OLy
49     DO i=1-OLx,sNx+OLx
50    
51     uVel (i,j,k,bi,bj) = 0. _d 0
52     vVel (i,j,k,bi,bj) = 0. _d 0
53     wVel (i,j,k,bi,bj) = 0. _d 0
54     theta(i,j,k,bi,bj) = 0. _d 0
55     salt (i,j,k,bi,bj) = 0. _d 0
56    
57     gU(i,j,k,bi,bj) = 0. _d 0
58     gV(i,j,k,bi,bj) = 0. _d 0
59 jmc 1.26 #ifdef USE_OLD_EXTERNAL_FORCING
60 jmc 1.24 gT(i,j,k,bi,bj) = 0. _d 0
61     gS(i,j,k,bi,bj) = 0. _d 0
62 jmc 1.26 #endif
63 jmc 1.15 #ifdef ALLOW_ADAMSBASHFORTH_3
64 jmc 1.24 guNm(i,j,k,bi,bj,1) = 0. _d 0
65     gvNm(i,j,k,bi,bj,1) = 0. _d 0
66     gtNm(i,j,k,bi,bj,1) = 0. _d 0
67     gsNm(i,j,k,bi,bj,1) = 0. _d 0
68     guNm(i,j,k,bi,bj,2) = 0. _d 0
69     gvNm(i,j,k,bi,bj,2) = 0. _d 0
70     gtNm(i,j,k,bi,bj,2) = 0. _d 0
71     gsNm(i,j,k,bi,bj,2) = 0. _d 0
72 jmc 1.15 #else
73 jmc 1.24 guNm1(i,j,k,bi,bj) = 0. _d 0
74     gvNm1(i,j,k,bi,bj) = 0. _d 0
75     gtNm1(i,j,k,bi,bj) = 0. _d 0
76     gsNm1(i,j,k,bi,bj) = 0. _d 0
77 jmc 1.15 #endif
78 jmc 1.24 totPhiHyd (i,j,k,bi,bj) = 0. _d 0
79     rhoInSitu (i,j,k,bi,bj) = 0. _d 0
80 jmc 1.14 IVDConvCount(i,j,k,bi,bj) = 0. _d 0
81 adcroft 1.2 ENDDO
82     ENDDO
83     ENDDO
84    
85     C- 2D arrays
86 jmc 1.24 DO J=1-OLy,sNy+OLy
87     DO I=1-OLx,sNx+OLx
88     etaN(i,j,bi,bj) = 0. _d 0
89     etaH(i,j,bi,bj) = 0. _d 0
90     phiHydLow(i,j,bi,bj) = 0. _d 0
91     hMixLayer(i,j,bi,bj) = 0. _d 0
92 adcroft 1.2 ENDDO
93     ENDDO
94    
95     ENDDO
96     ENDDO
97    
98     RETURN
99     END
  ViewVC Help
Powered by ViewVC 1.1.22