Changeset 2160

Timestamp:

Sep 23, 2019, 6:02:21 PM (5 years ago)

Author:

mvals

Message:

Mars GCM:
Set adpatable parameters for the rocket dust storm scheme (parameters included in callkeys.h, and adaptable according to the callphys.def with the function "call getin" in conf_phys.F):

• ti_injection, tf_injection (by default: ti_injection=10. and tf_injection=12., impacted files: compute_dtau_mod.F90, vdifc_mod.F)
• coeff_detrainment (by default: coeff_detrainment=0., impacted files: rocketduststorm_mod.F90)

MV

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/callkeys.h (modified) (2 diffs)
• libf/phymars/compute_dtau_mod.F90 (modified) (4 diffs)
• libf/phymars/conf_phys.F (modified) (2 diffs)
• libf/phymars/rocketduststorm_mod.F90 (modified) (5 diffs)
• libf/phymars/vdifc_mod.F (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -2737 +2737 @@
 == 12/09/2019 == MV
 - Update of compute_dtau_mod.F90: dtau is calculated in function of the dust opacity given by the dust scenario the day after
+
+== 23/09/2019 == MV
+Set adpatable parameters for the rocket dust storm scheme (parameters included in callkeys.h, and adaptable according to the callphys.def with the function "call getin" in conf_phys.F):
+- ti_injection, tf_injection (by default: ti_injection=10. and tf_injection=12., impacted files: compute_dtau_mod.F90, vdifc_mod.F)
+- coeff_detrainment (by default: coeff_detrainment=0., impacted files: rocketduststorm_mod.F90)
+

trunk/LMDZ.MARS/libf/phymars/callkeys.h

@@ -22 +22 @@
       COMMON/callkeys_r/topdustref,semi,alphan,euveff,                  &
      &   tke_heat_flux,dustrefir,fixed_euv_value,CLFfixval,             &
-     &   coeff_injection
+     &   coeff_injection,ti_injection,tf_injection,coeff_detrainment
      
       LOGICAL callrad,calldifv,calladj,callcond,callsoil,               &
@@ -44 +44 @@
       real tke_heat_flux
       real coeff_injection ! dust injection scheme coefficient
+      real ti_injection ! local time of beginning injection
+      real tf_injection ! local time of end injection
+      real coeff_detrainment ! rocket dust detrainment coefficient
       real CLFfixval
 

trunk/LMDZ.MARS/libf/phymars/compute_dtau_mod.F90

@@ -3 +3 @@
         IMPLICIT NONE
 
-        REAL,PARAMETER :: ti_injection=10/24. ! time of beginning injection
-        REAL,PARAMETER :: tf_injection=18/24. ! time of end injection
-        REAL,PARAMETER :: t_scenario=14/24.   ! time of tauref_scenario
+        include "callkeys.h"
+
+        REAL,SAVE :: ti_injection_sol ! time of beginning injection
+        REAL,SAVE :: tf_injection_sol ! time of end injection
+        REAL,PARAMETER :: t_scenario_sol=14/24.   ! time of tauref_scenario
 
         CONTAINS
@@ -52 +54 @@
                 ENDDO
                 nb_daystep=(daysec/dtphys)
+                ! Local time in sol fraction
+                ti_injection_sol=ti_injection/24.
+                tf_injection_sol=tf_injection/24.
                 firstcall=.FALSE.
         ENDIF
         
-        ! 1. Obtain tauref_scenario from dust scenario
+        ! 1. Obtain tauref_scenario from dust scenario at zday+1
         zday_scenario=zday-modulo(zday,1.) ! integer value of the day: the scenario opacity is measured at 14:00
         zday_scenario=zday_scenario+1      ! opacity of the dust scenario is read the day after
@@ -66 +71 @@
         ! 2. Compute dtau
         DO ig=1,ngrid
-         IF ((local_time(ig).ge.t_scenario).and.                        &
-       &          (local_time_prev(ig).lt.(t_scenario)))THEN
+         IF ((local_time(ig).ge.t_scenario_sol).and.                    &
+       &          (local_time_prev(ig).lt.(t_scenario_sol)))THEN
                  dtau(ig)=tauref_scenario(ig)-tauref(ig)
          ENDIF
@@ -77 +82 @@
          dustliftday(ig)=coeff_injection*                               &
        &                (dtau(ig)*pplev(ig,1)/odpref)                   &
-       &                /(daysec*(tf_injection-ti_injection))
+       &                /(daysec*(tf_injection_sol-ti_injection_sol))
         ENDDO 
          

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

@@ -292 +292 @@
          call getin("rdstorm",rdstorm)
          write(*,*)" rdstorm = ",rdstorm
+! rocket dust storm detrainment coefficient        
+        coeff_detrainment=0. ! default value
+        call getin("coeff_detrainment",coeff_detrainment)
+        write(*,*)" coeff_detrainment = ",coeff_detrainment
 
          write(*,*)"rad.transfer is computed every iradia",
@@ -379 +383 @@
         call getin("coeff_injection",coeff_injection)
         write(*,*)" coeff_in,jection = ",coeff_injection
+! timing for dust injection        
+        ti_injection=10. ! default value
+        tf_injection=12. ! default value
+        call getin("ti_injection",ti_injection)
+        write(*,*)" ti_injection = ",ti_injection
+        call getin("tf_injection",tf_injection)
+        write(*,*)" tf_injection = ",tf_injection
 
 ! free evolving dust

trunk/LMDZ.MARS/libf/phymars/rocketduststorm_mod.F90

@@ -39 +39 @@
       IMPLICIT NONE
 
+      include "callkeys.h"
+
 !--------------------------------------------------------
 ! Input Variables
@@ -126 +128 @@
 
       LOGICAL storm(ngrid)    ! true when there is a dust storm (if the opacity is high): trigger the rocket dust storm scheme
-      REAL coefdetrain(ngrid,nlayer)  ! coefficient for detrainment : % of stormdust detrained
+      REAL detrain(ngrid,nlayer)  ! coefficient for detrainment : % of stormdust detrained
       INTEGER scheme(ngrid)   ! triggered scheme
             
@@ -197 +199 @@
       deltahr(:,:)=0.
       scheme(:)=0
-      coefdetrain(:,:)=1.
+      detrain(:,:)=1.
 
       !! no update for the stormdust tracer and temperature fields 
@@ -410 +412 @@
                           wmin) .or. (zq_dust_mass(ig,l) .gt.  &
                                      10000.*zq_stormdust_mass(ig,l))) THEN
