1 |
|
1778 |
SUBROUTINE massbarxy(masse,massebxy) |
2 |
|
|
! |
3 |
|
|
!------------------------------------------------------------------------------- |
4 |
|
|
! Authors: P. Le Van , Fr. Hourdin. |
5 |
|
|
!------------------------------------------------------------------------------- |
6 |
|
|
! Purpose: Compute air mass mean along X and Y in each cell. |
7 |
|
|
! See iniconst for more details. |
8 |
|
|
IMPLICIT NONE |
9 |
|
|
include "dimensions.h" |
10 |
|
|
include "paramet.h" |
11 |
|
|
include "comgeom.h" |
12 |
|
|
!=============================================================================== |
13 |
|
|
! Arguments: |
14 |
|
|
REAL, INTENT(IN) :: masse (ip1jmp1,llm) |
15 |
|
|
REAL, INTENT(OUT) :: massebxy(ip1jm ,llm) |
16 |
|
|
!=============================================================================== |
17 |
|
|
! Local variables: |
18 |
|
|
INTEGER :: ij, l |
19 |
|
|
!=============================================================================== |
20 |
✓✓ |
71120 |
DO l=1,llm |
21 |
✓✓ |
73225152 |
DO ij=1,ip1jm-1 |
22 |
|
|
massebxy(ij,l)=masse(ij ,l)*alpha2(ij ) + & |
23 |
|
|
+ masse(ij+1 ,l)*alpha3(ij+1 ) + & |
24 |
|
|
+ masse(ij+iip1,l)*alpha1(ij+iip1) + & |
25 |
|
73225152 |
+ masse(ij+iip2,l)*alpha4(ij+iip2) |
26 |
|
|
END DO |
27 |
✓✓ |
2220722 |
DO ij=iip1,ip1jm,iip1; massebxy(ij,l)=massebxy(ij-iim,l); END DO |
28 |
|
|
END DO |
29 |
|
|
|
30 |
|
1778 |
END SUBROUTINE massbarxy |