55 |
DOUBLE PRECISION zglayd(kdlon) |
DOUBLE PRECISION zglayd(kdlon) |
56 |
DOUBLE PRECISION zglayu(kdlon) |
DOUBLE PRECISION zglayu(kdlon) |
57 |
DOUBLE PRECISION ztt(kdlon, ntra) |
DOUBLE PRECISION ztt(kdlon, ntra) |
|
DOUBLE PRECISION ztt1(kdlon, ntra) |
|
|
DOUBLE PRECISION ztt2(kdlon, ntra) |
|
58 |
DOUBLE PRECISION zuu(kdlon, nua) |
DOUBLE PRECISION zuu(kdlon, nua) |
59 |
|
|
60 |
INTEGER jk, jl, ja, im12, ind, inu, ixu, jg |
INTEGER jk, jl, ja, im12, ind, inu, ixu, jg |
65 |
|
|
66 |
DOUBLE PRECISION wg1(2) |
DOUBLE PRECISION wg1(2) |
67 |
SAVE wg1 |
SAVE wg1 |
68 |
DATA (wg1(jk), jk=1, 2)/1.0, 1.0/ |
DATA (wg1(jk), jk=1, 2)/1d0, 1d0/ |
69 |
! ----------------------------------------------------------------------- |
! ----------------------------------------------------------------------- |
70 |
|
|
71 |
! * 1. INITIALIZATION |
! * 1. INITIALIZATION |
90 |
DO ja = 1, ntra |
DO ja = 1, ntra |
91 |
DO jl = 1, kdlon |
DO jl = 1, kdlon |
92 |
ztt(jl, ja) = 1.0 |
ztt(jl, ja) = 1.0 |
|
ztt1(jl, ja) = 1.0 |
|
|
ztt2(jl, ja) = 1.0 |
|
93 |
END DO |
END DO |
94 |
END DO |
END DO |
95 |
|
|
182 |
pcntrb(jl, jk+1, jk) = zglayu(jl) |
pcntrb(jl, jk+1, jk) = zglayu(jl) |
183 |
pcntrb(jl, jk, jk) = 0.0 |
pcntrb(jl, jk, jk) = 0.0 |
184 |
END DO |
END DO |
|
|
|
185 |
END DO |
END DO |
186 |
|
|
187 |
DO jk = 1, kflev |
DO jk = 1, kflev |
194 |
END DO |
END DO |
195 |
END DO |
END DO |
196 |
|
|
|
RETURN |
|
|
|
|
197 |
END SUBROUTINE lwvn |
END SUBROUTINE lwvn |
198 |
|
|
199 |
end module lwvn_m |
end module lwvn_m |