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

Last change on this file since 250 was 226, checked in by aslmd, 14 years ago

MESOSCALE/LMDZ.MARS.new
--> modified to impact last changes

MESOSCALE/LMD_MM_MARS/makemeso
MESOSCALE/LMD_MM_MARS/SRC/WRFV2/call_meso_physiq?.inc
MESOSCALE/LMD_MM_MARS/SRC/WRFV2/call_meso_inifis?.inc
MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
--> modified to get rid of ecri_phys

and make changes related to meso_physiq and meso_inifis

LMDZ.MARS/libf/phymars
--> see LMDZ.MARS/README

15/07/2011 == AS

  • Modified the mesoscale part so that the previous change by EM does not imply an error in the mesoscale case. More development is needed though to get the "varying z0" capability in the mesoscale model.
  • Worked on versions of meso_physiq and meso_inifis as close as possible to physiq and inifis for more continuity in the process of impacting changes (and even possibly to reach a common version of physiq and inifis).

    The main point is to make the mesoscale significant specific parts

    coded into include files in meso_inc so that meso_physiq and meso_inifis looks very close to physiq and inifis.

    This is completely transparent for GCM users who does not need the

    contents of meso_inc.
  • Slight cosmetic changes to physiq.f and inifis.F --- some of them e.g. to prepare convergence between meso_physiq and physiq
File size: 1.4 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         ENDDO   
9
10!         write (*,*) 'PHYS HFX cp zdts', hfx(100), zflubid(100),
11!     .       capcal(100),
12!     .       zdtsdif(100)
13!         write (*,*) 'PHYS UST', ust(100)
14
15!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
16!!! LES LES
17       IF (flag_LES) THEN       
18
19         write (*,*) '************************************************'
20         write (*,*) '** LES mode: the difv part is only used to'
21         write (*,*) '**  provide HFX and UST to the dynamics'
22         write (*,*) '** NB: - dudif, dvdif, dhdif, dqdif are set to 0'
23         write (*,*) '**     - tsurf is updated'     
24         write (*,*) '************************************************'
25
26         DO ig=1,ngrid
27          DO l=1,nlayer
28            zdvdif(ig,l) = 0.
29            zdudif(ig,l) = 0.
30            zdhdif(ig,l) = 0.
31            DO iq=1, nq
32              zdqdif(ig,l,iq) = 0.
33              zdqsdif(ig,iq) = 0. !! sortir de la boucle
34            ENDDO
35          ENDDO
36         ENDDO
37      ENDIF
38!!! LES LES       
39!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Note: See TracBrowser for help on using the repository browser.