Changeset 607

Timestamp:

Mar 30, 2012, 3:45:49 PM (13 years ago)

Author:

emillour

Message:

Mars GCM:

Updated model so that 1) reference pressure for opacity is now 610 Pa (and
not 700 Pa anymore) and 2) the new MY24-MY30 dust scenarios (input files
dust_MY24.nc, dust_MY25.nc, ..., dust_MY30.nc) can be used
(when setting iaervar=24,25,...,30 in callphys.def): changed
"readtesassim.F90" into "read_dust_scenario.F90" and apadpted aeropacity.F.
one can still use the "old" MY24-MY25-MY26 (input files dust_tes_MY24.nc,
dust_tes_MY25.nc and dust_tes_MY26.nc) scenarios by setting
iaervar=124, 125 or 126 in callphys.def.

EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/aeropacity.F (modified) (10 diffs)
• libf/phymars/read_dust_scenario.F90 (added)
• libf/phymars/readtesassim.F90 (deleted)

trunk/LMDZ.MARS/README

@@ -1578 +1578 @@
 >> Update vertical grid resolution for diffusion.
  
+== 30/03/12 == EM
+>> Updated model so that 1) reference pressure for opacity is now 610 Pa (and
+   not 700 Pa anymore) and 2) the new MY24-MY30 dust scenarios (input files
+   dust_MY24.nc, dust_MY25.nc, ..., dust_MY30.nc) can be used
+   (when setting iaervar=24,25,...,30 in callphys.def): changed
+   "readtesassim.F90" into "read_dust_scenario.F90" and apadpted aeropacity.F.
+   one can still use the "old" MY24-MY25-MY26 (input files dust_tes_MY24.nc,
+   dust_tes_MY25.nc and dust_tes_MY26.nc) scenarios by setting
+   iaervar=124, 125 or 126 in callphys.def.
+   

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

@@ -19 +19 @@
 c       - dustopacity transferred from physiq.F to callradite.F,
 c           and renamed into aeropacity.F;
+c   update E. Millour, march 2012:
+c         - reference pressure is now set to 610Pa (not 700Pa)
 c   
 c   input:
@@ -37 +39 @@
 c   output:
 c   -------
-c   tauref       Prescribed mean column optical depth at 700 Pa 
+c   tauref       Prescribed mean column optical depth at 610 Pa 
 c   tau          Column total visible dust optical depth at each point
 c   aerosol      aerosol(ig,l,1) is the dust optical
@@ -114 +116 @@
       INTEGER, PARAMETER :: cstdustlevel = 7
 
-      LOGICAL firstcall
-      DATA firstcall/.true./
-      SAVE firstcall
+      LOGICAL,SAVE :: firstcall=.true.
 
 ! indexes of water ice and dust tracers:
       INTEGER,SAVE :: nqdust(nqmx) ! to store the indexes of dust tracers
       INTEGER,SAVE :: i_ice=0  ! water ice
+      real,parameter :: odpref=610. ! DOD reference pressure (Pa)
       CHARACTER(LEN=20) :: txt ! to temporarly store text
       CHARACTER(LEN=1) :: txt2 ! to temporarly store text
@@ -192 +193 @@
       END IF
 
-c     Vertical column optical depth at 700.Pa 
-c     ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+c     Vertical column optical depth at "odpref" Pa 
+c     ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
       IF(iaervar.eq.1) THEN 
          do ig=1, ngridmx
@@ -224 +225 @@
      &    (taueq-tauS)*0.5*(1+tanh((45+lati(ig)*180./pi)*6/60))
         end do
-      ELSE IF ((iaervar.eq.4).or.
-     &        ((iaervar.ge.24).and.(iaervar.le.26)))
-     &     THEN  ! << "TES assimilated dust scenarios >>
-        call readtesassim(ngrid,nlayer,zday,pplev,tauref)
-
       ELSE IF (iaervar.eq.5) THEN   ! << Escalier Scenario>>
 c         tauref(1) = 0.2
@@ -240 +236 @@
           tauref(ig) = tauref(1)
         end do
+      ELSE IF ((iaervar.ge.24).and.(iaervar.le.30))
+     &     THEN  ! << MY... dust scenarios >>
+        call read_dust_scenario(ngrid,nlayer,zday,pplev,tauref)
+      ELSE IF ((iaervar.eq.4).or.
+     &         ((iaervar.ge.124).and.(iaervar.le.126))) THEN
+       ! "old" TES assimation dust scenario (values at 700Pa in files!)
+        call read_dust_scenario(ngrid,nlayer,zday,pplev,tauref)
       ELSE
         stop 'problem with iaervar in aeropacity.F'
@@ -284 +287 @@
             DO ig=1,ngrid
 c             Typical mixing ratio profile 
-              if(pplay(ig,l).gt.700.
+              if(pplay(ig,l).gt.odpref
      $                        /(988.**(topdust(ig)/70.))) then
-                zp=(700./pplay(ig,l))**(70./topdust(ig))
+                zp=(odpref/pplay(ig,l))**(70./topdust(ig))
                  expfactor=max(exp(0.007*(1.-max(zp,1.))),1.e-3)
               else    
@@ -302 +305 @@
           DO l=1,nlayer
              DO ig=1,ngrid
-                zp=700./pplay(ig,l)
-                aerosol(ig,l,1)= tauref(ig)/700. *
+                zp=odpref/pplay(ig,l)
+                aerosol(ig,l,1)= tauref(ig)/odpref *
      s           (pplev(ig,l)-pplev(ig,l+1))
      s           *max( exp(.03*(1.-max(zp,1.))) , 1.E-3 )
@@ -409 +412 @@
 c Rescaling each layer to reproduce the choosen (or assimilated)
 c   dust extinction opacity at visible reference wavelength, which
-c   is originally scaled to an equivalent 700Pa pressure surface.
+c   is originally scaled to an equivalent odpref Pa pressure surface.
 c -----------------------------------------------------------------
 
@@ -427 +430 @@
       DO ig=1,ngrid
           tauscaling(ig) = tauref(ig) * 
-     &                     pplev(ig,1) / 700.E0 / taudusttmp(ig)
+     &                     pplev(ig,1) / odpref / taudusttmp(ig)
 c          tauscaling(ig) = 1.e-4
       ENDDO