Changeset 2161 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Sep 24, 2019, 4:18:50 PM (5 years ago)

Author:

emillour

Message:

Mars GCM:
Reactivate output of density scaled opacity
EM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• aeropacity_mod.F (modified) (2 diffs)
• physiq_mod.F (modified) (1 diff)

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

@@ -384 +384 @@
      &          pq(ig,cstdustlevel,igcm_dust_mass) *
      &          ( pplev(ig,l) - pplev(ig,l+1) ) / g
-                ! DENSITY SCALED OPACITY IN INFRARED:
-!                dsodust(ig,l) =
-!     &          (  0.75 * QREFir3d(ig,cstdustlevel,iaer) /
-!     &          ( rho_dust * reffrad(ig,cstdustlevel,iaer) )  ) *
-!     &          pq(ig,cstdustlevel,igcm_dust_mass)
+                ! DENSITY SCALED OPACITY IN THE INFRARED:
+                dsodust(ig,l) =
+     &          (  0.75 * QREFir3d(ig,cstdustlevel,iaer) /
+     &          ( rho_dust * reffrad(ig,cstdustlevel,iaer) )  ) *
+     &          pq(ig,cstdustlevel,igcm_dust_mass)
               ENDDO
             ELSE
@@ -397 +397 @@
      &          pq(ig,l,igcm_dust_mass) *
      &          ( pplev(ig,l) - pplev(ig,l+1) ) / g
-                ! DENSITY SCALED OPACITY IN INFRARED:
-!                dsodust(ig,l) =
-!     &          (  0.75 * QREFir3d(ig,l,iaer) /
-!     &          ( rho_dust * reffrad(ig,l,iaer) )  ) *
-!     &          pq(ig,l,igcm_dust_mass)
+                ! DENSITY SCALED OPACITY IN THE INFRARED:
+                dsodust(ig,l) =
+     &          (  0.75 * QREFir3d(ig,l,iaer) /
+     &          ( rho_dust * reffrad(ig,l,iaer) )  ) *
+     &          pq(ig,l,igcm_dust_mass)
               ENDDO
             ENDIF

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

@@ -2798 +2798 @@
      &                        'part/kg',3,ndust)
              call WRITEDIAGFI(ngrid,'dsodust',
-     &                       'density scaled optcial depth',
+     &                       'density scaled optical depth in the IR',
      &                       'm2.kg-1',3,dsodust)
              call WRITEDIAGFI(ngrid,'dso',
-     &                       'density scaled optcial depth',
+     &                       'density scaled optical depth in the IR',
      &                       'm2.kg-1',3,dsodust+dsords)
            else