Changeset 121 in lmdz_wrf

Timestamp:

Jul 28, 2014, 11:45:08 AM (10 years ago)

Author:

lfita

Message:

Adding checks of 'fraca' at the right place....

File:

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

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

@@ -936 +936 @@
     enddo
 !!jyg fin
+      lfname='thermcell_main before LCL'
+      lvarname = 'rhobarz0'
+      CALL check_var(lfname, lvarname, rhobarz0, ngrid, largest, .FALSE.)
+      lvarname = 'p_con'
+      CALL check_var(lfname, lvarname, p_con, ngrid, largest, .FALSE.)
+      lvarname = 'fraca'
+      CALL check_var3D(lfname, lvarname, fraca, ngrid, nlay+1, largest, .FALSE.)
+      lvarname = 'fraca0'
+      CALL check_var(lfname, lvarname, fraca0, ngrid, largest, .FALSE.)
+      lvarname = 'w_conv'
+      CALL check_var(lfname, lvarname, w_conv, ngrid, largest, .FALSE.)
+      lvarname = 'interp'
+      CALL check_var(lfname, lvarname, interp, ngrid, largest, .FALSE.)
+      lvarname = 'zlcl'
+      CALL check_var(lfname, lvarname, zlcl, ngrid, largest, .FALSE.)
 
 !------------Calcul des propriétés du thermique au LCL 
@@ -1010 +1025 @@
      if (pbl_tke_max0(ig).ge.20.) pbl_tke_max0(ig)=20.
    else 
+     IF (fraca0(ig) /= fraca0(ig) .OR. ABS(fraca0(ig)) > largest*10.e5) THEN
+       PRINT *,'  Lluis wrong fraca0(ig): ',fraca0(ig),' at : ',ig
+       PRINT *,'    klcl(ig): ', klcl(ig),' klcl(ig)+1: ',klcl(ig)+1,                &
+         ' fraca(ig,klcl(ig)): ',fraca(ig,klcl(ig)),' fraca(ig,klcl(ig)+1): ',       &
+         fraca(ig,klcl(ig)+1), ' interp(ig): ',interp(ig)
+     END IF
      IF (w_conv(ig) /= w_conv(ig) .OR. ABS(w_conv(ig)) > largest*10.e5) THEN
        PRINT *,'  Lluis wrong w_conv(ig): ',w_conv(ig),' at : ',ig