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

Timestamp:

Feb 15, 2012, 5:55:59 PM (13 years ago)

Author:

acolaitis

Message:

Most up-to-date thermals parameters. Optimized for high resolution runs.

File:

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

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

@@ -113 +113 @@
       REAL zdz,zbuoy(ngridmx,nlayermx),zw2m
       LOGICAL activecell(ngridmx),activetmp(ngridmx)
-      REAL a1,b1,ae,be,ad,bd,fdfu,b1inv,a1inv
+      REAL a1,b1,ae,be,ad,bd,fdfu,b1inv,a1inv,omega
       INTEGER tic
 
@@ -378 +378 @@
       ad = 0.00048088  ; bd = -0.6697
       fdfu = -1.3
+      a1inv=a1
+      b1inv=b1
+      omega=0.
+
+! One good config for 34/35 levels
+!      a1inv=a1
+!      b1inv=b1
+!      be=1.1*be
 
 ! Best configuration for 222 levels:
+
+      omega=0.06
       b1=0.
-      a1=0.7*a1
+      a1=1.
       a1inv=0.25*a1
       b1inv=0.0002
 
+! Best config for norad 222 levels:
+
+!       omega=0.06
+!       a1=1.
+!       b1=0.
+!       a1inv=a1
+!       be=1.1*be
+!       ad = 0.0004
+!       b1inv=0.00035
+
+! Trying stuff :
+
+!       omega=0.06
+!       a1=1.
+!       b1=0.
+!       be=be
+!       a1inv=0.25*a1
+!       b1inv=0.0002    
+!       ad=1.1*ad
+!       ae=1.*ae
 ! --------------------------------------------------------------------------
 ! --------------------------------------------------------------------------
@@ -503 +533 @@
                 if (((a1*zbuoy(ig,l)/w_est(ig,l)-b1) .gt. 0.) .and. (w_est(ig,l) .ne. 0.)) then
                 w_est(ig,l+1)=Max(0.0001,w_est(ig,l)+2.*zdz*a1*zbuoy(ig,l)-2.*zdz*w_est(ig,l)*b1 &
-     & -2.*zdz*w_est(ig,l)*ae*(a1*zbuoy(ig,l)/w_est(ig,l)-b1)**be)
+     & -2.*(1.-omega)*zdz*w_est(ig,l)*ae*(a1*zbuoy(ig,l)/w_est(ig,l)-b1)**be)
                 else
                 w_est(ig,l+1)=Max(0.0001,w_est(ig,l)+2.*zdz*a1inv*zbuoy(ig,l)-2.*zdz*w_est(ig,l)*b1inv)
@@ -598 +628 @@
 !---------------------------------------------------------------------------
 
-      DO tic=0,5  ! internal convergence loop
+      DO tic=0,0  ! internal convergence loop
       activetmp(:)=activecell(:) .and. f_star(:,l+1)>1.e-10
       do ig=1,ngridmx
@@ -619 +649 @@
            if (((a1*zbuoy(ig,l)/zw2(ig,l)-b1) .gt. 0.) .and. (zw2(ig,l) .ne. 0.) ) then
            zw2(ig,l+1)=Max(0.,zw2(ig,l)+2.*zdz*a1*zbuoy(ig,l)-         &
-     & 2.*zdz*zw2(ig,l)*b1-2.*zdz*zw2(ig,l)*ae*(a1*zbuoy(ig,l)/zw2(ig,l)-b1)**be)
+     & 2.*zdz*zw2(ig,l)*b1-2.*(1.-omega)*zdz*zw2(ig,l)*ae*(a1*zbuoy(ig,l)/zw2(ig,l)-b1)**be)
            else
            zw2(ig,l+1)=Max(0.,zw2(ig,l)+2.*zdz*a1inv*zbuoy(ig,l)-2.*zdz*zw2(ig,l)*b1inv)