Changeset 677 for trunk/LMDZ.MARS/libf


Ignore:
Timestamp:
May 29, 2012, 11:50:48 AM (13 years ago)
Author:
emillour
Message:

Mars GCM:

  • Added the possibility to use the "cold" (iaervar=6) or "warm" (iaervar=7) dust scenarios (derived from the 7 MY24 to MY30 scenarios) in aeropacity.F and read_dust_scenario.F

EM

Location:
trunk/LMDZ.MARS/libf/phymars
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/phymars/aeropacity.F

    r607 r677  
    236236          tauref(ig) = tauref(1)
    237237        end do
     238      ELSE IF ((iaervar.ge.6).and.(iaervar.le.7)) THEN
     239      ! cold or warm synthetic scenarios
     240        call read_dust_scenario(ngrid,nlayer,zday,pplev,tauref)
    238241      ELSE IF ((iaervar.ge.24).and.(iaervar.le.30))
    239242     &     THEN  ! << MY... dust scenarios >>
  • trunk/LMDZ.MARS/libf/phymars/inifis.F

    r633 r677  
    179179         write(*,*) "Dust scenario: 1=constant dust (read from startfi",
    180180     &   " or set as tauvis); 2=Viking scenario; =3 MGS scenario,",
    181      &   "=4 Mars Year 24 from TES assimilation, ",
    182      &   "=24,25 or 26 :Mars Year 24,25 or 26 from TES assimilation"
     181     &   "=6 cold (low dust) scenario; =7 warm (high dust) scenario ",
     182     &   "=24,25 ... 30 :Mars Year 24, ... or 30 from TES assimilation"
    183183         iaervar=3 ! default value
    184184         call getin("iaervar",iaervar)
  • trunk/LMDZ.MARS/libf/phymars/read_dust_scenario.F90

    r607 r677  
    4949   ! assimilated dust file: (NB: iaervar is a common in "callkeys.h")
    5050   ! iaervar=4 means read dust_tes.nc file
     51   ! iaervar=6 means read dust_cold.nc file
     52   ! iaervar=7 means read dust_warm.nc file
    5153   ! iaervar=24 means read dust_MY24.nc file
    5254   ! iaervar=25 means read dust_MY25.nc file
     
    5557   ! NB: 4: old TES assimilated MY24 dust scenarios (at 700Pa ref pressure!)
    5658     filename="dust_tes.nc"
     59   else if (iaervar.eq.6) then
     60     filename="dust_cold.nc"
     61   else if (iaervar.eq.7) then
     62     filename="dust_warm.nc"
    5763   else if (iaervar.eq.24) then
    5864     filename="dust_MY24.nc"
Note: See TracChangeset for help on using the changeset viewer.