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

Timestamp:

Mar 2, 2012, 4:46:26 PM (13 years ago)

Author:

acolaitis

Message:

Yamada4 is now ON when using thermals. Advised setup with thermals: modified z2sig (one more level in PBL) and at least 72 timesteps per day (for now...).

File:

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

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

@@ -410 +410 @@
        b1inv=0.00025
 
-
-!      b1=0.0007
 !      omega=0.06
 ! Trying stuff :
@@ -456 +454 @@
             if (ztv(ig,1)>=(ztv(ig,2))) then
                alim_star(ig,1)=MAX((ztv(ig,1)-ztv(ig,2)),0.)  &
-!     &                       *sqrt(zlev(ig,2))
+!     &                       *log(1.+zlev(ig,2))
+     &                       *sqrt(zlev(ig,2))
+!     &                       *sqrt(sqrt(zlev(ig,2)))
 !     &                       /sqrt(zlev(ig,2))
-      &                       *zlev(ig,2)
+!      &                       *zlev(ig,2)
 !      &                     *exp(-zlev(ig,2)/1000.)
                lalim(ig)=2
@@ -470 +470 @@
            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))
+!     &                       *log(1.+zlev(ig,l+1))
+     &                       *sqrt(zlev(ig,l+1))
+!     &                       *sqrt(sqrt(zlev(ig,l+1)))
 !     &                       /sqrt(zlev(ig,l+1))
-      &                       *zlev(ig,l+1)
+!      &                       *zlev(ig,l+1)
 !      &                     *exp(-zlev(ig,l+1)/1000.)
                 lalim(ig)=l+1