Changeset 3567


Ignore:
Timestamp:
02/07/18 14:22:26 (6 years ago)
Author:
tlurton
Message:

Update Th. Lurton
7 Feb. 2018
Changed input file paths to fetch new versions of the files.
For SOLAR, AER_TROP_EMISSIONS, GHG, NITROGEN and OZONE.

Location:
TOOLS/CMIP6_FORCING
Files:
14 edited

Legend:

Unmodified
Added
Removed
  • TOOLS/CMIP6_FORCING/AER_TROP_EMISSIONS/README

    r3398 r3567  
    3030Update: 02/10/2017 
    3131report NOx emissions as NO for both FF and BB emissions 
     32 
     33Update: 05/02/2018 
     34New version of regrid.sh 
     35Diffs: * 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 
    3238 
    3339********************************************************** 
  • TOOLS/CMIP6_FORCING/AER_TROP_EMISSIONS/REGRID/regrid.sh

    r3398 r3567  
     1#--updated on 5/2/2018 with new paths, ThL 
    12#--updated on 4/5/2017 with improved Sheng & Zwiers algorithm, ThL 
    23#--corrected some interpolation preprocessing (compared to v4) 
     
    1516 
    1617#--input directory for anthropogenic (non-BB) emissions 
    17 dirinPNNL="/prodigfs/project/input4MIPs/PNNL-JGCRI/emissions/CMIP/CEDS-2017-05-18/mon/" 
     18dirinPNNL="/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/PNNL-JGCRI/CEDS-2017-05-18/atmos/mon/" 
    1819##--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/ 
    2021#--input directory for anthropogenic (BB) emissions 
    21 dirinVUA="/prodigfs/project/input4MIPs/VUA/emissions/CMIP/VUA-CMIP-BB4CMIP6-1-2/mon/" 
     22dirinVUA="/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/VUA/VUA-CMIP-BB4CMIP6-1-2/atmos/mon/" 
    2223 
    2324#--LMDz grid information 
     
    7172fi 
    7273 
    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.nc 
     74#--input file PNNL (fossil fuel emissions)       Underscores in the old species-directory are now hyphens. 
     75filename=${dirinPNNL}${species}-em-anthro/gn/v20170519/${species}-em-anthro_input4MIPs_emissions_CMIP_CEDS-2017-05-18_gn_${year1}01-${year2}12.nc 
    7576 
    7677#--output files 
     
    266267fi 
    267268 
    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" 
     270filename=${dirinVUA}${species}"/gn/v20161213/"${species}"-em-biomassburning_input4MIPs_emissions_CMIP_VUA-CMIP-BB4CMIP6-1-2_gn_${year1}01-${year2}12.nc" 
    270271 
    271272#--output files 
  • TOOLS/CMIP6_FORCING/GHG/README

    r3386 r3567  
     1Script updated on 5 Feb. 2018 
     2Th. 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 
    17GHG data for CMIP6 are initially coming from  
    28ftp://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 
    16#--v3 of my script to extract global mean and put to IGCM format 
    27 
     
    1116 
    1217#--dir for annual mean data 
    13 dirin1='/prodigfs/project/input4MIPs/UoM/GHGConcentrations/CMIP/UoM-CMIP-1-2-0/yr/' 
    14 dirin2='gr1-GMNHSH/v20160902' 
     18dirin1='/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/UoM/UoM-CMIP-1-2-0/atmos/yr/' 
     19dirin2='gr1-GMNHSH/v20160830' 
     20#--dirin2 versioning date is different from previously 
    1521 
    1622dirout='./OUTPUT' 
     
    2531#--N2O we leave unit to ppb 
    2632GHG='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' 
     33name='nitrous-oxide' 
     34file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 
    2935rm -f tmp.txt ${dirout}/${GHG}_CMIP6_0000_2014.txt 
    30 dirin=${dirin1}'mole_fraction_of_'${name}'_in_air/'${dirin2} 
     36dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 
    3137cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt 
    3238paste years.txt tmp.txt > ${dirout}/${GHG}_CMIP6_0000_2014.txt 
     
    3541GHG='CH4' 
    3642name='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} 
     43file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 
     44dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 
    3945rm -f tmp.txt ${dirout}/${GHG}_CMIP6_0000_2014.txt 
    4046cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt 
     
    4349#--CO2 here we leave unit in ppm 
    4450GHG='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} 
     51name='carbon-dioxide' 
     52file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 
     53dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 
    4854rm -f tmp.txt ${dirout}/${GHG}_CMIP6_0000_2014.txt 
    4955cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt 
     
    5460GHG='CFC11eq' 
    5561name='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} 
     62file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 
     63dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 
    5864rm -f tmp.txt ${dirout}/${GHG}_0000_2014.txt 
    5965cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt 
     
    6369GHG='CFC12eq' 
    6470name='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} 
     71file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 
     72dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 
    6773rm -f tmp.txt ${dirout}/${GHG}_0000_2014.txt ${dirout}/${GHG}_CMIP6_0000_2014.txt 
    6874cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt 
     
    7581GHG='CFC11' 
    7682name='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} 
     83file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 
     84dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 
    7985rm -f tmp.txt ${dirout}/${GHG}_0000_2014.txt 
    8086cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt 
     
    8591GHG='CFC12' 
    8692name='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} 
     93file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 
     94dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 
    8995rm -f tmp.txt ${dirout}/${GHG}_0000_2014.txt 
    9096cdo outputf,%12.4e,3 -selvar,'mole_fraction_of_'${name}'_in_air' ${dirin}/${file} | awk '{print $1}' > tmp.txt 
     
    94100GHG='HFC134aeq' 
    95101name='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} 
     102file='mole-fraction-of-'${name}'-in-air_input4MIPs_GHGConcentrations_CMIP_UoM-CMIP-1-2-0_gr1-GMNHSH_0000-2014.nc' 
     103dirin=${dirin1}'mole-fraction-of-'${name}'-in-air/'${dirin2} 
    98104rm -f tmp.txt ${dirout}/${GHG}_0000_2014.txt ${dirout}/CFC11eq_CMIP6_0000_2014.txt 
    99105cdo 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  
     1Update 05 Feb. 2018 
     2Th. Lurton 
     3** Changed paths to new version 
     4 
     5 
    16Here we only add drynhx+wetnhx and drynoy+wetnoy into nhx and noy deposition fluxes (kg N/m2/s). 
    27There 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" 
     1dirin1="/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/NCAR/NCAR-CCMI-2-0/atmos/mon" 
    22dirin2="gn/v20161207" 
     3dirinC1="/prodigfs/esgf/mirror/input4MIPs/CMIP6/CMIP/NCAR/NCAR-CCMI-2-0/atmos/monC" 
     4 
    35dirout="./OUTPUT" 
    46if [ ! -d ${dirout} ] ; then mkdir ${dirout} ; fi 
     
    68#--1850 climatological N deposition  
    79# 
    8 drynhx=${dirin1}"/drynhx/"${dirin2}"/drynhx_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012-clim.nc" 
    9 wetnhx=${dirin1}"/wetnhx/"${dirin2}"/wetnhx_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012-clim.nc" 
    10 drynoy=${dirin1}"/drynoy/"${dirin2}"/drynoy_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012-clim.nc" 
    11 wetnoy=${dirin1}"/wetnoy/"${dirin2}"/wetnoy_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012-clim.nc" 
     10drynhx=${dirinC1}"/drynhx/"${dirin2}"/drynhx_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012.nc" 
     11wetnhx=${dirinC1}"/wetnhx/"${dirin2}"/wetnhx_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012.nc" 
     12drynoy=${dirinC1}"/drynoy/"${dirin2}"/drynoy_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012.nc" 
     13wetnoy=${dirinC1}"/wetnoy/"${dirin2}"/wetnoy_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-185012.nc" 
    1214# 
    13 cdo merge $drynhx $wetnhx $dirout/CMIP6_fluxnhx_merged_clim_1850.nc 
    14 cdo merge $drynoy $wetnoy $dirout/CMIP6_fluxnoy_merged_clim_1850.nc 
    15 cdo expr,'fluxnhx=drynhx+wetnhx' $dirout/CMIP6_fluxnhx_merged_clim_1850.nc $dirout/CMIP6_fluxnhx_clim_1850.nc 
    16 cdo expr,'fluxnoy=drynoy+wetnoy' $dirout/CMIP6_fluxnoy_merged_clim_1850.nc $dirout/CMIP6_fluxnoy_clim_1850.nc 
    17 rm -f $dirout/CMIP6_fluxnhx_merged_clim_1850.nc $dirout/CMIP6_fluxnoy_merged_clim_1850.nc 
     15cdo merge ${drynhx} ${wetnhx} ${dirout}/CMIP6_fluxnhx_merged_clim_1850.nc 
     16cdo merge ${drynoy} ${wetnoy} ${dirout}/CMIP6_fluxnoy_merged_clim_1850.nc 
     17cdo expr,'fluxnhx=drynhx+wetnhx' ${dirout}/CMIP6_fluxnhx_merged_clim_1850.nc ${dirout}/CMIP6_fluxnhx_clim_1850.nc 
     18cdo expr,'fluxnoy=drynoy+wetnoy' ${dirout}/CMIP6_fluxnoy_merged_clim_1850.nc ${dirout}/CMIP6_fluxnoy_clim_1850.nc 
     19rm -f ${dirout}/CMIP6_fluxnhx_merged_clim_1850.nc ${dirout}/CMIP6_fluxnoy_merged_clim_1850.nc 
    1820# 
    1921# 
     
    2527wetnoy=${dirin1}"/wetnoy/"${dirin2}"/wetnoy_input4MIPs_surfaceFluxes_CMIP_NCAR-CCMI-2-0_gn_185001-201412.nc" 
    2628# 
    27 cdo merge $drynhx $wetnhx $dirout/CMIP6_fluxnhx_merged_annual_1850-2014.nc 
    28 cdo merge $drynoy $wetnoy $dirout/CMIP6_fluxnoy_merged_annual_1850-2014.nc 
    29 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.nc 
     29cdo merge ${drynhx} ${wetnhx} ${dirout}/CMIP6_fluxnhx_merged_annual_1850-2014.nc 
     30cdo merge ${drynoy} ${wetnoy} ${dirout}/CMIP6_fluxnoy_merged_annual_1850-2014.nc 
     31cdo splityear -expr,'fluxnhx=drynhx+wetnhx' ${dirout}/CMIP6_fluxnhx_merged_annual_1850-2014.nc ${dirout}/CMIP6_fluxnhx_annual_1850-2014.nc 
     32cdo splityear -expr,'fluxnoy=drynoy+wetnoy' ${dirout}/CMIP6_fluxnoy_merged_annual_1850-2014.nc ${dirout}/CMIP6_fluxnoy_annual_1850-2014.nc 
     33rm -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  
    11#!/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  
     2module load nco/4.4.8 
     3 
     4ia="1"    #=== Indices of chosen Amip  datasets (see below AMIP_VERSION_NAMES/AMIP_VERSION_FOLDERS) 
     5io="1"    #=== Indices of chosen Ozone datasets (see below OZON_VERSION_NAMES/OZON_VERSION_FOLDERS) 
    96 
    107#=============================================================================== 
     
    1310if   [ ${HOSTNAME:0:3} = 'ada'    ]; then machine='ada';    work=$WORKDIR 
    1411elif [ ${HOSTNAME:0:5} = 'curie'  ]; then machine='curie';  work=$CCCWORKDIR 
    15 elif [ ${HOSTNAME:0:6} = 'ciclad' ]; then machine='ciclad'; work=/data/$USER 
     12elif [ ${HOSTNAME:0:6} = 'ciclad' ]; then machine='ciclad'; work=/home/$USER/tmp 
    1613else echo "Not set up for this machine yet, sorry."; exit; fi 
    1714export machine=$machine 
     
    2118#=============================================================================== 
    2219local=$PWD 
    23 echo $machine 
    2420case $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' ;; 
    3023  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'     ;; 
    4027esac 
    41  
    42 if [ ! -d ${DATAOU} ] ; then mkdir -p ${DATAOU} ; fi 
     28DATAOU=/home/tlurton/IGCM     ##### !!!!!! PROVISOIRE !!!!!!! 
    4329 
    4430#--- 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' 
     31AMIP_VERSION_NAMES=('PCMDI-AMIP-1-1-3') 
     32AMIP_VERSION_FOLDERS=('gn/v20161031') 
     33AMIP_FOLDER_IN=$DATAIN'/PCMDI/${amip_version_name}/${W}/mon/${V}/${amip_version_folder}' 
     34amip_tname_in=('ocean' 'seaIce') 
     35amip_vname_in=('tos' 'siconc') 
     36amip_fname_in='${V}_input4MIPs_SSTsAndSeaIce_CMIP_${amip_version_name}_${amip_version_folder%/*}_187001-20????.nc' 
    5037 
    5138#--- INPUT DATA: OZONE 
    5239OZON_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} 
     40OZON_VERSION_FOLDERS=('gn/v20160711') 
     41OZON_VERSION_FOLDERS_CLIM=('gn/v20160830') 
     42OZON_FOLDER_IN=$DATAIN/'UReading/${ozon_version_name}/atmos/mon${clim}/${V}/${ozon_version_folder}' 
    5543ozon_vname_in='vmro3' 
    56 ozon_fname_in='${V}_input4MIPs_ozone_CMIP_${ozon_version_name}_${ozon_version_folder%/*}*_${Yi}01-${Yf}12.nc' 
     44ozon_fname_in='${V}_input4MIPs_ozone_CMIP_${ozon_version_name}_${ozon_version_folder%/*}_${Yi}01-${Yf}12.nc' 
    5745 
    5846#--- VERSIONS OF INPUT FILES EFFECTIVELY CHOSEN 
    5947((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]} 
    6248 
    6349#--- 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' 
     50DATAOU=${DATAOU%/*}/ANTICHAMBRE; if [ ! -d $DATAOU ]; then mkdir -p $DATAOU; fi 
     51AMIP_FOLDER_OU=$DATAOU/'ATM/LIMIT/AMIP.${amip_version_folder##*/}/${suf}' 
     52OZON_FOLDER_OU=$DATAOU/'ATM/OZONE/UReading.${ozon_version_folder##*/}' 
     53amip_vname_ou=('tosbcs' 'sicbcs'); amip_fname_ou='amip${bc}_${V}_360x180_${Y}.nc' 
    6754ozon_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 
     55Y_amip_beg=1870 ; Y_amip_end=2017    #--- FOR AMIP  FILES (INTERANNUAL) 
     56Y_ozon_beg=1850 ; Y_ozon_end=2014    #--- FOR OZONE FILES (INTERANNUAL) 
     57Y_clim_beg=1979 ; Y_clim_end=2008    #--- FOR CLIMATOLOGIES (COMMON PERIOD FOR O3 AND SIC/SST) 
     58Y_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) 
     62recomp='n' 
     63do_amip='n' 
     64do_ozon='y' 
    7265 
    7366if [ "$do_amip" = "y" ]; then 
     67  amip_version_name=${AMIP_VERSION_NAMES[$ia]}; amip_version_folder=${AMIP_VERSION_FOLDERS[$ia]} 
    7468  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)" 
    7569  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)" 
    7670fi 
    7771if [ "$do_ozon" = "y" ]; then 
     72  ozon_version_name=${OZON_VERSION_NAMES[$io]}; ozon_version_folder=${OZON_VERSION_FOLDERS[$io]} 
    7873  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)" 
    7974fi 
     
    231226 
    232227#================================================================================== 
    233 function tim_len { 
    234   local n 
    235   n=`ncdump -h $1 | grep 'UNLIMITED'`; n=${n#*// (}; n=${n% currently*}; echo $n 
    236 } 
    237  
    238 #================================================================================== 
    239 function dim_len { 
    240   local n 
    241   if [ "$2" = "T" ]; then tim_len $1; else n=`ncdump -h $1 | grep $(dim_nam $1 $2)' = '`; n=${n#* = }; n=${n% ;*}; echo $n; fi 
    242 } 
    243  
    244 #================================================================================== 
    245228function test_var { 
    246229# Input arguments: 
     
    257240#    Rename dimensions, dimensional variables and possibly bounds_<variable>. 
    258241# 
    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 
    261243  shift 
    262244  while [ $# -ge 1 ]; do 
    263245    nam_ou=${1##*:}; nam_in=$(dim_nam $f ${1%:*}); shift 
    264246    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" 
    267248      nam_in=bounds_$nam_in; nam_ou=bounds_$nam_ou 
    268249      test_var $f $nam_in 
    269       if [ $? -eq 0 ]; then out="$out -v $nam_in,$nam_ou"; fi 
     250      if [ $? -eq 0 ]; then out="$out -d $nam_in,$nam_ou -v $nam_in,$nam_ou"; fi 
    270251    fi 
    271252    ((ia++)) 
     
    282263 
    283264#=== BUILD YEARLY INTERANNUAL SIC/SST FILES WITH 12 RECORDS 
    284 for is in $ia; do                                                       #--- Loop on Amip versions 
     265for is in $ia; do                                                            #--- Loop on Amip versions 
    285266  amip_version_name=${AMIP_VERSION_NAMES[$is]} 
    286267  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 
    309288 
    310289        #--- RECOMPUTE OR NOT 
     
    312291 
    313292        #--- 2D FIELDS EXTRACTION 
    314         extract $v_in,$v_ou $AMIP_FOLDER_IN/$amip_fname_in,$f_ou $Yb,$Ye,0 $Y 0 
    315  
    316         progress_bar $((Y-Yb+1)) $((Ye-Yb+1)) 50 
    317      done 
    318    done 
     293        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 
    319298  done 
    320299done 
    321300 
    322301#=== BUILD SIC/SST 12 RECORDS CLIMATOLOGIES 
    323 for is in $ia; do                                                       #--- Loop on Amip versions 
     302for is in $ia; do                                                            #--- Loop on Amip versions 
    324303  amip_version_name=${AMIP_VERSION_NAMES[$is]} 
    325304  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 
    344327    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/null 
    346328  done 
    347329done 
     
    358340 
    359341#=== BUILD YEARLY INTERANNUAL OZONE FILES WITH 12 RECORDS 
    360 Yb=${Y_ozon%%-*}; Ye=${Y_ozon##*-}; dY=${Ye#*:}; Ye=${Ye%:*} 
    361342for is in $io; do                                                 #--- Loop on Ozone versions 
    362343  ozon_version_name=${OZON_VERSION_NAMES[$is]} 
    363344  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 
    366348  v_in=${ozon_vname_in}                                           #--- Input  variable 
    367   v_ou=${ozon_vname_ou}                                           #--- Output variable 
    368  
    369   #--- INPUT/OUTPUT FILE NAME 
    370349  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]}) 
    371350  f_ou=$(V=$v_ou; Y=\$Y; eval echo $ozon_fname_ou) 
    372351 
    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) 
    381355  if [ ! -d $d_o2 ]; then mkdir -p $d_o2; fi; f_o2=$d_o2/$f_ou 
    382356  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 {${Yb}..${Ye}}); do      #--- Loop on years 
     357 
     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 
    386360  
    387361    #--- RECOMPUTE OR NOT 
     
    389363 
    390364    #--- 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 
    393366 
    394367    #--- RENAME DIMENSIONS ACCORDING TO $coords VARIABLE 
    395368    fo3=$(Y=$Y; eval echo $f_o3); fo2=$(Y=$Y; eval echo $f_o2) 
    396     if [ $Y -eq $Yb ]; then args_ren=$(renam_dims $fo3 $coords); fi 
     369    if [ $Y -eq $Y_ozon_beg ]; then args_ren=$(renam_dims $fo3 $coords); fi 
    397370    if [ "$args_ren" != "" ]; then ncrename $args_ren $fo3 2>&1 > /dev/null ; fi 
    398371 
     
    405378    if [ $? -eq 0 ]; then ncks -x -v $v $fo2 $tmp ; mv $tmp $fo2; fi 
    406379 
    407     progress_bar $((Y-Yb+1)) $((Ye-Yb+1)) $dY 
     380    progress_bar $((Y-Y_ozon_beg+1)) $((Y_ozon_end-Y_ozon_beg+1)) 50 
    408381 
    409382  done 
     
    416389  v_ou=${ozon_vname_ou}                                           #--- Output variable 
    417390 
    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 
    428423      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 
    455425  done 
    456426done 
  • 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 
    15% MATLAB run script to integrate spectral solar irradiance (SSI) data 
    26% recommended and provided by SPARC/SOLARIS-HEPPA for use in CMIP6 to 
     
    5155% placed in dir_in); adapt to respective file that is to be considered 
    5256%filename_in='solarforcing_ref_mon_3.1.nc'; 
    53 filename_in='/prodigfs/project/input4MIPs/SOLAR/transient/solarforcing_ref_day_3.1.nc'; 
     57filename_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'; 
    5458 
    5559% 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 
    15% MATLAB run script to integrate spectral solar irradiance (SSI) data 
    26% recommended and provided by SPARC/SOLARIS-HEPPA for use in CMIP6 to 
     
    5155% placed in dir_in); adapt to respective file that is to be considered 
    5256%filename_in='solarforcing_ref_mon_3.1.nc'; 
    53 filename_in='/prodigfs/project/input4MIPs/SOLAR/transient/solarforcing_ref_day_3.1.nc'; 
     57filename_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' 
    5458 
    5559% 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 
    15% MATLAB run script to integrate spectral solar irradiance (SSI) data 
    26% recommended and provided by SPARC/SOLARIS-HEPPA for use in CMIP6 to 
     
    5054%filename_in='solarforcing_ref_mon_3.1.nc'; 
    5155%filename_in='solarforcing_ref_day_3.1.nc'; 
    52 filename_in='/prodigfs/project/input4MIPs/SOLAR/solarforcing_picontrol_fx_3.1.nc'; 
     56filename_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' 
    5357 
    5458% 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 
    15% MATLAB run script to integrate spectral solar irradiance (SSI) data 
    26% recommended and provided by SPARC/SOLARIS-HEPPA for use in CMIP6 to 
     
    5054%filename_in='solarforcing_ref_mon_3.1.nc'; 
    5155%filename_in='solarforcing_ref_day_3.1.nc'; 
    52 filename_in='/prodigfs/project/input4MIPs/SOLAR/solarforcing_picontrol_fx_3.1.nc'; 
     56filename_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' 
    5357 
    5458% 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 
    15% MATLAB run script to integrate spectral solar irradiance (SSI) data 
    26% recommended and provided by SPARC/SOLARIS-HEPPA for use in CMIP6 to 
     
    5155% placed in dir_in); adapt to respective file that is to be considered 
    5256%filename_in='solarforcing_ext_mon_3.1.nc'; 
    53 filename_in='/prodigfs/project/input4MIPs/SOLAR/sensitivity/solarforcing_ext_day_3.1.nc'; 
     57filename_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' 
    5458 
    5559% 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 
    15% MATLAB run script to integrate spectral solar irradiance (SSI) data 
    26% recommended and provided by SPARC/SOLARIS-HEPPA for use in CMIP6 to 
     
    5155% placed in dir_in); adapt to respective file that is to be considered 
    5256%filename_in='solarforcing_ext_mon_3.1.nc'; 
    53 filename_in='/prodigfs/project/input4MIPs/SOLAR/sensitivity/solarforcing_ext_day_3.1.nc'; 
     57filename_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' 
    5458 
    5559% lines 53-54 only used to match input and output filenames 
  • TOOLS/CMIP6_FORCING/SOLAR/README

    r3382 r3567  
     1Update 05/02/2018 
     2Thibaut Lurton 
     3Scripts 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 
    17Solar forcing dataset from CMIP6  
    28See http://www.wcrp-climate.org/wgcm-cmip/wgcm-cmip6 
Note: See TracChangeset for help on using the changeset viewer.