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

Timestamp:

Nov 3, 2011, 2:28:57 PM (14 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

Files:

• libf/phymars/physiq.F (modified) (4 diffs)
• libf/phymars/readtesassim.F90.modif (added)
• makegcm_mpi (modified) (2 diffs)

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        

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/makegcm_mpi

@@ -177 +177 @@
 ##   set opt_link=" -Mfree -L/usr/local/pgi/linux86/lib -lpgf90 -lpgftnrtl -lpghpf -lpghpf2 -L$NCDFLIB -lnetcdf -Wl,-Bstatic"
 #   set opt_link=" -Mlarge_arrays -mcmodel=medium -Mfree -L/usr/local/pgi/linux86/lib -lpgf90 -lpgftnrtl -lpghpf -lpghpf2 -L$NCDFLIB -lnetcdf"
-   set optim90=" -O2 -I. -Munroll -Mcache_align" 
-   set optimtru90=" -O2 -I. -Munroll -Mcache_align" 
+   set optim90="-O2 -I. -Munroll -Mcache_align" # -Mbounds" 
+   set optimtru90="-O2 -I. -Munroll -Mcache_align" # -Mbounds" 
+   #set optim90=""#-Mvect=noaltcode -Msmartalloc -Mprefetch=distance:8 -Mfprelaxed" #-fast -I."
+   #set optimtru90=""#-Mvect=noaltcode -Msmartalloc -Mprefetch=distance:8 -Mfprelaxed"#-fast -I."
    #set optim90="-O3 -I. -tp=nehalem-64 -fastsse -Mscalarsse -Mvect=sse -Mipa=fast"
    #set optimtru90="-O3 -I. -tp=nehalem-64 -fastsse -Mscalarsse -Mvect=sse -Mipa=fast"
-   set opt_link=" -L$NCDFLIB -lnetcdf -Bstatic"
+   set opt_link=" -L$NCDFLIB -lnetcdf" # -Bstatic"
    set mod_loc_dir=""
    set mod_suffix=mod
@@ -362 +364 @@
            set optim90="$optim90"" -G1 "
         else if $LINUX then
-           set optim90="-g -Mbounds "
+           set optim90="-g -Mbounds -Ktrap=fp"
 #           set optim90="-mcmodel=medium"
         else