Changeset 120 in lmdz_wrf

Timestamp:

Jul 28, 2014, 11:25:22 AM (10 years ago)

Author:

lfita

Message:

Adding checks for fraca

File:

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

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

@@ -741 +741 @@
       lvarname = 'zw2'
       CALL check_var3D(lfname, lvarname, zw2, ngrid, nlay+1, largest, .FALSE.)
+      lvarname = 'rhobarz'
+      CALL check_var3D(lfname, lvarname, rhobarz, ngrid, nlay, largest, .FALSE.)
+      lvarname = 'fraca'
+      CALL check_var3D(lfname, lvarname, fraca, ngrid, nlay+1, largest, .FALSE.)
 
 !IM 090508  
@@ -979 +983 @@
      fraca0(ig)=fraca(ig,klcl(ig))+(fraca(ig,klcl(ig)+1) &
  &             -fraca(ig,klcl(ig)))*interp(ig)
+     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
+
      w0(ig)=zw2(ig,klcl(ig))+(zw2(ig,klcl(ig)+1) &
  &         -zw2(ig,klcl(ig)))*interp(ig)