Changeset 74 in lmdz_wrf

Timestamp:

Jul 24, 2014, 3:05:59 PM (10 years ago)

Author:

lfita

Message:

Getting deeper in the checkings...

Location:

branches/LMDZ_WRFmeas/WRFV3/lmdz

Files:

• thermcell_dq.F90 (modified) (4 diffs)
• thermcell_main.F90 (modified) (2 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_dq.F90

@@ -39 +39 @@
 
 
+! Lluis
+      INTEGER                                            :: llp
+      CHARACTER(LEN=50)                                  :: lvarname, lfname
+      REAL                                               :: largest
+
+      llp = 734
+      lfname = 'physiq'
+      largest = 10.e5
+
 ! Old explicite scheme
+      lfname='before old scheme'
+      lvarname = 'dq'
+      CALL check_var3D(lfname, lvarname, dq, klon, klev, largest, .FALSE.)
+
       if (impl==-1) then
          call thermcell_dq_o(ngrid,nlay,ptimestep,fm,entr,  &
@@ -45 +58 @@
          return
       endif
+      lfname='after old scheme'
+      lvarname = 'dq'
+      CALL check_var3D(lfname, lvarname, dq, klon, klev, largest, .FALSE.)
 
 ! Calcul du critere CFL pour l'advection dans la subsidence
@@ -128 +144 @@
       enddo
    endif
+      lfname='before tendencies _dq'
+      lvarname = 'dq'
+      CALL check_var3D(lfname, lvarname, dq, klon, klev, largest, .FALSE.)
 
 ! Tendencies
@@ -136 +155 @@
          enddo
       enddo
+      lfname='after tendencies _dq'
+      lvarname = 'q'
+      CALL check_var3D(lfname, lvarname, q, klon, klev, largest, .FALSE.)
+      lvarname = 'dq'
+      CALL check_var3D(lfname, lvarname, dq, klon, klev, largest, .FALSE.)
+      lvarname = 'qold'
+      CALL check_var3D(lfname, lvarname, qold, klon, klev, largest, .FALSE.)
 
 return

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

@@ -641 +641 @@
 !   calcul du transport vertical
 !------------------------------------------------------------------
+      lfname='before transport_vertical'
+      lvarname = 'zdthladj'
+      CALL check_var3D(lfname, lvarname, zdthladj, klon, klev, largest, .FALSE.)
 
       call thermcell_dq(ngrid,nlay,dqimpl,ptimestep,fm0,entr0,masse,  &
@@ -646 +649 @@
       call thermcell_dq(ngrid,nlay,dqimpl,ptimestep,fm0,entr0,masse,  &
      &                   po,pdoadj,zoa,lev_out)
+      lfname='after transport_vertical'
+      lvarname = 'zdthladj'
+      CALL check_var3D(lfname, lvarname, zdthladj, klon, klev, largest, .FALSE.)
 
 !------------------------------------------------------------------