Changeset 758 for trunk/LMDZ.MARS/libf

Timestamp:

Aug 8, 2012, 9:00:35 AM (12 years ago)

Author:

emillour

Message:

Mars GCM:

Some minor changes and adaptations for the MCDv5 runs:

• removed the 'nice' instruction from run0 script
• update dissipation factors in inidissip.F (fac_mid=3,fac_up=30,startalt=70)
• add some outputs in physiq.F
• put example def files used for MCDv5 simulations in deftank: callphys.def.MCD5 run.def.64x48x49.MCD5 traceur.def.MCD5

EM

Location:

trunk/LMDZ.MARS/libf

Files:

• dyn3d/inidissip.F (modified) (1 diff)
• phymars/physiq.F (modified) (6 diffs)

trunk/LMDZ.MARS/libf/dyn3d/inidissip.F

@@ -51 +51 @@
 c         FF 2004       
          if (pseudoalt(llm).lt.160.) then 
-             fac_mid=2  ! coeff  pour dissipation aux basses/moyennes altitudes
+             fac_mid=3  ! coeff  pour dissipation aux basses/moyennes altitudes
              fac_up=30 ! coeff multiplicateur pour dissipation hautes altitudes
              startalt=70. ! altitude en Km de la transition mid / up
              delta=20.! Intervalle (km) pour le changement mid / up
          else ! thermosphere model
-             fac_mid=2 ! coeff pour dissipation aux basses/moyennes altitudes
-             fac_up=25 ! coeff multiplicateur pour dissipation hautes altitudes
+             fac_mid=3 ! coeff pour dissipation aux basses/moyennes altitudes
+             fac_up=30 ! coeff multiplicateur pour dissipation hautes altitudes
 c            startalt: 95 OK for MY24
-             startalt=95. ! altitude en Km de la transition mid / up
+             startalt=70. ! altitude en Km de la transition mid / up
              delta=30.! Intervalle (km) pour le changement mid /up
          end if

trunk/LMDZ.MARS/libf/phymars/physiq.F

@@ -1585 +1585 @@
         call wstats(ngrid,"co2ice","CO2 ice cover",
      &                "kg.m-2",2,co2ice)
+        call wstats(ngrid,"tauref","reference dod at 610 Pa","NU",
+     &                2,tauref)
         call wstats(ngrid,"fluxsurf_lw",
      &                "Thermal IR radiative flux to surface","W.m-2",2,
@@ -1605 +1607 @@
         call wstats(ngrid,"rho","Atmospheric density","kg/m3",3,rho)
         call wstats(ngrid,"pressure","Pressure","Pa",3,pplay)
-c          call wstats(ngrid,"q2",
-c    &                "Boundary layer eddy kinetic energy",
-c    &                "m2.s-2",3,q2)
-c          call wstats(ngrid,"emis","Surface emissivity","w.m-1",2,
-c    &                emis)
+          call wstats(ngrid,"q2",
+     &                "Boundary layer eddy kinetic energy",
+     &                "m2.s-2",3,q2)
+          call wstats(ngrid,"emis","Surface emissivity","w.m-1",2,
+     &                emis)
 c          call wstats(ngrid,"ssurf","Surface stress","N.m-2",
 c    &                2,zstress)
@@ -1616 +1618 @@
 c          call wstats(ngrid,"lw_htrt","lw heat.rate",
 c    &                 "W.m-2",3,zdtlw)
+
+          if (calltherm) then
+            call wstats(ngrid,"zmax_th","Height of thermals",
+     &                "m",2,zmax_th)
+            call wstats(ngrid,"hfmax_th","Max thermals heat flux",
+     &                "K.m/s",2,hfmax_th)
+            call wstats(ngrid,"wstar",
+     &                "Max vertical velocity in thermals",
+     &                "m/s",2,wstar)
+          endif
 
            if (tracer) then
@@ -1637 +1649 @@
      &                    "surface h2o_ice","kg/m2",
      &                    2,qsurf(1,igcm_h2o_ice))
-c               call wstats(ngrid,'albedo',
-c     &                         'albedo',
-c     &                         '',2,albedo(1:ngridmx,1))
+               call wstats(ngrid,'albedo',
+     &                         'albedo',
+     &                         '',2,albedo(1:ngridmx,1))
                call wstats(ngrid,"mtot",
      &                    "total mass of water vapor","kg/m2",
@@ -1847 +1859 @@
 c        call WRITEDIAGFI(ngrid,"q2","q2","kg.m-3",3,q2)
 c        call WRITEDIAGFI(ngrid,'Teta','T potentielle','K',3,zh)
-         call WRITEDIAGFI(ngrid,"pressure","Pressure","Pa",3,pplay)
+         call WRITEDIAGFI(ngrid,"pressure","Pressure","Pa",3,zplay)
 c        call WRITEDIAGFI(ngrid,"ssurf","Surface stress","N.m-2",2,
 c    &                  zstress)
@@ -1980 +1992 @@
 c        ----------------------------------------------------------
 
-c        call WRITEDIAGFI(ngridmx,'tauref',
-c    &                    'Dust ref opt depth','NU',2,tauref)
+        call WRITEDIAGFI(ngridmx,'tauref',
+     &                    'Dust ref opt depth','NU',2,tauref)
 
          if (tracer.and.(dustbin.ne.0)) then