Changeset 2037 for LMDZ5/trunk/libf/dyn3d

Timestamp:

May 6, 2014, 4:56:20 PM (11 years ago)

Author:

lguez

Message:

PVteta computed by PVtheta was not used. Also there were a couple of
problems with PVtheta:

-- PVtheta calls tetalevel in phylmd, and interpolates at

Earth-specific values of potential temperature.

-- PVtheta calls tourabs, and the computation of rot in tourabs should

be modified (it is not correct when there is a zoom).

-- Even when there is no zoom, the computation of rot in tourabs

should probably be modified: directly combine vcov and ucov, as in
tourpot, instead of dividing by cv and cu.

File:

• LMDZ5/trunk/libf/dyn3d/calfis.F (modified) (3 diffs)

LMDZ5/trunk/libf/dyn3d/calfis.F

@@ -163 +163 @@
       REAL unskap, pksurcp
 c
-cIM diagnostique PVteta, Amip2
-      INTEGER,PARAMETER :: ntetaSTD=3
-      REAL,SAVE :: rtetaSTD(ntetaSTD)=(/350.,380.,405./) ! Earth-specific, beware !!
-      REAL PVteta(ngridmx,ntetaSTD)
-c
       REAL flxwfi(ngridmx,llm)  ! Flux de masse verticale sur la grille physiq
 c
@@ -431 +426 @@
 
       ENDDO
-c
-      if (planet_type=="earth") then
-#ifdef CPP_PHYS
-! PVtheta calls tetalevel, which is in the physics
-cIM calcul PV a teta=350, 380, 405K
-      CALL PVtheta(ngridmx,llm,pucov,pvcov,pteta,
-     $           ztfi,zplay,zplev,
-     $           ntetaSTD,rtetaSTD,PVteta)
-#endif
-      endif
 c
 c On change de grille, dynamique vers physiq, pour le flux de masse verticale
@@ -491 +476 @@
      .             zdqfi,
      .             zdpsrf,
-cIM diagnostique PVteta, Amip2          
-     .             pducov,
-     .             PVteta)
+     .             pducov)
 
       else if ( planet_type=="generic" ) then