/[MITgcm]/MITgcm/model/inc/DYNVARS.h
ViewVC logotype

Annotation of /MITgcm/model/inc/DYNVARS.h

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


Revision 1.16 - (hide annotations) (download)
Wed Sep 19 13:50:38 2001 UTC (22 years, 8 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint40
Changes since 1.15: +12 -4 lines
File MIME type: text/plain
add variable etaH for non-linear free-surface + Crank-Nickelson

1 jmc 1.16 C $Header: /u/gcmpack/models/MITgcmUV/model/inc/DYNVARS.h,v 1.15 2001/08/13 18:07:35 heimbach Exp $
2     C $Name: $
3 cnh 1.1 C
4     C /==========================================================\
5     C | DYNVARS.h |
6     C | o Dynamical model variables (common block DYNVARS_R) |
7     C |==========================================================|
8     C | The value and two levels of time tendency are held for |
9     C | each prognostic variable. |
10     C \==========================================================/
11     C
12 jmc 1.12 C etaN - free-surface r-anomaly (r unit) at current time level
13 cnh 1.1 C uVel - zonal velocity (m/s, i=1 held at western face)
14     C vVel - meridional velocity (m/s, j=1 held at southern face)
15     C theta - potential temperature (oC, held at pressure/tracer point)
16     C salt - salinity (ppt, held at pressure/tracer point)
17 cnh 1.4 C gX, gXNM1 - Time tendencies at current and prvious time levels.
18     C uVelD - D grid zonal velocity
19     C vVelD - D grid meridional velocity
20 jmc 1.16 C etaH - surface r-anomaly, advanced in time consistently
21     C with 2.D flow divergence (Exact-Conservation):
22     C etaH^n+1 = etaH^n - delta_t*Div.(H^n U^n)
23     C note: a) used with "exactConserv" but strictly necessary for NonLinFreeSurf
24     C b) same as etaN but not necessarely at the same time, e.g.:
25     C implicDiv2DFlow=0 => etaH=etaN ; =1 => etaH=etaNm1 ;
26 adcroft 1.9
27 adcroft 1.13 COMMON /DYNVARS_R/
28 jmc 1.16 & etaN, etaH,
29 adcroft 1.13 & uVel,vVel,wVel,theta,salt,
30 adcroft 1.8 & gu,gv,gt,gs,guNm1,gvNm1,gtNm1,gsNm1
31 adcroft 1.13 _RL etaN (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
32 jmc 1.16 _RL etaH (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
33 cnh 1.5 _RL uVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
34     _RL vVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
35 adcroft 1.10 _RL wVel (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
36 cnh 1.5 _RL theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
37     _RL salt (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
38     _RL gu(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
39     _RL gv(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
40     _RL gt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
41     _RL gs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
42     _RL guNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
43     _RL gvNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
44     _RL gtNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
45     _RL gsNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
46 heimbach 1.15
47 cnh 1.7 #ifdef INCLUDE_CD_CODE
48 cnh 1.4 COMMON /DYNVARS_CD/
49     & uVelD, vVelD,
50 jmc 1.16 & etaNm1,
51 cnh 1.4 & uNM1, vNM1,
52     & guCD, gvCD
53 cnh 1.5 _RL uVeld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
54     _RL vVeld (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
55 jmc 1.16 _RL etaNm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
56 cnh 1.5 _RL uNm1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
57     _RL vNm1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
58     _RL guCD (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
59     _RL gvCD (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
60 cnh 1.4 #endif
61 heimbach 1.15
62 adcroft 1.14 #ifdef ALLOW_NONHYDROSTATIC
63     COMMON /DYNVARS_NH/ phi_nh
64     _RL phi_nh(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
65     #endif /* ALLOW_NONHYDROSTATIC */
66 heimbach 1.15
67     #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
68     COMMON /DYNVARS_DIFFKR/
69     & diffKr
70     _RL diffKr (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
71     #endif
72     #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_KAPGM_CONTROL))
73     COMMON /DYNVARS_KAPGM/
74     & kapgm
75     _RL kapgm (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
76     #endif

  ViewVC Help
Powered by ViewVC 1.1.22