前言
近期,开始处理CMIP6模式数据,由于不同的模式数据的分辨率之间各不相同,为了方面后续分析,这里对于数据进行插值处理。而每个文件单独进行插值,无疑会浪费时间。因此,通过shell脚本结合cdo 插值命令可以实现高效、快速的多nc文件的插值。
运行下面的脚本之前,需要保证你已经成功安装了cdo软件。这里插值的分辨率为2.5x2.5°
#!/bin/bash
cd /Users/CMIP6_His_nc/
output_folder="/Users/CMIP6_His_nc/interp_2.5/"
if [ ! -d "$output_folder" ]; then
mkdir -p "$output_folder"
fi
for file in *1997-2014.nc
do
filename=$(basename -- "$file")
filename_no_ext="${filename%.*}"
echo "Interpolating $filename"
cdo remapcon,r144x73 "$file" "${output_folder}${filename_no_ext}_interp_2.5x2.5.nc"
done
echo "All files processed."
在这里,推荐使用remapcon 插值方法,特别是对于降水数据。
将以上脚本保存为sh文件,在Linux 命令段,输入 sh 文件名.sh
即可运行