Adds 3d postprocessing script

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
 ```
 
 

postprocessing/pp_3d.sh

+#!/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