Commit 31031beb authored by Eric Duminil's avatar Eric Duminil
Browse files

renames

parent ed5875ed
......@@ -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);
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment