Changeset 1242 for trunk

Timestamp:

May 10, 2014, 10:03:34 AM (11 years ago)

Author:

aslmd

Message:

LMDZ.MARS. MESOSCALE. Corrected problems affecting LES runs with new physics (sensibFlux calculated in vdifc is just used as an input to thermals, it does not have the right sign and apparently is not W/m2). Further reduction or factorisation of mesoscale includes. Also made ustar and tstar available in turb_mod from the calculations in vdifc.

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/phymars/meso_inc/meso_inc_les.F (deleted)
• libf/phymars/physiq.F (modified) (4 diffs)
• libf/phymars/turb_mod.F90 (modified) (2 diffs)
• libf/phymars/vdifc.F (modified) (5 diffs)

trunk/LMDZ.MARS/README

@@ -2085 +2085 @@
 - Warnings in conf_phys modified accordingly.
 
+== 10/05/2014 == AS
+- Further reduction (or 'factorisation') of mesoscale includes
+- ustar and tstar are calculated in vdifc. so no need to recalculate
+  those elsewhere. use turb_mod module to store it. 

trunk/LMDZ.MARS/libf/phymars/physiq.F

@@ -838 +838 @@
      &        zdqdif,zdqsdif,wstar,zcdv,zcdh,hfmax_th,sensibFlux)
 
-#ifdef MESOSCALE
-#include "meso_inc/meso_inc_les.F"
-#endif
 
           DO ig=1,ngrid
@@ -871 +868 @@
            ENDDO
           end if ! of if (tracer)
+         ELSE
+           write (*,*) '******************************************'
+           write (*,*) '** LES mode: the difv part is only used to'
+           write (*,*) '**  - provide HFX and UST to the dynamics'
+           write (*,*) '**  - update TSURF'
+           write (*,*) '******************************************'
+           !! Specific treatment for lifting in turbulent-resolving mode (AC)
+           IF (lifting .and. doubleq) THEN
+             !! lifted dust is injected in the first layer. 
+             !! Sedimentation must be called after turbulent mixing, i.e. on next step, after WRF. 
+             !! => lifted dust is not incremented before the sedimentation step.
+             zdqdif(1:ngrid,1,1:nq)=0.
+             zdqdif(1:ngrid,1,igcm_dust_number) = 
+     .                  -zdqsdif(1:ngrid,igcm_dust_number)
+             zdqdif(1:ngrid,1,igcm_dust_mass) = 
+     .                  -zdqsdif(1:ngrid,igcm_dust_mass)
+             zdqdif(1:ngrid,2:nlayermx,1:nq) = 0.
+             DO iq=1, nq
+              IF ((iq .ne. igcm_dust_mass)
+     &        .and. (iq .ne. igcm_dust_number)) THEN
+                zdqsdif(:,iq)=0.
+              ENDIF
+             ENDDO
+           ELSE
+             zdqdif(1:ngrid,1:nlayer,1:nq) = 0.
+             zdqsdif(1:ngrid,1:nq) = 0.
+           ENDIF
          ENDIF
-
       ELSE    
          DO ig=1,ngrid
             zdtsurf(ig)=zdtsurf(ig)+
-     s      (fluxrad(ig)+fluxgrd(ig))/capcal(ig)
+     s        (fluxrad(ig)+fluxgrd(ig))/capcal(ig)
          ENDDO
          IF (turb_resolved) THEN
@@ -1483 +1506 @@
      & g,zzlay,zzlev,zu,zv,wstar,hfmax_th,zmax_th,tsurf,zh,z_out,n_out,
      & T_out,u_out,ustar,tstar,L_mo,vhf,vvv)
+                   ! pourquoi ustar recalcule ici? fait dans vdifc.
 
 #ifndef MESOSCALE
@@ -1909 +1933 @@
       !state  real  RICE      ikj   misc  1  -  h  "RICE"      "ICE RADIUS"                      "m"
       comm_RICE(1:ngrid,1:nlayer) = rice(1:ngrid,1:nlayer)
