Changeset 2090

Timestamp:

Feb 5, 2019, 6:39:39 PM (6 years ago)

Author:

mvals

Message:

Mars GCM:

• follow-up of the last change in rocketduststorm_mod.F90

MV

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/rocketduststorm_mod.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/README

@@ -2694 +2694 @@
 == 24/01/2019 == MV
 - follow-up and correction of the last change in rocketduststorm_mod.F90: coefdetrain has to be a vector
+
+== 05/02/2019 == MV
+- follow-up of the last change in rocketduststorm_mod.F90

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

@@ -118 +118 @@
       REAL dqdet_stormdust_number(ngrid,nlayer) ! tendancy of detrainement (stormdust number)
 
-      REAL masse(nlayer)     ! mass of atmosphere (kg/m2)
+      REAL masse_col(nlayer)     ! mass of atmosphere (kg/m2)
       REAL zq(ngrid,nlayer,nq)   ! updated tracers
       
-      REAL w(nlayer)          ! air mass flux (calculated with the vertical wind velocity profile) used as input in Van Leer (kgair/m2)
-      REAL wqmass(nlayer+1)   ! tracer (dust_mass) mass flux in Van Leer (kg/m2)
-      REAL wqnumber(nlayer+1) ! tracer (dust_number) mass flux in Van Leer (kg/m2)
+      REAL w(ngrid,nlayer)          ! air mass flux (calculated with the vertical wind velocity profile) used as input in Van Leer (kgair/m2)
+      REAL wqmass(ngrid,nlayer+1)   ! tracer (dust_mass) mass flux in Van Leer (kg/m2)
+      REAL wqnumber(ngrid,nlayer+1) ! tracer (dust_number) mass flux in Van Leer (kg/m2)
 
       LOGICAL storm(ngrid)    ! true when there is a dust storm (if the opacity is high): trigger the rocket dust storm scheme
@@ -190 +190 @@
       dqdet_stormdust_number(:,:)=0.
       wrad(:,:)=0.
+      w(:,:)=0.
+      wqmass(:,:)=0.
+      wqnumber(:,:)=0.
       lapserate(:,:)=0.
       deltahr(:,:)=0.
@@ -329 +332 @@
             !! Mass of atmosphere in the layer
             DO l=1,nlayer
-               masse(l)=(pplev(ig,l)-pplev(ig,l+1))/g
+               masse_col(l)=(pplev(ig,l)-pplev(ig,l+1))/g
             ENDDO
             
             !! Mass flux in kg/m2
             DO l=1,nlayer
-               w(l)=wrad(ig,l)*(pplev(ig,l)/(r*ztlev(l)))*ptimestep
+               w(ig,l)=wrad(ig,l)*(pplev(ig,l)/(r*ztlev(l)))*ptimestep
             ENDDO
             
@@ -344 +347 @@
             
             !! Van Leer scheme
-            wqmass(:)=0.
-            wqnumber(:)=0.
             CALL vl_storm(nlayer,zq_vl_col,2.,   &
-                  masse,w,wqmass)
+                  masse_col,w(ig,:),wqmass(ig,:))
             CALL vl_storm(nlayer,zn_vl_col,2.,  &
-                  masse,w,wqnumber)
+                  masse_col,w(ig,:),wqnumber(ig,:))
             !! Mass mixing ratio after transport      
             mr_stormdust_mass(ig,:) = zq_vl_col(:)