#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
finally:
print ('done')
if tempfolder_full_path == None:
f.remove_dir(f.Tfolder)
