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

Diff of /trunk/dyn3d/massbar.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 91 by guez, Wed Mar 26 17:18:58 2014 UTC
# Line 1  Line 1 
1  SUBROUTINE massbar(  masse, massebx, masseby )  module massbar_m
2    
3    ! From LMDZ4/libf/dyn3d/massbar.F,v 1.1.1.1 2004/05/19 12:53:05    IMPLICIT NONE
4    
5    !  contains
   ! **********************************************************************  
   !  
   !  Calcule les moyennes en x et  y de la masse d'air dans chaque maille.  
   ! **********************************************************************  
   !    Auteurs : P. Le Van , Fr. Hourdin  .  
   !   ..........  
   !  
   !  ..  masse                 est  un argum. d'entree  pour le s-pg ...  
   !  ..  massebx,masseby      sont des argum. de sortie pour le s-pg ...  
   !      
   !  
   !     IMPLICIT NONE  
   !  
   use dimens_m  
   use paramet_m  
   use comconst  
   use comgeom  
   !  
   REAL, intent(in):: masse( ip1jmp1,llm )  
   real massebx( ip1jmp1,llm )  , masseby(   ip1jm,llm )  
   !  
   !  
   ! Méthode pour calculer massebx et masseby. A chaque point scalaire  
   ! P(i, j) sont affectés quatre coefficients d'aire.  
   
   ! 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  
   ! "massbar.txt".  
   
   
   DO       l = 1 , llm  
      !  
      DO  ij = 1, ip1jmp1 - 1  
         massebx(ij,l) =  masse( ij, l) * alpha1p2( ij  )     +  &  
              masse(ij+1, l) * alpha3p4(ij+1 )  
      ENDDO  
   
      !    .... correction pour massebx( iip1,j) .....  
      !    ...    massebx(iip1,j)= massebx(1,j) ...  
      !  
      !DIR$ IVDEP  
      DO  ij = iip1, ip1jmp1, iip1  
         massebx( ij,l ) = massebx( ij - iim,l )  
      ENDDO  
   
   
      DO  ij = 1,ip1jm  
         masseby( ij,l ) = masse(  ij   , l ) * alpha2p3(   ij    )  + &  
              masse(ij+iip1, l ) * alpha1p4( ij+iip1 )  
      ENDDO  
6    
7    end DO    SUBROUTINE massbar(masse, massebx, masseby)
8    
9  END SUBROUTINE massbar      ! From LMDZ4/libf/dyn3d/massbar.F, version 1.1.1.1, 2004/05/19 12:53:05
10        ! Authors: P. Le Van, F. Hourdin.
11    
12        ! Calcule les moyennes en x et y de la masse d'air dans chaque
13        ! maille. Cf. "inigeom.txt" et "massbar.txt".
14    
15        USE dimens_m, ONLY: iim, llm
16        USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1
17        USE comgeom, ONLY: alpha1p2, alpha1p4, alpha2p3, alpha3p4
18    
19        REAL, intent(in):: masse(ip1jmp1,llm)
20        real, intent(out):: massebx(ip1jmp1,llm), masseby(ip1jm,llm)
21    
22        ! Local:
23        integer l, ij
24    
25        !--------------------------------------------------------------
26    
27        DO l = 1, llm
28           DO ij = 1, ip1jmp1 - 1
29              massebx(ij,l) = masse(ij, l) * alpha1p2(ij) + &
30                   masse(ij+1, l) * alpha3p4(ij+1)
31           ENDDO
32    
33           ! correction pour massebx(iip1,j)
34           ! massebx(iip1,j)= massebx(1,j)
35           DO ij = iip1, ip1jmp1, iip1
36              massebx(ij,l) = massebx(ij - iim,l)
37           ENDDO
38    
39           DO ij = 1,ip1jm
40              masseby(ij,l) = masse(ij, l) * alpha2p3(ij) + &
41                   masse(ij+iip1, l) * alpha1p4(ij+iip1)
42           ENDDO
43        end DO
44    
45      END SUBROUTINE massbar
46    
47    end module massbar_m
Legend:
Removed from v.82  
changed lines
  Added in v.91
  ViewVC Help
Powered by ViewVC 1.1.21