#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()
final = fi.add_attributes(result, glob_attr, override=True)
path = f.get_file_path(f.Ofolder, f.get_basename(img)+'_completed.nc')
fi.save_ds(path, final)
#f.remove_folder_contents(f.Tfolder)
a += 1
