Changeset 4051

Timestamp:

10/11/18 13:52:39 (6 years ago)

Author:

tlurton

Message:

Corrected script to generate sflx files for scenarios (bug on 2080's decade).

File:

• TOOLS/CMIP6_FORCING/SCENARIOS/AER_TROP_EMISSIONS/REGRID/regrid_scen.sh (modified) (16 diffs)

TOOLS/CMIP6_FORCING/SCENARIOS/AER_TROP_EMISSIONS/REGRID/regrid_scen.sh

@@ +1 @@
+#--corrected 2080 bug 9/2018, ThL
 #--reworked 7/2018 for scenario emissions, ThL
 #--updated on 5/2/2018 with new paths, ThL
@@ -16 +17 @@
 fileINCAex="/home/oboucher/CMIP6/AER_EMISSIONS/INCAfile/sflx_lmdz_phy_1997.nc"
 
-#--input directory for emissions
+#--input directory for anthropogenic (non-BB) emissions
 dirin="/prodigfs/project/input4MIPs/CMIP6/ScenarioMIP/IAMC/"
 
@@ -28 +29 @@
 
 #--output directory
-dirout="/data/"${USER}"/CMIP6/AEROSOL/ScenarioMIP/${scen}/v1/"
+dirout="/data/"${USER}"/CMIP6/AEROSOL/ScenarioMIP/${scen}/v2_correcte/"
 if [ ! -d ${dirout} ] ; then mkdir -p ${dirout} ; fi
 
@@ -40 +41 @@
 if [ ${scen} == "ssp585" ] ; then prefix="REMIND-MAGPIE" ; fi
 
-#--Unique file for the whole period, so fixed variables:
+#--loop on years
+for year in {2015..2100}
+do
+
+if [ ${year} -ge "2015" -a ${year} -lt "2020" ]
+then
 year1=2015
+year2=2020
+fi
+if [ ${year} -ge "2020" -a ${year} -lt "2030" ]
+then
+year1=2020
+year2=2030
+fi
+if [ ${year} -ge "2030" -a ${year} -lt "2040" ]
+then
+year1=2030
+year2=2040
+fi
+if [ ${year} -ge "2040" -a ${year} -lt "2050" ]
+then
+year1=2040
+year2=2050
+fi
+if [ ${year} -ge "2050" -a ${year} -lt "2060" ]
+then
+year1=2050
+year2=2060
+fi
+if [ ${year} -ge "2060" -a ${year} -lt "2070" ]
+then
+year1=2060
+year2=2070
+fi
+if [ ${year} -ge "2070" -a ${year} -lt "2080" ]
+then
+year1=2070
+year2=2080
+fi
+if [ ${year} -ge "2080" -a ${year} -lt "2090" ]
+then
+year1=2080
+year2=2090
+fi
+if [ ${year} -ge "2090" -a ${year} -lt "2100" ]
+then
+year1=2090
 year2=2100
-
-#--Loop on species (anthro)
+fi
+if [ ${year} == 2100 ]
+then
+year1=2100
+year2=2100
+fi
+#--weights for linear interpolation
+if [ ${year} == 2015 ]
+then
+pond1="1"
+pond2="0"
+fi
+if [ ${year} == 2016 ]
+then
+pond1="0.8"
+pond2="0.2"
+fi
+if [ ${year} == 2017 ]
+then
+pond1="0.6"
+pond2="0.4"
+fi
+if [ ${year} == 2018 ]
+then
+pond1="0.4"
+pond2="0.6"
+fi
+if [ ${year} == 2019 ]
+then
+pond1="0.2"
+pond2="0.8"
+fi
+if [ ${year} == 2020 -o ${year} == 2030 -o ${year} == 2040 -o ${year} == 2050 -o ${year} == 2060 -o ${year} == 2070 -o ${year} == 2080 -o ${year} == 2090 ]
+then
+pond1="1"
+pond2="0"
+fi
+if [ ${year} == 2021 -o ${year} == 2031 -o ${year} == 2041 -o ${year} == 2051 -o ${year} == 2061 -o ${year} == 2071 -o ${year} == 2081 -o ${year} == 2091 ]
+then
+pond1="0.9"
+pond2="0.1"
+fi
+if [ ${year} == 2022 -o ${year} == 2032 -o ${year} == 2042 -o ${year} == 2052 -o ${year} == 2062 -o ${year} == 2072 -o ${year} == 2082 -o ${year} == 2092 ]
+then
+pond1="0.8"
+pond2="0.2"
+fi
+if [ ${year} == 2023 -o ${year} == 2033 -o ${year} == 2043 -o ${year} == 2053 -o ${year} == 2063 -o ${year} == 2073 -o ${year} == 2083 -o ${year} == 2093 ]
+then
+pond1="0.7"
+pond2="0.3"
+fi
+if [ ${year} == 2024 -o ${year} == 2034 -o ${year} == 2044 -o ${year} == 2054 -o ${year} == 2064 -o ${year} == 2074 -o ${year} == 2084 -o ${year} == 2094 ]
+then
+pond1="0.6"
+pond2="0.4"
+fi
+if [ ${year} == 2025 -o ${year} == 2035 -o ${year} == 2045 -o ${year} == 2055 -o ${year} == 2065 -o ${year} == 2075 -o ${year} == 2085 -o ${year} == 2095 ]
+then
+pond1="0.5"
+pond2="0.5"
+fi
+if [ ${year} == 2026 -o ${year} == 2036 -o ${year} == 2046 -o ${year} == 2056 -o ${year} == 2066 -o ${year} == 2076 -o ${year} == 2086 -o ${year} == 2096 ]
+then
+pond1="0.4"
+pond2="0.6"
+fi
+if [ ${year} == 2027 -o ${year} == 2037 -o ${year} == 2047 -o ${year} == 2057 -o ${year} == 2067 -o ${year} == 2077 -o ${year} == 2087 -o ${year} == 2097 ]
+then
+pond1="0.3"
+pond2="0.7"
+fi
+if [ ${year} == 2028 -o ${year} == 2038 -o ${year} == 2048 -o ${year} == 2058 -o ${year} == 2068 -o ${year} == 2078 -o ${year} == 2088 -o ${year} == 2098 ]
+then
+pond1="0.2"
+pond2="0.8"
+fi
+if [ ${year} == 2029 -o ${year} == 2039 -o ${year} == 2049 -o ${year} == 2059 -o ${year} == 2069 -o ${year} == 2079 -o ${year} == 2089 -o ${year} == 2099 ]
+then
+pond1="0.1"
+pond2="0.9"
+fi
+if [ ${year} == 2100 ]
+then
+pond1="0"
+pond2="1"
+fi
+
+#--species 
 for species in "BC" "NOx" "OC" "SO2" "NH3"
 do
-
-echo "... "${year}" : Dealing with "${species}" anthro input file..."
 
 if [ ${species} = "BC"  ] ; then speciesinca="bc"  ; fi
@@ -62 +193 @@
 if [ ${species} = "NH3" ] ; then speciesUp="NH3" ; fi
 
-#--input file anthro
-filename=${dirin}IAMC-${prefix}-${scen}-1-1/atmos/mon/${species}_em_anthro/gn/v20180628/${species}-em-anthro_input4MIPs_emissions_ScenarioMIP_IAMC-${prefix}-${scen}-1-1_gn_${year1}01-${year2}12.nc
-
-#--input file summed over sectors
-filesum=${dirout}${species}-em-anthro_input4MIPs_emissions_ScenarioMIP_IAMC-${prefix}-${scen}-1-1_gn_${year1}01-${year2}12_sum.nc
-#--reworked input file
-filetmp=${dirout}${species}-em-anthro_input4MIPs_emissions_ScenarioMIP_IAMC-${prefix}-${scen}-1-1_gn_${year1}01-${year2}12_tmp.nc
-#--time-interpolated input file
-fileintanthro=${dirout}${species}-em-anthro_input4MIPs_emissions_ScenarioMIP_IAMC-${prefix}-${scen}-1-1_gn_${year1}01-${year2}12_interp.nc
-
-if [ -f ${fileintanthro} ]
-then
-   echo "...Nothing to do!"
-else
-
-rm -f ${filesum} ${filetmp} ${fileintanthro}
-
-#--re-formatting the input file
-#--summing over sectors
-rm -f sumsectors.jnl
-cat << eod > sumsectors.jnl
-set memory/size=100
-use "${filename}"
-save/clobber/file="${filesum}" ${species}_EM_ANTHRO[k=@sum]
-eod
-#--run ferret script
-ferret << eod 
-go sumsectors.jnl
-exit
-eod
-
-#--reformatting input file for IDL to be able to read
-ncks -3 ${filesum} ${filetmp}
-#--cleaning up
-rm -f ${filesum}
-
-#--IDL code for time interpolation, to fill in 9-year gaps in input file
-rm -f time_interp_over_gaps.pro
-cat << eod >> time_interp_over_gaps.pro
-pro time_interp_over_gaps
-
-; Les dimensions des fichiers
-nlon=720
-nlat=360
-nt=1032
-
-une_annee=[31,30,29,31,30,31,30,31,31,30,31,30]
-les_mois = [31,28,31,30,31,30,31,31,30,31,30,31]
-
-jours_in = [15, 45, 74, 105, 135, 166, 196, 227, 258, 288, 319, 349, 1840, 1870, $
-    1899, 1930, 1960, 1991, 2021, 2052, 2083, 2113, 2144, 2174, 5490, 5520, $
-    5549, 5580, 5610, 5641, 5671, 5702, 5733, 5763, 5794, 5824, 9140, 9170, $
-    9199, 9230, 9260, 9291, 9321, 9352, 9383, 9413, 9444, 9474, 12790, 12820, $
-    12849, 12880, 12910, 12941, 12971, 13002, 13033, 13063, 13094, 13124, $
-    16440, 16470, 16499, 16530, 16560, 16591, 16621, 16652, 16683, 16713, $
-    16744, 16774, 20090, 20120, 20149, 20180, 20210, 20241, 20271, 20302, $
-    20333, 20363, 20394, 20424, 23740, 23770, 23799, 23830, 23860, 23891, $
-    23921, 23952, 23983, 24013, 24044, 24074, 27390, 27420, 27449, 27480, $ 
-    27510, 27541, 27571, 27602, 27633, 27663, 27694, 27724, 31040, 31070, $ 
-    31099, 31130, 31160, 31191, 31221, 31252, 31283, 31313, 31344, 31374]
-
-jours_out=make_array(1032)
-jours_out(0) = 15
-for cpt=1,1031 do jours_out(cpt) = jours_out(cpt-1) + une_annee(cpt mod 12)
-
-idin = ncdf_open('${filetmp}')
-varid = ncdf_varid(idin,'${speciesUp}_EM_ANTHRO')
-lonid = ncdf_varid(idin,'LON')
-latid = ncdf_varid(idin,'LAT')
-;lonbndid = ncdf_varid(idin,'LON_bnds')
-;latbndid = ncdf_varid(idin,'LAT_bnds')
-timbndid = ncdf_varid(idin,'TIME_bnds')
-ncdf_varget, idin, varid, emiss_init
-ncdf_attget, idin, varid, "units", emiss_unit
-ncdf_attget, idin, varid, "long_name", emiss_longname
-ncdf_varget, idin, lonid, lon
-ncdf_varget, idin, latid, lat
-;ncdf_varget, idin, lonbndid, lonbnds
-;ncdf_varget, idin, latbndid, latbnds
-ncdf_varget, idin, timbndid, timbnds_in
-ncdf_close, idin
-
-timbnds_out = make_array(2,1032)
-timbnds_out(*,0) = [0,31]
-timbnds_out(*,1) = [31,59]
-for cpt=2,1031 do timbnds_out(0,cpt) = timbnds_out(0,cpt-1) + les_mois((cpt-1) mod 12)
-for cpt=1,1030 do timbnds_out(1,cpt) = timbnds_out(0,cpt+1)
-timbnds_out(1,1031) = 31390
-
-emiss_interpol=make_array(nlon,nlat,1032,value=0.)
-
-indices_out_mois = make_array(86,12,value=0)
-for m=0,11 do begin
-   for n=0,85 do indices_out_mois(n,m) = 12*n + m
-endfor
-
-jours_out_mois = make_array(86,12,value=0)
-for n=0,85 do begin
-   for m=0,11 do jours_out_mois(n,m) = jours_out(indices_out_mois(n,m))
-endfor
-
-jours_in_mois = make_array(n_elements(jours_in)/12,12,value=0)
-for m=0,11 do begin
-   for n=0,n_elements(jours_in)/12-1 do jours_in_mois(n,m) = jours_in(n*12 + m)
-endfor
-
-indices_in_mois = make_array(n_elements(jours_in)/12,12,value=0)
-for m=0,11 do begin
-   for n=0,n_elements(jours_in)/12-1 do indices_in_mois(n,m) = where(jours_in eq jours_in_mois(n,m))
-endfor
-
-; Interpolation mois par mois...
-for i=0,nlon-1 do begin
-   for j=0,nlat-1 do begin
-      for m=0,11 do begin
-         emiss_interpol(i,j,indices_out_mois(*,m)) = interpol(emiss_init(i,j,indices_in_mois(*,m)),jours_in_mois(*,m),jours_out_mois(*,m))
-      endfor
-   endfor
-endfor
-
-idout = ncdf_create('${fileintanthro}', /clobber)
-dimlonid = ncdf_dimdef(idout,'lon',nlon)
-dimlatid = ncdf_dimdef(idout,'lat',nlat)
-dimtimid = ncdf_dimdef(idout,'time',1032)
-dimbndid = ncdf_dimdef(idout,'bound',2)
-varlonid = ncdf_vardef(idout, 'lon', [dimlonid], /double)
-varlatid = ncdf_vardef(idout, 'lat', [dimlatid], /double)
-vartimid = ncdf_vardef(idout, 'time', [dimtimid], /double)
-;varlonbndsid = ncdf_vardef(idout, 'lon_bnds', [dimbndid,dimlonid], /double)
-;varlatbndsid = ncdf_vardef(idout, 'lat_bnds', [dimbndid,dimlatid], /double)
-vartimbndsid = ncdf_vardef(idout, 'time_bnds', [dimbndid,dimtimid], /double)
-varemissid = ncdf_vardef(idout, '${species}_em_anthro', [dimlonid,dimlatid,dimtimid], /float)
-ncdf_control, idout, /endef
-ncdf_varput, idout, varlonid, lon
-ncdf_varput, idout, varlatid, lat
-;ncdf_varput, idout, varlonbndsid, lonbnds
-;ncdf_varput, idout, varlatbndsid, latbnds
-ncdf_varput, idout, vartimbndsid, timbnds_out
-ncdf_varput, idout, varemissid, emiss_interpol
-ncdf_varput, idout, vartimid, jours_out
-ncdf_control, idout, /redef
-ncdf_attput, idout, varlonid, "units", "degrees_east"
-ncdf_attput, idout, varlonid, "long_name", "longitude"
-ncdf_attput, idout, varlonid, "axis", "X"
-ncdf_attput, idout, varlonid, "bounds", "lon_bnds"
-ncdf_attput, idout, varlonid, "modulo", 360.
-ncdf_attput, idout, varlonid, "realtopology", "circular"
-ncdf_attput, idout, varlonid, "standard_name", "longitude"
-ncdf_attput, idout, varlonid, "topology", "circular"
-ncdf_attput, idout, varlatid, "units", "degrees_north"
-ncdf_attput, idout, varlatid, "long_name", "latitude"
-ncdf_attput, idout, varlatid, "axis", "Y"
-ncdf_attput, idout, varlatid, "bounds", "lat_bnds"
-ncdf_attput, idout, varlatid, "realtopology", "linear"
-ncdf_attput, idout, varlatid, "standard_name", "latitude"
-ncdf_attput, idout, varemissid, "units", string(emiss_unit)
-ncdf_attput, idout, varemissid, "_FillValue", 1.e+20
-ncdf_attput, idout, varemissid, "cell_methods", "time: mean"
-ncdf_attput, idout, varemissid, "long_name", string(emiss_longname)
-ncdf_attput, idout, varemissid, "missing_value", 1.e+20
-ncdf_attput, idout, vartimid, "units", "days since 2015-01-01 0:0:0"
-ncdf_attput, idout, vartimid, "long_name", "time"
-ncdf_attput, idout, vartimid, "calendar", '365_day'
-ncdf_attput, idout, vartimid, "axis", "T"
-ncdf_attput, idout, vartimid, "bounds", "time_bnds"
-ncdf_attput, idout, vartimid, "realtopology", "linear"
-ncdf_attput, idout, vartimid, "standard_name", "time"
-ncdf_close, idout
-end
-eod
-#
-#--calling IDL
-#
-/opt/idl-6.4/idl/bin/idl << eod
-.r time_interp_over_gaps.pro
-time_interp_over_gaps
-exit
-eod
-#
-
-#--adding lat_bnds and lon_bnds variables
-ncks -A -v lat_bnds,lon_bnds ${filename} ${fileintanthro}
-
-#--cleaning up
-rm -f ${filetmp}
-
-#--end test whether file tmp already exists
-fi
-
-
-
-#--now dealing with the BB input file
-
-echo "... "${year}" : Dealing with "${species}" openburning input file..."
-
-#--input file BB
-filename=${dirin}IAMC-${prefix}-${scen}-1-1/atmos/mon/${species}_em_openburning/gn/v20180628/${species}-em-openburning_input4MIPs_emissions_ScenarioMIP_IAMC-${prefix}-${scen}-1-1_gn_${year1}01-${year2}12.nc
-#--reworked input file
-filetmp=${dirout}${species}-em-openburning_input4MIPs_emissions_ScenarioMIP_IAMC-${prefix}-${scen}-1-1_gn_${year1}01-${year2}12_tmp.nc
-#--time-interpolated input file
-fileintbb=${dirout}${species}-em-openburning_input4MIPs_emissions_ScenarioMIP_IAMC-${prefix}-${scen}-1-1_gn_${year1}01-${year2}12_interp.nc
-
-if [ -f ${fileintbb} ]
-then
-   echo "...Nothing to do!"
-else
-
-rm -f ${filetmp} ${fileintbb}
-
-#--no need to sum over sectors, as none in BB files
-
-#--reformatting input file for IDL to be able to read
-ncks -3 ${filename} ${filetmp}
-
-#--IDL code for time interpolation, to fill in 9-year gaps in input file
-rm -f time_interp_over_gaps.pro
-cat << eod >> time_interp_over_gaps.pro
-pro time_interp_over_gaps
-
-; Les dimensions des fichiers
-nlon=720
-nlat=360
-nt=1032
-
-une_annee=[31,30,29,31,30,31,30,31,31,30,31,30]
-les_mois = [31,28,31,30,31,30,31,31,30,31,30,31]
-
-jours_in = [15, 45, 74, 105, 135, 166, 196, 227, 258, 288, 319, 349, 1840, 1870, $
-    1899, 1930, 1960, 1991, 2021, 2052, 2083, 2113, 2144, 2174, 5490, 5520, $
-    5549, 5580, 5610, 5641, 5671, 5702, 5733, 5763, 5794, 5824, 9140, 9170, $
-    9199, 9230, 9260, 9291, 9321, 9352, 9383, 9413, 9444, 9474, 12790, 12820, $
-    12849, 12880, 12910, 12941, 12971, 13002, 13033, 13063, 13094, 13124, $
-    16440, 16470, 16499, 16530, 16560, 16591, 16621, 16652, 16683, 16713, $
-    16744, 16774, 20090, 20120, 20149, 20180, 20210, 20241, 20271, 20302, $
-    20333, 20363, 20394, 20424, 23740, 23770, 23799, 23830, 23860, 23891, $
-    23921, 23952, 23983, 24013, 24044, 24074, 27390, 27420, 27449, 27480, $ 
-    27510, 27541, 27571, 27602, 27633, 27663, 27694, 27724, 31040, 31070, $ 
-    31099, 31130, 31160, 31191, 31221, 31252, 31283, 31313, 31344, 31374]
-
-jours_out=make_array(1032)
-jours_out(0) = 15
-for cpt=1,1031 do jours_out(cpt) = jours_out(cpt-1) + une_annee(cpt mod 12)
-
-idin = ncdf_open('${filetmp}')
-varid = ncdf_varid(idin,'${species}_em_openburning')      ; Here, small caps
-lonid = ncdf_varid(idin,'lon')
-latid = ncdf_varid(idin,'lat')
-;lonbndid = ncdf_varid(idin,'lon_bnds')
-;latbndid = ncdf_varid(idin,'lat_bnds')
-timbndid = ncdf_varid(idin,'time_bnds')
-ncdf_varget, idin, varid, emiss_init
-ncdf_attget, idin, varid, "units", emiss_unit
-ncdf_attget, idin, varid, "long_name", emiss_longname
-ncdf_varget, idin, lonid, lon
-ncdf_varget, idin, latid, lat
-;ncdf_varget, idin, lonbndid, lonbnds
-;ncdf_varget, idin, latbndid, latbnds
-ncdf_varget, idin, timbndid, timbnds_in
-ncdf_close, idin
-
-timbnds_out = make_array(2,1032)
-timbnds_out(*,0) = [0,31]
-timbnds_out(*,1) = [31,59]
-for cpt=2,1031 do timbnds_out(0,cpt) = timbnds_out(0,cpt-1) + une_annee((cpt-1) mod 12)
-for cpt=1,1030 do timbnds_out(1,cpt) = timbnds_out(0,cpt+1)
-timbnds_out(1,1031) = 31390
-
-emiss_interpol=make_array(nlon,nlat,1032,value=0.)
-
-indices_out_mois = make_array(86,12,value=0)
-for m=0,11 do begin
-   for n=0,85 do indices_out_mois(n,m) = 12*n + m
-endfor
-
-jours_out_mois = make_array(86,12,value=0)
-for n=0,85 do begin
-   for m=0,11 do jours_out_mois(n,m) = jours_out(indices_out_mois(n,m))
-endfor
-
-jours_in_mois = make_array(n_elements(jours_in)/12,12,value=0)
-for m=0,11 do begin
-   for n=0,n_elements(jours_in)/12-1 do jours_in_mois(n,m) = jours_in(n*12 + m)
-endfor
-
-indices_in_mois = make_array(n_elements(jours_in)/12,12,value=0)
-for m=0,11 do begin
-   for n=0,n_elements(jours_in)/12-1 do indices_in_mois(n,m) = where(jours_in eq jours_in_mois(n,m))
-endfor
-
-; Interpolation mois par mois...
-for i=0,nlon-1 do begin
-   for j=0,nlat-1 do begin
-      for m=0,11 do begin
-         emiss_interpol(i,j,indices_out_mois(*,m)) = interpol(emiss_init(i,j,indices_in_mois(*,m)),jours_in_mois(*,m),jours_out_mois(*,m))
-      endfor
-   endfor
-endfor
-
-idout = ncdf_create('${fileintbb}', /clobber)
-dimlonid = ncdf_dimdef(idout,'lon',nlon)
-dimlatid = ncdf_dimdef(idout,'lat',nlat)
-dimtimid = ncdf_dimdef(idout,'time',1032)
-dimbndid = ncdf_dimdef(idout,'bound',2)
-varlonid = ncdf_vardef(idout, 'lon', [dimlonid], /double)
-varlatid = ncdf_vardef(idout, 'lat', [dimlatid], /double)
-vartimid = ncdf_vardef(idout, 'time', [dimtimid], /double)
-;varlonbndsid = ncdf_vardef(idout, 'lon_bnds', [dimbndid,dimlonid], /double)
-;varlatbndsid = ncdf_vardef(idout, 'lat_bnds', [dimbndid,dimlatid], /double)
-vartimbndsid = ncdf_vardef(idout, 'time_bnds', [dimbndid,dimtimid], /double)
-varemissid = ncdf_vardef(idout, '${species}_em_openburning', [dimlonid,dimlatid,dimtimid], /float)
-ncdf_control, idout, /endef
-ncdf_varput, idout, varlonid, lon
-ncdf_varput, idout, varlatid, lat
-;ncdf_varput, idout, varlonbndsid, lonbnds
-;ncdf_varput, idout, varlatbndsid, latbnds
-ncdf_varput, idout, vartimbndsid, timbnds_out
-ncdf_varput, idout, varemissid, emiss_interpol
-ncdf_varput, idout, vartimid, jours_out
-ncdf_control, idout, /redef
-ncdf_attput, idout, varlonid, "units", "degrees_east"
-ncdf_attput, idout, varlonid, "long_name", "longitude"
-ncdf_attput, idout, varlonid, "axis", "X"
-ncdf_attput, idout, varlonid, "bounds", "lon_bnds"
-ncdf_attput, idout, varlonid, "modulo", 360.
-ncdf_attput, idout, varlonid, "realtopology", "circular"
-ncdf_attput, idout, varlonid, "standard_name", "longitude"
-ncdf_attput, idout, varlonid, "topology", "circular"
-ncdf_attput, idout, varlatid, "units", "degrees_north"
-ncdf_attput, idout, varlatid, "long_name", "latitude"
-ncdf_attput, idout, varlatid, "axis", "Y"
-ncdf_attput, idout, varlatid, "bounds", "lat_bnds"
-ncdf_attput, idout, varlatid, "realtopology", "linear"
-ncdf_attput, idout, varlatid, "standard_name", "latitude"
-ncdf_attput, idout, varemissid, "units", string(emiss_unit)
-ncdf_attput, idout, varemissid, "_FillValue", 1.e+20
-ncdf_attput, idout, varemissid, "cell_methods", "time: mean"
-ncdf_attput, idout, varemissid, "long_name", string(emiss_longname)
-ncdf_attput, idout, varemissid, "missing_value", 1.e+20
-ncdf_attput, idout, vartimid, "units", "days since 2015-01-01 0:0:0"
-ncdf_attput, idout, vartimid, "long_name", "time"
-ncdf_attput, idout, vartimid, "calendar", '365_day'
-ncdf_attput, idout, vartimid, "axis", "T"
-ncdf_attput, idout, vartimid, "bounds", "time_bnds"
-ncdf_attput, idout, vartimid, "realtopology", "linear"
-ncdf_attput, idout, vartimid, "standard_name", "time"
-ncdf_close, idout
-end
-eod
-#
-#--calling IDL
-#
-/opt/idl-6.4/idl/bin/idl << eod
-.r time_interp_over_gaps.pro
-time_interp_over_gaps
-exit
-eod
-#
-
-#--adding lat_bnds and lon_bnds variables
-ncks -A -v lat_bnds,lon_bnds ${filename} ${fileintbb}
-
-#--cleaning up
-rm -f ${filetmp}
-
-#--end test if interpolated file already exists
-fi
-
-#--end loop on species
-done
-
-
-
-#--anthro and BB input files fileintanthro and fileintbb are now ready, for the whole period and for the current scenario
-#--new double loop on species and years
-
-#--loop on years
-for year in {2015..2100}
-do
-
-#--loop on species
-for species in "BC" "NOx" "OC" "SO2" "NH3"
-do
-
-if [ ${species} = "BC"  ] ; then speciesinca="bc"  ; fi
-if [ ${species} = "NOx" ] ; then speciesinca="nox" ; fi
-if [ ${species} = "OC"  ] ; then speciesinca="pom" ; fi
-if [ ${species} = "SO2" ] ; then speciesinca="so2" ; fi
-if [ ${species} = "NH3" ] ; then speciesinca="nh3" ; fi
-
-#--dealing with anthro
-
-#--input file
-fileintanthro=${dirout}${species}-em-anthro_input4MIPs_emissions_ScenarioMIP_IAMC-${prefix}-${scen}-1-1_gn_${year1}01-${year2}12_interp.nc
+#--Dealing with anthro emissions file...
+
+echo "... "${year}" : Dealing with "${species}" anthro input file..."
+
+#--input file anthro (fossil fuel emissions)
+filename=${dirin}IAMC-${prefix}-${scen}-1-1/atmos/mon/${species}_em_anthro/gn/v20180628/${species}-em-anthro_input4MIPs_emissions_ScenarioMIP_IAMC-${prefix}-${scen}-1-1_gn_201501-210012.nc
+
+#--two temporary files for interpolation
+filetmp1=${dirout}tmp1.nc
+filetmp2=${dirout}tmp2.nc
+
 #--output files
 filenameout1=${dirout}flux_${speciesinca}_${year}.nc
@@ -461 +210 @@
 filenameout3=${dirout}flux_vector_${speciesinca}_${year}.nc
 
-echo ${fileintanthro} ${filenameout1} ${filenameout2} ${filenameout3}
-rm -f ${filenameout1} ${filenameout2} ${filenameout3}
-
-#--extracting the correct year
-rm -f select_yr.jnl
-cat << eod > select_yr.jnl
-set memory/size=100
-use "${fileintanthro}"
-set region/t=16-jan-${year}:16-dec-${year}
-save/clobber/file="${filenameout1}" ${species}_EM_ANTHRO
+echo ${filename} ${filenameout1} ${filenameout2} ${filenameout3}
+rm -f ${filetmp1} ${filetmp2} ${filenameout1} ${filenameout2} ${filenameout3}
+
+#--unfortunately idl not happy with input netcdf files so need to ferretize files
+#--I also sum over sectors and I extract the correct years as well
+#--first year (start of 5-Y or 10-Y period, for later interpolation)
+rm -f rewrite.jnl
+cat << eod > rewrite.jnl
+use "${filename}"
+set region/t=16-jan-${year1}:16-dec-${year1}
+save/clobber/file="${filetmp1}" ${speciesUp}_EM_ANTHRO[k=@sum]
+eod
+#--run ferret script
+ferret << eod
+go rewrite.jnl
+exit
+eod
+#--second year (end of 5-Y or 10-Y period, for later interpolation)
+rm -f rewrite.jnl
+cat << eod > rewrite.jnl
+use "${filename}"
+set region/t=16-jan-${year2}:16-dec-${year2}
+save/clobber/file="${filetmp2}" ${speciesUp}_EM_ANTHRO[k=@sum]
 eod
 #--run ferret script
 ferret << eod 
-go select_yr.jnl
+go rewrite.jnl
 exit
 eod
 
-#--replace undef with 0
+#--performing time interpolation
+cdo add -mulc,${pond1} -selname,${speciesUp}_EM_ANTHRO ${filetmp1} -mulc,${pond2} -selname,${speciesUp}_EM_ANTHRO ${filetmp2} ${filenameout1}
+
 cdo setmisstoc,0.0 ${filenameout1} ${filenameout1b}
 
@@ -624 +388 @@
 end
 eod
+
 #
 #--calling IDL
@@ -637 +402 @@
 #
 
-#--cleaning the mess
-rm -f ferret*
-rm -f regrid.pro
-
-
-#--now dealing with BB
-
-#--input file
-fileintbb=${dirout}${species}-em-openburning_input4MIPs_emissions_ScenarioMIP_IAMC-${prefix}-${scen}-1-1_gn_${year1}01-${year2}12_interp.nc
+
+#--now dealing with BB sources
+echo "... "${year}" : Dealing with "${species}" openburning input file..."
+
+filename=${dirin}IAMC-${prefix}-${scen}-1-1/atmos/mon/${species}_em_openburning/gn/v20180628/${species}-em-openburning_input4MIPs_emissions_ScenarioMIP_IAMC-${prefix}-${scen}-1-1_gn_201501-210012.nc
+
+if [ ${species} = "BC"  ] ; then speciesinca="bc"  ; fi
+if [ ${species} = "NOx" ] ; then speciesinca="nox" ; fi
+if [ ${species} = "OC"  ] ; then speciesinca="pom" ; fi
+if [ ${species} = "SO2" ] ; then speciesinca="so2" ; fi
+if [ ${species} = "NH3" ] ; then speciesinca="nh3" ; fi
+
+#--two temporary files for interpolation
+filetmp1=${dirout}tmp1.nc
+filetmp2=${dirout}tmp2.nc
+
 #--output files
 filenameout1=${dirout}flux_${speciesinca}bb_${year}.nc
@@ -652 +424 @@
 filenameout3=${dirout}flux_vector_${speciesinca}bb_${year}.nc
 
-echo ${fileintbb} ${filenameout1} ${filenameout1b} ${filenameout2} ${filenameout3}
-rm -f ${filenameout1} ${filenameout1b} ${filenameout2} ${filenameout3}
+echo ${filename} ${filenameout1} ${filenameout1b} ${filenameout2} ${filenameout3}
+rm -f ${filetmp1} ${filetmp2} ${filenameout1} ${filenameout1b} ${filenameout2} ${filenameout3}
 
 #--unfortunately idl not happy with input netcdf files so need to ferretize files
-#--I extract the correct year as well; no need to sum up sectors, as none in BB files
+#--I also sum over sectors and I extract the correct years as well
+#--first year (start of 5-Y or 10-Y period, for later interpolation)
 rm -f rewrite.jnl
 cat << eod > rewrite.jnl
-set memory/size=100
-use "${fileintbb}"
-set region/t=16-jan-${year}:16-dec-${year}
-save/clobber/file="${filenameout1}" ${species}_EM_OPENBURNING
+use "${filename}"
+set region/t=16-jan-${year1}:16-dec-${year1}
+save/clobber/file="${filetmp1}" ${speciesUp}_EM_OPENBURNING[k=@sum]
 eod
 #--run ferret script
@@ -669 +441 @@
 exit
 eod
-
-#--replace undef with 0 
+#--second year (end of 5-Y or 10-Y period, for later interpolation)
+rm -f rewrite.jnl
+cat << eod > rewrite.jnl
+use "${filename}"
+set region/t=16-jan-${year2}:16-dec-${year2}
+save/clobber/file="${filetmp2}" ${speciesUp}_EM_OPENBURNING[k=@sum]
+eod
+#--run ferret script
+ferret << eod 
+go rewrite.jnl
+exit
+eod
+
+#--performing time interpolation
+cdo add -mulc,${pond1} -selname,${speciesUp}_EM_OPENBURNING ${filetmp1} -mulc,${pond2} -selname,${speciesUp}_EM_OPENBURNING ${filetmp2} ${filenameout1}
+
+#--replace undef with 0
 cdo setmisstoc,0.0 ${filenameout1} ${filenameout1b}
 
@@ -676 +463 @@
 #--OC to POM conversion factor
 #--as ferret returns NOX, treat NOx NOX inconsistency in names by converting to upper case
-if [ ${species} != "OC" ] ; then
+if [ ${species} != "OC"  ] ; then
 echo cdo remapcon,${gridfile} -chname,`echo ${species}_EM_OPENBURNING | awk '{print toupper($0)}'`,flux ${filenameout1b} ${filenameout2}
 cdo remapcon,${gridfile} -chname,`echo ${species}_EM_OPENBURNING | awk '{print toupper($0)}'`,flux ${filenameout1b} ${filenameout2}
@@ -816 +603 @@
 end
 eod
+
 #
 #--calling IDL
@@ -832 +620 @@
 done
 
+#--cleaning up
+rm -f ferret*
+
 #--unfortunately idl use capital letters for variable names so need to change to small letters for now 
 rm -f ${dirout}flux_vector_h2s_${year}.nc ${dirout}flux_vector_so4_${year}.nc
@@ -845 +636 @@
 rm -f ${dirout}sflx_lmdz_cmip6_${year}.nc
 #--merging all files into a single one
-#--anthro NOx is NO2 so 30/46 scaling factor to change to NO 
-#--openburning NOx is NO so no change in unit
+#--PNNL NOx is NO2 so 30/46 scaling factor to change to NO 
+#--VUA NOx is NO so no change in unit
 cdo merge -expr,'flx_bc=FLX_BC' ${dirout}flux_vector_bc_${year}.nc -expr,'flx_bbbc=FLX_BBBC' ${dirout}flux_vector_bcbb_${year}.nc -expr,'flx_pom=FLX_POM' ${dirout}flux_vector_pom_${year}.nc -expr,'flx_bbpom=FLX_BBPOM' ${dirout}flux_vector_pombb_${year}.nc -expr,'flx_nox=30.*FLX_NOX/46.' ${dirout}flux_vector_nox_${year}.nc -expr,'flx_bbnox=FLX_BBNOX;' ${dirout}flux_vector_noxbb_${year}.nc -expr,'flx_so2=FLX_SO2' ${dirout}flux_vector_so2_${year}.nc -expr,'flx_bbso2=FLX_BBSO2' ${dirout}flux_vector_so2bb_${year}.nc -expr,'flx_nh3=FLX_NH3' ${dirout}flux_vector_nh3_${year}.nc -expr,'flx_bbnh3=FLX_BBNH3' ${dirout}flux_vector_nh3bb_${year}.nc  ${dirout}flux_vector_h2s_${year}.nc ${dirout}flux_vector_so4_${year}.nc -selname,conc_dms ${fileINCAex} ${dirout}sflx_lmdz_cmip6_${year}.nc
 
@@ -853 +644 @@
 
 #--cleaning up 
-#rm -f ${dirout}flux*_${year}.nc
-#--to be uncommented in final version
+rm -f ${dirout}flux*_${year}.nc
+#--to be uncommented in final script
 
 #--end loop on years
 done
-
-rm -f ${fileintanthro} ${fileintbb}
 
 #--end loop on scenarios
@@ -867 +656 @@
 rm -f ferret*
 rm -f regrid.pro
-rm -f time_interp_over_gaps.pro sumsectors.jnl select_yr.jnl rewrite.jnl
+rm -f rewrite.jnl