Changeset 117 in lmdz_wrf

Timestamp:

Jul 28, 2014, 11:00:13 AM (10 years ago)

Author:

lfita

Message:

Adding specific check for w_conv

File:

• branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90 (modified) (1 diff)

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

@@ -983 +983 @@
      w_conv(ig)=w_ls(ig,klcl(ig))+(w_ls(ig,klcl(ig)+1) &
  &             -w_ls(ig,klcl(ig)))*interp(ig)
+     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
+       PRINT *,'    klcl(ig): ', klcl(ig),' klcl(ig)+1: ',klcl(ig)+1,                &
+         ,' w_ls(ig,klcl(ig)): ',w_ls(ig,klcl(ig)),' w_ls(ig,klcl(ig)+1): ',         &
+         w_ls(ig,klcl(ig)+1), ' interp(ig): ',interp(ig)
+     END IF
      therm_tke_max0(ig)=therm_tke_max(ig,klcl(ig)) &
  &                     +(therm_tke_max(ig,klcl(ig)+1)-therm_tke_max(ig,klcl(ig)))*interp(ig)