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

Timestamp:

Feb 27, 2012, 10:44:32 AM (13 years ago)

Author:

acolaitis

Message:

27/02/12 == AC

Continuation of thermals setting, comparisons with mesoscale results on Case C
Added possibility to call gcm (or 1d) with constant prescribed sensible heat flux, in the spirit of direct comparisons with LES

... This is directly comparable to the variable tke_heat_flux in namelist.input
... Requires the use of yamada4.F from terrestrial GCM (mixes more, seem more numerically stable)
... Usually requires high timesteps (>1000) to avoids crashes. Best approach is to compare
height of first model level z1 and teta_1 between LES and 1D, and increase the timesteps until results
between the two models are comparable (might require a slitghly different tke_heat_flux between the two models
due to difference in vertical diffusion schemes and subgrid effects)
... Syntax for use is to add "tke_heat_flux = ..." in callphys.def

Corrected some stuff with tke transport in thermals (which is anyway deactivated but one day maybe...)

File:

• trunk/LMDZ.MARS/libf/phymars/thermcell_main_mars.F90 (modified) (8 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,omega
+      REAL a1,b1,ae,be,ad,bd,fdfu,b1inv,a1inv,omega,adalim
       INTEGER tic
 
@@ -381 +381 @@
       b1inv=b1
       omega=0.
+      adalim=0.
 
 ! One good config for 34/35 levels
@@ -389 +390 @@
 ! Best configuration for 222 levels:
 
-      omega=0.06
-      b1=0.
-      a1=1.
-      a1inv=0.25*a1
-      b1inv=0.0002
+!      omega=0.06
+!      b1=0.
+!      a1=1.
+!      a1inv=0.25*a1
+!      b1inv=0.0002
+!!
+!!
+!!      ae=0.9*ae
 
 ! 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
+! with yamada4 and alim at sqrt(zlev)
+
+       omega=0.06
+       a1=1.
+       b1=0.
+       a1inv=a1
+       be=1.1*be
+       ad = 0.0004
+       b1inv=0.00035
+       adalim=0.
+
+      b1inv=0.00025
 
 ! Trying stuff :
+
+!      omega=0.04
+!!      b1=0.
+!      a1=1.
+!      a1inv=a1
+!      b1inv=0.0005689
+!!      be=1.1*be
+!!      ae=0.96*ae
+
 
 !       omega=0.06
@@ -431 +448 @@
             if (ztv(ig,1)>=(ztv(ig,2))) then
                alim_star(ig,1)=MAX((ztv(ig,1)-ztv(ig,2)),0.)  &
-     &                       *sqrt(zlev(ig,2))
+!     &                       *sqrt(zlev(ig,2))
 !     &                       /sqrt(zlev(ig,2))
-!      &                       *zlev(ig,2)
+      &                       *zlev(ig,2)
                lalim(ig)=2
                alim_star_tot(ig)=alim_star_tot(ig)+alim_star(ig,1)
@@ -444 +461 @@
            if (ztv(ig,l)>(ztv(ig,l+1)+0.) .and. ztv(ig,1)>=ztv(ig,l) .and. (alim_star(ig,l-1) .ne. 0.)) then ! .and. (zlev(ig,l+1) .lt. 1000.)) then
                alim_star(ig,l)=MAX((ztv(ig,l)-ztv(ig,l+1)),0.)  &
-     &                       *sqrt(zlev(ig,l+1))
-!      &                       *zlev(ig,2)
+!     &                       *sqrt(zlev(ig,l+1))
+!     &                       /sqrt(zlev(ig,l+1))
+      &                       *zlev(ig,l+1)
                 lalim(ig)=l+1
                alim_star_tot(ig)=alim_star_tot(ig)+alim_star(ig,l)
@@ -562 +580 @@
           if(zbuoy(ig,l) .gt. 0.) then
              if(l .lt. lalim(ig)) then
-                detr_star(ig,l)=0.
+
+!                detr_star(ig,l)=0.
+                 detr_star(ig,l) = f_star(ig,l)*zdz*              &
+            &  adalim
              else
 
@@ -628 +649 @@
 !---------------------------------------------------------------------------
 
-      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
@@ -674 +695 @@
           if(zbuoy(ig,l) .gt. 0.) then
              if(l .lt. lalim(ig)) then
-                detr_star(ig,l)=0.
+
+!                detr_star(ig,l)=0.
+                 detr_star(ig,l) = f_star(ig,l)*zdz*              &
+            &  adalim
+
              else
                  detr_star(ig,l) = f_star(ig,l)*zdz*              &