diff --git a/tree-water-demand/src/main/java/de/hft/stuttgart/water/TreeWaterDemand.java b/tree-water-demand/src/main/java/de/hft/stuttgart/water/TreeWaterDemand.java index 9798d78316f0dfded3b60bb68094c24b9dc68ac7..5d5dae70cfd784b019d8ae28f3cb82ad8e25a402 100644 --- a/tree-water-demand/src/main/java/de/hft/stuttgart/water/TreeWaterDemand.java +++ b/tree-water-demand/src/main/java/de/hft/stuttgart/water/TreeWaterDemand.java @@ -101,8 +101,8 @@ public class TreeWaterDemand { int dayOfYear = counter.getDayOfYear(); double hourOfDay = counter.getHour() + 0.5; - double extraterrestrialRadiation = extraterrestrialRadiation(dayOfYear, hourOfDay, lz, lm); // ra - double netLongwaveRadiation = rnl(temp, actualVaporPressure, incomingRadiation, height, extraterrestrialRadiation); // rnl + double extraterrestrialRadiation = ra(dayOfYear, hourOfDay, lz, lm); + double netLongwaveRadiation = rnl(temp, actualVaporPressure, incomingRadiation, height, extraterrestrialRadiation); double netRadiation = rn(incomingRadiation, netLongwaveRadiation); double et0 = et0(temp, actualVaporPressure, incomingRadiation, height, windSpeed, humidity, netRadiation); double et0u = et0u(svf, advection, et0); @@ -190,12 +190,14 @@ public class TreeWaterDemand { } } + // Net radiation // unit is apparently [MJ / (m².h)] private static double rn(double shortWaveRadiation, double rnl) { // 1 - 0.23: albedo effect of grass return (1 - 0.23) * shortWaveRadiation - rnl; } + // NetLongwaveRadiation // unit is apparently [MJ / (m².h)] private static double rnl(double temperature, double actualVaporPressure, double shortWaveRadiation, double height, double extraterrestrialRadiation){ @@ -214,7 +216,8 @@ public class TreeWaterDemand { * (1.35 * rsRsoQuotient - 0.35); } - private static double extraterrestrialRadiation(int dayOfYear, double hourOfDay, double lz, double lm) { + // ExtraterrestrialRadiation + private static double ra(int dayOfYear, double hourOfDay, double lz, double lm) { double dr = 1 + 0.033 * Math.cos(2 * Math.PI * dayOfYear / 365d); double b = (2 * Math.PI * (dayOfYear - 81)) / 364; double sc = 0.1645 * Math.sin(2 * b) - 0.1255 * Math.cos(b) - 0.025 * Math.sin(b);