source: trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/scal_wind.F @ 1647

Last change on this file since 1647 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 scal_wind(xus,xvs,xu,xv)
2c=======================================================================
3c
4c
5c   Subject:
6c   ------
7c On passe  les variable xus, xvs  aux points de vent u et v (xu et xv)
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 xu(iip1,jjp1,llm),xv(iip1,jjm,llm)
23      REAL xus(iip1,jjp1,llm), xvs (iip1,jjp1,llm)
24
25c   Local:
26c   ------
27
28      INTEGER i,j,l
29
30c-----------------------------------------------------------------------
31
32c   transport zonal:
33c   ----------------
34      DO l=1,llm
35        Do j=1,jjp1
36              DO i=1,iim
37            xu(i,j,l)=0.5*(xus(i,j,l)+xus(i+1,j,l))
38              ENDDO
39          xu(iip1,j,l)=xu(1,j,l)
40            ENDDO
41      ENDDO
42
43
44c   Transport meridien:
45c   -------------------
46      DO l=1,llm
47         DO j=1,jjm
48           do i=1 ,iip1
49                 xv(i,j,l)=.5*(xvs(i,j,l)+xvs(i,j+1,l))
50           end do
51             ENDDO
52          ENDDO
53
54      RETURN
55      END
Note: See TracBrowser for help on using the repository browser.