/[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.7 - (show 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 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
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 phiHydLow(I,J,bi,bj)=0. _d 0
96 ENDDO
97 ENDDO
98
99 ENDDO
100 ENDDO
101
102 RETURN
103 END

  ViewVC Help
Powered by ViewVC 1.1.22