Changeset 2414

Timestamp:

Oct 6, 2020, 4:13:20 PM (4 years ago)

Author:

abierjon

Message:

Mars GCM:
Correction of a bad unit for qdusttotal0 and qdusttotal1 (diagnostic outputs of the stormdust scheme)
+ some corrections of comments about aerosol and tauref in aeropacity

AB

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/aeropacity_mod.F (modified) (13 diffs)
• libf/phymars/physiq_mod.F (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -3155 +3155 @@
 Cleanup around "aeropacity" to prepare future evolutions; added module
 dust_scaling_mod to handle computation of tauscaling.
+
+== 06/10/2020 == AB
+Correction of a bad unit for qdusttotal0 and qdusttotal1 (diagnostic outputs of the stormdust scheme)
++ some corrections of comments about aerosol and tauref in aeropacity

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

@@ -371 +371 @@
 c             avoid unrealistic values due to constant lifting:
               DO ig=1,ngrid
-              ! OPTICAL DEPTH for the computation of tauref, 
-              ! which is to be compared with tauref_scenario
+              ! OPTICAL DEPTH used in the radiative transfer
               ! => visible wavelength
                 aerosol(ig,l,iaer) = 
@@ -380 +379 @@
      &          ( pplev(ig,l) - pplev(ig,l+1) ) / g
               ! DENSITY SCALED OPACITY :
-              ! GCM output to be compared with observations
+              ! Diagnostic output to be compared with observations
               ! => infrared wavelength
                 dsodust(ig,l) =
@@ -389 +388 @@
             ELSE
               DO ig=1,ngrid
-              ! OPTICAL DEPTH for the computation of tauref, 
-              ! which is to be compared with tauref_scenario
+              ! OPTICAL DEPTH used in the radiative transfer
               ! => visible wavelength 
                 aerosol(ig,l,iaer) =
@@ -398 +396 @@
      &          ( pplev(ig,l) - pplev(ig,l+1) ) / g
               ! DENSITY SCALED OPACITY :
-              ! GCM output to be compared with observations
+              ! Diagnostic output to be compared with observations
               ! => infrared wavelength
                 dsodust(ig,l) =
@@ -491 +489 @@
 c           avoid unrealistic values due to constant lifting:
                DO ig=1,ngrid
-               ! OPTICAL DEPTH for the computation of tauref, 
-               ! which is to be compared with tauref_scenario
+               ! OPTICAL DEPTH used in the radiative transfer
                ! => visible wavelength
                  aerosol(ig,l,iaer) = 
@@ -500 +497 @@
      &           ( pplev(ig,l) - pplev(ig,l+1) ) / g
                ! DENSITY SCALED OPACITY :
-               ! GCM output to be compared with observations
+               ! Diagnostic output to be compared with observations
                ! => infrared wavelength
                  dsords(ig,l) =
@@ -509 +506 @@
              ELSE
                DO ig=1,ngrid
-               ! OPTICAL DEPTH for the computation of tauref, 
-               ! which is to be compared with tauref_scenario
+               ! OPTICAL DEPTH used in the radiative transfer
                ! => visible wavelength
                  aerosol(ig,l,iaer) =
@@ -518 +514 @@
      &           ( pplev(ig,l) - pplev(ig,l+1) ) / g
                ! DENSITY SCALED OPACITY :
-               ! GCM output to be compared with observations
+               ! Diagnostic output to be compared with observations
                ! => infrared wavelength
                  dsords(ig,l) =
@@ -543 +539 @@
 c           avoid unrealistic values due to constant lifting:
                DO ig=1,ngrid
+               ! OPTICAL DEPTH used in the radiative transfer
+               ! => visible wavelength
                   aerosol(ig,l,iaer) = 
      &           (  0.75 * QREFvis3d(ig,cstdustlevel,iaer) /
@@ -549 +547 @@
      &           ( pplev(ig,l) - pplev(ig,l+1) ) / g
                ! DENSITY SCALED OPACITY :
-               ! GCM output to be compared with observations
+               ! Diagnostic output to be compared with observations
                ! => infrared wavelength
                  dsotop(ig,l) =
@@ -558 +556 @@
              ELSE
                DO ig=1,ngrid
-               ! OPTICAL DEPTH for the computation of tauref, 
-               ! which is to be compared with tauref_scenario
+               ! OPTICAL DEPTH used in the radiative transfer
                ! => visible wavelength
                  aerosol(ig,l,iaer) =
@@ -567 +564 @@
      &           ( pplev(ig,l) - pplev(ig,l+1) ) / g
                ! DENSITY SCALED OPACITY :
-               ! GCM output to be compared with observations
+               ! Diagnostic output to be compared with observations
                ! => infrared wavelength
                  dsotop(ig,l) =
@@ -748 +745 @@
 #endif
 c      MV19: tauref must ALWAYS contain the opacity of all dust tracers
+       ! GCM DUST OPTICAL DEPTH tauref is to be compared
+       ! with the observation CDOD tauref_scenario
+       ! => visible wavelength
        IF (name_iaer(iaerdust(iaer)).eq."dust_doubleq") THEN
               tauref(ig) = tauref(ig) +

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

@@ -417 +417 @@
       REAL rdsndust(ngrid,nlayer) ! true n stormdust (kg/kg)
       REAL rdsqdust(ngrid,nlayer) ! true q stormdust (kg/kg)
-      REAL wspeed(ngrid,nlayer+1) ! vertical velocity stormdust tracer 
+      REAL wspeed(ngrid,nlayer+1) ! vertical velocity stormdust tracer
+      REAL wtop(ngrid,nlayer+1) ! vertical velocity topdust tracer
+       
       REAL dsodust(ngrid,nlayer) ! density scaled opacity for background dust
       REAL dsords(ngrid,nlayer) ! density scaled opacity for stormdust
       REAL dsotop(ngrid,nlayer) ! density scaled opacity for topdust
-      REAL wtop(ngrid,nlayer+1) ! vertical velocity topdust tracer
 
       REAL nccnco2(ngrid,nlayer)   ! true n ccnco2 (kg/kg)
@@ -1178 +1179 @@
          
          call writediagfi(ngrid,'qdusttotal0','q sum before rds',
-     &           'kg/kg ',2,qdusttotal0)
+     &           'kg/m2 ',2,qdusttotal0)
          call writediagfi(ngrid,'qdusttotal1','q sum after rds',
-     &           'kg/kg ',2,qdusttotal1)
+     &           'kg/m2 ',2,qdusttotal1)
 
       ENDIF ! end of if(rdstorm)