/[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.10 - (hide annotations) (download)
Thu Oct 9 04:19:18 2003 UTC (20 years, 7 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint51o_pre, checkpoint51l_post, checkpoint51n_post, checkpoint51n_pre, checkpoint51l_pre, checkpoint51i_post, checkpoint51o_post, checkpoint51m_post
Branch point for: tg2-branch, checkpoint51n_branch
Changes since 1.9: +3 -1 lines
 o first check-in for the "branch-genmake2" merge
 o verification suite as run on shelley (gcc 3.2.2):

Wed Oct  8 23:42:29 EDT 2003
                T           S           U           V
G D M    c        m  s        m  s        m  s        m  s
E p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
N n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
2 d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .

OPTFILE=NONE

Y Y Y Y 13 16 16 16  0 16 16 16 16 16 16 16 16 13 12  0  0 pass  adjustment.128x64x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16  0  0 16 16  0  0 pass  adjustment.cs-32x32x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16 22  0 16 16 22  0 pass  adjust_nlfs.cs-32x32x1
Y Y Y Y -- 13 13 16 16 13 13 13 13 16 16 16 16 16 16 16 16 N/O   advect_cs
Y Y Y Y -- 22 16 16 16 16 16 16 13 16 16 16 16 16 16 16 16 N/O   advect_xy
Y Y Y Y -- 13 16 13 16 16 16 16 16 16 16 22 16 16 16 16 16 N/O   advect_xz
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  aim.5l_cs
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 16 16 16 16 13 16 pass  aim.5l_Equatorial_Channel
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 13 16 16 13 13 16 pass  aim.5l_LatLon
Y Y Y Y 13 16 16 16 16 16 16 16 16 16 13 12 13 13 16 13 16 pass  exp0
Y Y Y Y 14 16 16 16 16 16 16 16 22 16 16 16 13 16 16 22 16 pass  exp1
Y Y Y Y 13 13 16 13 16 16 16 16 16 13 13 16 16 13 13 13 13 pass  exp2
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  exp4
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 22 16 pass  exp5
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  front_relax
Y Y Y Y 14 16 16 13 13 16 16 13 13 16 13 13 16 12 13 13 16 pass  global_ocean.90x40x15
Y Y Y Y 10 16 16 13 13 16 13 16 16 13 13 13 13 16 16 13 16 FAIL  global_ocean.cs32x15
Y Y Y Y  6 11 12 13 13 12 13 16 13  9  9  9  9 10  9  9 11 FAIL  global_ocean_pressure
Y Y Y Y 14 16 16 13 16 16 16 13 13 13 13 13 16 12 16 13 16 pass  global_with_exf
Y Y Y Y 14 16 16 16 16 16 16 16 16 11 13 22 13 16 16  9 16 pass  hs94.128x64x5
Y Y Y Y 13 16 16 16 16 16 16 16 16 11 16 16 16 13 16 22 13 pass  hs94.1x64x5
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 13 13 16 16 22 13 pass  hs94.cs-32x32x5
Y Y Y Y 10 10 16 13 13 16 16 16 22 16 13 13 13 13 13 22 13 FAIL  ideal_2D_oce
Y Y Y Y  8 16 16 16 16 16 16 16 16 13 13  8 16 16 16 16 16 FAIL  internal_wave
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 13 22 13 13 13 22 16 pass  inverted_barometer
Y Y Y Y 12 16 16 16 16 16 16 16 16 16 13 12 13 13 13 13 13 FAIL  lab_sea
Y Y Y Y 11 16 16 16 16 16 16 16 13 13 13 12 13 16 13 12 13 FAIL  natl_box
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  plume_on_slope
Y Y Y Y 13 16 16 16 16 13 16 16 16 16 16 16 16 13 16 16 16 pass  solid-body.cs-32x32x1

1 edhill 1.10 C $Header: /u/u3/gcmpack/MITgcm/model/src/ini_dynvars.F,v 1.9.2.1 2003/10/02 18:10:45 edhill Exp $
2 mlosch 1.7 C $Name: $
3 adcroft 1.2
4 edhill 1.10 #include "PACKAGES_CONFIG.h"
5 adcroft 1.2 #include "CPP_OPTIONS.h"
6    
7 cnh 1.5 CBOP
8     C !ROUTINE: INI_DYNVARS
9     C !INTERFACE:
10 adcroft 1.2 SUBROUTINE INI_DYNVARS( myThid )
11 cnh 1.5 C !DESCRIPTION: \bv
12     C *==========================================================*
13     C | SUBROUTINE INI_DYNVARS
14     C | o Initialise to zero all DYNVARS.h arrays
15     C *==========================================================*
16     C | Sets all the State variables to zero.
17     C | Sets all the Gs (arrays used for multi-level time-stepping)
18     C | to zero.
19     C *==========================================================*
20     C \ev
21    
22     C !USES:
23 adcroft 1.2 IMPLICIT NONE
24     C === Global variables ===
25     #include "SIZE.h"
26     #include "EEPARAMS.h"
27     #include "PARAMS.h"
28     #include "DYNVARS.h"
29     #ifdef ALLOW_NONHYDROSTATIC
30     #include "GW.h"
31     #endif
32    
33 cnh 1.5 C !INPUT/OUTPUT PARAMETERS:
34 adcroft 1.2 C == Routine arguments ==
35     C myThid - Number of this instance of INI_UVEL
36     INTEGER myThid
37    
38 cnh 1.5 C !LOCAL VARIABLES:
39 adcroft 1.2 C == Local variables ==
40     C bi,bj - Loop counters
41     C I,J,K
42     INTEGER bi, bj
43     INTEGER I, J, K
44 cnh 1.5 CEOP
45 adcroft 1.2
46     C-- Over all tiles
47     DO bj = myByLo(myThid), myByHi(myThid)
48     DO bi = myBxLo(myThid), myBxHi(myThid)
49    
50     C- 3D arrays
51     DO K=1,Nr
52     DO J=1-Oly,sNy+Oly
53     DO I=1-Olx,sNx+Olx
54    
55 heimbach 1.3 uVel(I,J,K,bi,bj)=0. _d 0
56     vVel(I,J,K,bi,bj)=0. _d 0
57     wVel(I,J,K,bi,bj)=0. _d 0
58     theta(I,J,K,bi,bj)=0. _d 0
59     salt(I,J,K,bi,bj)=0. _d 0
60    
61     gU(I,J,K,bi,bj)=0. _d 0
62     gV(I,J,K,bi,bj)=0. _d 0
63     gT(I,J,K,bi,bj)=0. _d 0
64     gS(I,J,K,bi,bj)=0. _d 0
65     gUnm1(I,J,K,bi,bj)=0. _d 0
66     gVnm1(I,J,K,bi,bj)=0. _d 0
67     gTnm1(I,J,K,bi,bj)=0. _d 0
68     gSnm1(I,J,K,bi,bj)=0. _d 0
69 adcroft 1.2
70     #ifdef INCLUDE_CD_CODE
71 heimbach 1.3 uNM1(I,J,K,bi,bj)=0. _d 0
72     vNM1(I,J,K,bi,bj)=0. _d 0
73     uVeld(I,J,K,bi,bj)=0. _d 0
74     vVeld(I,J,K,bi,bj)=0. _d 0
75 adcroft 1.2 #endif
76    
77     #ifdef ALLOW_NONHYDROSTATIC
78 edhill 1.10 ceh3 needs an IF ( useNONHYDROSTATIC ) THEN
79 heimbach 1.3 gW(I,J,K,bi,bj)=0. _d 0
80     gWnm1(I,J,K,bi,bj)=0. _d 0
81 adcroft 1.2 #endif
82    
83 jmc 1.8 totPhiHyd(I,J,K,bi,bj)=0. _d 0
84 adcroft 1.2 ENDDO
85     ENDDO
86     ENDDO
87    
88     C- 2D arrays
89     DO J=1-Oly,sNy+Oly
90     DO I=1-Olx,sNx+Olx
91 heimbach 1.3 etaN(I,J,bi,bj)=0. _d 0
92 jmc 1.4 etaH(I,J,bi,bj)=0. _d 0
93     #ifdef INCLUDE_CD_CODE
94 heimbach 1.3 etaNm1(I,J,bi,bj)=0. _d 0
95 jmc 1.4 #endif
96 mlosch 1.7 phiHydLow(I,J,bi,bj)=0. _d 0
97 adcroft 1.2 ENDDO
98     ENDDO
99    
100     ENDDO
101     ENDDO
102    
103     RETURN
104     END

  ViewVC Help
Powered by ViewVC 1.1.22