Changeset 575 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Mar 12, 2012, 2:10:08 PM (14 years ago)

Author:

emillour

Message:

Mars GCM:

correction in nirco2abs: io and ico2 should be declared as integers
Update sponge: remove posibility of specifying 'hsponge', all modes

apply to the last upper "nsponge" layers (default nsponge=3) and
sponge time scale doubles with every level, starting from top.

EM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• nirco2abs.F (modified) (1 diff)
• physiq.F (modified) (1 diff)

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

@@ -80 +80 @@
 c   ---------------------
       logical,save :: firstcall=.true.
-      real,save :: ico2=0 ! index of "co2" tracer
-      real,save :: io=0 ! index of "o" tracer
+      integer,save :: ico2=0 ! index of "co2" tracer
+      integer,save :: io=0 ! index of "o" tracer
 c     p0noonlte is a pressure below which non LTE effects are significant.
 c     REAL p0nonlte

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

@@ -1521 +1521 @@
            call wstats(ngrid,"ps","Surface pressure","Pa",2,ps)
            call wstats(ngrid,"tsurf","Surface temperature","K",2,tsurf)
-c           call wstats(ngrid,"co2ice","CO2 ice cover",
-c     &                "kg.m-2",2,co2ice)
-c           call wstats(ngrid,"fluxsurf_lw",
-c     &                "Thermal IR radiative flux to surface","W.m-2",2,
-c     &                fluxsurf_lw)
-c           call wstats(ngrid,"fluxsurf_sw",
-c     &                "Solar radiative flux to surface","W.m-2",2,
-c     &                fluxsurf_sw_tot)
-c           call wstats(ngrid,"fluxtop_lw",
-c     &                "Thermal IR radiative flux to space","W.m-2",2,
-c     &                fluxtop_lw)
-c           call wstats(ngrid,"fluxtop_sw",
-c     &                "Solar radiative flux to space","W.m-2",2,
-c     &                fluxtop_sw_tot)
-c           call wstats(ngrid,"temp","Atmospheric temperature","K",3,zt)
-c           call wstats(ngrid,"u","Zonal (East-West) wind","m.s-1",3,zu)
-c           call wstats(ngrid,"v","Meridional (North-South) wind",
-c     &                "m.s-1",3,zv)
+           call wstats(ngrid,"co2ice","CO2 ice cover",
+     &                "kg.m-2",2,co2ice)
+           call wstats(ngrid,"fluxsurf_lw",
+     &                "Thermal IR radiative flux to surface","W.m-2",2,
+     &                fluxsurf_lw)
+           call wstats(ngrid,"fluxsurf_sw",
+     &                "Solar radiative flux to surface","W.m-2",2,
+     &                fluxsurf_sw_tot)
+           call wstats(ngrid,"fluxtop_lw",
+     &                "Thermal IR radiative flux to space","W.m-2",2,
+     &                fluxtop_lw)
+           call wstats(ngrid,"fluxtop_sw",
+     &                "Solar radiative flux to space","W.m-2",2,
+     &                fluxtop_sw_tot)
+           call wstats(ngrid,"temp","Atmospheric temperature","K",3,zt)
+           call wstats(ngrid,"u","Zonal (East-West) wind","m.s-1",3,zu)
+           call wstats(ngrid,"v","Meridional (North-South) wind",
+     &                "m.s-1",3,zv)
 c           call wstats(ngrid,"w","Vertical (down-up) wind",
 c     &                "m.s-1",3,pw)
-c           call wstats(ngrid,"rho","Atmospheric density","none",3,rho)
+           call wstats(ngrid,"rho","Atmospheric density","kg/m3",3,rho)
 c           call wstats(ngrid,"pressure","Pressure","Pa",3,pplay)
 c          call wstats(ngrid,"q2",