update rad

ftomei · ftomei · commit e5d311fe1e23 · 2026-02-24T16:59:19.000+01:00
diff --git a/agrolib/interpolation/spatialControl.cpp b/agrolib/interpolation/spatialControl.cpp
@@ -69,7 +69,10 @@ float findThreshold(meteoVariable myVar, Crit3DMeteoSettings* meteoSettings,
         threshold += distWeight + stdDev * (nrStdDev + 1.f);
     }
     else if (myVar == atmTransmissivity)
-        threshold = MAXVALUE(stdDev * nrStdDev, 0.25f);
+    {
+        distWeight = minDistance / 5000.f;
+        threshold = std::max(stdDev * std::max(nrStdDev, distWeight), 0.3f);
+    }
     else
         threshold = stdDev * nrStdDev;
 
diff --git a/agrolib/solarRadiation/solarRadiation.cpp b/agrolib/solarRadiation/solarRadiation.cpp
@@ -735,7 +735,8 @@ bool computeRadiationRsun(Crit3DRadiationSettings* radSettings, float temperatur
         clearSkyTransmissivity = radSettings->getClearSky();
 
         // noon
-        Crit3DTime noonTime = myTime;
+        Crit3DTime noonTime;
+        noonTime.date = myTime.date;
         noonTime.time = 12*3600;
         if (radSettings->gisSettings->isUTC)
         {

Original file line number	Diff line number	Diff line change
`@@ -735,7 +735,8 @@ bool computeRadiationRsun(Crit3DRadiationSettings* radSettings, float temperatur`
`735`	`735`	`clearSkyTransmissivity = radSettings->getClearSky();`
`736`	`736`
`737`	`737`	`// noon`
`738`		`- Crit3DTime noonTime = myTime;`
	`738`	`+ Crit3DTime noonTime;`
	`739`	`+ noonTime.date = myTime.date;`
`739`	`740`	`noonTime.time = 12*3600;`
`740`	`741`	`if (radSettings->gisSettings->isUTC)`
`741`	`742`	`{`