Changeset 284 for trunk/LMDZ.MARS/libf/phymars/meso_inc

Timestamp:

Sep 7, 2011, 3:20:37 PM (13 years ago)

Author:

acolaitis

Message:

--- AC 07/09/2011 ---

• Added new flag for the Richardson-based surface layer :

callrichsl, can be changed in callphys.def

One should always use the thermals model when using this surface layer model.
Somes cases (weakly unstable with low winds), when not using thermals, won't be well represented by the
Richardson surface layer. This stands for Mesoscale and Gcm but not for LES model.

Correct configs :

callrichsl = .true.
calltherm = .true.

callrichsl = .false.
calltherm = .false.

callrichsl = .false.
calltherm = .true.

Previously unstable config :

callrichsl = .true.
calltherm = .false.

• To be able to run without thermals and with the new surface layer, a modification has been made to

physiq.F to account for gustiness in GCM and MESOSCALE for negative Richardson, so that :

callrichsl = .true.
calltherm = .false.

can now be used without problems, but is not recommended.

• Consequently, callrichsl = .false. is now the default configuration for thermals.

We recall the available options in callphys.def for thermals :

outptherm = BOOLEAN (.false. by default) : outputs thermals related quantities (lots of diagfi)
nsplit_thermals = INTEGER (50 by default in gcm, 2 in mesoscale) : subtimestep for thermals model.

It is recommended to use at least 40 in the gcm, and at least 2 in the mesoscale.
The user can lower these values but should check it's log for anomalies or errors regarding
tracer transport in the thermals, or "granulosity" in the outputs for wmax, lmax and hfmax.

---

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

@@ -1 @@
-         DO ig=1,ngrid
-          !! sensible heat flux in W/m2
-!          hfx(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)
+         if (callrichsl .eq. .false.) then
+       
+             DO ig=1,ngrid
+!! sensible heat flux in W/m2
+
+             hfx(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)
+
+!! u star in similarity theory in m/s
+             ust(ig) = 0.4
+     .               * sqrt( pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) )
+     .               / log( 1.E+0 + zzlay(ig,1)/z0_default )
+             ENDDO
+
+         else
+
+            DO ig=1,ngrid
 
 ! New SL parametrization, correct formulation for hfx :
 
-          hfx(ig) = (pplay(ig,1)/(r*pt(ig,1)))*cpp
-     &    *sqrt((pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)))
-     &    *zcdh(ig)*(tsurf(ig)-zh(ig,1))
+            hfx(ig) = (pplay(ig,1)/(r*pt(ig,1)))*cpp
+     &        *sqrt(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)
+     &        + (1.0*wmax_th(ig))**2)
+     &        *zcdh(ig)*(tsurf(ig)-zh(ig,1))
 
-          !! u star in similarity theory in m/s
-!          ust(ig) = 0.4
-!     .               * sqrt( pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) )
-!     .               / log( 1.E+0 + zzlay(ig,1)/z0_default )
 
 ! New SL parametrization, ust is more accurately computed in vdif_cd :
-        ust(ig) = sqrt(zcdv(ig)*(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)))
+            ust(ig) = sqrt(zcdv(ig)*
+     &   (pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) + (1.0*wmax_th(ig))**2)
+     &                     )
 
-         ENDDO   
+            ENDDO   
+
+         endif  !of if callrichsl
+
 !         write (*,*) 'PHYS HFX cp zdts', hfx(100), zflubid(100), 
 !     .       capcal(100),