#!/bin/bash

##NOTE:
#This is a prelimnary happy path, that still encounters some problems.

#1. load data:
#cate ds copy esacci.SOILMOISTURE.day.L3S.SSMV.multi-sensor.multi-platform.COMBINED.03-2.r1 --name SOILM_2006_2009_region_7_48_10_52 --time '2006-01-01,2009-12-31' --region '7,48,10,52'
cate ds copy esacci.CLOUD.mon.L3C.CLD_PRODUCTS.multi-sensor.multi-platform.ATSR2-AATSR.2-0.r1 --name CLOUD_2006_2009_region_7_48_10_52 --time '2006-01-01,2009-12-31' --region '7,48,10,52'

#2. open workspace:
cate ws new

#3. open resource:
cate res open cloud local.CLOUD_2006_2009_region_7_48_10_52 
#cate res open soilm local.SOILM_2006_2009_region_7_48_10_52

#3.b. lta:
cate res set cloud_lta long_term_average ds=@cloud var="cfc"
#cate res set soilm_lta long_term_average ds=@soilm var="???"

#4. Select timeseries:
cate res set cloud_point tseries_point ds=@cloud_lta point=8.78,50.1 var="cfc"
#cate res set soilm_point tseries_point ds=@soilm_lta point=8.78,50.1 var="???"

#5.. propagation of uncertainties?


#6. load stationdata (pointdata) #'read_text' or 'read_csv'

#cate res set res_1 read_text file="/data/sherzog/cate-0.9.0.dev4/UC02/produkt_klima_Tageswerte_20050715_20151231_07341_output.txt"
cate res set station read_csv file="/data/sherzog/cate-0.9.0.dev4/UC02/produkt_klima_Tageswerte_20050715_20151231_07341_output2.txt" delimiter=";" index_col="time"


#6.b. subset_temporal?

#7. filter the stationdata (detection of outliers)
cate res set station_ds from_dataframe df=@station
cate res set station_outl detect_outliers ds=@station_ds var=" precipitation" threshold_low=0.00 threshold_high=100.00 quantiles=False
#cate res set res_7 detect_outliers ds=@res_3 var=" precipitation"

#7.b. lta on station-data:
cate res set station_lta long_term_average ds=@station_outl var="precipitation"
 

#8. Plot the timeseries of precipitation 
cate res set res_8 plot ds=@station_ds var="precipitation" file="/data/sherzog/cate-0.9.0.dev4/UC02/station_ds_precipitation.jpg"
cate res set res_9 plot ds=@station_lta var="precipitation" file="/data/sherzog/cate-0.9.0.dev4/UC02/station_lta_precipitation.jpg"

#9. Plot the timeseries of the ODP datasets
cate res set cld_point_plot plot ds=@cloud_point var="cfc" file="/data/sherzog/cate-0.9.0.dev4/UC02/cld_point.jpg"

cate ws save
cate ws close
cate ws exit -y