85 |
REAL zc(klon, nsoilmx, nbsrf), zd(klon, nsoilmx, nbsrf) |
REAL zc(klon, nsoilmx, nbsrf), zd(klon, nsoilmx, nbsrf) |
86 |
REAL lambda |
REAL lambda |
87 |
SAVE dz1, dz2, zc, zd, lambda |
SAVE dz1, dz2, zc, zd, lambda |
88 |
LOGICAL firstcall, firstsurf(nbsrf) |
LOGICAL firstsurf(nbsrf) |
89 |
SAVE firstcall, firstsurf |
SAVE firstsurf |
90 |
REAL isol, isno, iice |
REAL isol, isno, iice |
91 |
SAVE isol, isno, iice |
SAVE isol, isno, iice |
92 |
|
|
|
DATA firstcall/.TRUE./ |
|
93 |
DATA firstsurf/.TRUE., .TRUE., .TRUE., .TRUE./ |
DATA firstsurf/.TRUE., .TRUE., .TRUE., .TRUE./ |
94 |
|
|
95 |
DATA isol, isno, iice/2000., 2000., 2000./ |
DATA isol, isno, iice/2000., 2000., 2000./ |
130 |
STOP 1 |
STOP 1 |
131 |
END IF |
END IF |
132 |
|
|
|
|
|
|
! $$$ IF (firstcall) THEN |
|
133 |
IF (firstsurf(indice)) THEN |
IF (firstsurf(indice)) THEN |
134 |
|
|
135 |
! ----------------------------------------------------------------------- |
! ----------------------------------------------------------------------- |
172 |
END DO |
END DO |
173 |
! PB |
! PB |
174 |
firstsurf(indice) = .FALSE. |
firstsurf(indice) = .FALSE. |
|
! $$$ firstcall =.false. |
|
175 |
|
|
176 |
! Initialisations: |
! Initialisations: |
177 |
! ---------------- |
! ---------------- |
178 |
|
|
179 |
ELSE !--not firstcall |
ELSE |
180 |
! ----------------------------------------------------------------------- |
! ----------------------------------------------------------------------- |
181 |
! Computation of the soil temperatures using the Cgrd and Dgrd |
! Computation of the soil temperatures using the Cgrd and Dgrd |
182 |
! coefficient computed at the previous time-step: |
! coefficient computed at the previous time-step: |
196 |
END DO |
END DO |
197 |
END DO |
END DO |
198 |
|
|
199 |
END IF !--not firstcall |
END IF |
200 |
! ----------------------------------------------------------------------- |
! ----------------------------------------------------------------------- |
201 |
! Computation of the Cgrd and Dgrd coefficient for the next step: |
! Computation of the Cgrd and Dgrd coefficient for the next step: |
202 |
! --------------------------------------------------------------- |
! --------------------------------------------------------------- |