Changeset 277

Timestamp:

Sep 5, 2011, 12:05:47 PM (13 years ago)

Author:

acolaitis

Message:

MESOSCALE: comments and minor modif for use of thermals.

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• calltherm_mars.F90 (modified) (1 diff)
• physiq.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/calltherm_mars.F90

@@ -95 +95 @@
 
 #ifdef MESOSCALE
+         !! valid for timesteps < 200s 
          nsplit_thermals=2
          r_aspect_thermals=0.7

trunk/LMDZ.MARS/libf/phymars/physiq.F

@@ -788 +788 @@
         endif
 
+        lmax_th_out(:)=real(lmax_th(:))
+
         else   !of if calltherm
         lmax_th(:)=0
         wmax_th(:)=0.
+        lmax_th_out(:)=0.
         end if
 
@@ -1631 +1634 @@
 !        endif
 
-        lmax_th_out(:)=real(lmax_th(:))
-
         call WRITEDIAGFI(ngridmx,'lmax_th',
      &              'hauteur du thermique','K',
@@ -1725 +1726 @@
          if(calltherm) then
 
-        lmax_th_out(:)=real(lmax_th(:))
-
         call WRITEDIAGFI(ngridmx,'lmax_th',
      &              'hauteur du thermique','point',
@@ -1755 +1754 @@
          call WRITEDIAGFI(ngrid,"tsurf","Surface temperature","K",0,
      &                  tsurf)
-!         call WRITEDIAGFI(ngrid,"u","u wind","m/s",1,zu)
-!         call WRITEDIAGFI(ngrid,"v","v wind","m/s",1,zv)
+         call WRITEDIAGFI(ngrid,"u","u wind","m/s",1,zu)
+         call WRITEDIAGFI(ngrid,"v","v wind","m/s",1,zv)
 
          call WRITEDIAGFI(ngrid,"pplay","Pressure","Pa",1,zplay)