source: trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F @ 260

Last change on this file since 260 was 256, checked in by acolaitis, 14 years ago

--- AC 03/08/2011 ---

M 255 libf/phymars/physiq.F
<> Modified to interface new SL parametrization

M 255 libf/phymars/vdif_cd.F
<> New SL parametrization based on a bulk Richardson Monin-Obukhov theory formulation

Stability functions are taken from D.E. England et al. (95)
Similarity functions coefficients based on Dyer and Hicks (70).
Includes thermal roughness length computation, heat and momentum drag coefficient computation
Can be used to output hydrodynamic-related SL quantities (bulk Richardson, turbulent Prandtl number estimation, Reynolds number)

M 255 libf/phymars/thermcell_dqupdown.F90
<> Minor modification to suit picky compilers

M 255 libf/phymars/vdifc.F
<> Now takes into account sub-grid gustiness, evaluated from thermals activity (it's proxy being the maximum vertical velocity)

M 255 libf/phymars/meso_inc/meso_inc_les.F
<> Minor modification : u* is now taken from the new vdifc and not recomputed from a simple law
File size: 1.5 KB
Line 
1         DO ig=1,ngrid
2          !! sensible heat flux in W/m2
3          hfx(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)
4          !! u star in similarity theory in m/s
5!          ust(ig) = 0.4
6!     .               * sqrt( pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) )
7!     .               / log( 1.E+0 + zzlay(ig,1)/z0_default )
8
9! New SL parametrization, ust is more accurately computed in vdif_cd :
10          ust(ig) = zcdv(ig)
11
12         ENDDO   
13!         write (*,*) 'PHYS HFX cp zdts', hfx(100), zflubid(100),
14!     .       capcal(100),
15!     .       zdtsdif(100)
16!         write (*,*) 'PHYS UST', ust(100)
17
18!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
19!!! LES LES
20       IF (flag_LES) THEN       
21
22         write (*,*) '************************************************'
23         write (*,*) '** LES mode: the difv part is only used to'
24         write (*,*) '**  provide HFX and UST to the dynamics'
25         write (*,*) '** NB: - dudif, dvdif, dhdif, dqdif are set to 0'
26         write (*,*) '**     - tsurf is updated'     
27         write (*,*) '************************************************'
28
29         DO ig=1,ngrid
30          DO l=1,nlayer
31            zdvdif(ig,l) = 0.
32            zdudif(ig,l) = 0.
33            zdhdif(ig,l) = 0.
34            DO iq=1, nq
35              zdqdif(ig,l,iq) = 0.
36              zdqsdif(ig,iq) = 0. !! sortir de la boucle
37            ENDDO
38          ENDDO
39         ENDDO
40      ENDIF
41!!! LES LES       
42!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Note: See TracBrowser for help on using the repository browser.