From c6977495c1ec6f07abd1afdffbda8fcb48fff293 Mon Sep 17 00:00:00 2001
From: Aiko Voigt <aiko.voigt@univie.ac.at>
Date: Tue, 7 May 2024 16:53:10 +0200
Subject: [PATCH] Updates Readme
---
README.md | 28 +++++++++++++++++++++++++++-
1 file changed, 27 insertions(+), 1 deletion(-)
diff --git a/README.md b/README.md
index 47a2e3c..46c337e 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`.
--
GitLab