1 |
jmc |
1.1 |
C $Header: $ |
2 |
|
|
C $Name: $ |
3 |
|
|
|
4 |
|
|
#ifdef ALLOW_AIM |
5 |
|
|
|
6 |
|
|
C-- COMMON /PHYCON/: Physical constants (initial. in INPHYS) |
7 |
|
|
C P0 = reference pressure |
8 |
|
|
C GG = gravity accel. |
9 |
|
|
C RD = gas constant for dry air |
10 |
|
|
C CP = specific heat at constant pressure |
11 |
|
|
C ALHC = latent heat of condensation |
12 |
|
|
C SBC = Stefan-Boltzmann constant |
13 |
|
|
COMMON /PHYCON/ P0, GG, RD, CP, ALHC, SBC |
14 |
|
|
_RL P0, GG, RD, CP, ALHC, SBC |
15 |
|
|
|
16 |
|
|
C-- COMMON /FSIGLT/: Functions of sigma and latitude (initial. in INPHYS) |
17 |
|
|
C (execpt SLAT & CLAT init. in aim_dyn2aim.F) |
18 |
|
|
C SIG = full-level sigma |
19 |
|
|
C SIGL = logarithm of full-level sigma |
20 |
|
|
C SIGH = half-level sigma |
21 |
|
|
C DSIG = layer depth in sigma |
22 |
|
|
C POUT = norm. pressure level [p/p0] for post-processing |
23 |
|
|
C GRDSIG = g/(d_sigma p0) : to convert fluxes of u,v,q into d(u,v,q)/dt |
24 |
|
|
C GRDSCP = g/(d_sigma p0 c_p): to convert energy fluxes into dT/dt |
25 |
|
|
C WVI = weights for vertical interpolation |
26 |
|
|
c SLAT = sin(lat) |
27 |
|
|
c CLAT = cos(lat) |
28 |
|
|
COMMON /FSIGLT/ |
29 |
|
|
& SIG, SIGL, SIGH, DSIG, GRDSIG, GRDSCP, WVI |
30 |
|
|
c & SIG, SIGL, SIGH, DSIG, POUT, GRDSIG, GRDSCP, WVI, SLAT, CLAT |
31 |
|
|
_RL SIG(NLEV) |
32 |
|
|
_RL SIGL(NLEV) |
33 |
|
|
_RL SIGH(0:NLEV) |
34 |
|
|
_RL DSIG(NLEV) |
35 |
|
|
c _RL POUT(NLEV) |
36 |
|
|
_RL GRDSIG(NLEV) |
37 |
|
|
_RL GRDSCP(NLEV) |
38 |
|
|
_RL WVI(NLEV,2) |
39 |
|
|
c _RL SLAT(NLAT) |
40 |
|
|
c _RL CLAT(NLAT) |
41 |
|
|
|
42 |
|
|
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
43 |
|
|
#endif /* ALLOW_AIM */ |