/[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.6 - (show annotations) (download)
Thu Nov 8 20:57:51 2001 UTC (22 years, 6 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint44e_post, checkpoint46g_pre, checkpoint46f_post, checkpoint44f_post, checkpoint46b_post, checkpoint43a-release1mods, chkpt44d_post, checkpoint44e_pre, checkpoint46d_pre, release1-branch_tutorials, checkpoint45d_post, chkpt44a_post, checkpoint44h_pre, checkpoint46a_post, chkpt44c_pre, checkpoint45a_post, checkpoint44g_post, checkpoint46e_pre, checkpoint45b_post, checkpoint46b_pre, release1-branch-end, release1_final_v1, checkpoint46c_pre, checkpoint46, checkpoint44b_post, checkpoint46a_pre, checkpoint45c_post, checkpoint44h_post, chkpt44a_pre, checkpoint46c_post, checkpoint46e_post, checkpoint44b_pre, checkpoint44, checkpoint45, chkpt44c_post, checkpoint44f_pre, checkpoint46d_post, release1-branch_branchpoint
Branch point for: release1_final, release1-branch
Changes since 1.5: +2 -2 lines 
Preparing adjoint of Held-Suarez:
- bugfix for storing in absence of CD code
- adding EXACT_CONSERV to AD list
- new routine ini_autodiff to add TAMC-specific initialisations
- adding Shapiro filter to AD list
1 C $Header: /u/gcmpack/models/MITgcmUV/model/src/ini_dynvars.F,v 1.5 2001/09/26 18:09:15 cnh Exp $
2 C $Name: checkpoint43 $
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 #ifdef ALLOW_NONHYDROSTATIC
29 #include "GW.h"
30 #endif
31
32 C !INPUT/OUTPUT PARAMETERS:
33 C == Routine arguments ==
34 C myThid - Number of this instance of INI_UVEL
35 INTEGER myThid
36
37 C !LOCAL VARIABLES:
38 C == Local variables ==
39 C bi,bj - Loop counters
40 C I,J,K
41 INTEGER bi, bj
42 INTEGER I, J, K
43 CEOP
44
45 C-- Over all tiles
46 DO bj = myByLo(myThid), myByHi(myThid)
47 DO bi = myBxLo(myThid), myBxHi(myThid)
48
49 C- 3D arrays
50 DO K=1,Nr
51 DO J=1-Oly,sNy+Oly
52 DO I=1-Olx,sNx+Olx
53
54 uVel(I,J,K,bi,bj)=0. _d 0
55 vVel(I,J,K,bi,bj)=0. _d 0
56 wVel(I,J,K,bi,bj)=0. _d 0
57 theta(I,J,K,bi,bj)=0. _d 0
58 salt(I,J,K,bi,bj)=0. _d 0
59
60 gU(I,J,K,bi,bj)=0. _d 0
61 gV(I,J,K,bi,bj)=0. _d 0
62 gT(I,J,K,bi,bj)=0. _d 0
63 gS(I,J,K,bi,bj)=0. _d 0
64 gUnm1(I,J,K,bi,bj)=0. _d 0
65 gVnm1(I,J,K,bi,bj)=0. _d 0
66 gTnm1(I,J,K,bi,bj)=0. _d 0
67 gSnm1(I,J,K,bi,bj)=0. _d 0
68
69 #ifdef INCLUDE_CD_CODE
70 guCD(I,J,K,bi,bj)=0. _d 0
71 gvCD(I,J,K,bi,bj)=0. _d 0
72 uNM1(I,J,K,bi,bj)=0. _d 0
73 vNM1(I,J,K,bi,bj)=0. _d 0
74 uVeld(I,J,K,bi,bj)=0. _d 0
75 vVeld(I,J,K,bi,bj)=0. _d 0
76 #endif
77
78 #ifdef ALLOW_NONHYDROSTATIC
79 gW(I,J,K,bi,bj)=0. _d 0
80 gWnm1(I,J,K,bi,bj)=0. _d 0
81 #endif
82
83 ENDDO
84 ENDDO
85 ENDDO
86
87 C- 2D arrays
88 DO J=1-Oly,sNy+Oly
89 DO I=1-Olx,sNx+Olx
90 etaN(I,J,bi,bj)=0. _d 0
91 etaH(I,J,bi,bj)=0. _d 0
92 #ifdef INCLUDE_CD_CODE
93 etaNm1(I,J,bi,bj)=0. _d 0
94 #endif
95 ENDDO
96 ENDDO
97
98 ENDDO
99 ENDDO
100
101 RETURN
102 END
  ViewVC Help
Powered by ViewVC 1.1.22