Changeset 310 for trunk/LMDZ.MARS/libf

Timestamp:

Sep 29, 2011, 3:59:01 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE:
-- Code in storm scenario corresponds to 'OMEGA' reference case [Julien Faure]
-- Easier settings for dust lifting without the need to recompile [see below]
-- A few modifications to plot and dust-devil detection PYTHON routines

29/09/11 == AS

--> To easily explore sensitivity to lifting thresholds: in dustlift.F, ustar_seuil=sqrt(stress/rho)

and alpha_lift[dust_mass] can be prescribed through an external stress.def parameter file.
--- alpha_lift[dust_number] is computed from alpha_lift[dust_mass] as in initracer.F
--- ustar_seuil is more user-friendly than stress because direct comparison with ustar from model

--> For the moment this is MESOSCALE only, but potentially useful to everyone

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• dustlift.F (modified) (1 diff)
• vdifc.F (modified) (2 diffs)

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

@@ -50 +50 @@
 !!!! AS: you have to compile with -DMESOSCALE to do so
       REAL alpha
+      REAL r0_lift
       INTEGER ierr
+      REAL ulim
         OPEN(99,file='stress.def',status='old',form='formatted'
      .   ,iostat=ierr)
         !!! no file => default values
         IF(ierr.EQ.0) THEN
-          READ(99,*) stress_seuil
+          READ(99,*) ulim !ulim = sqrt(stress_seuil/rho) avec rho = 0.02. 
+                          !prendre ulim = 1.061 m/s pour avoir stress_seuil = 0.0225 
           READ(99,*) alpha
+          stress_seuil = 0.02 * ulim * ulim
           write(*,*) 'USER-DEFINED threshold: ', stress_seuil, alpha
           CLOSE(99)
-          alpha_lift(1:nq) = alpha
+          alpha_lift(igcm_dust_mass) = alpha
+          r0_lift = radius(igcm_dust_mass) / ref_r0
+          alpha_lift(igcm_dust_number)=r3n_q*
+     &                        alpha_lift(igcm_dust_mass)/r0_lift**3
+          write(*,*) 'set dust number: ', alpha_lift(igcm_dust_number)
         ENDIF
 #endif

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

@@ -494 +494 @@
 c       Dust lifting:
         if (lifting) then
+#ifndef MESOSCALE
            if (doubleq.AND.submicron) then
              do ig=1,ngrid
@@ -522 +523 @@
      &                    pdqsdif)
            endif !doubleq.AND.submicron
+#else
+            call dustlift(ngrid,nlay,nq,rho,zcdh_true,zcdh,co2ice,
+     &                    pdqsdif)
+#endif
         else
            pdqsdif(1:ngrid,1:nq) = 0.