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 |
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 |
19 |
USE paramet_m, ONLY : iip1, iip2, ip1jm, ip1jmp1, jjp1, llmp1 |
USE paramet_m, ONLY : iip1, iip2, ip1jm, llmp1 |
20 |
use qminimum_m, only: qminimum |
use qminimum_m, only: qminimum |
21 |
|
|
22 |
REAL vcovm1(ip1jm, llm), ucovm1((iim + 1) * (jjm + 1), llm) |
REAL vcovm1(ip1jm, llm), ucovm1((iim + 1) * (jjm + 1), llm) |
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 |