source: LMDZ.3.3/trunk/libf/dyn3d/flumass_scal.F @ 242

Last change on this file since 242 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.5 KB
Line 
1      SUBROUTINE flumass_scal(pbaru,pbarv,w,v)
2c%W%%G%
3c=======================================================================
4c
5c   Author:    Frederic Hourdin      original: 17/02/93
6c   -------
7c
8c   Subject:
9c   ------
10c   On ramene les flux de masse aux points scalaires.
11c
12c   Method:
13c   --------
14c
15c   Interface:
16c   ----------
17c
18c      Input:
19c      ------
20c
21c      Output:
22c      -------
23c
24c=======================================================================
25      IMPLICIT NONE
26c-----------------------------------------------------------------------
27c   Declararations:
28c   ---------------
29
30#include "dimensions.h"
31#include "paramet.h"
32#include "comgeom.h"
33
34c   Arguments:
35c   ----------
36
37      REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm),w(ip1jmp1,llm)
38      REAL v(ip1jmp1,llm,3)
39
40c   Local:
41c   ------
42
43      INTEGER ij,l
44
45c-----------------------------------------------------------------------
46
47c   transport zonal:
48c   ----------------
49      DO l=1,llm
50         DO ij=2,ip1jmp1
51            v(ij,l,1)=.5*(pbaru(ij,l)+pbaru(ij-1,l))
52         ENDDO
53      ENDDO
54      CALL SCOPY(jjp1*llm,v(iip1,1,1),iip1,v(1,1,1),iip1)
55
56c   Transport meridien:
57c   -------------------
58      DO l=1,llm
59         DO ij=iip2,ip1jm
60            v(ij,l,2)=.5*(pbarv(ij,l)+pbarv(ij-iip1,l))
61         ENDDO
62         DO ij=1,iip1
63            v(ij,l,2)=0.
64            v(ip1jm+ij,l,2)=0.
65         ENDDO
66      ENDDO
67
68c   Transport vertical:
69c   -------------------
70      DO l=1,llm-1
71         DO ij=1,ip1jmp1
72            v(ij,l,3)=.5*(w(ij,l)+w(ij,l+1))
73         ENDDO
74      ENDDO
75      DO ij=1,ip1jmp1
76         v(ij,llm,3)=.5*w(ij,llm)
77      ENDDO
78
79      RETURN
80      END
Note: See TracBrowser for help on using the repository browser.