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