From 66245cdb7249dd05aaef13d4fbcfd891c034e782 Mon Sep 17 00:00:00 2001 From: Aiko Voigt <aiko.voigt@univie.ac.at> Date: Tue, 7 May 2024 17:21:55 +0200 Subject: [PATCH] Adds 3d postprocessing script --- README.md | 6 +++--- postprocessing/pp_3d.sh | 48 +++++++++++++++++++++++++++++++++++++++++ 2 files changed, 51 insertions(+), 3 deletions(-) create mode 100644 postprocessing/pp_3d.sh diff --git a/README.md b/README.md index 46c337e..44b2d20 100644 --- a/README.md +++ b/README.md @@ -109,9 +109,9 @@ cdo -s -P 16 -remapcon,r180x90 -setgrid,icon_grid_G.nc -ap2pl,$plev file.nc file 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 +cdo -O -P 16 mergetime slabctr_atm_3d_ml_????????T000000Z.nc slabctr_atm_3d_ml_197901-201701.nc +cdo -O -P 16 -ap2pl,$plev slabctr_atm_3d_ml_197901-201701.nc slabctr_atm_3d_ml_197901-201701.plev.nc +cdo -O -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 ``` diff --git a/postprocessing/pp_3d.sh b/postprocessing/pp_3d.sh new file mode 100644 index 0000000..32e0929 --- /dev/null +++ b/postprocessing/pp_3d.sh @@ -0,0 +1,48 @@ +#!/bin/bash +# +#SBATCH -J pp_3d +#SBATCH -N 1 +#SBATCH --ntasks-per-node=1 +#SBATCH --ntasks-per-core=48 +#SBATCH --mail-type=BEGIN # first have to state the type of event to occur +#SBATCH --mail-user=aiko.voigt@univie.ac.at # and then your email address +#SBATCH --time=16:00:00 + +module load cdo --auto + +# interpolate to these pressure levels +plev="100000,95000,90000,85000,80000,75000,70000,65000,60000,55000,50000,45000,40000,35000,30000,25000,20000,15000,10000,5000,10" + +echo "Working on slabctr:" +cd /gpfs/data/fs72044/avoigt_teach/experiments/s2024/slabctr +cdo -P 48 mergetime slabctr_atm_3d_ml_????????T000000Z.nc slabctr_atm_3d_ml_197901-201701.nc +cdo -P 48 -ap2pl,$plev slabctr_atm_3d_ml_197901-201701.nc slabctr_atm_3d_ml_197901-201701.plev.nc +cdo -s -P 48 -remapcon,r180x90 -setgrid,icon_grid_G.nc slabctr_atm_3d_ml_197901-201701.plev.nc slabctr_atm_3d_ml_197901-201701.plev.r180x90.nc +cp slabctr_atm_3d_ml_197901-201701.plev.r180x90.nc /gpfs/data/fs72044/avoigt_teach/experiments/s2024/output4students/slabctr +echo "slabctr finished!" + +echo "Working on slab4x:" +cd /gpfs/data/fs72044/avoigt_teach/experiments/s2024/slab4x +cdo -P 48 mergetime slab4x_atm_3d_ml_????????T000000Z.nc slab4x_atm_3d_ml_199002-201701.nc +cdo -P 48 -ap2pl,$plev slab4x_atm_3d_ml_197901-201701.nc slab4x_atm_3d_ml_199002-201701.plev.nc +cdo -s -P 48 -remapcon,r180x90 -setgrid,icon_grid_G.nc slab4x_atm_3d_ml_199002-201701.plev.nc slab4x_atm_3d_ml_199002-201701.plev.r180x90.nc +cp slab4x_atm_3d_ml_199002-201701.plev.r180x90.nc /gpfs/data/fs72044/avoigt_teach/experiments/s2024/output4students/slab4x +echo "slab4x finished!" + +echo "Working on slab4x-sun:" +cd /gpfs/data/fs72044/avoigt_teach/experiments/s2024/slab4x-sun +cdo -P 48 mergetime slab4x-sun_atm_3d_ml_????????T000000Z.nc slab4x-sun_atm_3d_ml_199002-201701.nc +cdo -P 48 -ap2pl,$plev slab4x-sun_atm_3d_ml_197901-201701.nc slab4x-sun_atm_3d_ml_199002-201701.plev.nc +cdo -s -P 48 -remapcon,r180x90 -setgrid,icon_grid_G.nc slab4x-sun_atm_3d_ml_199002-201701.plev.nc slab4x-sun_atm_3d_ml_199002-201701.plev.r180x90.nc +cp slab4x-sun_atm_3d_ml_199002-201701.plev.r180x90.nc /gpfs/data/fs72044/avoigt_teach/experiments/s2024/output4students/slab4x-sun +echo "slab4x-sun finished!" + +echo "Working on slab4x-vap:" +cd /gpfs/data/fs72044/avoigt_teach/experiments/s2024/slab4x-vap +cdo -P 48 mergetime slab4x-vap_atm_3d_ml_????????T000000Z.nc slab4x-vap_atm_3d_ml_199002-201701.nc +cdo -P 48 -ap2pl,$plev slab4x-sun_atm_3d_ml_197901-201701.nc slab4x-vap_atm_3d_ml_199002-201701.plev.nc +cdo -s -P 48 -remapcon,r180x90 -setgrid,icon_grid_G.nc slab4x-vap_atm_3d_ml_199002-201701.plev.nc slab4x-vap_atm_3d_ml_199002-201701.plev.r180x90.nc +cp slab4x-vap_atm_3d_ml_199002-201701.plev.r180x90.nc /gpfs/data/fs72044/avoigt_teach/experiments/s2024/output4students/slab4x-vap +echo "slab4x-vap finished!" + +cd /home/fs72044/avoigt_teach/climlab_s2024/msc-climmodlab-s2024/postprocessing -- GitLab