Joel Fiddes
Joel Fiddes
https://github.com/ArcticSnow/TopoPyScale/discussions/83#discussion-5231360
some more I think related info on this issue https://github.com/ecmwf/cfgrib/issues/110 Basically seems safer to use parallel =False with mf_opendataset otherwise there is a chance of conflict between threads doing "stuff"...
 
confirmed for the Alps (Kandersteg) This is Sw direct 2015-06-01 12:00. Values 0 (blue) - 529 (yellow) Wm-2. North face has higher SWdir than South face
I still dont understand l.149 in topo_param.py: ds['aspect'] = (["y", "x"], np.deg2rad(aspect - 180)) why do we subtract 180 deg? doesnt this reverses the aspect?
with GST results I just made I get cold N and warm S - just rerunning case above to see how SW is affected.....
.... with the 180 removed, that is.
removed the 180 ...looks a bit better :) this is Jan 31 2015 12:00 with similar color scale as above. S now has 400+ Wm2 and N 0-ish much better!
prob not steep north in January - I think its realistic. Better than before :)