Changeset 790 for trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F

Timestamp:

Sep 20, 2012, 5:50:17 PM (12 years ago)

Author:

acolaitis

Message:

MESOSCALE. Fixed lifting for LES runs (put lifting after newsedim otherwise lifting is stupidly counteracted). Transferred LES def files to SIMU DEF. Added def files for high resolution monster LES. Fixed a small bug in makemeso for debug.

File:

• trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F

@@ -55 +55 @@
             zdudif(ig,l) = 0.
             zdhdif(ig,l) = 0.
-            DO iq=1, nq
-              zdqdif(ig,l,iq) = 0.
-              zdqsdif(ig,iq) = 0. !! sortir de la boucle
-            ENDDO 
           ENDDO
          ENDDO
+         IF (lifting .and. doubleq) THEN 
+         !! lifted dust is injected in the first layer. 
+         !! Sedimentation must be called after turbulent mixing, i.e. on next step, after WRF. 
+         !! => lifted dust is not incremented before the sedimentation step.
+            zdqdif(:,1,:)=0.
+            zdqdif(:,1,igcm_dust_number)=-zdqsdif(:,igcm_dust_number)
+            zdqdif(:,1,igcm_dust_mass)=-zdqsdif(:,igcm_dust_mass)
+            zdqdif(:,2:nlayermx,:)=0.
+            DO iq=1, nq
+              IF ((iq .ne. igcm_dust_mass)
+     &        .and. (iq .ne. igcm_dust_number)) THEN
+                zdqsdif(:,iq)=0.
+              ENDIF
+            ENDDO
+         ELSE
+            DO ig=1,ngrid
+             DO l=1,nlayer
+              DO iq=1, nq
+                zdqdif(ig,l,iq) = 0.
+                zdqsdif(ig,iq) = 0. !! sortir de la boucle
+              ENDDO 
+             ENDDO
+            ENDDO
+         ENDIF
       ENDIF
 !!! LES LES