Changeset 4030

Timestamp:

Jan 21, 2026, 5:22:25 PM (11 hours ago)

Author:

emillour

Message:

Mars PCM:

• add possibility to read MY37 dust scenario and the possibility to read in a scenario provided in IR extinction opacity (at 9.3um) identified as such by being "cdodext" in the input file (whereas "cdod" is for IR absorption).
• add possibility to load EUV daily E10.7 input for MY37.

EM

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/aeronomars/param_read_e107.F (modified) (1 diff)
• libf/phymars/aeropacity_mod.F (modified) (1 diff)
• libf/phymars/read_dust_scenario_mod.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -5009 +5009 @@
 == 12/01/2026 == JBC
 Deletion of unused legacy code in "initracer.F90".
+
+== 21/01/2026 == EM
+- add possibility to read MY37 dust scenario and the possibility to read in
+  a scenario provided in IR extinction opacity (at 9.3um) identified as such
+  by being "cdodext" in the input file (whereas "cdod" is for IR absorption).
+- add possibility to load EUV daily E10.7 input for MY37.

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

@@ -255 +255 @@
 
       !Tabulated values for E10.7
-      if((solvaryear.ge.23).and.(solvaryear.le.36)) then
+      if((solvaryear.ge.23).and.(solvaryear.le.37)) 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 36"
+         write(*,*)"solvaryear must be between 24 and 37"
          stop
       endif

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

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

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

@@ -84 +84 @@
    else if (iaervar.eq.8) then
      filename="dust_clim.nc"
-   else if ((iaervar.ge.24).and.(iaervar.le.36)) then
+   else if ((iaervar.ge.24).and.(iaervar.le.37)) then
      write(filename,'(a7,i2,a3)')"dust_MY",iaervar,".nc"
    ! 124,125,126: old TES assimilated dust scenarios (at 700Pa ref pressure!)
@@ -139 +139 @@
    ! "dustop" if loading visible extinction opacity
    ! "cdod" if loading IR absorption opacity
+   ! "cdodext" if loading IR extinction opacity
    ierr=nf90_inq_varid(nid,"dustop",nvarid)
    if (ierr.eq.nf90_noerr) then
@@ -149 +150 @@
      ENDIF
    else
-     ! did not find "dustop" , look for "cdod"
+     ! did not find "dustop" , look for "cdod" or "cdodext"
      IRscenario = .true.
      ierr=nf90_inq_varid(nid,"cdod",nvarid)
-     ierr=nf90_get_var(nid,nvarid,tautes)
-     IF (ierr .NE. nf90_noerr) THEN
+     if (ierr.eq.nf90_noerr) then
+      ierr=nf90_get_var(nid,nvarid,tautes)
+      IF (ierr .NE. nf90_noerr) THEN
         PRINT*, "Error: read_dust_scenario <cdod> not found"
         write(*,*)trim(nf90_strerror(ierr))
         call abort_physic("read_dust_scenario","cdod not found",1)
-     ENDIF
+      ENDIF
+     else ! did not find "cdod" , look for "cdodext"
+      ierr=nf90_inq_varid(nid,"cdodext",nvarid)
+      ierr=nf90_get_var(nid,nvarid,tautes)
+      IF (ierr .NE. nf90_noerr) THEN
+        PRINT*, "Error: read_dust_scenario <cdodext> not found"
+        write(*,*)trim(nf90_strerror(ierr))
+        call abort_physic("read_dust_scenario","cdodext not found",1)
+      ELSE
+        ! convert IR extinction to IR absorption
+        tautes(:,:,:)=tautes(:,:,:)/1.3
+      ENDIF
+     endif
    endif