knutfrode / openwind Goto Github PK

Current usage is:

>>> s = Nansat(SAR_filename)
>>> w = Nansat(wind_filename)
>>> sarwind = openwind.SARwind(s, w)  # SARwind is a function, returning a SAR wind object
>>> sarwind.plot()

New usage will be similar:

>>> sarwind = openwind.SARwind(s, w) # SARwind is here a class

where
s is either a SAR filename or Nansat object
w (optional) is either a wind filename or Nansat object, or a constant direction.

If w is not given, SAR wind will be calculated with a second call:

>>> sarwind = openwind.SARwind(SAR_filename)
(>>> sarwind.reproject(some domain) )
>>> w = Nansat(wind_filename)
>>> sarwind.calculate_wind(w)
>>> sarwind.plot()

The latter two-step process allows the possibility to reproject the SARwind object before wind calculation, taking proper care that look direction is available also after reprojection (currently not possible).
If w is not given here either, download of NCEP GFS wind is attempted.