[GRASS-SVN] r32370 - grass-addons/gipe/i.eb.h_SEBAL95

[svn_grass at osgeo.org]

Author: ychemin
Date: 2008-07-30 07:44:56 -0400 (Wed, 30 Jul 2008)
New Revision: 32370

Modified:
   grass-addons/gipe/i.eb.h_SEBAL95/zom_0.c
Log:
bug fix Z0m

Modified: grass-addons/gipe/i.eb.h_SEBAL95/zom_0.c
===================================================================
--- grass-addons/gipe/i.eb.h_SEBAL95/zom_0.c	2008-07-30 11:39:33 UTC (rev 32369)
+++ grass-addons/gipe/i.eb.h_SEBAL95/zom_0.c	2008-07-30 11:44:56 UTC (rev 32370)
@@ -5,15 +5,15 @@
 double zom_0(double ndvi, double ndvi_max)
 {
 	double a, b, zom;
-	double hv_ndvimax=1.5; /* crop vegetation height (m) */
-	double hv_desert=0.002; /* desert base vegetation height (m) */
+//	double hv_ndvimax=1.5; /* crop vegetation height (m) */
+//	double hv_desert=0.002; /* desert base vegetation height (m) */
 
-	a = (log(hv_desert)-((log(hv_ndvimax/7)-log(hv_desert))/(ndvi_max-0.02)*0.02));
-	b = ((log(hv_ndvimax/7)-log(hv_desert))/(ndvi_max-0.02)) * ndvi;
-	zom = exp(a+b); 
+//	a = (log(hv_desert)-((log(hv_ndvimax/7)-log(hv_desert))/(ndvi_max-0.02)*0.02));
+//	b = ((log(hv_ndvimax/7)-log(hv_desert))/(ndvi_max-0.02)) * ndvi;
+//	zom = exp(a+b); 
 	
 //	printf("****************zom = %5.3f\n", zom);
-	
+	zom = exp(3.3219*ndvi-3.9939);
 	return (zom);
 }