Changeset 1452 for trunk/LMDZ.VENUS/libf/phyvenus/ini_histmth.h

Timestamp:

Jun 16, 2015, 1:26:46 PM (9 years ago)

Author:

slebonnois

Message:

SL: converging photochemistry/high atm, with bug corrections and rho output

File:

• trunk/LMDZ.VENUS/libf/phyvenus/ini_histmth.h (modified) (5 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/ini_histmth.h

@@ -101 +101 @@
 c    .                "ave(X)", zsto,zout)
 c
-c plusieurs traceurs
-          if (iflag_trac.eq.1) THEN
+         CALL histdef(nid_mth, "mmean", "Mean molecular mass", "g/mol",
+     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+
+         CALL histdef(nid_mth, "rho", "Air density [mass/Vol]", "kg/m3",
+     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+
+
+c plusieurs traceurs  !!!outputs in [vmr]
+       IF (iflag_trac.eq.1) THEN
             DO iq=1,nqmax
              IF (iq.LE.99) THEN
           WRITE(str2,'(i2.2)') iq
-          CALL histdef(nid_mth, tname(iq), ttext(iq), "ppm",
+          CALL histdef(nid_mth, tname(iq), ttext(iq), "mol/mol",
      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
      .                "ave(X)", zsto,zout)
@@ -114 +123 @@
              ENDIF 
             ENDDO
-          endif
+       ENDIF
+
+       IF (callthermos .and. ok_chem) THEN
+          CALL histdef(nid_mth, "d_qmoldifCO2", "Dif molec" , "kg/kg",
+     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+          CALL histdef(nid_mth, "d_qmoldifO3p", "Dif molec" , "kg/kg",
+     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+          CALL histdef(nid_mth, "d_qmoldifN2", "Dif molec" , "kg/kg",
+     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+       ENDIF
 c
          CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2",
@@ -192 +213 @@
      .                "ave(X)", zsto,zout)
 c
-         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
-     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
-     .                "ave(X)", zsto,zout)
+         CALL histdef(nid_mth,"dtswrNLTE", "SWNLTE radiation dT", 
+     .                "K/s", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+c
+         CALL histdef(nid_mth,"dtswrDCrisp","SWDCrisp radiation dT",
+     .                "K/s", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+c
+         CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s",
+     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+c
+         CALL histdef(nid_mth,"dtlwrNLTE", "LWNLTE radiation dT",
+     .                "K/s", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+c
+         CALL histdef(nid_mth,"dtlwrLTE","LW_LTE radiation dT",
+     .                "K/s", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+c
+         CALL histdef(nid_mth, "dteuv", "UV radiation dT", "K/s",
+     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+         CALL histdef(nid_mth, "dtcond", "Therm conduction", "K/s",
+     .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+         CALL histdef(nid_mth, "dumolvis", "molec viscosity (u)"
+     .                ,"m/s2", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+         CALL histdef(nid_mth, "dvmolvis", "molec viscosity (v)", 
+     .                "m/s2", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+     .                "ave(X)", zsto,zout)
+
 c
 c        CALL histdef(nid_mth, "dtec", "Cinetic dissip dT", "K/s",
@@ -211 +262 @@
      .                iim,jj_nb,nhori, klev,1,klev,nvert, 32,
      .                "ave(X)", zsto,zout)
+
+c         CALL histdef(nid_mth, "co2_vmr", "volume mixture ratio co2", 
+c     .         "cm-3",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+c     .                "ave(X)", zsto,zout)
+
+c         CALL histdef(nid_mth, "o_vmr", "volume mixture ratio o",
+c     .                "cm-3", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+c     .                "ave(X)", zsto,zout)
+c         CALL histdef(nid_mth, "co_vmr", "volume mixture ratio co", 
+c     .         "cm-3",iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+c     .                "ave(X)", zsto,zout)
+
+c         CALL histdef(nid_mth, "n2_vmr", "volume mixture ratio n2",
+c     .                "cm-3", iim,jj_nb,nhori, klev,1,klev,nvert, 32,
+c     .                "ave(X)", zsto,zout)
 c
 c        CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2",
@@ -246 +312 @@
 
       ENDIF ! fin de test sur ok_mensuel
+