Changeset 3567
- Timestamp:
- 02/07/18 14:22:26 (6 years ago)
- Location:
- TOOLS/CMIP6_FORCING
- Files:
-
- 14 edited
Legend:
- Unmodified
- Added
- Removed
-
TOOLS/CMIP6_FORCING/AER_TROP_EMISSIONS/README
r3398 r3567 30 30 Update: 02/10/2017 31 31 report NOx emissions as NO for both FF and BB emissions 32 33 Update: 05/02/2018 34 New version of regrid.sh 35 Diffs: * updated paths to access new version of input data; 36 * also some relevant corrections: underscores become hyphens for PNNL species-directories; VUA species-directory change from e.g. "CO-em-biomassburning" to "CO". 37 32 38 33 39 ********************************************************** -
TOOLS/CMIP6_FORCING/AER_TROP_EMISSIONS/REGRID/regrid.sh
r3398 r3567 1 #--updated on 5/2/2018 with new paths, ThL 1 2 #--updated on 4/5/2017 with improved Sheng & Zwiers algorithm, ThL 2 3 #--corrected some interpolation preprocessing (compared to v4) … … 15 16 16 17 #--input directory for anthropogenic (non-BB) emissions 17 dirinPNNL="/prodigfs/ project/input4MIPs/PNNL-JGCRI/emissions/CMIP/CEDS-2017-05-18/mon/"18 dirinPNNL="/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/PNNL-JGCRI/CEDS-2017-05-18/atmos/mon/" 18 19 ##--Be careful if ever one needs the 3D (=AIR) data, the latest are in the different directory 19 ## dirinPNNL_AIR="/prodigfs/project/input4MIPs/PNNL-JGCRI/emissions/CMIP/mon/CEDS-2017-08-30/"20 ## /prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/PNNL-JGCRI/CEDS-2017-08-30/atmos/mon/ 20 21 #--input directory for anthropogenic (BB) emissions 21 dirinVUA="/prodigfs/ project/input4MIPs/VUA/emissions/CMIP/VUA-CMIP-BB4CMIP6-1-2/mon/"22 dirinVUA="/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/VUA/VUA-CMIP-BB4CMIP6-1-2/atmos/mon/" 22 23 23 24 #--LMDz grid information … … 71 72 fi 72 73 73 #--input file PNNL (fossil fuel emissions) 74 filename=${dirinPNNL}${species} _em_anthro/gn/v20170519/${species}-em-anthro_input4MIPs_emissions_CMIP_CEDS-2017-05-18_gn_${year1}01-${year2}12.nc74 #--input file PNNL (fossil fuel emissions) Underscores in the old species-directory are now hyphens. 75 filename=${dirinPNNL}${species}-em-anthro/gn/v20170519/${species}-em-anthro_input4MIPs_emissions_CMIP_CEDS-2017-05-18_gn_${year1}01-${year2}12.nc 75 76 76 77 #--output files … … 266 267 fi 267 268 268 #--now dealing with BB sources 269 filename=${dirinVUA}${species}" -em-biomassburning/gn/v20161213/"${species}"-em-biomassburning_input4MIPs_emissions_CMIP_VUA-CMIP-BB4CMIP6-1-2_gn_${year1}01-${year2}12.nc"269 #--now dealing with BB sources Species-directory was previously of the form "CO-em-biomassburning", it is now just "CO" 270 filename=${dirinVUA}${species}"/gn/v20161213/"${species}"-em-biomassburning_input4MIPs_emissions_CMIP_VUA-CMIP-BB4CMIP6-1-2_gn_${year1}01-${year2}12.nc" 270 271 271 272 #--output files -
TOOLS/CMIP6_FORCING/GHG/README
r3386 r3567 1 Script updated on 5 Feb. 2018 2 Th. Lurton 3 * Paths changed to retrieve new input data versions; 4 * Correction of hyphens/underscores in directory/file names in order to suit the new data version. 5 6 1 7 GHG data for CMIP6 are initially coming from 2 8 ftp://data.iac.ethz.ch/CMIP6/input4MIPs/UoM/GHGConc/CMIP -
TOOLS/CMIP6_FORCING/GHG/extract.sh
r3386 r3567 1 #--Update Th.Lurton 05/02/2018 for paths to new version 2 # ** species-directory take hyphens instead of underscores previously; 3 # ** likewise for species names within filenames; 4 # ** variable names stay with underscores. 5 1 6 #--v3 of my script to extract global mean and put to IGCM format 2 7 … … 11 16 12 17 #--dir for annual mean data 13 dirin1='/prodigfs/project/input4MIPs/UoM/GHGConcentrations/CMIP/UoM-CMIP-1-2-0/yr/' 14 dirin2='gr1-GMNHSH/v20160902' 18 dirin1='/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/UoM/UoM-CMIP-1-2-0/atmos/yr/' 19 dirin2='gr1-GMNHSH/v20160830' 20 #--dirin2 versioning date is different from previously 15 21 16 22 dirout='./OUTPUT' … … 25 31 #--N2O we leave unit to ppb 26 32 GHG='N2O' 27 name='nitrous _oxide'28 file='mole _fraction_of_'${name}'_in_air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc'33 name='nitrous-oxide' 34 file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 29 35 rm -f tmp.txt ${dirout}/${GHG}_CMIP6_0000_2014.txt 30 dirin=${dirin1}'mole _fraction_of_'${name}'_in_air/'${dirin2}36 dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 31 37 cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt 32 38 paste years.txt tmp.txt > ${dirout}/${GHG}_CMIP6_0000_2014.txt … … 35 41 GHG='CH4' 36 42 name='methane' 37 file='mole _fraction_of_'${name}'_in_air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc'38 dirin=${dirin1}'mole _fraction_of_'${name}'_in_air/'${dirin2}43 file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 44 dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 39 45 rm -f tmp.txt ${dirout}/${GHG}_CMIP6_0000_2014.txt 40 46 cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt … … 43 49 #--CO2 here we leave unit in ppm 44 50 GHG='CO2' 45 name='carbon _dioxide'46 file='mole _fraction_of_'${name}'_in_air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc'47 dirin=${dirin1}'mole _fraction_of_'${name}'_in_air/'${dirin2}51 name='carbon-dioxide' 52 file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 53 dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 48 54 rm -f tmp.txt ${dirout}/${GHG}_CMIP6_0000_2014.txt 49 55 cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt … … 54 60 GHG='CFC11eq' 55 61 name='cfc11eq' 56 file='mole _fraction_of_'${name}'_in_air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc'57 dirin=${dirin1}'mole _fraction_of_'${name}'_in_air/'${dirin2}62 file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 63 dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 58 64 rm -f tmp.txt ${dirout}/${GHG}_0000_2014.txt 59 65 cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt … … 63 69 GHG='CFC12eq' 64 70 name='cfc12eq' 65 file='mole _fraction_of_'${name}'_in_air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc'66 dirin=${dirin1}'mole _fraction_of_'${name}'_in_air/'${dirin2}71 file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 72 dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 67 73 rm -f tmp.txt ${dirout}/${GHG}_0000_2014.txt ${dirout}/${GHG}_CMIP6_0000_2014.txt 68 74 cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt … … 75 81 GHG='CFC11' 76 82 name='cfc11' 77 file='mole _fraction_of_'${name}'_in_air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc'78 dirin=${dirin1}'mole _fraction_of_'${name}'_in_air/'${dirin2}83 file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 84 dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 79 85 rm -f tmp.txt ${dirout}/${GHG}_0000_2014.txt 80 86 cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt … … 85 91 GHG='CFC12' 86 92 name='cfc12' 87 file='mole _fraction_of_'${name}'_in_air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc'88 dirin=${dirin1}'mole _fraction_of_'${name}'_in_air/'${dirin2}93 file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 94 dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 89 95 rm -f tmp.txt ${dirout}/${GHG}_0000_2014.txt 90 96 cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt … … 94 100 GHG='HFC134aeq' 95 101 name='hfc134aeq' 96 file='mole _fraction_of_'${name}'_in_air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc'97 dirin=${dirin1}'mole _fraction_of_'${name}'_in_air/'${dirin2}102 file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 103 dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 98 104 rm -f tmp.txt ${dirout}/${GHG}_0000_2014.txt ${dirout}/CFC11eq_CMIP6_0000_2014.txt 99 105 cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt -
TOOLS/CMIP6_FORCING/NITROGEN/README
r3388 r3567 1 Update 05 Feb. 2018 2 Th. Lurton 3 ** Changed paths to new version 4 5 1 6 Here we only add drynhx+wetnhx and drynoy+wetnoy into nhx and noy deposition fluxes (kg N/m2/s). 2 7 There is a 1850 climatological file (10-year average) and yearly files for 1850-2014. -
TOOLS/CMIP6_FORCING/NITROGEN/nitrogen.sh
r3387 r3567 1 dirin1="/prodigfs/ project/input4MIPs/UReading/surfaceFluxes/CMIP/NCAR-CCMI-2-0/mon"1 dirin1="/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/NCAR/NCAR-CCMI-2-0/atmos/mon" 2 2 dirin2="gn/v20161207" 3 dirinC1="/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/NCAR/NCAR-CCMI-2-0/atmos/monC" 4 3 5 dirout="./OUTPUT" 4 6 if [ ! -d ${dirout} ] ; then mkdir ${dirout} ; fi … … 6 8 #--1850 climatological N deposition 7 9 # 8 drynhx=${dirin 1}"/drynhx/"${dirin2}"/drynhx_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012-clim.nc"9 wetnhx=${dirin 1}"/wetnhx/"${dirin2}"/wetnhx_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012-clim.nc"10 drynoy=${dirin 1}"/drynoy/"${dirin2}"/drynoy_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012-clim.nc"11 wetnoy=${dirin 1}"/wetnoy/"${dirin2}"/wetnoy_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012-clim.nc"10 drynhx=${dirinC1}"/drynhx/"${dirin2}"/drynhx_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012.nc" 11 wetnhx=${dirinC1}"/wetnhx/"${dirin2}"/wetnhx_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012.nc" 12 drynoy=${dirinC1}"/drynoy/"${dirin2}"/drynoy_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012.nc" 13 wetnoy=${dirinC1}"/wetnoy/"${dirin2}"/wetnoy_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012.nc" 12 14 # 13 cdo merge $ drynhx $wetnhx $dirout/CMIP6_fluxnhx_merged_clim_1850.nc14 cdo merge $ drynoy $wetnoy $dirout/CMIP6_fluxnoy_merged_clim_1850.nc15 cdo expr,'fluxnhx=drynhx+wetnhx' $ dirout/CMIP6_fluxnhx_merged_clim_1850.nc $dirout/CMIP6_fluxnhx_clim_1850.nc16 cdo expr,'fluxnoy=drynoy+wetnoy' $ dirout/CMIP6_fluxnoy_merged_clim_1850.nc $dirout/CMIP6_fluxnoy_clim_1850.nc17 rm -f $ dirout/CMIP6_fluxnhx_merged_clim_1850.nc $dirout/CMIP6_fluxnoy_merged_clim_1850.nc15 cdo merge ${drynhx} ${wetnhx} ${dirout}/CMIP6_fluxnhx_merged_clim_1850.nc 16 cdo merge ${drynoy} ${wetnoy} ${dirout}/CMIP6_fluxnoy_merged_clim_1850.nc 17 cdo expr,'fluxnhx=drynhx+wetnhx' ${dirout}/CMIP6_fluxnhx_merged_clim_1850.nc ${dirout}/CMIP6_fluxnhx_clim_1850.nc 18 cdo expr,'fluxnoy=drynoy+wetnoy' ${dirout}/CMIP6_fluxnoy_merged_clim_1850.nc ${dirout}/CMIP6_fluxnoy_clim_1850.nc 19 rm -f ${dirout}/CMIP6_fluxnhx_merged_clim_1850.nc ${dirout}/CMIP6_fluxnoy_merged_clim_1850.nc 18 20 # 19 21 # … … 25 27 wetnoy=${dirin1}"/wetnoy/"${dirin2}"/wetnoy_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-201412.nc" 26 28 # 27 cdo merge $ drynhx $wetnhx $dirout/CMIP6_fluxnhx_merged_annual_1850-2014.nc28 cdo merge $ drynoy $wetnoy $dirout/CMIP6_fluxnoy_merged_annual_1850-2014.nc29 cdo splityear -expr,'fluxnhx=drynhx+wetnhx' $ dirout/CMIP6_fluxnhx_merged_annual_1850-2014.nc $dirout/CMIP6_fluxnhx_annual_30 cdo splityear -expr,'fluxnoy=drynoy+wetnoy' $ dirout/CMIP6_fluxnoy_merged_annual_1850-2014.nc $dirout/CMIP6_fluxnoy_annual_31 rm -f $ dirout/CMIP6_fluxnhx_merged_annual_1850-2014.nc $dirout/CMIP6_fluxnoy_merged_annual_1850-2014.nc29 cdo merge ${drynhx} ${wetnhx} ${dirout}/CMIP6_fluxnhx_merged_annual_1850-2014.nc 30 cdo merge ${drynoy} ${wetnoy} ${dirout}/CMIP6_fluxnoy_merged_annual_1850-2014.nc 31 cdo splityear -expr,'fluxnhx=drynhx+wetnhx' ${dirout}/CMIP6_fluxnhx_merged_annual_1850-2014.nc ${dirout}/CMIP6_fluxnhx_annual_1850-2014.nc 32 cdo splityear -expr,'fluxnoy=drynoy+wetnoy' ${dirout}/CMIP6_fluxnoy_merged_annual_1850-2014.nc ${dirout}/CMIP6_fluxnoy_annual_1850-2014.nc 33 rm -f ${dirout}/CMIP6_fluxnhx_merged_annual_1850-2014.nc ${dirout}/CMIP6_fluxnoy_merged_annual_1850-2014.nc -
TOOLS/CMIP6_FORCING/OZONE/clims_CMIP6.bash
r3389 r3567 1 1 #!/bin/bash 2 3 ia="3" #=== Indices of chosen Amip datasets (see below AMIP_VERSION_NAMES/AMIP_VERSION_FOLDERS) 4 io="1" #=== Indices of chosen Ozone datasets (see below OZON_VERSION_NAMES/OZON_VERSION_FOLDERS) 5 recomp='n' #=== Recompute files that are already present (y/n) 6 do_amip='n' #=== Compute sic/sst files 7 do_ozon='y' #=== Compute ozone files 8 2 module load nco/4.4.8 3 4 ia="1" #=== Indices of chosen Amip datasets (see below AMIP_VERSION_NAMES/AMIP_VERSION_FOLDERS) 5 io="1" #=== Indices of chosen Ozone datasets (see below OZON_VERSION_NAMES/OZON_VERSION_FOLDERS) 9 6 10 7 #=============================================================================== … … 13 10 if [ ${HOSTNAME:0:3} = 'ada' ]; then machine='ada'; work=$WORKDIR 14 11 elif [ ${HOSTNAME:0:5} = 'curie' ]; then machine='curie'; work=$CCCWORKDIR 15 elif [ ${HOSTNAME:0:6} = 'ciclad' ]; then machine='ciclad'; work=/ data/$USER12 elif [ ${HOSTNAME:0:6} = 'ciclad' ]; then machine='ciclad'; work=/home/$USER/tmp 16 13 else echo "Not set up for this machine yet, sorry."; exit; fi 17 14 export machine=$machine … … 21 18 #=============================================================================== 22 19 local=$PWD 23 echo $machine24 20 case $machine in 25 ada) work=$WORKDIR 26 DATAIN="/workgpfs/rech/dzt/"$USER"/input4MIPs" 27 DATAOU="/workgpfs/rech/dzt/"$USER"/IGCM" 28 module load netcdf 29 module load nco/4.4.6 ;; 21 ada) work=$WORKDIR; DATAIN="" 22 DATAOU='/workgpfs/rech/psl/rpsl035/IGCM' ;; 30 23 curie) work=$CCCWORKDIR; DATAIN="" 31 DATAIN="/ccc/work/rech/psl/"$USER"/input4MIPs" 32 DATAOU="/ccc/work/rech/psl/"$USER"/IGCM" 33 module load netcdf 34 module load nco ;; 35 ciclad) work=/data/$USER 36 DATAIN='/prodigfs/project/input4MIPs' 37 DATAOU='/prodigfs/ipslfs/igcmg/IGCM' 38 DATAOU='/data/'${USER}'/CMIP6/OZONE' 39 module load nco/4.4.8 ;; 24 DATAOU='/ccc/work/rech/psl/rpsl035/IGCM' ;; 25 ciclad) work=/data/$USER; DATAIN=/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP 26 DATAOU='/prodigfs/ipslfs/igcmg/IGCM' ;; 40 27 esac 41 42 if [ ! -d ${DATAOU} ] ; then mkdir -p ${DATAOU} ; fi 28 DATAOU=/home/tlurton/IGCM ##### !!!!!! PROVISOIRE !!!!!!! 43 29 44 30 #--- INPUT DATA: AMIP 45 AMIP_VERSION_NAMES=('PCMDI-AMIP-1-1-0' 'PCMDI-AMIP-1-1-1' 'PCMDI-AMIP-1-1-2') 46 AMIP_VERSION_FOLDERS=('gs1x1/v20160609' 'gn/v20161020' 'gn/v20170419') 47 AMIP_FOLDER_IN=$DATAIN/PCMDI/SSTsAndSeaIce/CMIP/\${amip_version_name}/mon/\${V}/\${amip_version_folder} 48 amip_vname_in=('tosbcs' 'siconcbcs') 49 amip_fname_in='${V}_input4MIPs_SSTsAndSeaIce_CMIP_${amip_version_name}_${amip_version_folder%/*}*_18????-20????.nc' 31 AMIP_VERSION_NAMES=('PCMDI-AMIP-1-1-3') 32 AMIP_VERSION_FOLDERS=('gn/v20161031') 33 AMIP_FOLDER_IN=$DATAIN'/PCMDI/${amip_version_name}/${W}/mon/${V}/${amip_version_folder}' 34 amip_tname_in=('ocean' 'seaIce') 35 amip_vname_in=('tos' 'siconc') 36 amip_fname_in='${V}_input4MIPs_SSTsAndSeaIce_CMIP_${amip_version_name}_${amip_version_folder%/*}_187001-20????.nc' 50 37 51 38 #--- INPUT DATA: OZONE 52 39 OZON_VERSION_NAMES=('UReading-CCMI-1-0') 53 OZON_VERSION_FOLDERS=('gr/v20160711'); if [ $machine = 'ciclad' ]; then OZON_VERSION_FOLDERS=('gr/v20160907'); fi 54 OZON_FOLDER_IN=$DATAIN/UReading/ozone/CMIP/\${ozon_version_name}/mon/\${V}/\${ozon_version_folder} 40 OZON_VERSION_FOLDERS=('gn/v20160711') 41 OZON_VERSION_FOLDERS_CLIM=('gn/v20160830') 42 OZON_FOLDER_IN=$DATAIN/'UReading/${ozon_version_name}/atmos/mon${clim}/${V}/${ozon_version_folder}' 55 43 ozon_vname_in='vmro3' 56 ozon_fname_in='${V}_input4MIPs_ozone_CMIP_${ozon_version_name}_${ozon_version_folder%/*} *_${Yi}01-${Yf}12.nc'44 ozon_fname_in='${V}_input4MIPs_ozone_CMIP_${ozon_version_name}_${ozon_version_folder%/*}_${Yi}01-${Yf}12.nc' 57 45 58 46 #--- VERSIONS OF INPUT FILES EFFECTIVELY CHOSEN 59 47 ((ia--)); ((io--)) 60 amip_version_name=${AMIP_VERSION_NAMES[$ia]}; amip_version_folder=${AMIP_VERSION_FOLDERS[$ia]}61 ozon_version_name=${OZON_VERSION_NAMES[$io]}; ozon_version_folder=${OZON_VERSION_FOLDERS[$io]}62 48 63 49 #--- OUTPUT DATA 64 AMIP_FOLDER_OU=$DATAOU/ATM/LIMIT/AMIP.\${amip_version_folder##*/}/original/\${nlon}x\${nlat}/BCS 65 OZON_FOLDER_OU=$DATAOU/ATM/OZONE/UReading.\${ozon_version_folder##*/}/original 66 amip_vname_ou=('tosbcs' 'sicbcs'); amip_fname_ou='amipbc_${V}_${nlon}x${nlat}_${Y}.nc' 50 DATAOU=${DATAOU%/*}/ANTICHAMBRE; if [ ! -d $DATAOU ]; then mkdir -p $DATAOU; fi 51 AMIP_FOLDER_OU=$DATAOU/'ATM/LIMIT/AMIP.${amip_version_folder##*/}/${suf}' 52 OZON_FOLDER_OU=$DATAOU/'ATM/OZONE/UReading.${ozon_version_folder##*/}' 53 amip_vname_ou=('tosbcs' 'sicbcs'); amip_fname_ou='amip${bc}_${V}_360x180_${Y}.nc' 67 54 ozon_vname_ou='tro3' ; ozon_fname_ou='${V}_${Y}.nc' 68 Y_amip=1870-2016 #--- FOR AMIP FILES (INTERANNUAL) 69 Y_ozon=1850-2014:50 #--- FOR OZONE FILES (INTERANNUAL) 70 Y_clim=1979-2008 #--- FOR CLIMATOLOGIES (COMMON PERIODS FOR O3 AND SIC/SST) 71 Y_pind=1870-1899 #--- FOR SST/SIC PRE-INDUSTRIAL CLIMATOLOGY 55 Y_amip_beg=1870 ; Y_amip_end=2017 #--- FOR AMIP FILES (INTERANNUAL) 56 Y_ozon_beg=1850 ; Y_ozon_end=2014 #--- FOR OZONE FILES (INTERANNUAL) 57 Y_clim_beg=1979 ; Y_clim_end=2008 #--- FOR CLIMATOLOGIES (COMMON PERIOD FOR O3 AND SIC/SST) 58 Y_pind_beg=1870 ; Y_pind_end=1899 #--- FOR SST/SIC PRE-INDUSTRIAL CLIMATOLOGY 59 # NB: CLIMATOLOGIES ARE PRODUCED FOR EACH FILE CONTAINING "-clim" SUFFIX FOUND IN THE FOLDER TREE, IN PARTICULAR FOR OZONE. 60 61 #--- RECOMPUTE FILES THAT ARE ALREADY PRESENT (y/n) 62 recomp='n' 63 do_amip='n' 64 do_ozon='y' 72 65 73 66 if [ "$do_amip" = "y" ]; then 67 amip_version_name=${AMIP_VERSION_NAMES[$ia]}; amip_version_folder=${AMIP_VERSION_FOLDERS[$ia]} 74 68 echo ">> Chosen SST files: $(V=${amip_vname_in[0]}; amip_version_name=$amip_version_name; amip_version_folder=$amip_version_folder; eval echo $AMIP_FOLDER_IN)" 75 69 echo ">> Chosen SIC files: $(V=${amip_vname_in[1]}; amip_version_name=$amip_version_name; amip_version_folder=$amip_version_folder; eval echo $AMIP_FOLDER_IN)" 76 70 fi 77 71 if [ "$do_ozon" = "y" ]; then 72 ozon_version_name=${OZON_VERSION_NAMES[$io]}; ozon_version_folder=${OZON_VERSION_FOLDERS[$io]} 78 73 echo ">> Chosen OZONE files: $(V=${ozon_vname_in}; ozon_version_name=$ozon_version_name; ozon_version_folder=$ozon_version_folder; eval echo $OZON_FOLDER_IN)" 79 74 fi … … 231 226 232 227 #================================================================================== 233 function tim_len {234 local n235 n=`ncdump -h $1 | grep 'UNLIMITED'`; n=${n#*// (}; n=${n% currently*}; echo $n236 }237 238 #==================================================================================239 function dim_len {240 local n241 if [ "$2" = "T" ]; then tim_len $1; else n=`ncdump -h $1 | grep $(dim_nam $1 $2)' = '`; n=${n#* = }; n=${n% ;*}; echo $n; fi242 }243 244 #==================================================================================245 228 function test_var { 246 229 # Input arguments: … … 257 240 # Rename dimensions, dimensional variables and possibly bounds_<variable>. 258 241 # 259 local f=$1 out="" nam_in nam_ou vOK 260 vOK=0; if [ $f = "v" ]; then vOK=1; shift; f=$1; fi 242 local f=$1 out="" nam_in nam_ou 261 243 shift 262 244 while [ $# -ge 1 ]; do 263 245 nam_ou=${1##*:}; nam_in=$(dim_nam $f ${1%:*}); shift 264 246 if [ "$nam_in" != "$nam_ou" ]; then 265 out="$out -d $nam_in,$nam_ou" 266 if [ $vOK -eq 1 ]; then out="$out -v $nam_in,$nam_ou"; fi 247 out="$out -d $nam_in,$nam_ou -v $nam_in,$nam_ou" 267 248 nam_in=bounds_$nam_in; nam_ou=bounds_$nam_ou 268 249 test_var $f $nam_in 269 if [ $? -eq 0 ]; then out="$out - v $nam_in,$nam_ou"; fi250 if [ $? -eq 0 ]; then out="$out -d $nam_in,$nam_ou -v $nam_in,$nam_ou"; fi 270 251 fi 271 252 ((ia++)) … … 282 263 283 264 #=== BUILD YEARLY INTERANNUAL SIC/SST FILES WITH 12 RECORDS 284 for is in $ia; do #--- Loop on Amip versions265 for is in $ia; do #--- Loop on Amip versions 285 266 amip_version_name=${AMIP_VERSION_NAMES[$is]} 286 267 amip_version_folder=${AMIP_VERSION_FOLDERS[$is]} 287 for ivar in $(eval echo {1..${#amip_vname_in[@]}}); do iv=$((ivar-1)) #--- Loop on variables 288 289 #--- VARIABLES NAMES (INPUT, OUTPUT AND FOR FILE NAMING PURPOSE) 290 v_in=${amip_vname_in[$iv]}; v_ou=${amip_vname_ou[$iv]} 291 var='sst'; if [ $ivar -eq 2 ]; then var='sic'; fi 292 293 #--- INPUT FOLDER AND FILE NAMES 294 d_in=$(amip_version_name=$amip_version_name; amip_version_folder=$amip_version_folder; V=$v_in; eval echo $AMIP_FOLDER_IN) 295 f_in=$d_in/$(V=$v_in; eval echo ${amip_fname_in}) 296 297 #--- OUTPUT FILE AND FOLDER NAMES 298 nlon=$(dim_len $f_in X); nlat=$(dim_len $f_in Y) #--- Dimensions lengths 299 d_ou=$(amip_version_name=$amip_version_name; amip_version_folder=$amip_version_folder; nlon=$nlon; nlat=$nlat; eval echo $AMIP_FOLDER_OU) 300 if [ ! -d $d_ou ]; then mkdir -p $d_ou; fi 301 302 #--- FILES NAMES (INPUT, OUTPUT) 303 f_ou=$d_ou/$(V=$var; Y=\$Y; nlon=$nlon; nlat=$nlat; eval echo ${amip_fname_ou}) 304 305 for per in ${Y_amip}; do 306 Yb=${per%%-*}; Ye=${per##*-} 307 echo ">> BUILDING 12 MONTHS $v_in FILES FOR THE ${per} PERIOD USING $amip_version_name DATASET (${amip_version_folder##*/})..." 308 for Y in $(eval echo {${Yb}..${Ye}); do #--- Loop on years 268 for suf in AVE BCS; do #--- Loop on var type 269 for ivar in $(eval echo {1..${#amip_vname_in[@]}}); do iv=$((ivar-1)) #--- Loop on variables 270 271 #--- VARIABLES NAMES (INPUT, OUTPUT AND FOR FILE NAMING PURPOSE) 272 v_in=${amip_vname_in[$iv]}; if [ "$suf" = "BCS" ]; then v_in=${v_in}bcs; fi 273 v_ou=${amip_vname_ou[$iv]}; if [ "$suf" = "BCS" ]; then v_ou=${v_ou}bcs; fi 274 var='sst'; if [ $ivar -eq 2 ]; then var='sic'; fi 275 276 #--- INPUT/OUTPUT FOLDER NAME 277 d_in=$(amip_version_name=$amip_version_name; amip_version_folder=$amip_version_folder; suf=$suf; V=$v_in; eval echo $AMIP_FOLDER_IN) 278 d_ou=$(amip_version_name=$amip_version_name; amip_version_folder=$amip_version_folder; suf=$suf; eval echo $AMIP_FOLDER_OU) 279 if [ ! -d $d_ou ]; then mkdir -p $d_ou; fi 280 281 #--- FILES NAMES (INPUT, OUTPUT) 282 bc=""; if [ "$suf" = "BCS" ]; then bc=bc; fi 283 f_in=$d_in/$(V=$v_in; eval echo ${amip_fname_in}) 284 f_ou=$d_ou/$(V=$var; Y=\$Y; bc=$bc; eval echo ${amip_fname_ou}) 285 286 echo ">> BUILDING 12 MONTHS $v_in FILES FOR THE ${Y_amip_beg}-${Y_amip_end} PERIOD USING $amip_version_name DATASET (${amip_version_folder##*/})..." 287 for Y in $(eval echo {${Y_amip_beg}..${Y_amip_end}}); do #--- Loop on years 309 288 310 289 #--- RECOMPUTE OR NOT … … 312 291 313 292 #--- 2D FIELDS EXTRACTION 314 extract $v_in,$v_ou $AMIP_FOLDER_IN/$amip_fname_in,$f_ou $Y b,$Ye,0 $Y 0315 316 progress_bar $((Y-Y b+1)) $((Ye-Yb+1)) 50317 done318 done293 extract $v_in,$v_ou $AMIP_FOLDER_IN/$amip_fname_in,$f_ou $Y_amip_beg,$Y_amip_end,0 $Y 0 294 295 progress_bar $((Y-Y_amip_beg+1)) $((Y_amip_end-Y_amip_beg+1)) 50 296 done 297 done 319 298 done 320 299 done 321 300 322 301 #=== BUILD SIC/SST 12 RECORDS CLIMATOLOGIES 323 for is in $ia; do #--- Loop on Amip versions302 for is in $ia; do #--- Loop on Amip versions 324 303 amip_version_name=${AMIP_VERSION_NAMES[$is]} 325 304 amip_version_folder=${AMIP_VERSION_FOLDERS[$is]} 326 for ivar in $(eval echo {1..${#amip_vname_in[@]}}); do iv=$((ivar-1)) #--- Loop on variables 327 328 #--- VARIABLES (OUTPUT, FOR FILE NAMING PURPOSE) AND OUTPUT FILE NAMES 329 v_ou=${amip_vname_ou[$iv]} 330 var='sst'; if [ $ivar -eq 2 ]; then var='sic'; fi 331 d_ou=$(amip_version_name=$amip_version_name; amip_version_folder=$amip_version_folder; eval echo $AMIP_FOLDER_OU) 332 f_ou=$(V=$var; Y=\$Y; eval echo ${amip_fname_ou}) 333 334 for per in ${Y_clim} ${Y_pind}; do #--- Loop on periods 335 Yb=${per%%-*}; Ye=${per##*-} 336 echo ">> BUILDING 12 MONTHS $v_ou CLIMATOLOGY ON PERIOD ${per} USING $amip_version_name DATASET (${amip_version_folder##*/})..." 337 338 #--- RECOMPUTE OR NOT 339 if [[ -f $d_ou/$(Y=${Yb}_${Ye}_clim; eval echo $f_ou) && $recomp = 'n' ]]; then continue; fi 340 341 #--- COMPUTE CLIMATOLOGY 342 make_clim $v_ou $d_ou/$f_ou ${Yb},${Ye} 0 343 305 for suf in AVE BCS; do #--- Loop on var type 306 for ivar in $(eval echo {1..${#amip_vname_in[@]}}); do iv=$((ivar-1)) #--- Loop on variables 307 308 #--- VARIABLES (OUTPUT, FOR FILE NAMING PURPOSE) AND OUTPUT FILE NAMES 309 bc=""; if [ "$suf" = "BCS" ]; then bc=bc; fi 310 v_ou=${amip_vname_ou[$iv]}; if [ "$suf" = "BCS" ]; then v_ou=${v_ou}bcs; fi 311 var='sst'; if [ $ivar -eq 2 ]; then var='sic'; fi 312 d_ou=$(amip_version_name=$amip_version_name; amip_version_folder=$amip_version_folder; suf=$suf; eval echo $AMIP_FOLDER_OU) 313 f_ou=$(V=$var; Y=\$Y; bc=$bc; eval echo ${amip_fname_ou}) 314 315 for per in ${Y_clim_beg},${Y_clim_end} ${Y_pind_beg},${Y_pind_end}; do #--- Loop on periods 316 echo ">> BUILDING 12 MONTHS $v_ou CLIMATOLOGIES ON PERIOD ${per%%,*}-${per##*,} for $suf DATA..." 317 318 #--- RECOMPUTE OR NOT 319 Yint=${per%%,*}_${per##*,}_clim 320 if [[ -f $d_ou/$(Y=$Yint; eval echo $f_ou) && $recomp = 'n' ]]; then continue; fi 321 322 #--- COMPUTE CLIMATOLOGY 323 make_clim $v_ou $d_ou/$f_ou ${per} 0 324 325 done 326 cd $d_ou; ln -sf $(Y=2002; eval echo $f_ou) $(Y=2002_2002_clim; eval echo $f_ou); cd - 2>&1 > /dev/null 344 327 done 345 cd $d_ou; ln -sf $(Y=2002; eval echo $f_ou) $(Y=2002_2002_clim; eval echo $f_ou); cd - 2>&1 > /dev/null346 328 done 347 329 done … … 358 340 359 341 #=== BUILD YEARLY INTERANNUAL OZONE FILES WITH 12 RECORDS 360 Yb=${Y_ozon%%-*}; Ye=${Y_ozon##*-}; dY=${Ye#*:}; Ye=${Ye%:*}361 342 for is in $io; do #--- Loop on Ozone versions 362 343 ozon_version_name=${OZON_VERSION_NAMES[$is]} 363 344 ozon_version_folder=${OZON_VERSION_FOLDERS[$is]} 364 365 #--- VARIABLES NAMES (INPUT, OUTPUT) 345 d_ou=$(ozon_version_name=$ozon_version_name; ozon_version_folder=$ozon_version_folder; eval echo $OZON_FOLDER_OU) #--- Output folder 346 if [ ! -d $d_ou ]; then mkdir -p $d_ou; fi 347 v_ou=${ozon_vname_ou} #--- Output variable 366 348 v_in=${ozon_vname_in} #--- Input variable 367 v_ou=${ozon_vname_ou} #--- Output variable368 369 #--- INPUT/OUTPUT FILE NAME370 349 f_in=$(V=$v_in; Yi=\${Yi}; Yf=\${Yf}; ozon_version_name=$ozon_version_name; ozon_version_folder=$ozon_version_folder; eval echo $OZON_FOLDER_IN/${ozon_fname_in[$is]}) 371 350 f_ou=$(V=$v_ou; Y=\$Y; eval echo $ozon_fname_ou) 372 351 373 #--- INPUT FILE HORIZONTAL RESOLUTION 374 f0=$(Yi=${Yb}; Yf=$((Yb+dY-1)); eval echo $f_in) #--- First file name 375 nlon=$(dim_len $f0 X); nlat=$(dim_len $f0 Y) #--- Dimensions lengths 376 377 #--- OUTPUT FOLDERS NAMES 378 d_ou=$(ozon_version_name=$ozon_version_name; ozon_version_folder=$ozon_version_folder; eval echo $OZON_FOLDER_OU) 379 d_o2=$d_ou/1x${nlat} #--- Output folder (2D) 380 d_o3=$d_ou/${nlon}x${nlat} #--- Output folder (3D) 352 nm=12 353 d_o2=$d_ou/2D_12Fields #--- Output folder (2D) 354 d_o3=$d_ou/3D_12Fields #--- Output folder (3D) 381 355 if [ ! -d $d_o2 ]; then mkdir -p $d_o2; fi; f_o2=$d_o2/$f_ou 382 356 if [ ! -d $d_o3 ]; then mkdir -p $d_o3; fi; f_o3=$d_o3/$f_ou 383 if [ ! -d $d_ou ]; then mkdir -p $d_ou; fi 384 385 for Y in $(eval echo {${Y b}..${Ye}}); do #--- Loop on years357 358 echo ">> BUILDING $nm MONTHS OZONE FILES FOR PERIOD ${Y_ozon_beg}-${Y_ozon_end} USING $OZON_VERSION_NAMES DATASET (${ozon_version_folder##*/})..." 359 for Y in $(eval echo {${Y_ozon_beg}..${Y_ozon_end}}); do #--- Loop on years 386 360 387 361 #--- RECOMPUTE OR NOT … … 389 363 390 364 #--- 3D FIELDS EXTRACTION 391 extract $v_in,$v_ou $f_in,$f_o3 $Yb,$Ye,$dY $Y 0 392 #$(((nm-12)/2)) 365 extract $v_in,$v_ou $f_in,$f_o3 $Y_ozon_beg,$Y_ozon_end,50 $Y 0 393 366 394 367 #--- RENAME DIMENSIONS ACCORDING TO $coords VARIABLE 395 368 fo3=$(Y=$Y; eval echo $f_o3); fo2=$(Y=$Y; eval echo $f_o2) 396 if [ $Y -eq $Y b]; then args_ren=$(renam_dims $fo3 $coords); fi369 if [ $Y -eq $Y_ozon_beg ]; then args_ren=$(renam_dims $fo3 $coords); fi 397 370 if [ "$args_ren" != "" ]; then ncrename $args_ren $fo3 2>&1 > /dev/null ; fi 398 371 … … 405 378 if [ $? -eq 0 ]; then ncks -x -v $v $fo2 $tmp ; mv $tmp $fo2; fi 406 379 407 progress_bar $((Y-Y b+1)) $((Ye-Yb+1)) $dY380 progress_bar $((Y-Y_ozon_beg+1)) $((Y_ozon_end-Y_ozon_beg+1)) 50 408 381 409 382 done … … 416 389 v_ou=${ozon_vname_ou} #--- Output variable 417 390 418 for per in ${Y_clim} ${Y_pind}; do #--- Loop on periods 419 Yb=${Y_clim%%-*}; Ye=${Y_clim##*-} 420 for rank in 3 2; do #--- Loop on rank 421 if [ $rank -eq 2 ]; then d_ou=$d_o2; else d_ou=$d_o3; fi 422 423 #--- CLIMATOLOGY ON GIVEN YEARS INTERVAL (Y_clim) 424 f_ou=$(V=$v_ou; Y=${Yb}_${Ye}_clim; eval echo $ozon_fname_ou) #--- Regular climatology name 425 if [[ ! -f $d_ou/$f_ou || $recomp = 'y' ]]; then 426 echo ">> BUILDING 12 MONTHS ${rank}D $v_ou CLIMATOLOGY ON PERIOD ${Y_clim}..." 427 make_clim $v_ou $d_ou/\$ozon_fname_ou ${Yb},${Ye} 0 391 for rank in 3 2; do #--- Loop on rank 392 d12=$d_ou/${rank}D_12Fields 393 394 #--- CLIMATOLOGY ON GIVEN YEARS INTERVAL (Y_clim_beg-Y_clim_end) 395 f_ou=$(V=$v_ou; Y=${Y_clim_beg}_${Y_clim_end}_clim; eval echo $ozon_fname_ou) #--- Regular climatology name 396 if [[ ! -f $d12/$f_ou || $recomp = 'y' ]]; then 397 echo ">> BUILDING 12 MONTHS ${rank}D $v_ou CLIMATOLOGY ON PERIOD ${Y_clim_beg}-${Y_clim_end}..." 398 make_clim $v_ou $d12/\$ozon_fname_ou ${Y_clim_beg},${Y_clim_end} 0 399 fi 400 401 #--- PRE-INDUSTRIAL CLIMATOLOGY USING SPECIAL FILE 402 ozon_version_folder=$(eval echo ${OZON_VERSION_FOLDERS_CLIM[$is]}) 403 404 f_in=$(V=$v_in; Yi=1850; Yf=1850; ozon_version_name=$ozon_version_name; ozon_version_folder=$ozon_version_folder; clim='C'; eval echo $OZON_FOLDER_IN/$ozon_fname_in); #f_in=${f_in%.nc*}-clim.nc 405 f_ou=$(V=$v_ou; Y=1850_1850_clim; eval echo $ozon_fname_ou) 406 if [[ ! -f $d12/$f_ou || $recomp = 'y' ]]; then 407 echo ">> BUILDING 12 MONTHS ${rank}D $v_ou PRE-INDUSTRIAL CLIMATOLOGY..." 408 if [ $rank -eq 3 ]; then f_ou0=$d12/$f_ou 409 410 #--- DUPLICATE FILE AND RENAME DIMENSIONS 411 cp $f_in $d12/$f_ou ; args_ren=$(renam_dims $d12/$f_ou $coords) 412 ncrename $args_ren -v $ozon_vname_in,$ozon_vname_ou $d12/$f_ou 2>&1 > /dev/null 413 414 else 415 #--- ZONAL MEAN 416 ncwa -a longitude $f_ou0 $tmp ; mv $tmp $d12/$f_ou 417 ncks -x -v longitude $d12/$f_ou $tmp ; mv $tmp $d12/$f_ou 418 419 #--- REMOVE USELESS VARIABLE 420 v="bounds_longitude" 421 test_var $d12/$f_ou $v 422 if [ $? -eq 0 ]; then ncks -x -v $v $d12/$f_ou $tmp ; mv $tmp $d12/$f_ou; fi 428 423 fi 429 430 #--- PRE-INDUSTRIAL CLIMATOLOGY USING SPECIAL FILE 431 #set -vx 432 f_in=$(V=$v_in; Yi=1850; Yf=1850; eval echo $OZON_FOLDER_IN/$ozon_fname_in); f_in=${f_in%.nc*}-clim.nc 433 f_ou=$(V=$v_ou; Y=1850_1850_clim; eval echo $ozon_fname_ou) 434 if [[ ! -f $d_ou/$f_ou || $recomp = 'y' ]]; then 435 echo ">> BUILDING 12 MONTHS ${rank}D $v_ou PRE-INDUSTRIAL CLIMATOLOGY..." 436 if [ $rank -eq 3 ]; then f_ou0=$d_ou/$f_ou 437 438 #--- DUPLICATE FILE AND RENAME DIMENSIONS 439 cp $f_in $d_ou/$f_ou ; args_ren=$(renam_dims v $d_ou/$f_ou $coords) 440 ncrename $args_ren -v $ozon_vname_in,$ozon_vname_ou $d_ou/$f_ou 2>&1 > /dev/null 441 442 else 443 #--- ZONAL MEAN 444 ncwa -a longitude $f_ou0 $tmp ; mv $tmp $d_ou/$f_ou 445 ncks -x -v longitude $d_ou/$f_ou $tmp ; mv $tmp $d_ou/$f_ou 446 447 #--- REMOVE USELESS VARIABLE 448 v="bounds_longitude" 449 test_var $d_ou/$f_ou $v 450 if [ $? -eq 0 ]; then ncks -x -v $v $d_ou/$f_ou $tmp ; mv $tmp $d_ou/$f_ou; fi 451 fi 452 fi 453 #set +vx 454 done 424 fi 455 425 done 456 426 done -
TOOLS/CMIP6_FORCING/SOLAR/CMIP6_historical_2bands.m
r3383 r3567 1 % Updated 05/02/2018, ThL 2 % * changed paths to retrieve new version of input data; 3 % * version considered is now 3.2. 4 1 5 % MATLAB run script to integrate spectral solar irradiance (SSI) data 2 6 % recommended and provided by SPARC/SOLARIS-HEPPA for use in CMIP6 to … … 51 55 % placed in dir_in); adapt to respective file that is to be considered 52 56 %filename_in='solarforcing_ref_mon_3.1.nc'; 53 filename_in='/prodigfs/ project/input4MIPs/SOLAR/transient/solarforcing_ref_day_3.1.nc';57 filename_in='/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/SOLARIS-HEPPA/SOLARIS-HEPPA-3-2/atmos/day/multiple/gz/v20170103/solarforcing-ref-day_input4MIPs_solar_CMIP_SOLARIS-HEPPA-3-2_gn_18500101-22991231.nc'; 54 58 55 59 % lines 53-54 only used to match input and output filenames -
TOOLS/CMIP6_FORCING/SOLAR/CMIP6_historical_6bands.m
r3383 r3567 1 % Updated 05/02/2018, ThL 2 % * changed paths to retrieve new version of input data; 3 % * version considered is now 3.2. 4 1 5 % MATLAB run script to integrate spectral solar irradiance (SSI) data 2 6 % recommended and provided by SPARC/SOLARIS-HEPPA for use in CMIP6 to … … 51 55 % placed in dir_in); adapt to respective file that is to be considered 52 56 %filename_in='solarforcing_ref_mon_3.1.nc'; 53 filename_in='/prodigfs/ project/input4MIPs/SOLAR/transient/solarforcing_ref_day_3.1.nc';57 filename_in='/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/SOLARIS-HEPPA/SOLARIS-HEPPA-3-2/atmos/day/multiple/gz/v20170103/solarforcing-ref-day_input4MIPs_solar_CMIP_SOLARIS-HEPPA-3-2_gn_18500101-22991231.nc' 54 58 55 59 % lines 53-54 only used to match input and output filenames -
TOOLS/CMIP6_FORCING/SOLAR/CMIP6_piControl_2bands.m
r3383 r3567 1 % Updated 05/02/2018, ThL 2 % * changed paths to retrieve new version of input data; 3 % * version considered is now 3.2. 4 1 5 % MATLAB run script to integrate spectral solar irradiance (SSI) data 2 6 % recommended and provided by SPARC/SOLARIS-HEPPA for use in CMIP6 to … … 50 54 %filename_in='solarforcing_ref_mon_3.1.nc'; 51 55 %filename_in='solarforcing_ref_day_3.1.nc'; 52 filename_in='/prodigfs/ project/input4MIPs/SOLAR/solarforcing_picontrol_fx_3.1.nc';56 filename_in='/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/SOLARIS-HEPPA/SOLARIS-HEPPA-3-2/atmos/fx/multiple/gz/v20170103/solarforcing-picontrol-fx_input4MIPs_solar_CMIP_SOLARIS-HEPPA-3-2_gn_18500101-18730128.nc' 53 57 54 58 % lines 53-54 only used to match input and output filenames -
TOOLS/CMIP6_FORCING/SOLAR/CMIP6_piControl_6bands.m
r3383 r3567 1 % Updated 05/02/2018, ThL 2 % * changed paths to retrieve new version of input data; 3 % * version considered is now 3.2. 4 1 5 % MATLAB run script to integrate spectral solar irradiance (SSI) data 2 6 % recommended and provided by SPARC/SOLARIS-HEPPA for use in CMIP6 to … … 50 54 %filename_in='solarforcing_ref_mon_3.1.nc'; 51 55 %filename_in='solarforcing_ref_day_3.1.nc'; 52 filename_in='/prodigfs/ project/input4MIPs/SOLAR/solarforcing_picontrol_fx_3.1.nc';56 filename_in='/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/SOLARIS-HEPPA/SOLARIS-HEPPA-3-2/atmos/fx/multiple/gz/v20170103/solarforcing-picontrol-fx_input4MIPs_solar_CMIP_SOLARIS-HEPPA-3-2_gn_18500101-18730128.nc' 53 57 54 58 % lines 53-54 only used to match input and output filenames -
TOOLS/CMIP6_FORCING/SOLAR/CMIP6_sensitivity_2bands.m
r3383 r3567 1 % Updated 05/02/2018, ThL 2 % * changed paths to retrieve new version of input data; 3 % * version considered is now 3.2. 4 1 5 % MATLAB run script to integrate spectral solar irradiance (SSI) data 2 6 % recommended and provided by SPARC/SOLARIS-HEPPA for use in CMIP6 to … … 51 55 % placed in dir_in); adapt to respective file that is to be considered 52 56 %filename_in='solarforcing_ext_mon_3.1.nc'; 53 filename_in='/prodigfs/ project/input4MIPs/SOLAR/sensitivity/solarforcing_ext_day_3.1.nc';57 filename_in='/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/SOLARIS-HEPPA/SOLARIS-HEPPA-3-2/atmos/day/multiple/gz/v20170103/solarforcing-ref-day_input4MIPs_solar_CMIP_SOLARIS-HEPPA-3-2_gn_18500101-22991231.nc' 54 58 55 59 % lines 53-54 only used to match input and output filenames -
TOOLS/CMIP6_FORCING/SOLAR/CMIP6_sensitivity_6bands.m
r3383 r3567 1 % Updated 05/02/2018, ThL 2 % * changed paths to retrieve new version of input data; 3 % * version considered is now 3.2. 4 1 5 % MATLAB run script to integrate spectral solar irradiance (SSI) data 2 6 % recommended and provided by SPARC/SOLARIS-HEPPA for use in CMIP6 to … … 51 55 % placed in dir_in); adapt to respective file that is to be considered 52 56 %filename_in='solarforcing_ext_mon_3.1.nc'; 53 filename_in='/prodigfs/ project/input4MIPs/SOLAR/sensitivity/solarforcing_ext_day_3.1.nc';57 filename_in='/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/SOLARIS-HEPPA/SOLARIS-HEPPA-3-2/atmos/day/multiple/gz/v20170103/solarforcing-ref-day_input4MIPs_solar_CMIP_SOLARIS-HEPPA-3-2_gn_18500101-22991231.nc' 54 58 55 59 % lines 53-54 only used to match input and output filenames -
TOOLS/CMIP6_FORCING/SOLAR/README
r3382 r3567 1 Update 05/02/2018 2 Thibaut Lurton 3 Scripts now have paths directed to version 3.2 of the data. 4 /prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/SOLARIS-HEPPA/SOLARIS-HEPPA-3-2/atmos/ 5 6 1 7 Solar forcing dataset from CMIP6 2 8 See http://www.wcrp-climate.org/wgcm-cmip/wgcm-cmip6
Note: See TracChangeset
for help on using the changeset viewer.