source: TOOLS/CMIP6_FORCING/AER_TROP_CONCENTRATIONS/dep_from_inca_DAMIP.sh @ 4006

Last change on this file since 4006 was 4006, checked in by tlurton, 6 years ago

Adding script to generate 2015-2020 dusts for DAMIP, starting from ssp245-AER experiment output.

  • Property svn:executable set to *
File size: 2.7 KB
Line 
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
9for year in {2015..2020}
10do
11
12# Chemins
13chemin=/ccc/store/cont003/gencmip6/lurtont/IGCM_OUT/LMDZORINCA/PROD/scenarioAER/AER-ssp245/CHM/Output/MO/
14workchemin=/ccc/work/cont003/gencmip6/lurtont/aerosols_from_inca/dep_for_DAMIP/
15
16if [ ! -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 ].
20prec=$((${year}-1))
21suiv=$((${year}+1))
22
23echo '***** Annee '${year}' : Moyenne glissante sur ans '${prec}', '${year}', et '${suiv}'.'
24
25# Adresses generiques des fichiers a traiter
26fileprec=${chemin}AER-ssp245_${prec}0101_${prec}1230_1M_inca_dep.nc
27fileyear=${chemin}AER-ssp245_${year}0101_${year}1230_1M_inca_dep.nc
28filesuiv=${chemin}AER-ssp245_${suiv}0101_${suiv}1230_1M_inca_dep.nc
29
30# Override de l'annee precedente pour 2015 :
31if [ ${year} = 2015 ]
32then
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
35fi
36
37# Fichiers temporaires avec ponderation (cas generique)
38pondprec=${workchemin}dust_ssp245_${prec}0101_${prec}1230_pond.nc
39pondyear=${workchemin}dust_ssp245_${year}0101_${year}1230_pond.nc
40pondsuiv=${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)
45tmp=${workchemin}temp.nc
46cdo selname,seddustci,wetdustci,drydustci ${fileprec} ${tmp}
47cdo mulc,.75 ${tmp} ${pondprec}
48rm -f ${tmp}
49cdo selname,seddustci,wetdustci,drydustci ${fileyear} ${tmp}
50cdo mulc,1.5 ${tmp} ${pondyear}
51rm -f ${tmp}
52cdo selname,seddustci,wetdustci,drydustci ${filesuiv} ${tmp}
53cdo mulc,.75 ${tmp} ${pondsuiv}
54rm -f ${tmp}
55
56pondtroisans=${workchemin}dep_${prec}0101_${suiv}1230_pond.nc
57
58cdo mergetime ${pondprec} ${pondyear} ${pondsuiv} ${pondtroisans}
59rm -f ${pondprec} ${pondyear} ${pondsuiv}
60
61fileout=${workchemin}DUST_INCA_ssp245_1M_${year}.nc
62
63cdo ymonmean ${pondtroisans} ${fileout}
64rm -f ${pondtroisans}
65
66# Fix du vecteur-temps
67ncks -A -v time_counter,time_centered,time_counter_bounds,time_centered_bounds ${fileyear} ${fileout}
68
69# Fin de boucle sur les années
70done
Note: See TracBrowser for help on using the repository browser.