7 |
SUBROUTINE ppm3d(igd, q, ps1, ps2, u, v, w, ndt, iord, jord, kord, nc, imr, & |
SUBROUTINE ppm3d(igd, q, ps1, ps2, u, v, w, ndt, iord, jord, kord, nc, imr, & |
8 |
jnp, j1, nlay, ap, bp, pt, ae, fill, umax) |
jnp, j1, nlay, ap, bp, pt, ae, fill, umax) |
9 |
|
|
10 |
! rajout de déclarations |
! rajout de d\'eclarations |
11 |
! integer Jmax,kmax,ndt0,nstep,k,j,i,ic,l,js,jn,imh,iad,jad,krd |
! integer Jmax,kmax,ndt0,nstep,k,j,i,ic,l,js,jn,imh,iad,jad,krd |
12 |
! integer iu,iiu,j2,jmr,js0,jt |
! integer iu,iiu,j2,jmr,js0,jt |
13 |
! real dtdy,dtdy5,rcap,iml,jn0,imjm,pi,dl,dp |
! real dtdy,dtdy5,rcap,iml,jn0,imjm,pi,dl,dp |
224 |
INTEGER igd, iord, jord, kord, nc, jnp, j1, nlay, ae |
INTEGER igd, iord, jord, kord, nc, jnp, j1, nlay, ae |
225 |
REAL q(imr, jnp, nlay, nc), ps1(imr, jnp), ps2(imr, jnp), & |
REAL q(imr, jnp, nlay, nc), ps1(imr, jnp), ps2(imr, jnp), & |
226 |
u(imr, jnp, nlay), v(imr, jnp, nlay), ap(nlay+1), bp(nlay+1), & |
u(imr, jnp, nlay), v(imr, jnp, nlay), ap(nlay+1), bp(nlay+1), & |
227 |
w(imr, jnp, nlay), ndt, val(nlay), umax |
w(imr, jnp, nlay), ndt, umax |
228 |
INTEGER imrd2 |
INTEGER imrd2 |
229 |
REAL pt |
REAL pt |
230 |
LOGICAL cross, fill |
LOGICAL cross, fill |
706 |
IF (j1/=2) THEN |
IF (j1/=2) THEN |
707 |
DO k = 1, nlay |
DO k = 1, nlay |
708 |
DO i = 1, imr |
DO i = 1, imr |
709 |
! j=1 c'est le pôle Sud, j=JNP c'est le pôle Nord |
! j=1 c'est le p\^ole Sud, j=JNP c'est le p\^ole Nord |
710 |
q(i, 2, k, ic) = q(i, 1, k, ic) |
q(i, 2, k, ic) = q(i, 1, k, ic) |
711 |
q(i, jmr, k, ic) = q(i, jmp, k, ic) |
q(i, jmr, k, ic) = q(i, jmp, k, ic) |
712 |
END DO |
END DO |