Changeset 103 in lmdz_wrf

Timestamp:

Jul 25, 2014, 4:14:28 PM (10 years ago)

Author:

lfita

Message:

Checking only on horizontal (inside big klev loop)

File:

• branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_flux2.F90 (modified) (5 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_flux2.F90

@@ -59 +59 @@
       CHARACTER(LEN=50)                                  :: lvarname, lfname
       REAL                                               :: largest
+      CHARACTER(LEN=4)                                   :: lS
 
       llp = 734
@@ -197 +198 @@
 
       do l=1,klev
-
+         WRITE(lS,'(I4)')l
          do ig=1,ngrid
             if (l.lt.lmax(ig)) then
@@ -209 +210 @@
          enddo
 
-      lfname='thermcell_flux2 after first loop'
+      lfname='thermcell_flux2 after first loop ' // lS
       lvarname = 'fm'
-      CALL check_var3D(lfname, lvarname, fm, ngrid, klev+1, largest, .FALSE.)
+      CALL check_var(lfname, lvarname, fm(:,l), ngrid, largest, .FALSE.)
       lvarname = 'entr'
-      CALL check_var3D(lfname, lvarname, entr, ngrid, klev, largest, .FALSE.)
+      CALL check_var(lfname, lvarname, entr(:,l), ngrid, largest, .FALSE.)
       lvarname = 'detr'
-      CALL check_var3D(lfname, lvarname, detr, ngrid, klev, largest, .FALSE.)
+      CALL check_var(lfname, lvarname, detr(:,l), ngrid, largest, .FALSE.)
 
 !-------------------------------------------------------------------------
@@ -257 +258 @@
       endif
 
-      lfname='thermcell_flux2 after fraca croissant'
+      lfname='thermcell_flux2 after fraca croissant ' //lS
       lvarname = 'fm'
-      CALL check_var3D(lfname, lvarname, fm, ngrid, klev+1, largest, .FALSE.)
+      CALL check_var(lfname, lvarname, fm(:,l), ngrid, klev+1, largest, .FALSE.)
+      lvarname = '1/zw2'
+      CALL check_var(lfname, lvarname, 1./zw2(:,l), ngrid, klev+1, largest, .FALSE.)
+      lvarname = '1/rhobarz'
+      CALL check_var(lfname, lvarname, 1./rhobarz(:,l), ngrid, klev+1, largest, .FALSE.)
       lvarname = 'entr'
-      CALL check_var3D(lfname, lvarname, entr, ngrid, klev, largest, .FALSE.)
+      CALL check_var(lfname, lvarname, entr(:,l), ngrid, klev, largest, .FALSE.)
       lvarname = 'detr'
-      CALL check_var3D(lfname, lvarname, detr, ngrid, klev, largest, .FALSE.)
+      CALL check_var(lfname, lvarname, detr(:,l), ngrid, klev, largest, .FALSE.)
 
       if (prt_level.ge.10) &
@@ -286 +291 @@
       endif
 
-      lfname='thermcell_flux2 after massa croissant'
+      lfname='thermcell_flux2 after massa croissant ' // lS
       lvarname = 'fm'
-      CALL check_var3D(lfname, lvarname, fm, ngrid, klev+1, largest, .FALSE.)
+      CALL check_var(lfname, lvarname, fm(:,l), ngrid, klev+1, largest, .FALSE.)
       lvarname = 'entr'
-      CALL check_var3D(lfname, lvarname, entr, ngrid, klev, largest, .FALSE.)
+      CALL check_var(lfname, lvarname, entr(:,l), ngrid, klev, largest, .FALSE.)
       lvarname = 'detr'
-      CALL check_var3D(lfname, lvarname, detr, ngrid, klev, largest, .FALSE.)
+      CALL check_var(lfname, lvarname, detr(:,l), ngrid, klev, largest, .FALSE.)
 
       if (prt_level.ge.10) &