1 | #Â ThL --- 2018 |
---|
2 | |
---|
3 | # Script permettant de diminuer d'un pourcentage voulu la TSI |
---|
4 | # a partir d'un fichier solarforcing_v3.1_piControl_ave_6bands_366days.nc |
---|
5 | # (recupere du compte commun, repertoire /v3.1corr/ ) |
---|
6 | |
---|
7 | # Repertoire de travail, A PERSONNALISER |
---|
8 | ici="/ccc/work/cont003/gencmip6/lurtont/MesInputs/" |
---|
9 | |
---|
10 | # Le fichier d'origine est suppose se trouver dans le rep de travail |
---|
11 | filein=${ici}"solarforcing_v3.1_piControl_ave_6bands_366days.nc" |
---|
12 | |
---|
13 | # Pourcentage de diminution voulu |
---|
14 | pctdim="4.1" |
---|
15 | # Apparaitra en suffixe du nom du fichier avec TSI reduite |
---|
16 | |
---|
17 | fileout=${ici}"solarforcing_v3.1_piControl_ave_6bands_366days_TSI-"${pctdim}"pct.nc" |
---|
18 | |
---|
19 | tsidelta=${ici}tsidelta.nc |
---|
20 | tsireduite=${ici}tsireduite.nc |
---|
21 | |
---|
22 | ssidelta=${ici}ssidelta.nc |
---|
23 | ssireduite=${ici}ssireduite.nc |
---|
24 | |
---|
25 | cdo mulc,-1. -mulc,${pctdim} -divc,100. -selvar,tsi ${filein} ${tsidelta} |
---|
26 | cdo mulc,-1. -mulc,${pctdim} -divc,100. -selvar,ssi ${filein} ${ssidelta} |
---|
27 | cdo add -selvar,"tsi" ${filein} -selvar,"tsi" ${tsidelta} ${tsireduite} |
---|
28 | cdo add -selvar,"ssi" ${filein} -selvar,"ssi" ${ssidelta} ${ssireduite} |
---|
29 | |
---|
30 | cp ${filein} ${fileout} |
---|
31 | ncks -A -v tsi ${tsireduite} ${fileout} |
---|
32 | ncks -A -v ssi ${ssireduite} ${fileout} |
---|
33 | |
---|
34 | # Menage |
---|
35 | rm -f tsi*.nc ssi*.nc |
---|