v486b - throw warning if gauges don't have required LW Incoming

Author: James Craig

src/Diagnostics.cpp

@@ -1045,7 +1045,7 @@ double CDiagnostic::CalculateDiagnostic(CTimeSeriesABC  *pTSMod,
 
     if (_type==DIAG_KLING_GUPTA_DEVIATION){Beta=1.0;} //remove penalty for difference in means
 
-    if (_type==DIAG_KGE_PRIME){Alpha/=Beta;}// Uses C.O.V. instead of std dev. from Kling et al. (2012) Runoff conditions in the upper Danube basin under an ensemble of climate change scenarios, Journal of Hydrology
+    if (_type==DIAG_KGE_PRIME){if (Beta!=0.0){Alpha/=Beta;}}// Uses C.O.V. instead of std dev. from Kling et al. (2012) Runoff conditions in the upper Danube basin under an ensemble of climate change scenarios, Journal of Hydrology
 
     if ((N>0) && ((ObsAvg!=0.0) || (Beta==1.0)) && (ObsStd!=0.0) && (ModStd!=0.0))
     {

src/Gauge.cpp

@@ -225,6 +225,7 @@ void CGauge::Initialize(const optStruct   &Options,
   WarnAboutForcing(Options.SW_radiation  ==SW_RAD_DATA,   F_SW_RADIA);
   WarnAboutForcing(Options.SW_radia_net  ==NETSWRAD_DATA, F_SW_RADIA_NET);
   WarnAboutForcing(Options.LW_radiation  ==LW_RAD_DATA,   F_LW_RADIA_NET);
+  WarnAboutForcing(Options.LW_incoming   ==LW_INC_DATA,   F_LW_INCOMING);
   WarnAboutForcing(Options.cloud_cover   ==CLOUDCOV_DATA, F_CLOUD_COVER);
   WarnAboutForcing(Options.ow_evaporation==PET_DATA,      F_OW_PET);
   WarnAboutForcing(Options.evaporation   ==PET_DATA,      F_PET);