====== a. Shell and CDO =======
Shown below is an example that interpolates one data file from the dpp0016 simulation onto a regular 1°x°1 global grid.
For the interpolation example shown here, a grid description text-file and a weight netcdf-file are needed.
The grid description text-file for an intented 1°x°1 globally grid is as follows:
# gridID 1
#
gridtype = lonlat
gridsize = 64800
xsize = 360
ysize = 180
xname = lon
xlongname = "longitude"
xunits = "degrees_east"
yname = lat
ylongname = "latitude"
yunits = "degrees_north"
xfirst = -179.5
xinc = 1
yfirst = -89.5
yinc = 1
Simply create a text file with any text editor, copy and paste the text above into the file and name the file, for example, "griddes_1x1.txt". It can be convenient to do this only once and keep it for later multiple applications in the users ''$HOME'' or ''$WORK''-directory. If the intended grid should be of a different resolution, say 0.1°x0.1°, then the grid discription would be :
# gridID 1
#
gridtype = lonlat
gridsize = 6480000
xsize = 3600
ysize = 1800
xname = lon
xlongname = "longitude"
xunits = "degrees_east"
yname = lat
ylongname = "latitude"
yunits = "degrees_north"
xfirst = -179.95
xinc = 0.1
yfirst = -89.95
yinc = 0.1
The weight file can be created using CDO. For this you need the grid description text-file shown above, the ICON grid-file for the native model grid (in this case R2B9 for the dpp0016 simulation) and also one dpp0016 data file is required. The ICON grid-file of the dpp0016 simulation can be found at /pool/data/ICON/grids/public/mpim/0015/icon_grid_0015_R02B09_G.nc. The general structure of the CDO command creating the weight file is:
cdo gendis,grid-description-text-file -setgrid,ICON_grid_file one_dpp0016_data weightfile
For example, storing the output weight file to the ''$SCRATCH''-directory:
cdo gendis,griddes_1x1.txt -setgrid,/pool/data/ICON/grids/public/mpim/0015/icon_grid_0015_R02B09_G.nc /work/mh0287/k203123/GIT/icon-aes-dyw_albW/experiments/dpp0016/dpp0016_atm_2d_ml_20200201T000000Z.nc $SCRATCH/weight_dpp0016_1x1.nc
Now, the interpolation can be done via:
cdo -remap,griddes_1x1.txt,$SCRATCH/weight_dpp0016_1x1.nc /work/mh0287/k203123/GIT/icon-aes-dyw_albW/experiments/dpp0016/dpp0016_atm_2d_ml_20200201T000000Z.nc $SCRATCH/dpp0016_atm_2d_ml_20200201T000000Z_1x1.nc