Index: trunk/LMDZ.MARS/libf/phymars/physiq_mod.F
===================================================================
--- trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 2246)
+++ trunk/LMDZ.MARS/libf/phymars/physiq_mod.F	(revision 2247)
@@ -2924,11 +2924,21 @@
              call WRITEDIAGFI(ngrid,'dustN','Dust number',
      &                        'part/kg',3,ndust)
-             call WRITEDIAGFI(ngrid,'dsodust',
-     &                        'density scaled optical depth in the IR',
-     &                        'm2.kg-1',3,dsodust)
-             call WRITEDIAGFI(ngrid,'dso',
-     &                        'density scaled optical depth in the IR',
-     &                        'm2.kg-1',3,dsodust+dsords+dsotop)
-	     
+             
+	     select case (trim(dustiropacity))
+              case ("tes")
+               call WRITEDIAGFI(ngrid,'dsodust',
+     &  'density scaled extinction opacity of std dust at 9.3um(TES)',
+     &         'm2.kg-1',3,dsodust)
+               call WRITEDIAGFI(ngrid,'dso',
+     &  'density scaled extinction opacity of all dust at 9.3um(TES)',
+     &         'm2.kg-1',3,dsodust+dsords+dsotop)
+              case ("mcs")
+               call WRITEDIAGFI(ngrid,'dsodust',
+     &  'density scaled extinction opacity of std dust at 21.6um(MCS)',
+     &         'm2.kg-1',3,dsodust)
+               call WRITEDIAGFI(ngrid,'dso',
+     &  'density scaled extinction opacity of all dust at 21.6um(MCS)',
+     &         'm2.kg-1',3,dsodust+dsords+dsotop)
+             end select
            else ! (doubleq=.false.)
              do iq=1,dustbin
@@ -2987,13 +2997,28 @@
              call WRITEDIAGFI(ngrid,'totaldustq','total dust mass',
      &                        'kg/kg',3,qdusttotal)
-             call WRITEDIAGFI(ngrid,'dsords',
-     &                       'density scaled opacity of stormdust',
-     &                       'm2.kg-1',3,dsords)
+             
+	     select case (trim(dustiropacity))
+              case ("tes")
+               call WRITEDIAGFI(ngrid,'dsords',
+     &  'density scaled extinction opacity of stormdust at 9.3um(TES)',
+     &         'm2.kg-1',3,dsords)
+              case ("mcs")
+               call WRITEDIAGFI(ngrid,'dsords',
+     &  'density scaled extinction opacity of stormdust at 21.6um(MCS)',
+     &         'm2.kg-1',3,dsords)
+             end select
            endif ! (rdstorm)
 
            if (slpwind) then
-	     call WRITEDIAGFI(ngrid,'dsotop',
-     &                       'density scaled opacity of topdust',
-     &                       'm2.kg-1',3,dsotop)
+	     select case (trim(dustiropacity))
+              case ("tes")
+               call WRITEDIAGFI(ngrid,'dsotop',
+     &  'density scaled extinction opacity of topdust at 9.3um(TES)',
+     &         'm2.kg-1',3,dsotop)
+              case ("mcs")
+               call WRITEDIAGFI(ngrid,'dsotop',
+     &  'density scaled extinction opacity of topdust at 21.6um(MCS)',
+     &         'm2.kg-1',3,dsotop)
+             end select
            endif ! (slpwind)