Changeset 190

Timestamp:

Jul 4, 2011, 11:04:57 AM (14 years ago)

Author:

acolaitis

Message:

04/07/2011 == AC

• Minor setting modification to thermals
• Added new flux optimization in thermcell_main_mars.F90, to run properly

in 3D

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (2 diffs)
• libf/phymars/calltherm_mars.F90 (modified) (1 diff)
• libf/phymars/thermcell_main_mars.F90 (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -764 +764 @@
   boundaries for soil temperature.
 
-== 17/06/2011 == AC
+== 01/07/2011 == AC
 - Added new settings for the Martian thermals from new LES observations
 - Revamped thermcell's module variables to allow it's removal
@@ -770 +770 @@
 - Switched from dynamic to static memory allocation for all thermals variable 
 to gain computation speed
+
+== 04/07/2011 == AC
+- Minor setting modification to thermals
+- Added new flux optimization in thermcell_main_mars.F90, to run properly
+in 3D

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

@@ -87 +87 @@
 !  ===================
  
-         r_aspect_thermals=3.
+         r_aspect_thermals=2.
          nsplit_thermals=40
          call getin("nsplit_thermals",nsplit_thermals)

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

@@ -865 +865 @@
 
          do ig=1,ngridmx
+
             if (fm(ig,l+1).lt.0.) then
-              print*,'fm1<0',l+1,lmax(ig),fm(ig,l+1)
-                ncorecfm1=ncorecfm1+1
+               if((l+1) .eq. lmax(ig)) then
+               detr(ig,l)=detr(ig,l)+fm(ig,l+1)
+               fm(ig,l+1)=0.
+               ncorecfm2=ncorecfm2+1
+               else
+          print*,'fm(l+1)<0 : ig, l+1,lmax :',ig,l+1,lmax(ig),fm(ig,l+1)
+               ncorecfm1=ncorecfm1+1
                fm(ig,l+1)=fm(ig,l)
                detr(ig,l)=entr(ig,l)
-            endif
+               endif
+            endif
+
          enddo