46 |
!* 0.2 LOCAL ARRAYS |
!* 0.2 LOCAL ARRAYS |
47 |
! ------------ |
! ------------ |
48 |
LOGICAL lifthigh |
LOGICAL lifthigh |
49 |
INTEGER klevm1, jl, ilevh, jk |
INTEGER jl, jk |
50 |
REAL zcons1, ztmst, zrtmst, zpi, zhgeo |
REAL zcons1, ztmst, zpi, zhgeo |
51 |
REAL zdelp, zslow, zsqua, zscav, zbet |
REAL zdelp, zslow, zsqua, zscav, zbet |
52 |
INTEGER iknub(klon), iknul(klon) |
INTEGER iknub(klon), iknul(klon) |
53 |
LOGICAL ll1(klon,klev+1) |
LOGICAL ll1(klon,klev+1) |
64 |
|
|
65 |
IF (nlon/=klon .OR. nlev/=klev) STOP |
IF (nlon/=klon .OR. nlev/=klev) STOP |
66 |
zcons1 = 1./rd |
zcons1 = 1./rd |
|
klevm1 = klev - 1 |
|
67 |
ztmst = ptsphy |
ztmst = ptsphy |
|
zrtmst = 1./ztmst |
|
68 |
zpi = acos(-1.) |
zpi = acos(-1.) |
69 |
|
|
70 |
DO 1001 jl = 1, klon |
DO 1001 jl = 1, klon |
73 |
pvlow(jl) = 0.0 |
pvlow(jl) = 0.0 |
74 |
iknub(jl) = klev |
iknub(jl) = klev |
75 |
iknul(jl) = klev |
iknul(jl) = klev |
|
ilevh = klev/3 |
|
76 |
ll1(jl,klev+1) = .FALSE. |
ll1(jl,klev+1) = .FALSE. |
77 |
DO 1000 jk = 1, klev |
DO 1000 jk = 1, klev |
78 |
pvom(jl,jk) = 0.0 |
pvom(jl,jk) = 0.0 |