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 |
|
|