Ignore:
Timestamp:
Mar 20, 2014, 10:57:19 AM (10 years ago)
Author:
Laurent Fairhead
Message:

Merged trunk changes r1920:1997 into testing branch

Location:
LMDZ5/branches/testing
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/branches/testing

  • LMDZ5/branches/testing/libf/dyn3dpar/calfis_p.F

    r1910 r1999  
    3838      Use Write_field_p
    3939      USE Times
    40       USE infotrac
    41       USE control_mod
     40      USE infotrac, ONLY: nqtot, niadv, tname
     41      USE control_mod, ONLY: planet_type, nsplit_phys
    4242
    4343      IMPLICIT NONE
     
    112112c    Arguments :
    113113c    -----------
    114       LOGICAL  lafin
    115 !      REAL heure
    116       REAL, intent(in):: jD_cur, jH_cur
    117       REAL pvcov(iip1,jjm,llm)
    118       REAL pucov(iip1,jjp1,llm)
    119       REAL pteta(iip1,jjp1,llm)
    120       REAL pmasse(iip1,jjp1,llm)
    121       REAL pq(iip1,jjp1,llm,nqtot)
    122       REAL pphis(iip1,jjp1)
    123       REAL pphi(iip1,jjp1,llm)
    124 c
    125       REAL pdvcov(iip1,jjm,llm)
    126       REAL pducov(iip1,jjp1,llm)
    127       REAL pdteta(iip1,jjp1,llm)
    128       REAL pdq(iip1,jjp1,llm,nqtot)
    129       REAL flxw(iip1,jjp1,llm)  ! Flux de masse verticale sur la grille dynamique
    130 c
    131       REAL pps(iip1,jjp1)
    132       REAL pp(iip1,jjp1,llmp1)
    133       REAL ppk(iip1,jjp1,llm)
    134 c
    135       REAL pdvfi(iip1,jjm,llm)
    136       REAL pdufi(iip1,jjp1,llm)
    137       REAL pdhfi(iip1,jjp1,llm)
    138       REAL pdqfi(iip1,jjp1,llm,nqtot)
    139       REAL pdpsfi(iip1,jjp1)
    140 
    141       INTEGER        longcles
    142       PARAMETER    ( longcles = 20 )
    143       REAL clesphy0( longcles )
     114      LOGICAL,INTENT(IN) ::  lafin ! .true. for the very last call to physics
     115      REAL,INTENT(IN) :: jD_cur, jH_cur
     116      REAL,INTENT(IN) :: pvcov(iip1,jjm,llm) ! covariant meridional velocity
     117      REAL,INTENT(IN) :: pucov(iip1,jjp1,llm) ! covariant zonal velocity
     118      REAL,INTENT(IN) :: pteta(iip1,jjp1,llm) ! potential temperature
     119      REAL,INTENT(IN) :: pmasse(iip1,jjp1,llm) ! mass in each cell ! not used
     120      REAL,INTENT(IN) :: pq(iip1,jjp1,llm,nqtot) ! tracers
     121      REAL,INTENT(IN) :: pphis(iip1,jjp1) ! surface geopotential
     122      REAL,INTENT(IN) :: pphi(iip1,jjp1,llm) ! geopotential
     123
     124      REAL,INTENT(IN) :: pdvcov(iip1,jjm,llm) ! dynamical tendency on vcov ! not used
     125      REAL,INTENT(IN) :: pducov(iip1,jjp1,llm) ! dynamical tendency on ucov
     126      REAL,INTENT(IN) :: pdteta(iip1,jjp1,llm) ! dynamical tendency on teta
     127      ! NB: pdteta is used only to compute pcvgt which is in fact not used...
     128      REAL,INTENT(IN) :: pdq(iip1,jjp1,llm,nqtot) ! dynamical tendency on tracers
     129      ! NB: pdq is only used to compute pcvgq which is in fact not used...
     130
     131      REAL,INTENT(IN) :: pps(iip1,jjp1) ! surface pressure (Pa)
     132      REAL,INTENT(IN) :: pp(iip1,jjp1,llmp1) ! pressure at mesh interfaces (Pa)
     133      REAL,INTENT(IN) :: ppk(iip1,jjp1,llm) ! Exner at mid-layer
     134      REAL,INTENT(IN) :: flxw(iip1,jjp1,llm)  ! Vertical mass flux on dynamics grid
     135
     136      ! tendencies (in */s) from the physics
     137      REAL,INTENT(OUT) :: pdvfi(iip1,jjm,llm) ! tendency on covariant meridional wind
     138      REAL,INTENT(OUT) :: pdufi(iip1,jjp1,llm) ! tendency on covariant zonal wind
     139      REAL,INTENT(OUT) :: pdhfi(iip1,jjp1,llm) ! tendency on potential temperature (K/s)
     140      REAL,INTENT(OUT) :: pdqfi(iip1,jjp1,llm,nqtot) ! tendency on tracers
     141      REAL,INTENT(OUT) :: pdpsfi(iip1,jjp1) ! tendency on surface pressure (Pa/s)
     142
     143      INTEGER,PARAMETER :: longcles = 20
     144      REAL,INTENT(IN) :: clesphy0( longcles ) ! unused
    144145
    145146#ifdef CPP_PHYS
     
    217218c
    218219cIM diagnostique PVteta, Amip2
    219       INTEGER ntetaSTD
    220       PARAMETER(ntetaSTD=3)
    221       REAL rtetaSTD(ntetaSTD)
    222       DATA rtetaSTD/350., 380., 405./ ! Earth-specific values, beware !!
     220      INTEGER,PARAMETER :: ntetaSTD=3
     221      REAL,SAVE :: rtetaSTD(ntetaSTD)=(/350.,380.,405./) ! Earth-specific, beware !!
    223222      REAL PVteta(klon,ntetaSTD)
    224223     
    225      
    226224      REAL SSUM
    227225
    228       LOGICAL firstcal, debut
    229       DATA firstcal/.true./
    230       SAVE firstcal,debut
     226      LOGICAL,SAVE :: firstcal=.true., debut=.true.
    231227c$OMP THREADPRIVATE(firstcal,debut)
    232228     
Note: See TracChangeset for help on using the changeset viewer.