Changeset 2921 for trunk/LMDZ.MARS

Timestamp:

Mar 25, 2023, 7:05:48 AM (20 months ago)

Author:

emillour

Message:

Mars PCM:
Add the possibility to use MY36 dust scenario and MY36 EUV inputs
EM

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/aeronomars/param_read_e107.F (modified) (1 diff)
• libf/phymars/aeropacity_mod.F (modified) (2 diffs)
• libf/phymars/read_dust_scenario_mod.F90 (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -3958 +3958 @@
 (removing co2ice and added islope dimension to some arrays).
 
+== 25/03/2023 == EM
+Add the possibility to use MY36 dust scenario and MY36 EUV inputs
+

trunk/LMDZ.MARS/libf/aeronomars/param_read_e107.F

@@ -252 +252 @@
 
       !Tabulated values for E10.7
-      if((solvaryear.ge.23).and.(solvaryear.le.35)) then
+      if((solvaryear.ge.23).and.(solvaryear.le.36)) then
         write(filename,fmt="(a,i2,a)")"e107_MY",solvaryear,".dat"
       else
          write(*,*)"param_read_e107: "
          write(*,*)"bad value for solvaryear in callphys.def"
-         write(*,*)"solvaryear must be between 24 and 35"
+         write(*,*)"solvaryear must be between 24 and 36"
          stop
       endif

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

@@ -311 +311 @@
         call read_dust_scenario(ngrid,nlayer,zday,pplev,
      &                          IRtoVIScoef,tau_pref_scenario)
-      ELSE IF ((iaervar.ge.24).and.(iaervar.le.35))
+      ELSE IF ((iaervar.ge.24).and.(iaervar.le.36))
      &     THEN  ! << MY... dust scenarios >>
         call read_dust_scenario(ngrid,nlayer,zday,pplev,
@@ -706 +706 @@
           call read_dust_scenario(ngrid,nlayer,zday,pplev,
      &                            IRtoVIScoef,tau_pref_scenario)
-        ELSE IF ((iaervar.ge.24).and.(iaervar.le.35))
+        ELSE IF ((iaervar.ge.24).and.(iaervar.le.36))
      &       THEN  ! << MY... dust scenarios >>
           call read_dust_scenario(ngrid,nlayer,zday,pplev,

trunk/LMDZ.MARS/libf/phymars/read_dust_scenario_mod.F90

@@ -83 +83 @@
    else if (iaervar.eq.8) then
      filename="dust_clim.nc"
-   else if ((iaervar.ge.24).and.(iaervar.le.35)) then
+   else if ((iaervar.ge.24).and.(iaervar.le.36)) then
      write(filename,'(a7,i2,a3)')"dust_MY",iaervar,".nc"
    ! 124,125,126: old TES assimilated dust scenarios (at 700Pa ref pressure!)