Changeset 2083 for trunk/LMDZ.MARS/libf/phymars/rocketduststorm_mod.F90

Timestamp:

Jan 24, 2019, 4:19:42 PM (6 years ago)

Author:

mvals

Message:

Mars GCM:

• follow-up and correction of the last change in rocketduststorm_mod.F90: coefdetrain has to be a vector

MV

File:

• trunk/LMDZ.MARS/libf/phymars/rocketduststorm_mod.F90 (modified) (4 diffs)

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

@@ -126 +126 @@
 
       LOGICAL storm(ngrid)    ! true when there is a dust storm (if the opacity is high): trigger the rocket dust storm scheme
-      REAL coefdetrain        ! coefficient for detrainment : % of stormdust detrained
+      REAL coefdetrain(ngrid,nlayer)  ! coefficient for detrainment : % of stormdust detrained
       INTEGER scheme(ngrid)   ! triggered scheme
             
@@ -193 +193 @@
       deltahr(:,:)=0.
       scheme(:)=0
-      
+      coefdetrain(:,:)=1.
+
       !! no update for the stormdust tracer and temperature fields 
       !! because previous callradite was for background dust
@@ -405 +406 @@
                           wmin) .or. (zq_dust_mass(ig,l) .gt.  &
                                      10000.*zq_stormdust_mass(ig,l))) THEN
-               coefdetrain=1.
+               coefdetrain(ig,l)=1.
             ELSE IF (max(abs(wrad(ig,l)),abs(wrad(ig,l+1)))   &
                                                        .le. wmax) THEN
-               coefdetrain=1.*(((1-coefmin)/(wmin-wmax)**2)*     &
+               coefdetrain(ig,l)=1.*(((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=coefmin
+               coefdetrain(ig,l)=coefmin
             ELSE
-               coefdetrain=coefmin
+               coefdetrain(ig,l)=coefmin
             ENDIF
           ENDDO ! DO l=1,nlayer
@@ -423 +424 @@
         DO ig=1,ngrid
           DO l=1,nlayer
-            dqdet_stormdust_mass(ig,l)=-coefdetrain*zq_stormdust_mass(ig,l) &
+            dqdet_stormdust_mass(ig,l)=-coefdetrain(ig,l)*zq_stormdust_mass(ig,l) &
                                                         /ptimestep
-            dqdet_stormdust_number(ig,l)=-coefdetrain*zq_stormdust_number(ig,l) &
+            dqdet_stormdust_number(ig,l)=-coefdetrain(ig,l)*zq_stormdust_number(ig,l) &
                                                         /ptimestep
           ENDDO ! DO l=1,nlayer