Changeset 2037 for LMDZ5/trunk/libf/phydev/physiq.F90

Timestamp:

May 6, 2014, 4:56:20 PM (10 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/phydev/physiq.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phydev/physiq.F90

@@ -8 +8 @@
      &            flxmass_w, &
      &            d_u, d_v, d_t, d_qx, d_ps &
-     &            , dudyn &
-     &            , PVteta)
+     &            , dudyn)
 
       USE dimphy, only : klon,klev
@@ -59 +58 @@
       real,intent(out) :: d_ps(klon) ! physics tendency on surface pressure
       real,intent(in) :: dudyn(iim+1,jjmp1,klev) ! Not used
-!FH! REAL PVteta(klon,nbteta)
-!      REAL PVteta(klon,1)
-      real,intent(in) :: PVteta(klon,3) ! Not used ; should match definition
-                                        ! in calfis.F
 
 integer,save :: itau=0 ! counter to count number of calls to physics