[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