[GRASS-SVN] r68581 - in grass/branches/releasebranch_7_0/imagery: i.eb.eta i.eb.evapfr i.evapo.mh i.evapo.pt
svn_grass at osgeo.org
svn_grass at osgeo.org
Thu Jun 2 12:26:48 PDT 2016
Author: neteler
Date: 2016-06-02 12:26:48 -0700 (Thu, 02 Jun 2016)
New Revision: 68581
Modified:
grass/branches/releasebranch_7_0/imagery/i.eb.eta/eta.c
grass/branches/releasebranch_7_0/imagery/i.eb.evapfr/i.eb.evapfr.html
grass/branches/releasebranch_7_0/imagery/i.eb.evapfr/main.c
grass/branches/releasebranch_7_0/imagery/i.evapo.mh/mh_eto.c
grass/branches/releasebranch_7_0/imagery/i.evapo.mh/mh_original.c
grass/branches/releasebranch_7_0/imagery/i.evapo.mh/mh_samani.c
grass/branches/releasebranch_7_0/imagery/i.evapo.pt/pt_daily_et.c
Log:
i.eb.* + i.evapo.*: fix for MegaJoules per meter square per day (backport of r68564, r68572, r68573, r68574, r68579, trac #3052); all sync'ed to trunk
Modified: grass/branches/releasebranch_7_0/imagery/i.eb.eta/eta.c
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.eb.eta/eta.c 2016-06-02 19:20:37 UTC (rev 68580)
+++ grass/branches/releasebranch_7_0/imagery/i.eb.eta/eta.c 2016-06-02 19:26:48 UTC (rev 68581)
@@ -6,7 +6,7 @@
double latent, t_celsius, result;
t_celsius = tempk - 273.15;
- latent = 86400 / ((2.501 - 0.002361 * t_celsius) * pow(10, 6));
+ latent = (24.0*60.0*60.0) / ((2.501 - 0.002361 * t_celsius) * pow(10, 6));
result = r_net_day * evap_fr * latent;
Modified: grass/branches/releasebranch_7_0/imagery/i.eb.evapfr/i.eb.evapfr.html
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.eb.evapfr/i.eb.evapfr.html 2016-06-02 19:20:37 UTC (rev 68580)
+++ grass/branches/releasebranch_7_0/imagery/i.eb.evapfr/i.eb.evapfr.html 2016-06-02 19:26:48 UTC (rev 68581)
@@ -1,37 +1,48 @@
<h2>DESCRIPTION</h2>
-<em>i.eb.evapfr</em> calculates the evaporative fraction after [1]. Main
-implementation in [3].
-It takes input of Net Radiation (see <i>r.sun</i>, <i>i.eb.netrad</i>), soil
-heat flux (see i.eb.soilheatflux) and sensible heat flux (see i.eb.hsebal01).
-A flag adds a root zone empirical soil moisture output from the article of
-Makin, Molden and Bastiaanssen (2001).
+<em>i.eb.evapfr</em> calculates the evaporative fraction after
+Bastiaanssen 1995. The main implementation follows Alexandridis et al.
+(2009). The module takes as input the net radiation (see <em>r.sun</em>,
+<em>i.eb.netrad</em>), soil heat flux (see <em>i.eb.soilheatflux</em>)
+and sensible heat flux (see <em>i.eb.hsebal01</em>). A flag adds a
+root zone empirical soil moisture output from the article of
+Bastiaanssen, et al. (2000).
+<!-- TODO
+<h2>EXAMPLE</h2>
+-->
+
<h2>SEE ALSO</h2>
<em>
-<a href="r.sun.html">r.sun</a>,
+<a href="i.eb.hsebal01.html">i.eb.hsebal01</a>,
+<a href="i.eb.netrad.html">i.eb.netrad</a>,
<a href="i.eb.soilheatflux.html">i.eb.soilheatflux</a>,
-<a href="i.eb.hsebal01.html">i.eb.hsebal01</a>
+<a href="r.sun.html">r.sun</a>
</em>
<h2>REFERENCES</h2>
-<p>[1] Bastiaanssen, W.G.M., 1995.
+<p>Bastiaanssen, W.G.M., 1995.
Estimation of Land surface parameters by remote sensing under clear-sky
conditions. PhD thesis, Wageningen University, Wageningen, The Netherlands.
(<a href="http://edepot.wur.nl/206553">PDF</a>)
-<p>[2] Chemin Y., Alexandridis T.A., 2001. Improving spatial resolution of ET
+<p>Bastiaanssen, W.G.M., Molden, D.J., Makin, I.W., 2000.
+Remote sensing for irrigated agriculture: examples from research and
+possible applications. Agricultural water management 46.2: 137-155.
+
+<p>Chemin Y., Alexandridis T.A., 2001. Improving spatial resolution of ET
seasonal for irrigated rice in Zhanghe, China. Asian Journal of Geoinformatics.
-5(1):3-11,2004.
+5(1):3-11.
-<p>[3] Alexandridis T.K., Cherif I., Chemin Y., Silleos N.G., Stavrinos E.,
-Zalidis G.C. Integrated methodology for estimating water use in Mediterranean
-agricultural areas. Remote Sensing. 2009, 1, 445-465.
+<p>Alexandridis T.K., Cherif I., Chemin Y., Silleos N.G., Stavrinos E.,
+Zalidis G.C., 2009.
+Integrated methodology for estimating water use in Mediterranean
+agricultural areas. Remote Sensing. 1, 445-465.
(<a href="http://www.mdpi.com/2072-4292/1/3/445">PDF</a>)
-<p>[4] Chemin, Y., 2012.
+<p>Chemin, Y., 2012.
A Distributed Benchmarking Framework for Actual ET Models,
in: Irmak, A. (Ed.), Evapotranspiration - Remote Sensing and Modeling. InTech.
(<a href="http://www.intechopen.com/books/evapotranspiration-remote-sensing-and-modeling/a-distributed-benchmarking-framework-for-actual-et-models">PDF</a>)
Modified: grass/branches/releasebranch_7_0/imagery/i.eb.evapfr/main.c
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.eb.evapfr/main.c 2016-06-02 19:20:37 UTC (rev 68580)
+++ grass/branches/releasebranch_7_0/imagery/i.eb.evapfr/main.c 2016-06-02 19:26:48 UTC (rev 68581)
@@ -50,8 +50,7 @@
G_add_keyword(_("evaporative fraction"));
G_add_keyword(_("SEBAL"));
module->description =
- _("Computes evaporative fraction (Bastiaanssen, 1995) and "
- "root zone soil moisture (Makin, Molden and Bastiaanssen, 2001).");
+ _("Computes evaporative fraction and root zone soil moisture.");
/* Define the different options */
input1 = G_define_standard_option(G_OPT_R_INPUT);
Modified: grass/branches/releasebranch_7_0/imagery/i.evapo.mh/mh_eto.c
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.evapo.mh/mh_eto.c 2016-06-02 19:20:37 UTC (rev 68580)
+++ grass/branches/releasebranch_7_0/imagery/i.evapo.mh/mh_eto.c 2016-06-02 19:26:48 UTC (rev 68581)
@@ -12,7 +12,7 @@
if (tavg > 100.0) {
tavg = tavg - 273.15; /*in case temperature is in Kelvin */
}
- ra = ra * (84600.0 * 1000.0); /*convert W -> MJ/d */
+ ra = ra * (24.0 * 60.0 * 60.0 / 1000.0); /*convert W -> MJ/d */
result =
0.0013 * 0.408 * ra * (tavg + 17.0) * pow((td - 0.0123 * p), 0.76);
return result;
Modified: grass/branches/releasebranch_7_0/imagery/i.evapo.mh/mh_original.c
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.evapo.mh/mh_original.c 2016-06-02 19:20:37 UTC (rev 68580)
+++ grass/branches/releasebranch_7_0/imagery/i.evapo.mh/mh_original.c 2016-06-02 19:26:48 UTC (rev 68581)
@@ -11,7 +11,7 @@
if (tavg > 100.0) {
tavg = tavg - 273.15; /*in case Temperature is in Kelvin */
}
- ra = ra * (84600.0 * 1000.0); /*convert W -> MJ/d */
+ ra = ra * (24.0 * 60.0 * 60.0 / 1000.0); /*convert W -> MJ/d */
result = 0.0023 * 0.408 * ra * (tavg + 17.8) * pow(td, 0.5);
return result;
}
Modified: grass/branches/releasebranch_7_0/imagery/i.evapo.mh/mh_samani.c
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.evapo.mh/mh_samani.c 2016-06-02 19:20:37 UTC (rev 68580)
+++ grass/branches/releasebranch_7_0/imagery/i.evapo.mh/mh_samani.c 2016-06-02 19:26:48 UTC (rev 68581)
@@ -11,7 +11,7 @@
if (tavg > 100.0) {
tavg = tavg - 273.15; /*in case Temperature is in Kelvin */
}
- ra = ra * (84600.0 * 1000.0); /* convert W -> MJ/d */
+ ra = ra * (24.0 * 60.0 * 60.0 / 1000.0); /* convert W -> MJ/d */
result =
0.0023 * 0.408 * ra * pow(td,
0.5) * ((tmax + tmin) / 2 + 17.8) / 2.45;
Modified: grass/branches/releasebranch_7_0/imagery/i.evapo.pt/pt_daily_et.c
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.evapo.pt/pt_daily_et.c 2016-06-02 19:20:37 UTC (rev 68580)
+++ grass/branches/releasebranch_7_0/imagery/i.evapo.pt/pt_daily_et.c 2016-06-02 19:26:48 UTC (rev 68581)
@@ -12,7 +12,7 @@
/*Latent Heat of vaporization (W/m2/d) */
t_celsius = tempka - 273.15;
- latentHv = 86400 / ((2.501 - 0.002361 * t_celsius) * pow(10, 6));
+ latentHv = (24.0*60.0*60.0) / ((2.501 - 0.002361 * t_celsius) * pow(10, 6));
/* Ratio of slope of saturation-vapour pressure Vs Temperature */
/* ghamma_pt = psychrometric constant */
More information about the grass-commit
mailing list