wi = WIND(dataset = wind_ds, output = wind_ds, drop_var_output = True, inter_layer = wind_extra_layers)
wi.get_direction_matrix()
wi.get_speed_matrix()
result = wi.return_ds(wi.out)
else:
#just prepare canvas
gen_amp = WIND(dataset=amp, output=None)
canvas = gen_amp.drop_var(gen_amp.ds, list(gen_amp.ds.data_vars))
result = gen_amp.return_ds(canvas)
if wav and sla:
#merging sigma and slant 10x10 products
print ('Calculating wave layers')
gen_wav = Wave(dataset=wav, output=None)
gen_sla = Wave(dataset=sla, output=None)
wave_ds = gen_wav.merge_datasets([gen_wav.ds, gen_sla.ds])
#calculating wave parameters
wa = Wave(dataset = wave_ds, output = result, inter_layer = wave_extra_layers)
wa.get_wave_arrays()
result = wa.return_ds(wa.out)
if obj:
#calculating objects
print ('Calculating object layers')
ob = Object(dataset = obj, output = result, inter_layer = object_extra_layers)
ob.get_objects_array()
result = ob.return_ds(ob.out)
#saving netcdf with results
fi = Object()
