98 |
logical leapf |
logical leapf |
99 |
real dt ! time step, in s |
real dt ! time step, in s |
100 |
|
|
101 |
REAL p3d(iim + 1, jjm + 1, llm+1) ! pressure at layer interfaces, in Pa |
REAL p3d(iim + 1, jjm + 1, llm + 1) ! pressure at layer interfaces, in Pa |
102 |
! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)", |
! ("p3d(i, j, l)" is at longitude "rlonv(i)", latitude "rlatu(j)", |
103 |
! for interface "l") |
! for interface "l") |
104 |
|
|
202 |
forall (l = 1: llm) |
forall (l = 1: llm) |
203 |
teta(:, 1, l) = SUM(aire_2d(:iim, 1) * teta(:iim, 1, l)) & |
teta(:, 1, l) = SUM(aire_2d(:iim, 1) * teta(:iim, 1, l)) & |
204 |
/ apoln |
/ apoln |
205 |
teta(:, jjm + 1, l) = SUM(aire_2d(:iim, jjm+1) & |
teta(:, jjm + 1, l) = SUM(aire_2d(:iim, jjm + 1) & |
206 |
* teta(:iim, jjm + 1, l)) / apols |
* teta(:iim, jjm + 1, l)) / apols |
207 |
END forall |
END forall |
208 |
END IF |
END IF |