Examine the hourly Rs values for a few days

[qjhart]

First I need a set of stations in the solar location. I'll use the ones from 20190524. I can get them by (in gdb17/cimis/qhart)

g.copy vect=et@20190524,et
# copy for later
v.out.ascii input=et output=~/stations.csv separator=','

sudo ln -s gdb17/cimis/qhart gdb17/solar/foo

Then back in the solar location (~gdb17/solar/quinn)

 g.copy vect=et@foo,et

Getting the points

Once you've done that, getting the points is easy.

for g in G Gi; do 
  for m in 20190524 20190525; do 
    echo ${g}_$m; 
    G=$(g.list -m separator=',' type=rast mapset=$m pattern=*-$g );   
    echo "station_id,x,y,z,$G" > ~/${g}_$m.csv; 
    r.what -v separator=','  map=$G points=et >> ~/${g}_$m.csv; 
 done;
done

[qjhart]

These data were calculated and are available via an Instant Radiation Example