Changeset 1236 for trunk/LMDZ.MARS/libf/phymars/meso_inc

Timestamp:

May 5, 2014, 11:38:51 AM (11 years ago)

Author:

aslmd

Message:

MESOSCALE. A necessary complement commit to r1234 where a upgraded interface making use of modules was proposed. Completed the new formulation for module_lmd_driver for newphys with improved interface with both ini/bdy conditions and physical parameterizations. Changed the Registry accordingly. Finished changes about I/O with the LMD physics (see LMDZ.MARS/README). Made all those changes compatible for old interface, and LES runs (checked with test cases), as well as old input files. Changed makemeso to account for full flexibility on changin nx ny ntracers nproc with newphys. Cleaned the now obsolete bits of code used in LMD physics shared with the GCM. ----- Everything is now ready to properly code both restart runs and nesting for mesoscale runs with new physics.

Location:

trunk/LMDZ.MARS/libf/phymars/meso_inc

Files:

• meso_inc_caps.F (deleted)
• meso_inc_inifisini.F (deleted)
• meso_inc_inifisinvar.F (deleted)
• meso_inc_inifisvar.F (deleted)
• meso_inc_les.F (modified) (4 diffs)
• meso_inc_lift_les.F (deleted)
• meso_inc_var.F (deleted)

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

@@ -2 +2 @@
        
              DO ig=1,ngrid
-!! sensible heat flux in W/m2
-
-             sensheat(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)
+!!! sensible heat flux in W/m2
+!
+!             sensheat(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)
 
 !! u star in similarity theory in m/s
@@ -16 +16 @@
             DO ig=1,ngrid
 
-! New SL parametrization, correct formulation for sensheat :
-
-            sensheat(ig) = (pplay(ig,1)/(r*pt(ig,1)))*cpp
-     &        *sqrt(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)
-     &        + (log(1.+0.7*wstar(ig) + 2.3*wstar(ig)**2))**2)
-     &        *zcdh(ig)*(tsurf(ig)-zh(ig,1))
+!! New SL parametrization, correct formulation for sensheat :
+!
+!            sensheat(ig) = (pplay(ig,1)/(r*pt(ig,1)))*cpp
+!     &        *sqrt(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1)
+!     &        + (log(1.+0.7*wstar(ig) + 2.3*wstar(ig)**2))**2)
+!     &        *zcdh(ig)*(tsurf(ig)-zh(ig,1))
 
 ! New SL parametrization, ustar is more accurately computed in vdif_cd :
@@ -40 +40 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !!! LES LES 
-       IF (flag_LES) THEN        
+       IF (turb_resolved) THEN        
 
          write (*,*) '************************************************' 
@@ -49 +49 @@
          write (*,*) '************************************************'
 
-         DO ig=1,ngrid
-          wstar(ig)=0.  !! no additional gustiness needed in surface layer (see vdifc.F)
-          DO l=1,nlayer
-            zdvdif(ig,l) = 0.
-            zdudif(ig,l) = 0.
-            zdhdif(ig,l) = 0.
-          ENDDO
-         ENDDO
          IF (lifting .and. doubleq) THEN 
          !! lifted dust is injected in the first layer.