/[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.7 - (hide annotations) (download)
Wed Sep 18 16:38:02 2002 UTC (21 years, 8 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint47e_post, checkpoint46l_post, checkpoint47c_post, checkpoint48e_post, checkpoint46l_pre, checkpoint48b_post, checkpoint48c_pre, checkpoint47d_pre, checkpoint47a_post, checkpoint48d_pre, checkpoint47i_post, checkpoint47d_post, checkpoint48d_post, checkpoint48f_post, checkpoint46j_pre, checkpoint47g_post, checkpoint46j_post, checkpoint46k_post, checkpoint48a_post, checkpoint47j_post, branch-exfmods-tag, checkpoint48c_post, checkpoint47b_post, checkpoint46h_pre, checkpoint46m_post, checkpoint46g_post, checkpoint47f_post, checkpoint46i_post, checkpoint47, checkpoint48, checkpoint46h_post, checkpoint47h_post
Branch point for: branch-exfmods-curt
Changes since 1.6: +3 -2 lines
o Include a new diagnostic variable phiHydLow for the ocean model
  - in z-coordinates, it is the bottom pressure anomaly
  - in p-coordinates, it is the sea surface elevation
  - in both cases, these variable have global drift, reflecting the mass
    drift in z-coordinates and the volume drift in p-coordinates
  - included time averaging for phiHydLow, be aware of the drift!
o depth-dependent computation of Bo_surf for pressure coordinates
  in the ocean (buoyancyRelation='OCEANICP')
  - requires a new routine (FIND_RHO_SCALAR) to compute density with only
    Theta, Salinity, and Pressure in the parameter list. This routine is
    presently contained in find_rho.F. This routine does not give the
    correct density for 'POLY3', which would be a z-dependent reference
    density.
o cleaned up find_rho
  - removed obsolete 'eqn' from the parameter list.
o added two new verification experiments: gop and goz
  (4x4 degree global ocean, 15 layers in pressure and height coordinates)

1 mlosch 1.7 C $Header: /u/gcmpack/MITgcm/model/src/ini_dynvars.F,v 1.6 2001/11/08 20:57:51 heimbach Exp $
2     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     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 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     #ifdef ALLOW_NONHYDROSTATIC
29     #include "GW.h"
30     #endif
31    
32 cnh 1.5 C !INPUT/OUTPUT PARAMETERS:
33 adcroft 1.2 C == Routine arguments ==
34     C myThid - Number of this instance of INI_UVEL
35     INTEGER myThid
36    
37 cnh 1.5 C !LOCAL VARIABLES:
38 adcroft 1.2 C == Local variables ==
39     C bi,bj - Loop counters
40     C I,J,K
41     INTEGER bi, bj
42     INTEGER I, J, K
43 cnh 1.5 CEOP
44 adcroft 1.2
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 heimbach 1.3 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 adcroft 1.2
69     #ifdef INCLUDE_CD_CODE
70 heimbach 1.3 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 adcroft 1.2 #endif
77    
78     #ifdef ALLOW_NONHYDROSTATIC
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     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 heimbach 1.3 etaN(I,J,bi,bj)=0. _d 0
91 jmc 1.4 etaH(I,J,bi,bj)=0. _d 0
92     #ifdef INCLUDE_CD_CODE
93 heimbach 1.3 etaNm1(I,J,bi,bj)=0. _d 0
94 jmc 1.4 #endif
95 mlosch 1.7 phiHydLow(I,J,bi,bj)=0. _d 0
96 adcroft 1.2 ENDDO
97     ENDDO
98    
99     ENDDO
100     ENDDO
101    
102     RETURN
103     END

  ViewVC Help
Powered by ViewVC 1.1.22