[GRASS-SVN] r73274 - grass-addons/grass7/imagery/i.rh

[svn_grass at osgeo.org]

Author: ychemin
Date: 2018-09-05 06:14:16 -0700 (Wed, 05 Sep 2018)
New Revision: 73274

Modified:
   grass-addons/grass7/imagery/i.rh/rh.c
Log:
fixed typos

Modified: grass-addons/grass7/imagery/i.rh/rh.c
===================================================================
--- grass-addons/grass7/imagery/i.rh/rh.c	2018-09-05 13:13:34 UTC (rev 73273)
+++ grass-addons/grass7/imagery/i.rh/rh.c	2018-09-05 13:14:16 UTC (rev 73274)
@@ -1,32 +1,14 @@
 #include<stdio.h>
 #include<math.h>
 
-double rh(double PW,double Pa,double Ta,double dem){
-	/*
-	https{//www.researchgate.net/publication/227247013_High-resolution_Surface_Relative_Humidity_Computation_Using_MODIS_Image_in_Peninsular_Malaysia/
-	PW <- MOD05_L2 product
-	Pa <- MOD07 product
-	Pa <- 1013.3-0.1038*dem
-	Ta <- MOD07 product
-	Ta <- -0.0065*dem+TaMOD07 (if dem>400m)
-	
-	https{//ladsweb.modaps.eosdis.nasa.gov/archive/allData/61/MOD05_L2/
-	https{//ladsweb.modaps.eosdis.nasa.gov/archive/allData/61/MOD07_L2/
-	*/
-	/*Specific Humidity*/
-	double q=0.001*(-0.0762*PW*PW+1.753*PW+12.405);
-	double ta=-0.0065*dem+Ta;
-	double a=17.2694*ta/(ta+238.3);
-	return(q*Pa/(380*exp(a)));
-}
-
 double esat(double tamean){
 	/*
 	esat{ saturated vapour pressure
 	tamean{ air temperature daily mean
 	*/
-	return(610.78*exp(17.2694*tamean/(tamean+238.3)));
+	return(0.61078*exp(17.2694*tamean/(tamean+237.3)));
 }
+
 double eact(double esat,double rh){
 	/*
 	eact{ actual vapour pressure
@@ -35,6 +17,7 @@
 	*/
 	return(0.01*esat*rh);
 }
+
 double eatm(double eact){
 	/*
 	eatm: atmospheric emissivity
@@ -42,3 +25,23 @@
 	*/
 	return(1.34-0.14*sqrt(eact));
 }
+
+double rh(double PW,double Pa,double Ta,double dem){
+	/*
+	https{//www.researchgate.net/publication/227247013_High-resolution_Surface_Relative_Humidity_Computation_Using_MODIS_Image_in_Peninsular_Malaysia/
+	PW <- MOD05_L2 product [mm]
+	Pa <- MOD07 product [kPa]
+	Pa <- 1013.3-0.1038*dem [kPa]
+	Ta <- MOD07 product [Celsius]
+	Ta <- -0.0065*dem+TaMOD07 (if dem>400m) [Celsius]
+	
+	https{//ladsweb.modaps.eosdis.nasa.gov/archive/allData/61/MOD05_L2/
+	https{//ladsweb.modaps.eosdis.nasa.gov/archive/allData/61/MOD07_L2/
+	*/
+	/*Specific Humidity*/
+	double q=0.001*(-0.0762*PW*PW+1.753*PW+12.405);
+	double ta=-0.0065*dem+Ta;
+	double a=17.2694*ta/(ta+237.3);
+	return(2.63224*q*Pa/exp(a));
+}
+