Changeset 3928

Timestamp:

06/21/18 09:35:33 (5 years ago)

Author:

tlurton

Message:

Script for scenario GHG extraction: updated for subdirectory organisation by scenario.

File:

• TOOLS/CMIP6_FORCING/GHG/extract_scenario_v.1.2.0.sh (modified) (10 diffs)

TOOLS/CMIP6_FORCING/GHG/extract_scenario_v.1.2.0.sh

@@ -33 +33 @@
 do
 
+rep=${scenario}'/'
+
 if [ ${scenario} == 'ssp119' ] ; then prefixe='IMAGE' ; fi
 if [ ${scenario} == 'ssp126' ] ; then prefixe='IMAGE' ; fi
@@ -51 +53 @@
 if [ ${scenario} == 'ssp585' ] ; then suffixe='' ; fi
 
+if [ ! -d ${dirout}${rep} ] ; then mkdir -p ${dirout}${rep} ; fi
 
 #--N2O we leave unit to ppb
@@ -60 +63 @@
 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
+paste years.txt tmp.txt > ${dirout}${rep}${GHG}_CMIP6_2015_2500.txt
 
 #--CH4 here we leave unit in ppb
@@ -70 +73 @@
 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
+paste years.txt tmp.txt > ${dirout}${rep}${GHG}_CMIP6_2015_2500.txt
 
 #--CO2 here we leave unit in ppm
@@ -80 +83 @@
 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
+paste years.txt tmp.txt > ${dirout}${rep}${GHG}_CMIP6_2015_2500.txt
 
 #--CFC11eq we leave in pptv
@@ -91 +94 @@
 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
+paste years.txt tmp.txt > ${dirout}${rep}${GHG}_2015_2500.txt
 
 #--CFC12eq we leave in pptv
@@ -101 +104 @@
 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
+paste years.txt tmp.txt > ${dirout}${rep}${GHG}_2015_2500.txt
 #--we will use CFC12eq for CMIP6
-cp ${dirout}/${GHG}_${scenario}_2015_2500.txt ${dirout}/${GHG}_CMIP6_${scenario}_2015_2500.txt
+cp ${dirout}${rep}${GHG}_2015_2500.txt ${dirout}${rep}${GHG}_CMIP6_2015_2500.txt
 
 #--CFC11 we leave in pptv
@@ -114 +117 @@
 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
+paste years.txt tmp.txt > ${dirout}${rep}${GHG}_2015_2500.txt
 
 #--CFC12 we leave in pptv
@@ -125 +128 @@
 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
+paste years.txt tmp.txt > ${dirout}${rep}${GHG}_2015_2500.txt
 
 #--HFC134aeq we leave in pptv
@@ -135 +138 @@
 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 
+paste years.txt tmp.txt > ${dirout}${rep}${GHG}_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
+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}${rep}CFC11eq_CMIP6_2015_2500.txt
 
 done