Changeset 108 in lmdz_wrf

Timestamp:

Jul 25, 2014, 5:17:25 PM (10 years ago)

Author:

lfita

Message:

Finding problems in rhobarz

File:

• branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90 (modified) (2 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

@@ -342 +342 @@
 !   --------------------------------------------------------------------
 !
-      lfname='before thermcell_env'
+      lfname='thermcell_main before thermcell_env'
       lvarname = 'pt'
       CALL check_var3D(lfname, lvarname, pt, ngrid, nlay, largest, .FALSE.)
       lvarname = 'pdtadj'
       CALL check_var3D(lfname, lvarname, pdtadj, ngrid, nlay, largest, .FALSE.)
+      lvarname = 'pplev'
+      CALL check_var3D(lfname, lvarname, pplev, ngrid, nlay, largest, .FALSE.)
 
       CALL thermcell_env(ngrid,nlay,po,pt,pu,pv,pplay,  &
      &           pplev,zo,zh,zl,ztv,zthl,zu,zv,zpspsk,zqsat,lev_out)
 
-      lfname='after thermcell_env'
+      lfname='thermcell_main after thermcell_env'
       lvarname = 'pt'
       CALL check_var3D(lfname, lvarname, pt, ngrid, nlay, largest, .FALSE.)
       lvarname = 'pdtadj'
       CALL check_var3D(lfname, lvarname, pdtadj, ngrid, nlay, largest, .FALSE.)
+      lvarname = 'pplev'
+      CALL check_var3D(lfname, lvarname, pplev, ngrid, nlay, largest, .FALSE.)
        
       if (prt_level.ge.1) print*,'thermcell_main apres thermcell_env'
@@ -413 +417 @@
          masse(:,l)=(pplev(:,l)-pplev(:,l+1))/RG
       enddo
+
+      lfname='thermcell_main after initialization'
+      lvarname = 'rhobarz'
+      CALL check_var3D(lfname, lvarname, rhobarz, ngrid, nlay, largest, .FALSE.)
+      lvarname = 'rho'
+      CALL check_var3D(lfname, lvarname, rho, ngrid, nlay, largest, .FALSE.)
+      lvarname = 'zpspsk'
+      CALL check_var3D(lfname, lvarname, zpspsk, ngrid, nlay, largest, .FALSE.)
+      lvarname = 'pplay'
+      CALL check_var3D(lfname, lvarname, pplay, ngrid, nlay, largest, .FALSE.)
+      lvarname = 'ztv'
+      CALL check_var3D(lfname, lvarname, ztv, ngrid, nlay, largest, .FALSE.)
 
       if (prt_level.ge.1) print*,'thermcell_main apres initialisation'