Changeset 341 for trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf

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.

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars

Files:

• physiq.F (modified) (4 diffs)
• readtesassim.F90.modif (added)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/physiq.F

@@ -178 +178 @@
       REAL sl_ra, sl_di0
       REAL sky
-      REAL hfx(ngridmx)    !! pour LES avec isfflx!=0
-      REAL ust(ngridmx)    !! pour LES avec isfflx!=0
+      REAL sensheat(ngridmx)    !! pour LES avec isfflx!=0
+      REAL ustar(ngridmx)    !! pour LES avec isfflx!=0
       LOGICAL flag_LES     !! pour LES avec isfflx!=0
       REAL qsurflast(ngridmx) !! pour diagnostics
@@ -925 +925 @@
          DO ig=1,ngrid
           !! 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 )
@@ -944 +944 @@
          write (*,*) '************************************************'
 
-!         DO ig=1,ngrid
-!          !! sensible heat flux in W/m2
-!          hfx(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)
-!          !! u star in similarity theory in m/s
-!          ust(ig) = 0.4
-!     .               * sqrt( pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) )
-!     .               / log( 1.E+0 + zzlay(ig,1)/z0 )
 !
           DO l=1,nlayer
@@ -966 +959 @@
          !write (*,*) 'GRD ',fluxgrd(igout)
          !write (*,*) 'dTs/dt ',capcal(igout)*zdtsurf(igout)
-         !write (*,*) 'HFX ', hfx(igout)
-         !write (*,*) 'UST ', ust(igout)
       ENDIF
 !!! LES LES