Changeset 1974 for trunk/LMDZ.MARS/libf/phymars/updatereffrad_mod.F

Timestamp:

Jul 18, 2018, 4:48:34 PM (6 years ago)

Author:

mvals

Message:

Mars GCM:
Integration of the detached dust layer parametrizations (rocket dust storm, slope wind lifting, CW, and dust injection scheme, DB).
Still experimental, default behaviour (rdstorm=.false., dustinjection=0) identical to previous revision.
NB: Updated newstart requires an updated "surface.nc" containing the "hmons" field.
EM+MV

File:

• trunk/LMDZ.MARS/libf/phymars/updatereffrad_mod.F (modified) (7 diffs)

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

@@ -6 +6 @@
       
       SUBROUTINE updatereffrad(ngrid,nlayer,
-     &                rdust,rice,nuice,
+     &                rdust,rstormdust,rice,nuice,
      &                reffrad,nueffrad,
      &                pq,tauscaling,tau,pplay)
@@ -14 +14 @@
      &                       igcm_h2o_ice, igcm_ccn_mass, radius,
      &                       igcm_ccn_number, nuice_ref, varian,
-     &                       ref_r0, igcm_dust_submicron
+     &                       ref_r0, igcm_dust_submicron,
+     &                       igcm_stormdust_mass,igcm_stormdust_number
        USE dimradmars_mod, only: nueffdust,naerkind,
      &            name_iaer,
      &            iaer_dust_conrath,iaer_dust_doubleq,
-     &            iaer_dust_submicron,iaer_h2o_ice
+     &            iaer_dust_submicron,iaer_h2o_ice,
+     &            iaer_stormdust_doubleq 
 
        IMPLICIT NONE
@@ -45 +47 @@
 
 c-----------------------------------------------------------------------
-c     Inputs:
+c     Inputs/outputs:
 c     ------
 
-      INTEGER ngrid,nlayer
+      INTEGER, INTENT(in) :: ngrid,nlayer
 c     Ice geometric mean radius (m)
-      REAL :: rice(ngrid,nlayer)
+      REAL, INTENT(out) :: rice(ngrid,nlayer)
 c     Estimated effective variance of the size distribution (n.u.)
-      REAL :: nuice(ngrid,nlayer)
+      REAL, INTENT(out) :: nuice(ngrid,nlayer)
 c     Tracer mass mixing ratio (kg/kg)
-      REAL pq(ngrid,nlayer,nqmx)
-      REAL rdust(ngrid,nlayer) ! Dust geometric mean radius (m)
-      
-      REAL pplay(ngrid,nlayer) ! altitude at the middle of the layers
-      REAL tau(ngrid,naerkind)
-
-
-c     Outputs:
-c     -------
-
+      REAL, INTENT(in) :: pq(ngrid,nlayer,nqmx)
+      REAL, INTENT(out) :: rdust(ngrid,nlayer) ! Dust geometric mean radius (m)
+      REAL, INTENT(out) :: rstormdust(ngrid,nlayer) ! Dust geometric mean radius (m) 
+      REAL, INTENT(in) :: pplay(ngrid,nlayer) ! altitude at the middle of the layers
+      REAL, INTENT(in) :: tau(ngrid,naerkind)
 c     Aerosol effective radius used for radiative transfer (meter)
-      REAL :: reffrad(ngrid,nlayer,naerkind)
+      REAL, INTENT(out) :: reffrad(ngrid,nlayer,naerkind)
 c     Aerosol effective variance used for radiative transfer (n.u.)
-      REAL :: nueffrad(ngrid,nlayer,naerkind)
-
+      REAL, INTENT(out) :: nueffrad(ngrid,nlayer,naerkind)
+      REAL, INTENT(in) :: tauscaling(ngrid)         ! Convertion factor for qccn and Nccn
+      
 c     Local variables:
 c     ---------------
@@ -85 +83 @@
       REAL Mo,No                       ! Mass and number of ccn
       REAL rhocloud(ngrid,nlayer)  ! Cloud density (kg.m-3)
-      REAL tauscaling(ngrid)         ! Convertion factor for qccn and Nccn
 
       LOGICAL,SAVE :: firstcall=.true.
@@ -114 +111 @@
           ENDDO
         ENDIF
+
+        ! updating radius of stormdust particles
+        IF (rdstorm.AND.active) THEN
+          DO l=1,nlayer
+            DO ig=1, ngrid
+              call updaterdust(pq(ig,l,igcm_stormdust_mass),
+     &                 pq(ig,l,igcm_stormdust_number),rstormdust(ig,l))
+              nueffdust(ig,l) = exp(varian**2.)-1.
+             ENDDO
+           ENDDO
+        ENDIF
         
 c       1.2 Water-ice particles
@@ -126 +134 @@
 
           IF (firstcall) THEN
-            IF (minval(tauscaling).lt.0) tauscaling(:) = 1.e-3 ! default value when non-read in startfi is -1
-            IF (freedust)                tauscaling(:) = 1.    ! if freedust, enforce no rescaling at all
+            !IF (minval(tauscaling).lt.0) tauscaling(:) = 1.e-3 ! default value when non-read in startfi is -1
+            !IF (freedust)                tauscaling(:) = 1.    ! if freedust, enforce no rescaling at all
             firstcall = .false.
           ENDIF
@@ -206 +214 @@
           ENDDO
 c==================================================================
+        CASE("stormdust_doubleq") aerkind! Two-moment scheme for
+c       stormdust; same distribution than normal dust
+c==================================================================
+          DO l=1,nlayer
+            DO ig=1,ngrid
+              reffrad(ig,l,iaer) = rstormdust(ig,l) * ref_r0
+              nueffrad(ig,l,iaer) = nueffdust(ig,l)
+            ENDDO
+          ENDDO
+c==================================================================
         END SELECT aerkind
       ENDDO ! iaer (loop on aerosol kind)