source: TOOLS/CMIP6_FORCING/GHG/extract_scenario.sh @ 4226

#--Update Th.Lurton 05/02/2018 for paths to new version
# ** species-directory take hyphens instead of underscores previously;
# ** likewise for species names within filenames;
# ** variable names stay with underscores.

#--v3 of my script to extract global mean and put to IGCM format

#--according to Malte Meinsheausen two options are possible 
#    CFC11=CFC11eq and CCFC12=CFC12
# or CFC11=HFC134a-eq * (0.0100/62.7)/(0.062/238) [values from AR5] and CFC12=CFC12eq 
# in the second case the HFC134a-eq is converted into CFC11eq using RF efficiencies from AR5 in Wm-2pptv-1 

#--11/10/2016 New format for the data, now global, NH and SH means are in a sector dimension
#--hence I parse the first value from the cdo outputf to get the global mean
#--outputf is used with format %12.4e to make sure there is a space between values and 3 values per line

scenario='ssp126'

#--dir for annual mean data
dirin1='/prodigfs/project/input4MIPs/CMIP6/ScenarioMIP/UoM/UoM-'${scenario}'-1-1-0/atmos/yr/multiple/'
dirin2='gr1/v20180503'

dirout='./OUTPUT'
if [ ! -d ${dirout} ] ; then mkdir ${dirout} ; fi

rm -f years.txt
for i in `seq -w 2015 2500`
do
echo 'Annee_'$i'=' >> years.txt
done

for scenario in 'ssp126'
do

#--N2O we leave unit to ppb
GHG='N2O'
name='nitrous-oxide'
namevar='nitrous_oxide'
file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${scenario}'-1-1-0_gr1-GMNHSH_2015-2500.nc'
rm -f tmp.txt ${dirout}/${GHG}_CMIP6_${scenario}_2015_2500.txt
dirin=${dirin1}${dirin2}
cdo outputf,%12.4e,3 -selvar,${namevar}'_GM' ${dirin}/${file} > tmp.txt
paste years.txt tmp.txt > ${dirout}/${GHG}_CMIP6_${scenario}_2015_2500.txt

#--CH4 here we leave unit in ppb
GHG='CH4'
name='methane'
namevar='methane'
file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${scenario}'-1-1-0_gr1-GMNHSH_2015-2500.nc'
dirin=${dirin1}${dirin2}
rm -f tmp.txt ${dirout}/${GHG}_CMIP6_${scenario}_2015_2500.txt
cdo outputf,%12.4e,3 -selvar,${namevar}'_GM' ${dirin}/${file} > tmp.txt
paste years.txt tmp.txt > ${dirout}/${GHG}_CMIP6_${scenario}_2015_2500.txt

#--CO2 here we leave unit in ppm
GHG='CO2'
name='carbon-dioxide'
namevar='carbon_dioxide'
file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${scenario}'-1-1-0_gr1-GMNHSH_2015-2500.nc'
dirin=${dirin1}${dirin2}
rm -f tmp.txt ${dirout}/${GHG}_CMIP6_${scenario}_2015_2500.txt
cdo outputf,%12.4e,3 -selvar,${namevar}'_GM' ${dirin}/${file} > tmp.txt
paste years.txt tmp.txt > ${dirout}/${GHG}_CMIP6_${scenario}_2015_2500.txt

#--CFC11eq we leave in en pptv
#--not used for CMIP6
GHG='CFC11eq'
name='cfc11eq'
file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${scenario}'-1-1-0_gr1-GMNHSH_2015-2500.nc'
dirin=${dirin1}${dirin2}
rm -f tmp.txt ${dirout}/${GHG}_${scenario}_2015_2500.txt
cdo outputf,%12.4e,3 -selvar,${name}'_GM' ${dirin}/${file} > tmp.txt
paste years.txt tmp.txt > ${dirout}/${GHG}_${scenario}_2015_2500.txt

#--CFC12eq we leave in en pptv
GHG='CFC12eq'
name='cfc12eq'
file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${scenario}'-1-1-0_gr1-GMNHSH_2015-2500.nc'
dirin=${dirin1}${dirin2}
rm -f tmp.txt ${dirout}/${GHG}_${scenario}_2015_2500.txt ${dirout}/${GHG}_CMIP6_${scenario}_2015_2500.txt
cdo outputf,%12.4e,3 -selvar,${name}'_GM' ${dirin}/${file} > tmp.txt
paste years.txt tmp.txt > ${dirout}/${GHG}_${scenario}_2015_2500.txt
#--we will use CFC12eq for CMIP6
cp ${dirout}/${GHG}_${scenario}_2015_2500.txt ${dirout}/${GHG}_CMIP6_${scenario}_2015_2500.txt

#--CFC11 we leave in en pptv
#--not used for CMIP6
GHG='CFC11'
name='cfc11'
file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${scenario}'-1-1-0_gr1-GMNHSH_2015-2500.nc'
dirin=${dirin1}${dirin2}
rm -f tmp.txt ${dirout}/${GHG}_${scenario}_2015_2500.txt
cdo outputf,%12.4e,3 -selvar,${name}'_GM' ${dirin}/${file} > tmp.txt
paste years.txt tmp.txt > ${dirout}/${GHG}_${scenario}_2015_2500.txt

#--CFC12 we leave in en pptv
#--not used for CMIP6
GHG='CFC12'
name='cfc12'
file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${scenario}'-1-1-0_gr1-GMNHSH_2015-2500.nc'
dirin=${dirin1}${dirin2}
rm -f tmp.txt ${dirout}/${GHG}_${scenario}_2015_2500.txt
cdo outputf,%12.4e,3 -selvar,${name}'_GM' ${dirin}/${file} > tmp.txt
paste years.txt tmp.txt > ${dirout}/${GHG}_${scenario}_2015_2500.txt

#--HFC134aeq we leave in en pptv
GHG='HFC134aeq'
name='hfc134aeq'
file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${scenario}'-1-1-0_gr1-GMNHSH_2015-2500.nc'
dirin=${dirin1}${dirin2}
rm -f tmp.txt ${dirout}/${GHG}_${scenario}_2015_2500.txt ${dirout}/CFC11eq_CMIP6_${scenario}_2015_2500.txt
cdo outputf,%12.4e,3 -selvar,${name}'_GM' ${dirin}/${file} > tmp.txt
paste years.txt tmp.txt > ${dirout}/${GHG}_${scenario}_2015_2500.txt 
#--HFC134a-eq converted to CFC11-eq for CMIP6
#--by multiplying by (0.0100/62.7)/(0.062/238)=0.6122344
cdo outputf,%12.4e,3 -mulc,0.6122344 -selvar,${name}'_GM' ${dirin}/${file} > tmp.txt
paste years.txt tmp.txt > ${dirout}/CFC11eq_CMIP6_${scenario}_2015_2500.txt

done
#--fin boucle scenario

#--clean up
rm -f tmp.txt years.txt