| 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 |
|---|
