Changeset 1356 for trunk/LMDZ.TITAN/libf/phytitan/physiq.F

Timestamp:

Oct 8, 2014, 9:26:28 AM (10 years ago)

Author:

slebonnois

Message:

SL: update to newstart/start2archive tools in Venus+Titan / additional diagnostics in radiative fluxes for Titan

File:

• trunk/LMDZ.TITAN/libf/phytitan/physiq.F (modified) (6 diffs)

trunk/LMDZ.TITAN/libf/phytitan/physiq.F

@@ -7 +7 @@
      .            paprs,pplay,ppk,pphi,pphis,presnivs,
      .            u,v,t,qx,
-     .            omega,
+     .            flxmw,
      .            d_u, d_v, d_t, d_qx, d_ps)
 
@@ -46 +46 @@
 c qx------input-R-mass mixing ratio traceurs (kg/kg) 
 c d_t_dyn-input-R-tendance dynamique pour "t" (K/s)
-c omega---input-R-vitesse verticale en Pa/s
+c flxmw---input-R-flux de masse vertical en kg/s
 c
 c d_u-----output-R-tendance physique de "u" (m/s/s)
@@ -133 +133 @@
       REAL d_t_dyn(klon,klev)
 
-      REAL omega(klon,klev)
+      REAL flxmw(klon,klev)
 
       REAL d_u(klon,klev)
@@ -146 +146 @@
       INTEGER,save :: itap        ! compteur pour la physique
       REAL delp(klon,klev)        ! epaisseur d'une couche
+      REAL omega(klon,klev)
       
       INTEGER igwd,idx(klon),itest(klon)
@@ -722 +723 @@
 c====================================================================
 c
+c Calcule de vitesse verticale a partir de flux de masse verticale
+      DO k = 1, klev
+       DO i = 1, klon
+        omega(i,k) = RG*flxmw(i,k) / airephy(i)
+       END DO
+      END DO
+
 c Ajouter le geopotentiel du sol:
 c
@@ -1316 +1324 @@
 c
 c A ADAPTER POUR VENUS!!!
-        CALL drag_noro(klon,klev,dtime,paprs,pplay,zphi,zn2,
+        CALL drag_noro(klon,klev,dtime,paprs,pplay,pphi,zn2,
      e                   zmea,zstd, zsig, zgam, zthe,zpic,zval,
      e                   igwd,idx,itest,