source: trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/wind_scal.F @ 1644

Last change on this file since 1644 was 1403, checked in by emillour, 10 years ago

All models: Reorganizing the physics/dynamics interface.

  • makelmdz and makelmdz_fcm scripts adapted to handle the new directory settings
  • misc: (replaces what was the "bibio" directory)
  • Should only contain extremely generic (and non physics or dynamics-specific) routines
  • Therefore moved initdynav.F90, initfluxsto.F, inithist.F, writedynav.F90, write_field.F90, writehist.F to "dyn3d_common"
  • dynlonlat_phylonlat: (new interface directory)
  • This directory contains routines relevent to physics/dynamics grid interactions, e.g. routines gr_dyn_fi or gr_fi_dyn and calfis
  • Moreover the dynlonlat_phylonlat contains directories "phy*" corresponding to each physics package "phy*" to be used. These subdirectories should only contain specific interfaces (e.g. iniphysiq) or main programs (e.g. newstart)
  • phy*/dyn1d: this subdirectory contains the 1D model using physics from phy*

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.