Changeset 113 in lmdz_wrf

Timestamp:

Jul 25, 2014, 6:08:11 PM (10 years ago)

Author:

lfita

Message:

It is using thermcellV1_plume!!!! :(

File:

• branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_plume.F90 (modified) (4 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_plume.F90

@@ -592 +592 @@
       REAL fact_gamma,fact_epsilon,fact_gamma2,fact_epsilon2
       REAL c2(ngrid,klev)
+
+! Lluis
+      INTEGER                                            :: llp
+      CHARACTER(LEN=50)                                  :: lvarname, lfname
+      REAL                                               :: largest
+
+      llp = 734
+      lfname = 'physiq'
+      largest = 10.e5
+
       Zsat=.false.
 ! Initialisation
@@ -660 +670 @@
                lalim(ig)=l+1
                alim_star_tot(ig)=alim_star_tot(ig)+alim_star(ig,l)
+            IF (zlev(ig,l+1) < 0. .OR. alim_star_tot(ig) == 0.) THEN
+              PRINT *,'Lluis ',ig,',',l,' negative zlev: ',zlev(ig,l+1),' null: ',alim_star_tot(ig)
+            END IF
             endif
          enddo
@@ -672 +685 @@
       alim_star_tot(:)=1.
 
+      lfname='thermcellV1_plume after computing alim_star'
+      lvarname = 'alim_star'
+      CALL check_var3D(lfname, lvarname, alim_star, ngrid, klev, largest, .FALSE.)
+      lvarname = 'ztv'
+      CALL check_var3D(lfname, lvarname, ztv, ngrid, klev, largest, .FALSE.)
+      lvarname = 'zlev'
+      CALL check_var3D(lfname, lvarname, zlev, ngrid, klev+1, largest, .FALSE.)
 
 
@@ -796 +816 @@
 
 
+      lfname='thermcellV1_plume after k=2'
+      lvarname = 'alim_star'
+      CALL check_var3D(lfname, lvarname, alim_star, ngrid, klev, largest, .FALSE.)
+      lvarname = 'entr_star'
+      CALL check_var3D(lfname, lvarname, entr_star, ngrid, klev, largest, .FALSE.)
+      lvarname = 'detr_star'
+      CALL check_var3D(lfname, lvarname, detr_star, ngrid, klev, largest, .FALSE.)
+      lvarname = 'f_star'
+      CALL check_var3D(lfname, lvarname, f_star, ngrid, klev, largest, .FALSE.)
+
 !----------------------------------------------------------------------------
 !calcul de la vitesse verticale en melangeant Tl et qt du thermique