[GRASS-SVN] r32287 - in grass-addons/gipe: i.eb.netrad i.eb.z0m0

[svn_grass at osgeo.org]

Author: ychemin
Date: 2008-07-25 06:33:32 -0400 (Fri, 25 Jul 2008)
New Revision: 32287

Modified:
   grass-addons/gipe/i.eb.netrad/r_net.c
   grass-addons/gipe/i.eb.z0m0/zom_0.c
Log:
bug fixing

Modified: grass-addons/gipe/i.eb.netrad/r_net.c
===================================================================
--- grass-addons/gipe/i.eb.netrad/r_net.c	2008-07-25 10:29:46 UTC (rev 32286)
+++ grass-addons/gipe/i.eb.netrad/r_net.c	2008-07-25 10:33:32 UTC (rev 32287)
@@ -20,9 +20,9 @@
 // 	printf("**rnet: bbalb = %5.3f\n\tndvi = %5.3f\n\ttempk = %5.3f\n\te0 = %5.3f\n\ttsw = %5.3f\n\tdoy = %f\n\tutc = %5.3f\n\tsunzangle = %5.3f\n\tdtair = %5.3f\n",bbalb,ndvi,tempk,e0,tsw,doy,utc,sunzangle,dtair);
 	
 	// Atmospheric emissivity (Bastiaanssen, 1995)
-	e_atm	=  1.08 * pow(-log10(tsw_for_e_atm),0.265) ;
+	e_atm	=  1.08 * pow(-log(tsw_for_e_atm),0.265) ;
 	// Atmospheric emissivity (Pawan, 2004)
-// 	e_atm	= 0.85 * pow(-log10(tsw),0.09);
+// 	e_atm	= 0.85 * pow(-log(tsw),0.09);
 // 	printf("rnet: e_atm = %5.3f\n",e_atm);
 
 	ds = 1.0 + 0.01672 * sin(2*PI*(doy-93.5)/365);

Modified: grass-addons/gipe/i.eb.z0m0/zom_0.c
===================================================================
--- grass-addons/gipe/i.eb.z0m0/zom_0.c	2008-07-25 10:29:46 UTC (rev 32286)
+++ grass-addons/gipe/i.eb.z0m0/zom_0.c	2008-07-25 10:33:32 UTC (rev 32287)
@@ -7,8 +7,8 @@
 	//double hv_ndvimax=1.5; /* crop vegetation height (m) */
 	//double hv_desert=0.002; /* desert base vegetation height (m) */
 	
-	a = (log10(hv_desert)-((log10(hv_ndvimax/7)-log10(hv_desert))/(ndvi_max-0.02)*0.02));
-	b = (log10(hv_ndvimax/7)-log10(hv_desert))/(ndvi_max-0.02)* ndvi;
+	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;
 	z0m = exp(a+b); 
 	
 	return (z0m);