1 |
C $Header: /u/gcmpack/MITgcm/pkg/aim_v23/com_physvar.h,v 1.4 2010/01/21 00:10:04 jmc Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
#ifdef ALLOW_AIM |
5 |
|
6 |
C-- Note: Variables which do not need to stay in common block (local var) |
7 |
C are declare locally in each S/R that use them (commented with "cL"); |
8 |
C Some variables are not needed at all (commented with single "c"). |
9 |
|
10 |
C-- COMMON /PHYGR1/ : Model variables on gaussian grid (updated in PHYPAR) |
11 |
C UG1 = u-wind |
12 |
C VG1 = v-wind |
13 |
C TG1 = abs. temperature |
14 |
C QG1 = specific humidity (g/kg) |
15 |
C PHIG1 = geopotential |
16 |
C PSLG1 = log. of surface pressure |
17 |
C VsurfSq = Square of surface wind speed (grid position = as T,Q) |
18 |
c COMMON /PHYGR1/ TG1, QG1, Vsurfsq |
19 |
c _RL UG1 (NGP,NLEV) |
20 |
c _RL VG1 (NGP,NLEV) |
21 |
cL _RL TG1 (NGP,NLEV) |
22 |
cL _RL QG1 (NGP,NLEV) |
23 |
c _RL PHIG1 (NGP,NLEV) |
24 |
c _RL PSLG1 (NGP) |
25 |
cL _RL VsurfSq(NGP) |
26 |
|
27 |
C-- COMMON /PHYGR2/ : Diagnosed upper-air variables (updated in PHYPAR) |
28 |
C SE = dry static energy <- replaced by Pot.Temp. |
29 |
C RH = relative humidity |
30 |
C QSAT = saturation specific humidity (g/kg) |
31 |
c COMMON /PHYGR2/ SE, RH, QSAT |
32 |
COMMON /PHYGR2/ RH |
33 |
cL _RL SE (NGP,NLEV) |
34 |
_RL RH (NGP,NLEV,MAX_NO_THREADS) |
35 |
cL _RL QSAT (NGP,NLEV) |
36 |
|
37 |
|
38 |
C-- COMMON /PHYGR3/ : Diagnosed surface variables (updated in PHYPAR) |
39 |
C PSG = surface pressure (normalized) |
40 |
C TS = surface temperature |
41 |
C TSKIN = skin temperature |
42 |
C U0 = near-surface u-wind |
43 |
C V0 = near-surface v-wind |
44 |
C T0 = near-surface air temperature |
45 |
C Q0 = near-surface specific humidity (g/kg) |
46 |
C CLOUDC = total cloud cover (fraction) |
47 |
C CLTOP = norm. pressure at cloud top |
48 |
c COMMON /PHYGR3/ PSG, TS, TSKIN, |
49 |
COMMON /PHYGR3/ TS, TSKIN, T0, Q0, CLOUDC, CLTOP |
50 |
cL _RL PSG (NGP) |
51 |
_RL TS (NGP,MAX_NO_THREADS) |
52 |
_RL TSKIN (NGP,MAX_NO_THREADS) |
53 |
c _RL U0(NGP), V0(NGP) |
54 |
_RL T0 (NGP,MAX_NO_THREADS) |
55 |
_RL Q0 (NGP,MAX_NO_THREADS) |
56 |
_RL CLTOP (NGP,MAX_NO_THREADS) |
57 |
_RL CLOUDC(NGP,MAX_NO_THREADS) |
58 |
|
59 |
C-- COMMON /PHYTEN/ : Physical param. tendencies (updated in PHYPAR) |
60 |
C TT_CNV = temperature tendency due to convection |
61 |
C QT_CNV = sp. humidity tendency due to convection |
62 |
C TT_LSC = temperature tendency due to large-scale condensation |
63 |
C QT_LSC = sp. humidity tendency due to large-scale condensation |
64 |
C TT_RSW = temperature tendency due to short-wave radiation |
65 |
C TT_RLW = temperature tendency due to long-wave radiation |
66 |
C UT_PBL = u-wind tendency due to PBL and diffusive processes |
67 |
C VT_PBL = v-wind tendency due to PBL and diffusive processes |
68 |
C TT_PBL = temperature tendency due to PBL and diffusive processes |
69 |
C QT_PBL = sp. humidity tendency due to PBL and diffusive processes |
70 |
COMMON /PHYTEN/ TT_CNV, QT_CNV, TT_LSC, QT_LSC, |
71 |
& TT_RSW, TT_RLW, TT_PBL, QT_PBL |
72 |
_RL TT_CNV (NGP,NLEV,MAX_NO_THREADS) |
73 |
_RL QT_CNV (NGP,NLEV,MAX_NO_THREADS) |
74 |
_RL TT_LSC (NGP,NLEV,MAX_NO_THREADS) |
75 |
_RL QT_LSC (NGP,NLEV,MAX_NO_THREADS) |
76 |
_RL TT_RSW (NGP,NLEV,MAX_NO_THREADS) |
77 |
_RL TT_RLW (NGP,NLEV,MAX_NO_THREADS) |
78 |
c _RL UT_PBL (NGP,NLEV,MAX_NO_THREADS) |
79 |
c _RL VT_PBL (NGP,NLEV,MAX_NO_THREADS) |
80 |
_RL TT_PBL (NGP,NLEV,MAX_NO_THREADS) |
81 |
_RL QT_PBL (NGP,NLEV,MAX_NO_THREADS) |
82 |
|
83 |
C-- COMMON /FLUXES/ : Surface and upper boundary fluxes (updated in PHYPAR) |
84 |
C PRECNV = convective precipitation [g/m2/s] |
85 |
C PRECLS = large-scale precipitation [g/m2/s] |
86 |
C EnPrec = energy of precipitation (snow, rain temp) [J/g] |
87 |
C CBMF = cloud-base mass flux |
88 |
C TSR = top-of-atm. shortwave radiation (downward) |
89 |
C SSR = surf. shortwave radiation (downward) (1:land, 2:sea, 3:sea-ice) |
90 |
C SLR = surface longwave radiation (upward) (1:land, 2:sea, 3:sea-ice) |
91 |
C OLR = outgoing longwave radiation (upward) |
92 |
C USTR = u-stress (1:land, 2:sea, 3:weighted average) |
93 |
C VSTR = v-stress (1:land, 2:sea, 3:weighted average) |
94 |
C SHF = sensible heat flux (1:land, 2:sea, 3:sea-ice) |
95 |
C EVAP = evaporation [g/m2/s] (1:land, 2:sea, 3:sea-ice) |
96 |
C DRAG = surface Drag term (= Cd*Rho*|V|) (1:land, 2:sea, 3:sea-ice) |
97 |
COMMON /FLUXES/ PRECNV, PRECLS, EnPrec, |
98 |
& CBMF, TSR, SSR, SLR, OLR, |
99 |
& SHF, EVAP, SPEED0, DRAG |
100 |
_RL PRECNV (NGP,MAX_NO_THREADS) |
101 |
_RL PRECLS (NGP,MAX_NO_THREADS) |
102 |
_RL EnPrec (NGP,MAX_NO_THREADS) |
103 |
_RL CBMF (NGP,MAX_NO_THREADS) |
104 |
_RL TSR (NGP,MAX_NO_THREADS) |
105 |
_RL SSR (NGP,0:3,MAX_NO_THREADS) |
106 |
_RL SLR (NGP,0:3,MAX_NO_THREADS) |
107 |
_RL OLR (NGP,MAX_NO_THREADS) |
108 |
c _RL USTR (NGP,3) |
109 |
c _RL VSTR (NGP,3) |
110 |
_RL SHF (NGP,0:3,MAX_NO_THREADS) |
111 |
_RL EVAP (NGP,0:3,MAX_NO_THREADS) |
112 |
_RL SPEED0 (NGP,MAX_NO_THREADS) |
113 |
_RL DRAG (NGP,0:3,MAX_NO_THREADS) |
114 |
|
115 |
#ifdef ALLOW_CLR_SKY_DIAG |
116 |
C TT_SWclr = temp. tendency due to clear-sky short-wave radiation |
117 |
C TSWclr = top-of-atm. clear-sky shortwave radiation (downward) |
118 |
C SSWclr = clear-sky surf. (net) shortwave radiation (downward) |
119 |
C TT_LWclr = temp. tendency due to clear-sky long-wave radiation |
120 |
C OLWclr = clear-sky outgoing longwave radiation (upward) |
121 |
C SLWclr = clear-sky surf. (net) longwave radiation (upward) |
122 |
COMMON /CLRSKYDIAG/ |
123 |
& TT_SWclr, TSWclr, SSWclr, |
124 |
& TT_LWclr, OLWclr, SLWclr |
125 |
_RL TT_SWclr(NGP,NLEV,MAX_NO_THREADS) |
126 |
_RL TSWclr (NGP,MAX_NO_THREADS) |
127 |
_RL SSWclr (NGP,MAX_NO_THREADS) |
128 |
_RL TT_LWclr(NGP,NLEV,MAX_NO_THREADS) |
129 |
_RL OLWclr (NGP,MAX_NO_THREADS) |
130 |
_RL SLWclr (NGP,MAX_NO_THREADS) |
131 |
#endif |
132 |
|
133 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
134 |
#endif /* ALLOW_AIM */ |