[GRASS-SVN] r68572 - grass/trunk/imagery/i.evapo.mh

[svn_grass at osgeo.org]

Author: ychemin
Date: 2016-06-02 00:20:31 -0700 (Thu, 02 Jun 2016)
New Revision: 68572

Modified:
   grass/trunk/imagery/i.evapo.mh/mh_eto.c
   grass/trunk/imagery/i.evapo.mh/mh_original.c
   grass/trunk/imagery/i.evapo.mh/mh_samani.c
Log:
corrected typo 84600 with actual seconds per day computation

Modified: grass/trunk/imagery/i.evapo.mh/mh_eto.c
===================================================================
--- grass/trunk/imagery/i.evapo.mh/mh_eto.c	2016-06-02 03:14:15 UTC (rev 68571)
+++ grass/trunk/imagery/i.evapo.mh/mh_eto.c	2016-06-02 07:20:31 UTC (rev 68572)
@@ -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/trunk/imagery/i.evapo.mh/mh_original.c
===================================================================
--- grass/trunk/imagery/i.evapo.mh/mh_original.c	2016-06-02 03:14:15 UTC (rev 68571)
+++ grass/trunk/imagery/i.evapo.mh/mh_original.c	2016-06-02 07:20:31 UTC (rev 68572)
@@ -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/trunk/imagery/i.evapo.mh/mh_samani.c
===================================================================
--- grass/trunk/imagery/i.evapo.mh/mh_samani.c	2016-06-02 03:14:15 UTC (rev 68571)
+++ grass/trunk/imagery/i.evapo.mh/mh_samani.c	2016-06-02 07:20:31 UTC (rev 68572)
@@ -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;