4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, dudyn, dteta, & |
SUBROUTINE integrd(vcovm1, ucovm1, tetam1, psm1, massem1, dv, du, dteta, & |
8 |
dp, vcov, ucov, teta, q, ps, masse, dt, leapf) |
dp, vcov, ucov, teta, q, ps, masse, dt, leapf) |
9 |
|
|
10 |
! From dyn3d/integrd.F, version 1.1.1.1, 2004/05/19 12:53:05 |
! From dyn3d/integrd.F, version 1.1.1.1, 2004/05/19 12:53:05 |
12 |
! Objet: incrémentation des tendances dynamiques |
! Objet: incrémentation des tendances dynamiques |
13 |
|
|
14 |
USE comgeom, ONLY : aire, aire_2d, apoln, apols |
USE comgeom, ONLY : aire, aire_2d, apoln, apols |
15 |
USE dimens_m, ONLY : iim, jjm, llm |
USE dimensions, ONLY : iim, jjm, llm |
16 |
USE disvert_m, ONLY : ap, bp |
USE disvert_m, ONLY : ap, bp |
17 |
use massdair_m, only: massdair |
use massdair_m, only: massdair |
18 |
use nr_util, only: assert |
use nr_util, only: assert |
23 |
REAL, intent(inout):: tetam1(iim + 1, jjm + 1, llm) |
REAL, intent(inout):: tetam1(iim + 1, jjm + 1, llm) |
24 |
REAL, intent(inout):: psm1((iim + 1) * (jjm + 1)) |
REAL, intent(inout):: psm1((iim + 1) * (jjm + 1)) |
25 |
real, intent(inout):: massem1(iim + 1, jjm + 1, llm) |
real, intent(inout):: massem1(iim + 1, jjm + 1, llm) |
26 |
REAL, intent(in):: dv(ip1jm, llm), dudyn((iim + 1) * (jjm + 1), llm) |
REAL, intent(in):: dv(ip1jm, llm), du((iim + 1) * (jjm + 1), llm) |
27 |
REAL, intent(in):: dteta(iim + 1, jjm + 1, llm), dp((iim + 1) * (jjm + 1)) |
REAL, intent(in):: dteta(iim + 1, jjm + 1, llm), dp((iim + 1) * (jjm + 1)) |
28 |
REAL, intent(inout):: vcov(ip1jm, llm), ucov((iim + 1) * (jjm + 1), llm) |
REAL, intent(inout):: vcov(ip1jm, llm), ucov((iim + 1) * (jjm + 1), llm) |
29 |
real, intent(inout):: teta(iim + 1, jjm + 1, llm) |
real, intent(inout):: teta(iim + 1, jjm + 1, llm) |
92 |
DO l = 1, llm |
DO l = 1, llm |
93 |
DO ij = iip2, ip1jm |
DO ij = iip2, ip1jm |
94 |
uscr(ij) = ucov(ij, l) |
uscr(ij) = ucov(ij, l) |
95 |
ucov(ij, l) = ucovm1(ij, l) + dt * dudyn(ij, l) |
ucov(ij, l) = ucovm1(ij, l) + dt * du(ij, l) |
96 |
END DO |
END DO |
97 |
|
|
98 |
DO ij = 1, ip1jm |
DO ij = 1, ip1jm |