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

Timestamp:

Jul 26, 2013, 12:54:55 PM (12 years ago)

Author:

emillour

Message:

Mars GCM:

• Upgrade on the thermospheric photochemical reaction rates. These are now read in from a file "chemthermos_reactionrates.def".

FGG + JYC

File:

• trunk/LMDZ.MARS/libf/phymars/physiq.F (modified) (3 diffs)

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

@@ -470 +470 @@
          if(callnlte.and.nltemodel.eq.2) call nlte_setup
          if(callnirco2.and.nircorr.eq.1) call NIR_leedat
+         if(thermochem) call chemthermos_readini
 
         IF (tracer.AND.water.AND.(ngridmx.NE.1)) THEN
@@ -1727 +1728 @@
                call wstats(ngrid,'albedo',
      &                         'albedo',
-     &                         '',2,albedo(1:ngridmx,1))
+     &                         '',2,albedo(1,1))
                call wstats(ngrid,"mtot",
      &                    "total mass of water vapor","kg/m2",
@@ -2044 +2045 @@
             CALL WRITEDIAGFI(ngridmx,'albedo',
      &                         'albedo',
-     &                         '',2,albedo(1:ngridmx,1))
+     &                         '',2,albedo(1,1))
               if (tifeedback) then
                  call WRITEDIAGSOIL(ngridmx,"soiltemp",