/[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 76 by guez, Fri Nov 15 18:45:49 2013 UTC revision 78 by guez, Wed Feb 5 17:51:07 2014 UTC
# Line 9  contains Line 9  contains
9      ! From LMDZ4/libf/dyn3d/massdair.F, version 1.1.1.1 2004/05/19 12:53:07      ! From LMDZ4/libf/dyn3d/massdair.F, version 1.1.1.1 2004/05/19 12:53:07
10    
11      ! Calcule la masse d'air dans chaque maille.      ! Calcule la masse d'air dans chaque maille.
12      ! Auteurs : P. Le Van , F. Hourdin.      ! Auteurs : P. Le Van, F. Hourdin
   
     ! Méthode pour calculer massebx et masseby. A chaque point  
     ! scalaire P(i, j) sont affectés quatre coefficients d'aires.  
   
     ! alpha1(i, j) calculé au point (i + 1/4, j - 1/4)  
     ! alpha2(i, j) calculé au point (i + 1/4, j + 1/4)  
     ! alpha3(i, j) calculé au point (i - 1/4, j + 1/4)  
     ! alpha4(i, j) calculé au point (i - 1/4, j - 1/4)  
   
     ! Avec alpha1(i, j) = aire(i + 1/4, j - 1/4)/ aire(i, j)  
   
     ! Pour plus de détails, voir sous-programme "iniconst" et  
     ! "massdair.txt".  
13    
14      USE comgeom, ONLY: airesurg      USE comgeom, ONLY: airesurg
15      USE dimens_m, ONLY: iim, llm      USE dimens_m, ONLY: iim, llm
# Line 36  contains Line 23  contains
23    
24      !----------------------------------------------------------      !----------------------------------------------------------
25    
26      DO l = 1 , llm      DO l = 1, llm
27         masse(:, l) = airesurg(:) * (p(:, l) - p(:, l + 1))         masse(:, l) = airesurg * (p(:, l) - p(:, l + 1))
28    
29         DO ij = 1, ip1jmp1, iip1         DO ij = 1, ip1jmp1, iip1
30            masse(ij + iim, l) = masse(ij, l)            masse(ij + iim, l) = masse(ij, l)
Legend:
Removed from v.76  
changed lines
  Added in v.78
  ViewVC Help
Powered by ViewVC 1.1.21