"""
Aina Johannessen. 
2017/2018 - Masterthesis

Using ECMWF analysis data to create thickness from geopotential height and get wind at 250
"""



#get wind at 250 hpa and height 500-1000 hpa:


#---PATHS-----------------------------------------------------------------------------------
path_gribs="../gribs/"
#path_figures="../figs/"
#path_norestore_files="/Users/ainajohannessen/Documents/Aina/skole/master/Masterthesis/data_norstore/"

#-------------------------------------------------------------------------------------------
dates=(20171127 20171128 20171129 20171130 20171201 20171202 20171203 20171204 20171205 20171206 20171207 20171208 20171209 20171210) #
times=( 0000 0600 1200 1800 ) # 1200 1800 )

for date in ${dates[@]}; do
	for time in ${times[@]}; do	
		path_name_wind="${path_gribs}param_uv_${date}_${time}.grib"
                #wind at 250 hpa
		wind="cdo -sellevel,25000 ${path_name_wind} ${path_gribs}wind250_${date}_${time}.grib"
		echo $wind
		$wind
		path_name_z="${path_gribs}param_z_${date}_${time}.grib"
                #Thickness=z_500 - z_1000 is thickness of layer between 1000 hpa and 500 hpa
		Thickness="cdo -R -sub -sellevel,50000 ${path_name_z} -sellevel,100000 ${path_name_z} ${path_gribs}Z_1000_500_${date}_${time}.grib"
		echo $Thickness
		$Thickness
	done
done

#path_name_z="${path_norestore_files}SCA/pl/z/param_z_${date}_${time}.grib"
#z_500="cdo sellevidx,9 ${path_name_z} ${path_gribs}z500.grib"
#Thicness="cdo -R -sub -sellevidx,9 ${path_name_z} -sellevidx,1 ${path_name_z} ${path_gribs}Z_1000_500.grib"
#cdo -R -sub -sellevidx,9 data_norstore/SCA/pl/z/param_z_20160920_0000.grib -sellevidx,1 ...
#  ..  data_norstore/SCA/pl/z/param_z_20160920_0000.grib test_Z_1000_500.grib