[4993] | 1 | #!/bin/bash |
---|
| 2 | set -vx |
---|
| 3 | |
---|
| 4 | # Script that prepares the PCMDI-AMIP-1-1-6 files to be read by ce0l in LMDZ. |
---|
| 5 | # These new files are stored in IGCM/ATM/LIMIT/AMIP.v20191121/original |
---|
| 6 | # and can be used to produce limit.nc files on model grid. |
---|
| 7 | # Done by J. Ghattas, 2020/03/03 |
---|
| 8 | |
---|
| 9 | dirold=original/360x180.orig |
---|
| 10 | dirnew=original/360x180.new |
---|
| 11 | |
---|
| 12 | # 0) Copy original files from ciclad Input4MIPS and save in $dirold |
---|
| 13 | #/bdd/input4MIPs/CMIP6/CMIP/PCMDI/PCMDI-AMIP-1-1-6/ocean/mon/tosbcs/gn/v20191121/tosbcs_input4MIPs_SSTsAndSeaIce_CMIP_PCMDI-AMIP-1-1-6_gn_187001-201812.nc |
---|
| 14 | #/bdd/input4MIPs/CMIP6/CMIP/PCMDI/PCMDI-AMIP-1-1-6/seaIce/mon/siconcbcs/gn/v20191121/siconcbcs_input4MIPs_SSTsAndSeaIce_CMIP_PCMDI-AMIP-1-1-6_gn_187001-201812.nc |
---|
| 15 | |
---|
| 16 | # 1) Split files to yearly files |
---|
| 17 | cd $dirold |
---|
[5011] | 18 | cdo splityear tosbcs_input4MIPs_SSTsAndSeaIce_CMIP_PCMDI-AMIP-1-1-6_gn_187001-201812.nc amipbc_sst_360x180_ |
---|
| 19 | cdo splityear siconcbcs_input4MIPs_SSTsAndSeaIce_CMIP_PCMDI-AMIP-1-1-6_gn_187001-201812.nc amipbc_sic_360x180_ |
---|
[4993] | 20 | |
---|
| 21 | |
---|
| 22 | # 2) Produce climatologie files |
---|
| 23 | ncea amipbc_sst_360x180_1979.nc amipbc_sst_360x180_198* amipbc_sst_360x180_199* amipbc_sst_360x180_200[012345678]* amipbc_sst_360x180_1979_2008_clim.nc |
---|
| 24 | ncea amipbc_sic_360x180_1979.nc amipbc_sic_360x180_198* amipbc_sic_360x180_199* amipbc_sic_360x180_200[012345678]* amipbc_sic_360x180_1979_2008_clim.nc |
---|
| 25 | |
---|
| 26 | # 3) Rename variable name for sea-ice concentration into sicbcs |
---|
| 27 | # and 4) The time axis is set fixed instead of unlimited |
---|
| 28 | mkdir -p $dirnew |
---|
| 29 | |
---|
| 30 | Date_Begin=1870 |
---|
| 31 | Date_End=2018 |
---|
| 32 | year=${Date_Begin} |
---|
| 33 | while [ ${year} -le ${Date_End} ] ; do |
---|
| 34 | echo year = $year |
---|
| 35 | filesic=amipbc_sic_360x180_${year}.nc |
---|
| 36 | filesst=amipbc_sst_360x180_${year}.nc |
---|
| 37 | |
---|
| 38 | echo "Rename variable sicbcs in file $dirold/$filesic " |
---|
[5011] | 39 | ncrename -v siconcbcs,sicbcs $dirold/${filesic} $dirnew/tmpsic.nc |
---|
[4993] | 40 | echo "Copy sst file to new folder $dirold/$filesst " |
---|
| 41 | cp $dirold/${filesst} $dirnew/tmpsst.nc |
---|
| 42 | |
---|
| 43 | # The time axis is set fixed instead of unlimited and the files are transformend into netcdf-4 format. |
---|
| 44 | # This is done so that the files can be read by XIOS. |
---|
| 45 | echo "Change time axis" |
---|
| 46 | ncks -4 --no_rec_dmn time $dirnew/tmpsic.nc -o $dirnew/${filesic} |
---|
| 47 | ncks -4 --no_rec_dmn time $dirnew/tmpsst.nc -o $dirnew/${filesst} |
---|
| 48 | |
---|
| 49 | let year=${year}+1 |
---|
| 50 | echo "" |
---|
| 51 | rm -f $dirnew/tmpsic.nc $dirnew/tmpsst.nc |
---|
| 52 | done |
---|
| 53 | |
---|
| 54 | # Do 3) and 4) also for the climatologie files |
---|
| 55 | filesic=amipbc_sic_360x180_1979_2008_clim.nc |
---|
| 56 | filesst=amipbc_sst_360x180_1979_2008_clim.nc |
---|
| 57 | |
---|
| 58 | echo "Rename variable sicbcs in file $dirold/$filesic " |
---|
[5011] | 59 | ncrename -v siconcbcs,sicbcs $dirold/${filesic} $dirnew/tmpsic.nc |
---|
[4993] | 60 | |
---|
| 61 | echo "Copy sst file to new folder $dirold/$filesst " |
---|
| 62 | cp $dirold/${filesst} $dirnew/tmpsst.nc |
---|
| 63 | |
---|
| 64 | echo "Change time axis" |
---|
| 65 | ncks -4 --no_rec_dmn time $dirnew/tmpsic.nc -o $dirnew/${filesic} |
---|
| 66 | ncks -4 --no_rec_dmn time $dirnew/tmpsst.nc -o $dirnew/${filesst} |
---|
| 67 | |
---|
| 68 | rm -f $dirnew/tmpsic.nc $dirnew/tmpsst.nc |
---|
| 69 | |
---|
| 70 | |
---|