#!/bin/sh

# some variables
insti=CNRM-CERFACS # institution name
workspa=/bdd/CFMIP/workspace/gcesana/CMIP6/$insti # workspace directory
mod=CNRM-ESM2-1 # model name
grid=gr # type of grid
rr=r1i1p1f2 # initialization type


exp=abrupt-4xCO2  # type of experiment (abrupt-4x-CO2, piControl, amip, amip4k)

# create and access working directory
mkdir $workspa
cd $workspa

# 2D variables
for var in "ts" "psl" "rsut" "rsutcs";
 do
echo $var
# extracting specific time steps (year 120 to 149)
time1=1440  # first index = 0 not 1 
time2=1799
period1=185001-199912
period2=197101-200012
pct=/bdd/CMIP6/CMIP/$insti/$mod/$exp/$rr/Amon/$var/$grid/latest; # location of the files in climserv
# ncks is a nco function that extracts specific variables from netcdf files (here extracting specific time)
ncks -C -d time,$time1,$time2 ${pct}/${var}_Amon_${mod}_${exp}_${rr}_${grid}_${period1}.nc ${var}_Amon_${mod}_${exp}_${rr}_${grid}_${period2}.nc
done

# 3D variables
for var in "ta" "zg";
 do
echo $var
pct=/bdd/CMIP6/CMIP/$insti/$mod/$exp/$rr/Amon/$var/$grid/latest;
# extracting specific time and one particular vertical level (700 mb)
ncks -C -d time,$time1,$time2 -C -d plev,3 ${pct}/${var}_Amon_${mod}_${exp}_${rr}_${grid}_${period1}.nc ${var}700_Amon_${mod}_${exp}_${rr}_${grid}_${period1}.nc
done

for var in "wap";
 do
echo $var
pct=/bdd/CMIP6/CMIP/$insti/$mod/$exp/$rr/Amon/$var/$grid/latest;
# extracting specific time and one particular level (500 mb)
ncks -C -d time,$time1,$time2 -C -d plev,5 $pct/${var}_Amon_${mod}_${exp}_${rr}_${grid}_${period1}.nc ${var}500_Amon_${mod}_${exp}_${rr}_${grid}_${period1}.nc
done


# Repeat for piControl
#piControl
exp=piControl

cd $workspa
for var in "ts" "psl" "rsut" "rsutcs";
 do
echo $var
time1=1440
time2=1799
period1=185001-234912
period2=197101-200012
pct=/bdd/CMIP6/CMIP/$insti/$mod/$exp/$rr/Amon/$var/$grid/latest;
ncks -C -d time,$time1,$time2 $pct/${var}_Amon_${mod}_${exp}_${rr}_${grid}_${period1}.nc ${var}_Amon_${mod}_${exp}_${rr}_${grid}_${period2}.nc
done

period1=185001-234912
for var in "ta" "zg";
 do
echo $var
pct=/bdd/CMIP6/CMIP/$insti/$mod/$exp/$rr/Amon/$var/$grid/latest;
ncks -C -d time,$time1,$time2 -C -d plev,3 $pct/${var}_Amon_${mod}_${exp}_${rr}_${grid}_${period1}.nc ${var}700_Amon_${mod}_${exp}_${rr}_${grid}_${period1}.nc
done

eriod1=185001-209912
for var in "wap";
 do
echo $var
pct=/bdd/CMIP6/CMIP/$insti/$mod/$exp/$rr/Amon/$var/$grid/latest;
ncks -C -d time,$time1,$time2 -C -d plev,5 $pct/${var}_Amon_${mod}_${exp}_${rr}_${grid}_${period1}.nc ${var}500_Amon_${mod}_${exp}_${rr}_${grid}_${period1}.nc
done