Index: trunk/LMDZ.MARS/libf/phymars/aeropacity_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/aeropacity_mod.F	(revision 4029)
+++ trunk/LMDZ.MARS/libf/phymars/aeropacity_mod.F	(revision 4030)
@@ -314,5 +314,5 @@
         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,
Index: trunk/LMDZ.MARS/libf/phymars/read_dust_scenario_mod.F90
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/read_dust_scenario_mod.F90	(revision 4029)
+++ trunk/LMDZ.MARS/libf/phymars/read_dust_scenario_mod.F90	(revision 4030)
@@ -84,5 +84,5 @@
    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,4 +139,5 @@
    ! "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,13 +150,26 @@
      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