[4006] | 1 | # Script prenant en entree les outputs de la simu AER-ssp245, et generant les depots |
---|
| 2 | # sur la periode 2015-2020, a etre utilises ensuite par les simus DAMIP. |
---|
| 3 | # Une moyenne glissante sur 3 ans est effectuée, avec pondération [0.25 ; 0.5 ; 0.25]. |
---|
| 4 | # Un fichier depots par an est genere au final. |
---|
| 5 | # Th. Lurton IPSL aout 2018 |
---|
| 6 | |
---|
| 7 | #!/bin/sh |
---|
| 8 | |
---|
| 9 | for year in {2015..2020} |
---|
| 10 | do |
---|
| 11 | |
---|
| 12 | # Chemins |
---|
| 13 | chemin=/ccc/store/cont003/gencmip6/lurtont/IGCM_OUT/LMDZORINCA/PROD/scenarioAER/AER-ssp245/CHM/Output/MO/ |
---|
| 14 | workchemin=/ccc/work/cont003/gencmip6/lurtont/aerosols_from_inca/dep_for_DAMIP/ |
---|
| 15 | |
---|
| 16 | if [ ! -d ${workchemin} ] ; then mkdir -p ${workchemin} ; fi |
---|
| 17 | |
---|
| 18 | # On effectue une moyenne glissante sur trois années, centrée en l'année d'intérêt, |
---|
| 19 | # et de poids [ 1/4 ; 1/2 ; 1/4 ]. |
---|
| 20 | prec=$((${year}-1)) |
---|
| 21 | suiv=$((${year}+1)) |
---|
| 22 | |
---|
| 23 | echo '***** Annee '${year}' : Moyenne glissante sur ans '${prec}', '${year}', et '${suiv}'.' |
---|
| 24 | |
---|
| 25 | # Adresses generiques des fichiers a traiter |
---|
| 26 | fileprec=${chemin}AER-ssp245_${prec}0101_${prec}1230_1M_inca_dep.nc |
---|
| 27 | fileyear=${chemin}AER-ssp245_${year}0101_${year}1230_1M_inca_dep.nc |
---|
| 28 | filesuiv=${chemin}AER-ssp245_${suiv}0101_${suiv}1230_1M_inca_dep.nc |
---|
| 29 | |
---|
| 30 | # Override de l'annee precedente pour 2015 : |
---|
| 31 | if [ ${year} = 2015 ] |
---|
| 32 | then |
---|
| 33 | echo '*** Comme on est en 2015, je prends 2014 de la simu historical.' |
---|
| 34 | fileprec=/ccc/store/cont003/gencmip6/lurtont/IGCM_OUT/LMDZORINCA/PROD/historicalAER/LOI6012-histAER-partE/CHM/Output/MO/LOI6012-histAER-partE_20140101_20141230_1M_inca_dep.nc |
---|
| 35 | fi |
---|
| 36 | |
---|
| 37 | # Fichiers temporaires avec ponderation (cas generique) |
---|
| 38 | pondprec=${workchemin}dust_ssp245_${prec}0101_${prec}1230_pond.nc |
---|
| 39 | pondyear=${workchemin}dust_ssp245_${year}0101_${year}1230_pond.nc |
---|
| 40 | pondsuiv=${workchemin}dust_ssp245_${suiv}0101_${suiv}1230_pond.nc |
---|
| 41 | |
---|
| 42 | # Ponderation des trois anees successives par [ 1/4 ; 1/2 ; 1/4 ] |
---|
| 43 | # et multiplication par un facteur 3 qui disparaitra lors de la moyenne par mois. |
---|
| 44 | # (je procede en deux etapes, l'enchainement cdo ne fonctionne pas) |
---|
| 45 | tmp=${workchemin}temp.nc |
---|
| 46 | cdo selname,seddustci,wetdustci,drydustci ${fileprec} ${tmp} |
---|
| 47 | cdo mulc,.75 ${tmp} ${pondprec} |
---|
| 48 | rm -f ${tmp} |
---|
| 49 | cdo selname,seddustci,wetdustci,drydustci ${fileyear} ${tmp} |
---|
| 50 | cdo mulc,1.5 ${tmp} ${pondyear} |
---|
| 51 | rm -f ${tmp} |
---|
| 52 | cdo selname,seddustci,wetdustci,drydustci ${filesuiv} ${tmp} |
---|
| 53 | cdo mulc,.75 ${tmp} ${pondsuiv} |
---|
| 54 | rm -f ${tmp} |
---|
| 55 | |
---|
| 56 | pondtroisans=${workchemin}dep_${prec}0101_${suiv}1230_pond.nc |
---|
| 57 | |
---|
| 58 | cdo mergetime ${pondprec} ${pondyear} ${pondsuiv} ${pondtroisans} |
---|
| 59 | rm -f ${pondprec} ${pondyear} ${pondsuiv} |
---|
| 60 | |
---|
| 61 | fileout=${workchemin}DUST_INCA_ssp245_1M_${year}.nc |
---|
| 62 | |
---|
| 63 | cdo ymonmean ${pondtroisans} ${fileout} |
---|
| 64 | rm -f ${pondtroisans} |
---|
| 65 | |
---|
| 66 | # Fix du vecteur-temps |
---|
| 67 | ncks -A -v time_counter,time_centered,time_counter_bounds,time_centered_bounds ${fileyear} ${fileout} |
---|
| 68 | |
---|
| 69 | # Fin de boucle sur les années |
---|
| 70 | done |
---|