source: trunk/LMDZ.MARS/libf/dyn3d/wind_scal.F @ 937

Last change on this file since 937 was 38, checked in by emillour, 14 years ago

Ajout du modè Martien (mon LMDZ.MARS.BETA, du 28/01/2011) dans le rértoire mars, pour pouvoir suivre plus facilement les modifs.
EM
File size: 1.2 KB
Line 
1      SUBROUTINE wind_scal(pbaru,pbarv,us,vs)
2c=======================================================================
3c
4c
5c   Subject:
6c   ------
7c   On ramene les flux de masse /vents  aux points scalaires.
8c
9c=======================================================================
10      IMPLICIT NONE
11c-----------------------------------------------------------------------
12c   Declararations:
13c   ---------------
14
15#include "dimensions.h"
16#include "paramet.h"
17#include "comgeom.h"
18
19c   Arguments:
20c   ----------
21
22      REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
23      REAL us(ip1jmp1,llm), vs (ip1jmp1,llm)
24
25c   Local:
26c   ------
27
28      INTEGER ij,l
29
30c-----------------------------------------------------------------------
31
32c   transport zonal:
33c   ----------------
34      DO l=1,llm
35         DO ij=2,ip1jmp1
36            us(ij,l)=.5*(pbaru(ij,l)+pbaru(ij-1,l))
37         ENDDO
38      ENDDO
39      CALL SCOPY(jjp1*llm,us(iip1,1),iip1,us(1,1),iip1)
40
41
42c   Transport meridien:
43c   -------------------
44      DO l=1,llm
45         DO ij=iip2,ip1jm
46            vs(ij,l)=.5*(pbarv(ij,l)+pbarv(ij-iip1,l))
47         ENDDO
48         DO ij=1,iip1
49            vs(ij,l)=0.
50            vs(ip1jm+ij,l)=0.
51         ENDDO
52      ENDDO
53
54      RETURN
55      END
Note: See TracBrowser for help on using the repository browser.