Changeset 341 for trunk/LMDZ.MARS/libf/phymars/meso_inc
- Timestamp:
- Nov 3, 2011, 2:28:57 PM (13 years ago)
- Location:
- trunk/LMDZ.MARS/libf/phymars/meso_inc
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F
r330 r341 4 4 !! sensible heat flux in W/m2 5 5 6 hfx(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)6 sensheat(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig) 7 7 8 8 !! u star in similarity theory in m/s 9 ust (ig) = 0.49 ustar(ig) = 0.4 10 10 . * sqrt( pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) ) 11 11 . / log( 1.E+0 + zzlay(ig,1)/z0_default ) … … 16 16 DO ig=1,ngrid 17 17 18 ! New SL parametrization, correct formulation for hfx:18 ! New SL parametrization, correct formulation for sensheat : 19 19 20 hfx(ig) = (pplay(ig,1)/(r*pt(ig,1)))*cpp20 sensheat(ig) = (pplay(ig,1)/(r*pt(ig,1)))*cpp 21 21 & *sqrt(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) 22 22 & + (1.2*wmax_th(ig))**2) 23 23 & *zcdh(ig)*(tsurf(ig)-zh(ig,1)) 24 24 25 ! New SL parametrization, ust is more accurately computed in vdif_cd :26 ust (ig) = sqrt(zcdv(ig)*25 ! New SL parametrization, ustar is more accurately computed in vdif_cd : 26 ustar(ig) = sqrt(zcdv(ig)* 27 27 & (pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) + (0.7*wmax_th(ig))**2) 28 28 & ) … … 32 32 endif !of if callrichsl 33 33 34 ! write (*,*) 'PHYS HFX cp zdts', hfx(100), zflubid(100),34 ! write (*,*) 'PHYS HFX cp zdts', sensheat(100), zflubid(100), 35 35 ! . capcal(100), 36 36 ! . zdtsdif(100) 37 ! write (*,*) 'PHYS UST', ust (100)37 ! write (*,*) 'PHYS UST', ustar(100) 38 38 39 39 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! -
trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_var.F
r308 r341 10 10 REAL output_tab2d(ngridmx,n2d) 11 11 REAL output_tab3d(ngridmx,nlayer,n3d) 12 REAL hfx(ngridmx) !! pour LES avec isfflx!=013 REAL ust (ngridmx) !! pour LES avec isfflx!=012 REAL sensheat(ngridmx) !! pour LES avec isfflx!=0 13 REAL ustar(ngridmx) !! pour LES avec isfflx!=0 14 14 LOGICAL flag_LES !! pour LES avec isfflx!=0 15 15 REAL qsurfice(ngridmx) !! pour diagnostics
Note: See TracChangeset
for help on using the changeset viewer.