* Nikos Alexandris <[email protected]> [2015-05-15 16:34:38 +0300]:
> * Vaclav Petras <[email protected]> [2015-05-15 08:41:17 -0400]: > > > On Fri, May 15, 2015 at 8:10 AM, Nikos Alexandris <[email protected]> > > wrote: > > > > > I have another use-case, hence another update to the script. Does > > > anyone care to have a look at the proposed diff (really a quick hack, no > > > extensive testing done but it works for me). > > > > > > > Looks good except for computing the slope and aspect on the fly. > > r.slope.aspect runs twice once for slope, once for aspect. With current > > code, it actually looks that there might be some override issues. Try > > harder to preserve current mechanism where one r.slope.aspect call is build > > according to missing inputs. You need to change the if conditions somehow. > > (Only reading the diff.) > > Known issue, thanks for reading Vaclav. Will rework another day on it. I am not getting easily away from this script. I want to check whether a strds already exists, and if it does, don't re-create it (overwrite it). Rather, register follow-up maps in the existing strds. Makes sense? This need rose from the attempt to loop r.sun.daily over different monthly Linke Turbitidy values. Something like the following: --%<--- for LINE in `cat monthly_linke_map_start_day_end_day.csv` ;do set -- $(echo ${LINE} |tr "|" " ") && r.sun.daily elevation=srtm4@PERMANENT aspect_value=0 slope_value=0.0 albedo=albedo@albedo linke=$1 start_day=$2 end_day=$3 glob_rad_basename=global_rad_zero_$4 insol_time_basename=insol_time_$4 beam_rad_basename=beam_rad_zero_$4 diff_rad_basename=diff_rad_zero_$4 --o -t ;done --->%-- Currently, each subsequent call of creating a strds, will overwrite if it exists. Hence, only the maps involved in the last run of the above loop will appear in the strds. I would like to add support for multiple linke maps (linke_map=), and then, necessarily, multiple start_day= and end_day= inputs. I don't have the time to try it now, so I just want to add a check if a strds exists. To do so, is something like the `check_stds` [0] exposed to the scripter/appropriate? Thanks, Nikos [0] <http://grass.osgeo.org/grass71/manuals/libpython/_modules/temporal/temporal_algebra.html#TemporalAlgebraParser.check_stds> _______________________________________________ grass-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-dev
