[GRASS-SVN] r73288 - grass/trunk/imagery/i.eb.hsebal01

[svn_grass at osgeo.org]

Author: ychemin
Date: 2018-09-07 14:36:50 -0700 (Fri, 07 Sep 2018)
New Revision: 73288

Modified:
   grass/trunk/imagery/i.eb.hsebal01/main.c
Log:
updated code for eact and z0m nullval in auto mode

Modified: grass/trunk/imagery/i.eb.hsebal01/main.c
===================================================================
--- grass/trunk/imagery/i.eb.hsebal01/main.c	2018-09-07 19:03:54 UTC (rev 73287)
+++ grass/trunk/imagery/i.eb.hsebal01/main.c	2018-09-07 21:36:50 UTC (rev 73288)
@@ -300,17 +300,25 @@
 	/*********************/
 	for (row = 0; row < nrows; row++) {
 	    DCELL d_t0dem;
+	    DCELL d_z0m;
+	    DCELL d_eact;
 
 	    G_percent(row, nrows, 2);
 	    Rast_get_d_row(infd_t0dem, inrast_t0dem, row);
+	    Rast_get_d_row(infd_z0m, inrast_z0m, row);
+	    Rast_get_d_row(infd_eact, inrast_eact, row);
 	    Rast_get_d_row(infd_Rn, inrast_Rn, row);
 	    Rast_get_d_row(infd_g0, inrast_g0, row);
 	    /*process the data */
 	    for (col = 0; col < ncols; col++) {
 		d_t0dem = ((DCELL *) inrast_t0dem)[col];
+		d_z0m = ((DCELL *) inrast_z0m)[col];
+		d_eact = ((DCELL *) inrast_eact)[col];
 		d_Rn = ((DCELL *) inrast_Rn)[col];
 		d_g0 = ((DCELL *) inrast_g0)[col];
 		if (Rast_is_d_null_value(&d_t0dem) ||
+		    Rast_is_d_null_value(&d_eact) ||
+		    Rast_is_d_null_value(&d_z0m) ||
 		    Rast_is_d_null_value(&d_Rn) ||
 		    Rast_is_d_null_value(&d_g0)) {
 		    /* do nothing */
@@ -431,6 +439,11 @@
 		/* do nothing */
 		d_Roh[row][col] = -999.9;
 		d_Rah[row][col] = -999.9;
+		if (row == rowDry && col == colDry) {	/*collect dry pix info */
+		    d_rah_dry = d_rah1;
+		    d_roh_dry = d_roh1;
+		    G_message("d_rah_dry=%f d_roh_dry=%f",d_rah_dry,d_roh_dry);
+		}
 	    }
 	    else {
 		d_u5 = (ustar / 0.41) * log(5 / d_z0m);