3 |
! |
! |
4 |
SUBROUTINE integrd |
SUBROUTINE integrd |
5 |
$ ( nq,vcovm1,ucovm1,tetam1,psm1,massem1, |
$ ( nq,vcovm1,ucovm1,tetam1,psm1,massem1, |
6 |
$ dv,du,dteta,dq,dp,vcov,ucov,teta,q,ps,masse,phis,finvmaold ) |
$ dv,du,dteta,dq,dp,vcov,ucov,teta,q,ps,masse,phis,finvmaold, |
7 |
|
$ leapf ) |
8 |
|
|
9 |
use dimens_m |
use dimens_m |
10 |
use paramet_m |
use paramet_m |
51 |
REAL dv(ip1jm,llm),du(ip1jmp1,llm) |
REAL dv(ip1jm,llm),du(ip1jmp1,llm) |
52 |
REAL dteta(ip1jmp1,llm),dp(ip1jmp1) |
REAL dteta(ip1jmp1,llm),dp(ip1jmp1) |
53 |
REAL dq(ip1jmp1,llm,nq), finvmaold(ip1jmp1,llm) |
REAL dq(ip1jmp1,llm,nq), finvmaold(ip1jmp1,llm) |
54 |
|
logical, intent(in):: leapf |
55 |
|
|
56 |
c Local: |
c Local: |
57 |
c ------ |
c ------ |
159 |
|
|
160 |
10 CONTINUE |
10 CONTINUE |
161 |
|
|
|
|
|
|
c |
|
|
c ....... integration de q ...... |
|
|
c |
|
|
c$$$ IF( iadv(1).NE.3.AND.iadv(2).NE.3 ) THEN |
|
|
c$$$c |
|
|
c$$$ IF( forward. OR . leapf ) THEN |
|
|
c$$$ DO iq = 1,2 |
|
|
c$$$ DO l = 1,llm |
|
|
c$$$ DO ij = 1,ip1jmp1 |
|
|
c$$$ q(ij,l,iq) = ( q(ij,l,iq)*finvmaold(ij,l) + dtvr *dq(ij,l,iq) )/ |
|
|
c$$$ $ finvmasse(ij,l) |
|
|
c$$$ ENDDO |
|
|
c$$$ ENDDO |
|
|
c$$$ ENDDO |
|
|
c$$$ ELSE |
|
|
c$$$ DO iq = 1,2 |
|
|
c$$$ DO l = 1,llm |
|
|
c$$$ DO ij = 1,ip1jmp1 |
|
|
c$$$ q( ij,l,iq ) = q( ij,l,iq ) * finvmaold(ij,l) / finvmasse(ij,l) |
|
|
c$$$ ENDDO |
|
|
c$$$ ENDDO |
|
|
c$$$ ENDDO |
|
|
c$$$ |
|
|
c$$$ END IF |
|
|
c$$$c |
|
|
c$$$ ENDIF |
|
|
|
|
162 |
DO l = 1, llm |
DO l = 1, llm |
163 |
DO ij = 1, ip1jmp1 |
DO ij = 1, ip1jmp1 |
164 |
deltap(ij,l) = p(ij,l) - p(ij,l+1) |
deltap(ij,l) = p(ij,l) - p(ij,l+1) |