source: TOOLS/CMIP6_FORCING/SCENARIOS/RCP_to_ssp_OXY.sh @ 4004

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

A few scripts to generate oxidant files and SST/SIC files for the scenarios.

  • Property svn:executable set to *
File size: 3.8 KB
Line 
1#--weights on RCPs to yield ssp oxidant files
2#--Th. Lurton, IPSL, 07/2018
3
4# ssp1-19 = 19/26*RCP26
5# ssp1-26 = RCP26
6# ssp2-45 = RCP45
7# ssp3-70 = .4*RCP85 + .6*RCP60
8# ssp4-34 = 8/19*RCP45 + 11/19*RCP26
9# ssp4-60 = .6*RCP60 + .4*RCP45
10# ssp5-34-over = 8/19*RCP45 + 11/19*RCP26
11# ssp5-85 = RCP85
12
13if [ ! -d ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp119 ] ; then mkdir -p ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp119 ; fi
14if [ ! -d ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp126 ] ; then mkdir -p ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp126 ; fi
15if [ ! -d ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp245 ] ; then mkdir -p ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp245 ; fi
16if [ ! -d ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp370 ] ; then mkdir -p ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp370 ; fi
17if [ ! -d ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp434 ] ; then mkdir -p ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp434 ; fi
18if [ ! -d ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp460 ] ; then mkdir -p ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp460 ; fi
19if [ ! -d ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp534-over ] ; then mkdir -p ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp534-over ; fi
20if [ ! -d ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp585 ] ; then mkdir -p ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp585 ; fi
21
22#--creating specific file "RCP0" for RCP26-->ssp119 weighting
23fileinRCP0=/ccc/work/cont003/gencmip6/lurtont/CMIP6_OXIDANTS/LOI_AR6_RCP0.nc
24rm -f ${fileinRCP0}
25ncrcat /ccc/work/cont003/gencmip6/lurtont/CMIP6_OXIDANTS/yearly_interp/LOI_AR6_RCP26_201[05].nc /ccc/work/cont003/gencmip6/lurtont/CMIP6_OXIDANTS/LOI_AR6_RCP26_2010-2015.nc
26cdo ymonmean /ccc/work/cont003/gencmip6/lurtont/CMIP6_OXIDANTS/LOI_AR6_RCP26_2010-2015.nc ${fileinRCP0}
27rm -f /ccc/work/cont003/gencmip6/lurtont/CMIP6_OXIDANTS/LOI_AR6_RCP26_2010-2015.nc
28
29for year in {2015..2100}
30do
31
32echo '... processing year '${year}'...'
33
34fileinRCP26=/ccc/work/cont003/gencmip6/lurtont/CMIP6_OXIDANTS/yearly_interp/LOI_AR6_RCP26_${year}.nc
35fileinRCP45=/ccc/work/cont003/gencmip6/lurtont/CMIP6_OXIDANTS/yearly_interp/LOI_AR6_RCP45_${year}.nc
36fileinRCP60=/ccc/work/cont003/gencmip6/lurtont/CMIP6_OXIDANTS/yearly_interp/LOI_AR6_RCP60_${year}.nc
37fileinRCP85=/ccc/work/cont003/gencmip6/lurtont/CMIP6_OXIDANTS/yearly_interp/LOI_AR6_RCP85sstA2_${year}.nc
38
39cdo add -mulc,19. -divc,26. -selname,H2O2,HNO3,NO2,NO3,O1D,O3,OH ${fileinRCP26} -mulc,7. -divc,26. -selname,H2O2,HNO3,NO2,NO3,O1D,O3,OH ${fileinRCP0} ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp119/oxidants_ssp119_${year}.nc
40cp ${fileinRCP26} ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp126/oxidants_ssp126_${year}.nc
41cp ${fileinRCP45} ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp245/oxidants_ssp245_${year}.nc
42cdo add -mulc,.4 -selname,H2O2,HNO3,NO2,NO3,O1D,O3,OH ${fileinRCP85} -mulc,.6 -selname,H2O2,HNO3,NO2,NO3,O1D,O3,OH ${fileinRCP60} ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp370/oxidants_ssp370_${year}.nc
43cdo add -mulc,8. -divc,19. -selname,H2O2,HNO3,NO2,NO3,O1D,O3,OH ${fileinRCP45} -mulc,11. -divc,19. -selname,H2O2,HNO3,NO2,NO3,O1D,O3,OH ${fileinRCP26} ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp434/oxidants_ssp434_${year}.nc
44cdo add -mulc,.6 -selname,H2O2,HNO3,NO2,NO3,O1D,O3,OH ${fileinRCP60} -mulc,.4 -selname,H2O2,HNO3,NO2,NO3,O1D,O3,OH ${fileinRCP45} ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp460/oxidants_ssp460_${year}.nc
45cdo add -mulc,8. -divc,19. -selname,H2O2,HNO3,NO2,NO3,O1D,O3,OH ${fileinRCP45} -mulc,11. -divc,19. -selname,H2O2,HNO3,NO2,NO3,O1D,O3,OH ${fileinRCP26} ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp534-over/oxidants_ssp534-over_${year}.nc
46cp ${fileinRCP85} ${GENCMIP6_CCCWORKDIR}/CMIP6_OXIDANTS/Scenarios/ssp585/oxidants_ssp585_${year}.nc
47
48#--end loop on years
49done
Note: See TracBrowser for help on using the repository browser.