Changeset 165 for trunk/LMDZ.MARS/libf

Timestamp:

Jun 17, 2011, 3:08:24 PM (13 years ago)

Author:

acolaitis

Message:

17/06/2011 == AC

Updates to thermal parameters

• Tuned aspect ratio of thermals to suit Buoyancy estimations from LES in CLOSURE relation
• Renormalization of alim_star after plume
• Removed alimentation mixing of estimated Teta in plume

  Minor change in makegcm_ifort
  

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• calltherm_interface.F90 (modified) (2 diffs)
• thermcell.F90 (modified) (1 diff)
• thermcell_main_mars.F90 (modified) (7 diffs)

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

@@ -103 +103 @@
 
             pq_therm(:,:,:)=0.
+            qtransport_thermals=.true.
             call getin("qtransport_thermals",qtransport_thermals)
             if(qtransport_thermals) then
@@ -119 +120 @@
             buoyancyEst(:,:)=0.
 
+       dtke_thermals=.false.
        call getin("dtke_thermals",dtke_thermals)
          if(dtke_thermals) then

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

@@ -72 +72 @@
 ! iflag_clos : USELESS ?                                  [useless ?]
 ! iflag_wake : USELESS ?                                  [useless ?]
-      PARAMETER (r_aspect_thermals = 3.)
+      PARAMETER (r_aspect_thermals = 1.5)
       PARAMETER (l_mix_thermals = 30.)
       PARAMETER (w2di_thermals = 1)

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

@@ -324 +324 @@
                alim_star(ig,1)=MAX((ztv(ig,1)-ztv(ig,2)),0.)  &
      &                       *sqrt(zlev(ig,2))
+!      &                       *zlev(ig,2)
                lalim(ig)=2
                alim_star_tot(ig)=alim_star_tot(ig)+alim_star(ig,1)
@@ -334 +335 @@
                alim_star(ig,l)=MAX((ztv(ig,l)-ztv(ig,l+1)),0.)  &
      &                       *sqrt(zlev(ig,l+1))
+!      &                       *zlev(ig,2)
                 lalim(ig)=l+1
                alim_star_tot(ig)=alim_star_tot(ig)+alim_star(ig,l)
@@ -348 +350 @@
 
       alim_star_tot(:)=1.
+      if(alim_star(1,1) .ne. 0.) then
+      print*, alim_star(:,:)
+      endif
 
 !------------------------------------------------------------------------------
@@ -403 +408 @@
              if(active(ig)) then
 
-                if(l .lt. lalim(ig)) then
-                ztva_est(ig,l)=(f_star(ig,l)*ztla(ig,l-1)+  &
-     &            alim_star(ig,l)*ztv(ig,l))  &
-     &            /(f_star(ig,l)+alim_star(ig,l))
-                else
+!                if(l .lt. lalim(ig)) then
+!                ztva_est(ig,l)=(f_star(ig,l)*ztla(ig,l-1)+  &
+!     &            alim_star(ig,l)*ztv(ig,l))  &
+!     &            /(f_star(ig,l)+alim_star(ig,l))
+!                else
                 ztva_est(ig,l)=ztla(ig,l-1)
-                endif
+!                endif
 
                 zdz=zlev(ig,l+1)-zlev(ig,l)
@@ -566 +571 @@
        enddo
 
+      do l=1,nlay
+         do ig=1,ngrid
+            if (alim_star_tot(ig) > 1.e-10 ) then
+               alim_star(ig,l)=alim_star(ig,l)/alim_star_tot(ig)
+            endif
+         enddo
+      enddo
 
 ! ===========================================================================
@@ -1154 +1166 @@
      &   fm(ig,k)*ztv(ig,k)+fm(ig,k)*ztva(ig,k)-fm(ig,k+1)*ztva(ig,k+1))
             if (ztv(ig,k) + ptimestep*zdthladj(ig,k) .le. 0.) then
-              print*,'q<0 in thermcell_dTeta up: qenv .. dq : ', ztv(ig,k),ptimestep*zdthladj(ig,k)
+              print*,'Teta<0 in thermcell_dTeta up: qenv .. dq : ', ztv(ig,k),ptimestep*zdthladj(ig,k)
+              if(ztv(ig,k) .gt. 0.) then
+              zdthladj(ig,k)=0.
+              endif
             endif
          enddo
@@ -1211 +1226 @@
             if (ztv(ig,k) + ptimestep*zdthladj_down(ig,k) .le. 0.) then
               print*,'q<0 in thermcell_dTeta down: qenv .. dq : ', ztv(ig,k),ptimestep*zdthladj_down(ig,k)
+              if(ztv(ig,k) .gt. 0.) then
+              zdthladj(ig,k)=0.
+              endif
             endif
          enddo