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

Last change on this file since 272 was 268, checked in by acolaitis, 14 years ago

--- AC 03/08/2011 ---
M 267 libf/phymars/physiq.F
<> Minor modification to pass Ch from vdifc to meso_inc_les

M 267 libf/phymars/surflayer_interpol.F
<> Major modification to the formulation of integrals

Now stable for most cases. Some cases with highly negative Monin Obukhov length
remain to be explored.

M 267 libf/phymars/vdif_cd.F
<> Added gustiness to the Richardson computation. Gustiness factor is for now of beta=1., after

several comparisons with LES aerodynamic conductances. May be subject to a minor change (+/- 0.1)
in the near future. (almost transparent for the user)

M 267 libf/phymars/vdifc.F
<> Minor modifications relative to variables.

M 267 libf/phymars/calltherm_mars.F90
<> Added a comment on a "sensitive" parameter that should not be changed without knowing the consequence !

M 267 libf/phymars/meso_inc/meso_inc_les.F
<> Changed the definition for HFX computation in the LES (to be discussed with Aymeric). New definition yields

very similar results to old one and follows a strict definition of what HFX should be.
File size: 1.9 KB
Line 
1         DO ig=1,ngrid
2          !! sensible heat flux in W/m2
3!          hfx(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)
4
5! New SL parametrization, correct formulation for hfx :
6
7          hfx(ig) = (pplay(ig,1)/(r*pt(ig,1)))*cpp
8     &    *sqrt((pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)))
9     &    *zcdh(ig)*(tsurf(ig)-zh(ig,1))
10
11          !! u star in similarity theory in m/s
12!          ust(ig) = 0.4
13!     .               * sqrt( pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) )
14!     .               / log( 1.E+0 + zzlay(ig,1)/z0_default )
15
16! New SL parametrization, ust is more accurately computed in vdif_cd :
17        ust(ig) = sqrt(zcdv(ig)*(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)))
18
19         ENDDO   
20!         write (*,*) 'PHYS HFX cp zdts', hfx(100), zflubid(100),
21!     .       capcal(100),
22!     .       zdtsdif(100)
23!         write (*,*) 'PHYS UST', ust(100)
24
25!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
26!!! LES LES
27       IF (flag_LES) THEN       
28
29         write (*,*) '************************************************'
30         write (*,*) '** LES mode: the difv part is only used to'
31         write (*,*) '**  provide HFX and UST to the dynamics'
32         write (*,*) '** NB: - dudif, dvdif, dhdif, dqdif are set to 0'
33         write (*,*) '**     - tsurf is updated'     
34         write (*,*) '************************************************'
35
36         DO ig=1,ngrid
37          wmax_th(ig)=0.  !! no additional gustiness needed in surface layer (see vdifc.F)
38          DO l=1,nlayer
39            zdvdif(ig,l) = 0.
40            zdudif(ig,l) = 0.
41            zdhdif(ig,l) = 0.
42            DO iq=1, nq
43              zdqdif(ig,l,iq) = 0.
44              zdqsdif(ig,iq) = 0. !! sortir de la boucle
45            ENDDO
46          ENDDO
47         ENDDO
48      ENDIF
49!!! LES LES       
50!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Note: See TracBrowser for help on using the repository browser.