/[lmdze]/trunk/libf/dyn3d/integrd.f90
ViewVC logotype

Diff of /trunk/libf/dyn3d/integrd.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 12 by guez, Mon Jul 21 16:05:07 2008 UTC
# Line 3  Line 3 
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
# Line 50  c   ---------- Line 51  c   ----------
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   ------
# Line 157  c Line 159  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)
Legend:
Removed from v.3  
changed lines
  Added in v.12
  ViewVC Help
Powered by ViewVC 1.1.21