Commit e5919908 authored by Eric Duminil's avatar Eric Duminil
Browse files

Bugfix: Eq. 37

parent 636e41dc
This source diff could not be displayed because it is too large. You can view the blob instead.
...@@ -206,7 +206,7 @@ public class TreeWaterDemand { ...@@ -206,7 +206,7 @@ public class TreeWaterDemand {
double extraterrestrialRadiation){ double extraterrestrialRadiation){
// divide by 24 for hourly // divide by 24 for hourly
double sigma = 4.903E-9 / 24; double sigma = 4.903E-9 / 24;
double rso = (0.75 + 2 * 10E-5 * altitude) * extraterrestrialRadiation; double rso = (0.75 + 2E-5 * altitude) * extraterrestrialRadiation;
double rsRsoQuotient = 0.5; double rsRsoQuotient = 0.5;
if (rso != 0) { if (rso != 0) {
rsRsoQuotient = rs / rso; rsRsoQuotient = rs / rso;
...@@ -219,6 +219,7 @@ public class TreeWaterDemand { ...@@ -219,6 +219,7 @@ public class TreeWaterDemand {
} }
// ExtraterrestrialRadiation // ExtraterrestrialRadiation
// unit is apparently [MJ / (m².h)]
private static double ra(int dayOfYear, double hourOfDay, double lz, double lm, double latitude) { private static double ra(int dayOfYear, double hourOfDay, double lz, double lm, double latitude) {
double dr = 1 + 0.033 * Math.cos(2 * Math.PI * dayOfYear / 365d); double dr = 1 + 0.033 * Math.cos(2 * Math.PI * dayOfYear / 365d);
double b = (2 * Math.PI * (dayOfYear - 81)) / 364; double b = (2 * Math.PI * (dayOfYear - 81)) / 364;
......
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