-               coefdetrain(ig,l)=1.
+               detrain(ig,l)=1.
             ELSE IF (max(abs(wrad(ig,l)),abs(wrad(ig,l+1)))   &
                                                        .le. wmax) THEN
-               coefdetrain(ig,l)=1.*(((1-coefmin)/(wmin-wmax)**2)*     &
-                   (max(abs(wrad(ig,l)),abs(wrad(ig,l+1)))-wmax)**2 &
-                                                               +coefmin)
+               detrain(ig,l)=coeff_detrainment*                 &
+                             (((1-coefmin)/(wmin-wmax)**2)*     &
+                             (max(abs(wrad(ig,l)),abs(wrad(ig,l+1)))-wmax)**2 &
+                              +coefmin)
             ELSE IF (max(abs(wrad(ig,l)),abs(wrad(ig,l+1))).gt. wmax ) THEN
-               coefdetrain(ig,l)=coefmin
+               detrain(ig,l)=coefmin
             ELSE
-               coefdetrain(ig,l)=coefmin
+               detrain(ig,l)=coefmin
             ENDIF
           ENDDO ! DO l=1,nlayer
@@ -428 +431 @@
         DO ig=1,ngrid
           DO l=1,nlayer
-            dqdet_stormdust_mass(ig,l)=-coefdetrain(ig,l)*zq_stormdust_mass(ig,l) &
+            dqdet_stormdust_mass(ig,l)=-detrain(ig,l)*zq_stormdust_mass(ig,l) &
                                                         /ptimestep
-            dqdet_stormdust_number(ig,l)=-coefdetrain(ig,l)*zq_stormdust_number(ig,l) &
+            dqdet_stormdust_number(ig,l)=-detrain(ig,l)*zq_stormdust_number(ig,l) &
                                                         /ptimestep
           ENDDO ! DO l=1,nlayer

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

@@ -23 +23 @@
       use watersat_mod, only: watersat
       use turb_mod, only: turb_resolved, ustar, tstar
-      use compute_dtau_mod, only: ti_injection,tf_injection
+      use compute_dtau_mod, only: ti_injection_sol,tf_injection_sol
       IMPLICIT NONE
 
@@ -717 +717 @@
               do ig=1,ngrid
                if(co2ice(ig).lt.1) then ! pas de soulevement si glace CO2
-                IF((ti_injection.LE.local_time(ig)).and.
-     &                  (local_time(ig).LE.tf_injection)) THEN
+                IF((ti_injection_sol.LE.local_time(ig)).and.
+     &                  (local_time(ig).LE.tf_injection_sol)) THEN
                   if (rdstorm) then !Rocket dust storm scheme
                         pdqsdif(ig,igcm_stormdust_mass) =