Changeset 94 in lmdz_wrf

Timestamp:

Jul 25, 2014, 2:41:55 PM (10 years ago)

Author:

lfita

Message:

NaNs? from zw2...

Location:

branches/LMDZ_WRFmeas/WRFV3/lmdz

Files:

• thermcell_main.F90 (modified) (1 diff)
• thermcell_plume.F90 (modified) (4 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

@@ -516 +516 @@
       lvarname = 'pdtadj'
       CALL check_var3D(lfname, lvarname, pdtadj, ngrid, nlay, largest, .FALSE.)
+      lvarname = 'zw2'
+      CALL check_var3D(lfname, lvarname, zw2, ngrid, nlay+1, largest, .FALSE.)
 
       if (prt_level.ge.1) print*,'apres thermcell_plume ',lev_out

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_plume.F90

@@ -90 +90 @@
       REAL zw2fact,expa
       Zsat=.false.
+
+! Lluis
+      INTEGER                                            :: llp
+      CHARACTER(LEN=50)                                  :: lvarname, lfname
+      REAL                                               :: largest
+
+      llp = 734
+      lfname = 'physiq'
+      largest = 10.e5
+
 ! Initialisation
       RLvCp = RLVTT/RCPD
@@ -193 +203 @@
 !
 
+      lfname='thermcell_plume after k=2'
+      lvarname = 'zw2'
+      CALL check_var3D(lfname, lvarname, zw2, ngrid, klev+1, 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.)
+      lvarname = 'pphi'
+      CALL check_var3D(lfname, lvarname, pphi, ngrid, klev, largest, .FALSE.)
+
+
 !==============================================================================
 !boucle de calcul de la vitesse verticale dans le thermique
@@ -427 +448 @@
       enddo
 !=========================================================================
+
+  lfname='thermcell_plume end of vertical loop'
+  lvarname = 'zw2'
+  CALL check_var3D(lfname, lvarname, zw2, ngrid, klev+1, 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.)
+  lvarname = 'pphi'
+  CALL check_var3D(lfname, lvarname, pphi, ngrid, klev, largest, .FALSE.)
+
 
 !on recalcule alim_star_tot
@@ -686 +718 @@
 
 !------------------------------------------------
-!AJAM:nouveau calcul de w²  
+!AJAM:nouveau calcul de w\B2  
 !------------------------------------------------
               zdz=zlev(ig,l+1)-zlev(ig,l)