# ThL --- 2018

# Script permettant de diminuer d'un pourcentage voulu la TSI
# a partir d'un fichier solarforcing_v3.1_piControl_ave_6bands_366days.nc
# (recupere du compte commun, repertoire /v3.1corr/ )

# Repertoire de travail, A PERSONNALISER
ici="/ccc/work/cont003/gencmip6/lurtont/MesInputs/"

# Le fichier d'origine est suppose se trouver dans le rep de travail
filein=${ici}"solarforcing_v3.1_piControl_ave_6bands_366days.nc"

# Pourcentage de diminution voulu
pctdim="4.1"
# Apparaitra en suffixe du nom du fichier avec TSI reduite

fileout=${ici}"solarforcing_v3.1_piControl_ave_6bands_366days_TSI-"${pctdim}"pct.nc"

tsidelta=${ici}tsidelta.nc
tsireduite=${ici}tsireduite.nc

ssidelta=${ici}ssidelta.nc
ssireduite=${ici}ssireduite.nc

cdo mulc,-1. -mulc,${pctdim} -divc,100. -selvar,tsi ${filein} ${tsidelta}
cdo mulc,-1. -mulc,${pctdim} -divc,100. -selvar,ssi ${filein} ${ssidelta}
cdo add -selvar,"tsi" ${filein} -selvar,"tsi" ${tsidelta} ${tsireduite}
cdo add -selvar,"ssi" ${filein} -selvar,"ssi" ${ssidelta} ${ssireduite}

cp ${filein} ${fileout}
ncks -A -v tsi ${tsireduite} ${fileout}
ncks -A -v ssi ${ssireduite} ${fileout}

# Menage
rm -f tsi*.nc ssi*.nc
