Changeset 341 for trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F

Timestamp:

Nov 3, 2011, 2:28:57 PM (13 years ago)

Author:

aslmd

Message:

MESOSCALE: tests pour faire marcher le modele en parallele sur la ferme. toujours infructueux... toutes les notes incluses et options explorees en commentaire. les options par defaut restent les memes en attendant. ajout de scripts pour compiler NETCDF et MPI. correction d un probleme de Registry et de makemeso pour les runs LES ancienne physique. ajout d un cas test LES phoenix.

File:

• trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F

@@ -4 +4 @@
 !! sensible heat flux in W/m2
 
-             hfx(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)
+             sensheat(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)
 
 !! u star in similarity theory in m/s
-             ust(ig) = 0.4
+             ustar(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 )
@@ -16 +16 @@
             DO ig=1,ngrid
 
-! New SL parametrization, correct formulation for hfx :
+! New SL parametrization, correct formulation for sensheat :
 
-            hfx(ig) = (pplay(ig,1)/(r*pt(ig,1)))*cpp
+            sensheat(ig) = (pplay(ig,1)/(r*pt(ig,1)))*cpp
      &        *sqrt(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)
      &        + (1.2*wmax_th(ig))**2)
      &        *zcdh(ig)*(tsurf(ig)-zh(ig,1))
 
-! New SL parametrization, ust is more accurately computed in vdif_cd :
-            ust(ig) = sqrt(zcdv(ig)*
+! New SL parametrization, ustar is more accurately computed in vdif_cd :
+            ustar(ig) = sqrt(zcdv(ig)*
      &   (pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) + (0.7*wmax_th(ig))**2)
      &                     )
@@ -32 +32 @@
          endif  !of if callrichsl
 
-!         write (*,*) 'PHYS HFX cp zdts', hfx(100), zflubid(100), 
+!         write (*,*) 'PHYS HFX cp zdts', sensheat(100), zflubid(100), 
 !     .       capcal(100), 
 !     .       zdtsdif(100)
-!         write (*,*) 'PHYS UST', ust(100) 
+!         write (*,*) 'PHYS UST', ustar(100) 
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!