source: trunk/LMDZ.COMMON/libf/dynphy_lonlat/phyvenus/wind_scal.F @ 3705

Last change on this file since 3705 was 1525, checked in by emillour, 10 years ago

All GCMs:
More on enforcing dynamics/physics separation: get rid of references to "control_mod" from physics packages.
EM
File size: 1.1 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
18c   Arguments:
19c   ----------
20
21      REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
22      REAL us(ip1jmp1,llm), vs (ip1jmp1,llm)
23
24c   Local:
25c   ------
26
27      INTEGER ij,l
28
29c-----------------------------------------------------------------------
30
31c   transport zonal:
32c   ----------------
33      DO l=1,llm
34         DO ij=2,ip1jmp1
35            us(ij,l)=.5*(pbaru(ij,l)+pbaru(ij-1,l))
36         ENDDO
37      ENDDO
38      CALL SCOPY(jjp1*llm,us(iip1,1),iip1,us(1,1),iip1)
39
40
41c   Transport meridien:
42c   -------------------
43      DO l=1,llm
44         DO ij=iip2,ip1jm
45            vs(ij,l)=.5*(pbarv(ij,l)+pbarv(ij-iip1,l))
46         ENDDO
47         DO ij=1,iip1
48            vs(ij,l)=0.
49            vs(ip1jm+ij,l)=0.
50         ENDDO
51      ENDDO
52
53      RETURN
54      END
Note: See TracBrowser for help on using the repository browser.