/[lmdze]/trunk/dyn3d/massdair.f
ViewVC logotype

Diff of /trunk/dyn3d/massdair.f

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

revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC revision 102 by guez, Tue Jul 15 13:43:24 2014 UTC
# Line 11  contains Line 11  contains
11      ! Calcule la masse d'air dans chaque maille.      ! Calcule la masse d'air dans chaque maille.
12      ! Authors: P. Le Van, F. Hourdin      ! Authors: P. Le Van, F. Hourdin
13    
14      USE comgeom, ONLY: airesurg      USE comgeom, ONLY: airesurg_2d
15      USE dimens_m, ONLY: iim, llm      USE dimens_m, ONLY: iim, jjm, llm
     USE paramet_m, ONLY: iip1, ip1jmp1, llmp1  
16    
17      REAL, intent(in):: p(ip1jmp1, llmp1) ! aux interfaces des llm couches      REAL, intent(in):: p(iim + 1, jjm + 1, llm + 1)
18      real, intent(out):: masse(ip1jmp1, llm)      ! aux interfaces des llm couches
19    
20      ! Variables locales      real, intent(out):: masse(iim + 1, jjm + 1, llm)
     INTEGER l, ij  
21    
22      !----------------------------------------------------------      ! Local:
23        INTEGER l
24    
25      DO l = 1, llm      !----------------------------------------------------------
        masse(:, l) = airesurg * (p(:, l) - p(:, l + 1))  
26    
27         DO ij = 1, ip1jmp1, iip1      forall (l = 1: llm) masse(:iim, :, l) = airesurg_2d(:iim, :) &
28            masse(ij + iim, l) = masse(ij, l)           * (p(:iim, :, l) - p(:iim, :, l + 1))
29         ENDDO      masse(iim + 1, :, :) = masse(1, :, :)
     end DO  
30    
31    END SUBROUTINE massdair    END SUBROUTINE massdair
32    

Legend:
Removed from v.82  
changed lines
  Added in v.102

  ViewVC Help
Powered by ViewVC 1.1.21