Changeset 593 for trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90

Timestamp:

Mar 20, 2012, 2:27:39 PM (13 years ago)

Author:

acolaitis

Message:

Minor change to internal parameters of the thermal: smoothed detrainment profiles, reduced the contribution of upward subsidence to entrainment drag, and increased pressure torque induced drag. Consequences : more precise prediction of zmax (it was overestimated by 5 to 10 percent during the late afternoon), no consequence on the thermal structure, better predictions for internal variables (much better w profile and detrainment rate profile). Tested on both high resolution and low resolution (33 level grid from me and 36 level grid from Forget)

File:

• trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90 (modified) (5 diffs)

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

@@ -401 +401 @@
 ! Best config for norad 222 levels:
 
-       omega=0.06
+!       omega=0.06
+       omega=0.04
        a1=1.
-       b1=0.
+!       b1=0.
+       b1=0.0001
        a1inv=a1
        be=1.1*be
@@ -409 +411 @@
        adalim=0.
        b1inv=0.00025
+!       b1inv = 0.0003
 
 !      omega=0.06
@@ -621 +624 @@
           else
           detr_star(ig,l)=f_star(ig,l)*zdz*                        &
-                &     bd*zbuoy(ig,l)/zw2m
+                &     MAX(ad,bd*zbuoy(ig,l)/zw2m)
 !             &  Max(0., 0.001 - 0.45*zbuoy(ig,l)/zw2m)
 !             &  Max(0., Min(0.001,0.0005 - 0.55*zbuoy(ig,l)/zw2m))
@@ -660 +663 @@
 !---------------------------------------------------------------------------
 
-      DO tic=0,0  ! internal convergence loop
+      DO tic=0,5  ! internal convergence loop
       activetmp(:)=activecell(:) .and. f_star(:,l+1)>1.e-10
       do ig=1,ngridmx
@@ -718 +721 @@
           else
           detr_star(ig,l)=f_star(ig,l)*zdz*                        &
-                &     bd*zbuoy(ig,l)/zw2m
+                &     MAX(ad,bd*zbuoy(ig,l)/zw2m)
 !             &  Max(0.,Min(0.001,0.0005 - 0.55*zbuoy(ig,l)/zw2m))