Changeset 72 in lmdz_wrf

Timestamp:

Jul 24, 2014, 2:48:30 PM (10 years ago)

Author:

lfita

Message:

Adding checkings

File:

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

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

@@ -251 +251 @@
 !
 
+! Lluis
+      INTEGER                                            :: llp
+      CHARACTER(LEN=50)                                  :: lvarname, lfname
+      REAL                                               :: largest
+
+      llp = 734
+      lfname = 'physiq'
+      largest = 10.e5
+
 !-----------------------------------------------------------------------
 !   initialisation:
@@ -329 +338 @@
 !   --------------------------------------------------------------------
 !
+      lfname='before thermcell_env'
+      lvarname = 'pt'
+      CALL check_var3D(lfname, lvarname, pt, klon, klev, largest, .FALSE.)
+      lvarname = 'pdtadj'
+      CALL check_var3D(lfname, lvarname, pdtadj, klon, klev, 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'
+      lvarname = 'pt'
+      CALL check_var3D(lfname, lvarname, pt, klon, klev, largest, .FALSE.)
+      lvarname = 'pdtadj'
+      CALL check_var3D(lfname, lvarname, pdtadj, klon, klev, largest, .FALSE.)
        
       if (prt_level.ge.1) print*,'thermcell_main apres thermcell_env'
@@ -486 +507 @@
 
       endif
+      lfname='after thermcell_plume'
+      lvarname = 'pt'
+      CALL check_var3D(lfname, lvarname, pt, klon, klev, largest, .FALSE.)
+      lvarname = 'pdtadj'
+      CALL check_var3D(lfname, lvarname, pdtadj, klon, klev, largest, .FALSE.)
 
       if (prt_level.ge.1) print*,'apres thermcell_plume ',lev_out
@@ -641 +667 @@
 !  calcul du transport vertical du moment horizontal
 !------------------------------------------------------------------
+      lfname='before thermcell_dv2'
+      lvarname = 'pt'
+      CALL check_var3D(lfname, lvarname, pt, klon, klev, largest, .FALSE.)
+      lvarname = 'pdtadj'
+      CALL check_var3D(lfname, lvarname, pdtadj, klon, klev, largest, .FALSE.)
 
 !IM 090508  
@@ -669 +700 @@
          enddo
       enddo
+      lfname='after thermcell_dv2'
+      lvarname = 'pt'
+      CALL check_var3D(lfname, lvarname, pt, klon, klev, largest, .FALSE.)
+      lvarname = 'pdtadj'
+      CALL check_var3D(lfname, lvarname, pdtadj, klon, klev, largest, .FALSE.)
 
       if (prt_level.ge.1) print*,'14 OK convect8'
@@ -933 +969 @@
   ENDIF ! iflag_trig_bl
 !  print *,'ENDIF  iflag_trig_bl'    !!jyg
+      lfname='before closure'
+      lvarname = 'pt'
+      CALL check_var3D(lfname, lvarname, pt, klon, klev, largest, .FALSE.)
+      lvarname = 'pdtadj'
+      CALL check_var3D(lfname, lvarname, pdtadj, klon, klev, largest, .FALSE.)
 
 !------------Closure------------------
@@ -1099 +1140 @@
 
       if (prt_level.ge.1) print*,'thermcell_main FIN  OK'
+
+      lfname='before leaving thermcell_main'
+      lvarname = 'pt'
+      CALL check_var3D(lfname, lvarname, pt, klon, klev, largest, .FALSE.)
+      lvarname = 'pdtadj'
+      CALL check_var3D(lfname, lvarname, pdtadj, klon, klev, largest, .FALSE.)
+
 
       return