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

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

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

Revision 1.27 - (show 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 C $Header: /u/gcmpack/MITgcm/model/src/ini_dynvars.F,v 1.26 2014/08/15 19:20:57 jmc Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: INI_DYNVARS
8 C !INTERFACE:
9 SUBROUTINE INI_DYNVARS( myThid )
10 C !DESCRIPTION: \bv
11 C *==========================================================*
12 C | SUBROUTINE INI_DYNVARS
13 C | o Initialise to zero all DYNVARS.h arrays
14 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 IMPLICIT NONE
23 C === Global variables ===
24 #include "SIZE.h"
25 #include "EEPARAMS.h"
26 #include "PARAMS.h"
27 #include "DYNVARS.h"
28
29 C !INPUT/OUTPUT PARAMETERS:
30 C == Routine arguments ==
31 C myThid :: my Thread Id. number
32 INTEGER myThid
33
34 C !LOCAL VARIABLES:
35 C == Local variables ==
36 C bi,bj :: tile indices
37 C i,j,k :: loop index
38 INTEGER bi, bj
39 INTEGER i, j, k
40 CEOP
41
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 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 #ifdef USE_OLD_EXTERNAL_FORCING
60 gT(i,j,k,bi,bj) = 0. _d 0
61 gS(i,j,k,bi,bj) = 0. _d 0
62 #endif
63 #ifdef ALLOW_ADAMSBASHFORTH_3
64 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 #else
73 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 #endif
78 totPhiHyd (i,j,k,bi,bj) = 0. _d 0
79 rhoInSitu (i,j,k,bi,bj) = 0. _d 0
80 IVDConvCount(i,j,k,bi,bj) = 0. _d 0
81 ENDDO
82 ENDDO
83 ENDDO
84
85 C- 2D arrays
86 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 ENDDO
93 ENDDO
94
95 ENDDO
96 ENDDO
97
98 RETURN
99 END
  ViewVC Help
Powered by ViewVC 1.1.22