Index: TOOLS/CMIP6_FORCING/GHG/extract_scenario_v.1.2.0.sh
===================================================================
--- TOOLS/CMIP6_FORCING/GHG/extract_scenario_v.1.2.0.sh	(revision 3926)
+++ TOOLS/CMIP6_FORCING/GHG/extract_scenario_v.1.2.0.sh	(revision 3926)
@@ -0,0 +1,148 @@
+#--Update Th.Lurton 18/06/2018 for v.1.2.0
+#
+#--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
+
+#--dir for annual mean data
+dirin1='/prodigfs/project/input4MIPs/CMIP6/ScenarioMIP/UoM/'
+dirin2='gr1-GMNHSH/v20180611/'
+
+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 'ssp119' 'ssp126' 'ssp245' 'ssp370' 'ssp434' 'ssp460' 'ssp534' 'ssp585'
+do
+
+if [ ${scenario} == 'ssp119' ] ; then prefixe='IMAGE' ; fi
+if [ ${scenario} == 'ssp126' ] ; then prefixe='IMAGE' ; fi
+if [ ${scenario} == 'ssp245' ] ; then prefixe='MESSAGE-GLOBIOM' ; fi
+if [ ${scenario} == 'ssp370' ] ; then prefixe='AIM' ; fi
+if [ ${scenario} == 'ssp434' ] ; then prefixe='GCAM4' ; fi
+if [ ${scenario} == 'ssp460' ] ; then prefixe='GCAM4' ; fi
+if [ ${scenario} == 'ssp534' ] ; then prefixe='REMIND-MAGPIE' ; fi
+if [ ${scenario} == 'ssp585' ] ; then prefixe='REMIND-MAGPIE' ; fi
+
+if [ ${scenario} == 'ssp119' ] ; then suffixe='' ; fi
+if [ ${scenario} == 'ssp126' ] ; then suffixe='' ; fi
+if [ ${scenario} == 'ssp245' ] ; then suffixe='' ; fi
+if [ ${scenario} == 'ssp370' ] ; then suffixe='' ; fi
+if [ ${scenario} == 'ssp434' ] ; then suffixe='' ; fi
+if [ ${scenario} == 'ssp460' ] ; then suffixe='' ; fi
+if [ ${scenario} == 'ssp534' ] ; then suffixe='over-' ; fi
+if [ ${scenario} == 'ssp585' ] ; then suffixe='' ; fi
+
+
+#--N2O we leave unit to ppb
+GHG='N2O'
+name='nitrous-oxide'
+namevar='nitrous_oxide'
+dirin=${dirin1}'UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0/atmos/yr/mole_fraction_of_'${namevar}'_in_air/'${dirin2}
+file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0_gr1-GMNHSH_2015-2500.nc'
+rm -f tmp.txt ${dirout}${GHG}_CMIP6_${scenario}_2015_2500.txt
+cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${namevar}'_in_air' ${dirin}${file} | awk '{print $1}' > 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'
+dirin=${dirin1}'UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0/atmos/yr/mole_fraction_of_'${namevar}'_in_air/'${dirin2}
+file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0_gr1-GMNHSH_2015-2500.nc'
+rm -f tmp.txt ${dirout}${GHG}_CMIP6_${scenario}_2015_2500.txt
+cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${namevar}'_in_air' ${dirin}${file} | awk '{print $1}' > 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'
+dirin=${dirin1}'UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0/atmos/yr/mole_fraction_of_'${namevar}'_in_air/'${dirin2}
+file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0_gr1-GMNHSH_2015-2500.nc'
+rm -f tmp.txt ${dirout}${GHG}_CMIP6_${scenario}_2015_2500.txt
+cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${namevar}'_in_air' ${dirin}${file} | awk '{print $1}' > tmp.txt
+paste years.txt tmp.txt > ${dirout}${GHG}_CMIP6_${scenario}_2015_2500.txt
+
+#--CFC11eq we leave in pptv
+#--not used for CMIP6
+GHG='CFC11eq'
+name='cfc11eq'
+namevar=${name}
+dirin=${dirin1}'UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0/atmos/yr/mole_fraction_of_'${namevar}'_in_air/'${dirin2}
+file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0_gr1-GMNHSH_2015-2500.nc'
+rm -f tmp.txt ${dirout}${GHG}_${scenario}_2015_2500.txt
+cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${namevar}'_in_air' ${dirin}${file} | awk '{print $1}' > tmp.txt
+paste years.txt tmp.txt > ${dirout}${GHG}_${scenario}_2015_2500.txt
+
+#--CFC12eq we leave in pptv
+GHG='CFC12eq'
+name='cfc12eq'
+namevar=${name}
+dirin=${dirin1}'UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0/atmos/yr/mole_fraction_of_'${namevar}'_in_air/'${dirin2}
+file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0_gr1-GMNHSH_2015-2500.nc'
+rm -f tmp.txt ${dirout}${GHG}_${scenario}_2015_2500.txt ${dirout}/${GHG}_CMIP6_${scenario}_2015_2500.txt
+cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${namevar}'_in_air' ${dirin}${file} | awk '{print $1}' > 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 pptv
+#--not used for CMIP6
+GHG='CFC11'
+name='cfc11'
+namevar=${name}
+dirin=${dirin1}'UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0/atmos/yr/mole_fraction_of_'${namevar}'_in_air/'${dirin2}
+file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0_gr1-GMNHSH_2015-2500.nc'
+rm -f tmp.txt ${dirout}${GHG}_${scenario}_2015_2500.txt
+cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${namevar}'_in_air' ${dirin}${file} | awk '{print $1}' > tmp.txt
+paste years.txt tmp.txt > ${dirout}${GHG}_${scenario}_2015_2500.txt
+
+#--CFC12 we leave in pptv
+#--not used for CMIP6
+GHG='CFC12'
+name='cfc12'
+namevar=${name}
+dirin=${dirin1}'UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0/atmos/yr/mole_fraction_of_'${namevar}'_in_air/'${dirin2}
+file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0_gr1-GMNHSH_2015-2500.nc'
+rm -f tmp.txt ${dirout}${GHG}_${scenario}_2015_2500.txt
+cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${namevar}'_in_air' ${dirin}${file} | awk '{print $1}' > tmp.txt
+paste years.txt tmp.txt > ${dirout}${GHG}_${scenario}_2015_2500.txt
+
+#--HFC134aeq we leave in pptv
+GHG='HFC134aeq'
+name='hfc134aeq'
+namevar=${name}
+dirin=${dirin1}'UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0/atmos/yr/mole_fraction_of_'${namevar}'_in_air/'${dirin2}
+file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_ScenarioMIP_UoM-'${prefixe}'-'${scenario}'-'${suffixe}'1-2-0_gr1-GMNHSH_2015-2500.nc'
+rm -f tmp.txt ${dirout}${GHG}_${scenario}_2015_2500.txt ${dirout}/CFC11eq_CMIP6_${scenario}_2015_2500.txt
+cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${namevar}'_in_air' ${dirin}${file} | awk '{print $1}' > 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,'mole_fraction_of_'${namevar}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt
+paste years.txt tmp.txt > ${dirout}CFC11eq_CMIP6_${scenario}_2015_2500.txt
+
+done
+#--fin boucle sur scenarios
+
+#--clean up
+rm -f tmp.txt years.txt
+