+
+   
+      !! calculate sensible heat flux in W/m2 for outputs
+      !! -- the one computed in vdifc is not the real one
+      !! -- vdifc must have been called
+      if (.not.callrichsl) then
+        sensibFlux(1:ngrid) = zflubid(1:ngrid)
+     .         - capcal(1:ngrid)*zdtsdif(1:ngrid)
+      else
+        sensibFlux(1:ngrid) = 
+     &   (pplay(1:ngrid,1)/(r*pt(1:ngrid,1)))*cpp
+     &   *sqrt(pu(1:ngrid,1)*pu(1:ngrid,1)+pv(1:ngrid,1)*pv(1:ngrid,1)
+     &         +(log(1.+0.7*wstar(1:ngrid) + 2.3*wstar(1:ngrid)**2))**2)
+     &   *zcdh(1:ngrid)*(tsurf(1:ngrid)-zh(1:ngrid,1))
+      endif
+          
 #else
 #ifndef MESOINI

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

@@ -6 +6 @@
   REAL,SAVE,ALLOCATABLE :: ustar(:)
   REAL,SAVE,ALLOCATABLE :: wstar(:)
+  REAL,SAVE,ALLOCATABLE :: tstar(:)
   REAL,SAVE,ALLOCATABLE :: hfmax_th(:)
   REAL,SAVE,ALLOCATABLE :: zmax_th(:)
@@ -25 +26 @@
     allocate(wstar(ngrid))
     allocate(ustar(ngrid))
+    allocate(tstar(ngrid))
     allocate(hfmax_th(ngrid))
     allocate(zmax_th(ngrid))

trunk/LMDZ.MARS/libf/phymars/vdifc.F

@@ -13 +13 @@
       use surfdat_h, only: watercaptag, frost_albedo_threshold, dryness
       USE comcstfi_h
-      use turb_mod, only: turb_resolved
+      use turb_mod, only: turb_resolved, ustar, tstar
       IMPLICIT NONE
 
@@ -84 +84 @@
 
       REAL :: pt(ngrid,nlay)
-      REAL ust(ngrid),tst(ngrid)
  
       INTEGER ilev,ig,ilay,nlev
@@ -335 +334 @@
      &     (log(1.+0.7*wstar(:) + 2.3*wstar(:)**2))**2)
 
-           ust(:)=sqrt(zcdv_true(:))*sqrt(zu2(:)+
+           ustar(:)=sqrt(zcdv_true(:))*sqrt(zu2(:)+
      &     (log(1.+0.7*wstar(:) + 2.3*wstar(:)**2))**2)
 
-           tst(:)=0.
+           tstar(:)=0.
            DO ig=1,ngrid
               IF (zcdh_true(ig) .ne. 0.) THEN      ! When Cd=Ch=0, u*=t*=0
-                 tst(ig)=(ph(ig,1)-ptsrf(ig))*zcdh(ig)/ust(ig)
+                 tstar(ig)=(ph(ig,1)-ptsrf(ig))*zcdh(ig)/ustar(ig)
               ENDIF
            ENDDO
@@ -348 +347 @@
            zcdv(:)=zcdv_true(:)*sqrt(zu2(:))     ! 1 / bulk aerodynamic momentum conductance 
            zcdh(:)=zcdh_true(:)*sqrt(zu2(:))     ! 1 / bulk aerodynamic heat conductance
-           ust(:)=sqrt(zcdv_true(:))*sqrt(zu2(:))
-           tst(:)=(ph(:,1)-ptsrf(:))*zcdh_true(:)/sqrt(zcdv_true(:))
+           ustar(:)=sqrt(zcdv_true(:))*sqrt(zu2(:))
+           tstar(:)=(ph(:,1)-ptsrf(:))*zcdh_true(:)/sqrt(zcdv_true(:))
         ENDIF
 
@@ -383 +382 @@
       pt(:,:)=ph(:,:)*ppopsk(:,:)
       CALL yamada4(ngrid,nlay,nq,ptimestep,g,r,pplev,pt
-     s   ,pzlev,pzlay,pu,pv,ph,pq,zcdv_true,pq2,zkv,zkh,zkq,ust
+     s   ,pzlev,pzlay,pu,pv,ph,pq,zcdv_true,pq2,zkv,zkh,zkq,ustar
      s   ,9)
       ENDIF