diff --git a/README.md b/README.md index 47a2e3c7a2ff25fda86ed4ac79965e3e0a9a38ac..46c337eee8b64c236452dccc3afdee261527e9a6 100644 --- a/README.md +++ b/README.md @@ -79,14 +79,40 @@ The simulation is AMIP-like because SSTs, sea ice etc. are prescribed as climato - slab4x-sun: 4xCO2 and sun dimmed by 3% via fsolrad=0.97 namelist switch, restarted from 1990-01-01 of slabctr - slab4x-vap: 4xCO2 and transparent stratospheric water vapor at klev=22 and above, restarted from 1990-01-01 of slabctr -Output is on native triangular grid and model levels. To interpolate to lat-lon-p, use +Output is on native triangular grid and model levels. To interpolate to lat-lon-p: ``` # 2d data: lat-lon interpolation cdo -P 16 -remapcon,r180x90 -setgrid,icon_grid_G.nc file.nc file.r180x90.nc +``` + +``` +# 2d daily-mean precip data interpolated to nearest neighbour for analysis of extreme precipitation +cdo -P 16 -remapnn,r360x90 -setgrid,icon_grid_G.nc -selvar,pr file_dm.nc file_dm.pr.remapnn-r360x180.nc +``` +For example: +``` +for file in slab4x-vap_atm_2d_daily_ml_????????T000000Z.nc; do echo $file; cdo -remapnn,r360x180 -setgrid,icon_grid_G.nc -selvar,pr $file ${file}.pr.remapnn-r360x180.nc; done + +cdo mergetime slab4x-vap_atm_2d_daily_ml_*pr.remapnn-r360x180.nc ../output4students/slab4x-vap/slab4x-vap_atm_2d_daily_ml_199001-201701.pr.remapnn-r360x180.nc + +# clean up if needed: files with one month each +rm slab4x-vap_atm_2d_daily_ml_*.pr.remapnn-r360x180.nc +``` + +``` # 3d data: lat-lon-p interpolation in one go plev="100000,95000,90000,85000,80000,75000,70000,65000,60000,55000,50000,45000,40000,35000,30000,25000,20000,15000,10000,5000,10" cdo -s -P 16 -remapcon,r180x90 -setgrid,icon_grid_G.nc -ap2pl,$plev file.nc file.r180x90.plev.nc ``` + +For example: +``` +cdo -P 16 mergetime slabctr_atm_3d_ml_????????T000000Z.nc slabctr_atm_3d_ml_197901-201701.nc +cdo -P 16 -ap2pl,$plev slabctr_atm_3d_ml_197901-201701.nc slabctr_atm_3d_ml_197901-201701.plev.nc +cdo -s -P 16 -remapcon,r180x90 -setgrid,icon_grid_G.nc slabctr_atm_3d_ml_197901-201701.plev.nc slabctr_atm_3d_ml_197901-201701.plev.r180x90.nc +``` + + Interpolated output is on VSC4 in `/gpfs/data/fs72044/avoigt_teach/experiments/s2024/output4students`.