Datasets: EC-Earth historical (1850-2005) and scenario (up to 2100) runs for CMIP5/IPCC AR5.
Temporal resolution:  3 hours, chunked month by month. 
Objective: to merge the chunked files and make daily means for analyzing the land-atmosphere interactions.

First, build a job script by:
nedit  cat.CMIP5.data.daily.job

#!/bin/bash
# setup cdo path
cdo='/usr/local/free/bin/cdo'
# setup data path of historical or senarios
cd /.../CMIP5/historic
# setup variables to process
vars="CP LSP EVA T2M SWL1"  # Precipitation, evaporation, 2M temperature and top layer soil moisture.
for var in ${vars};
   do
     ${cdo}  -r -f nc copy ${var}_*.nc ${var}_historic_1850-2005.nc
     ${cdo}  daymean -setcalendar,standard ${var}_historic_1850-2005.nc ${var}_historic_1850-2005_daymean.nc
     # rm ${var}_18*.nc
     # rm ${var}_19*.nc
     # rm ${var}_20*.nc
   done
     ${cdo}  add CP_*_daymean.nc LSP_*_daymean.nc PRCP_*_daymean.nc
for var in ${vars};
   do
     ${cdo}  sellonlatbox,lon1,lon2,lat1,lat2  ${var}_*_daymean.nc  ${var}_*_daymean_US.nc 
     ${cdo}  sellonlatbox,lon1,lon2,lat1,lat2  ${var}_*_daymean.nc  ${var}_*_daymean_Europe.nc 
     ${cdo}  sellonlatbox,lon1,lon2,lat1,lat2  ${var}_*_daymean.nc  ${var}_*_daymean_Africa.nc 
   done

chmod +x cat.CMIP5.data.daily.job
nohup ./cat.CMIP5.data.daily > foo.out 2> foo.err < /dev/null &





转载本文请联系原作者获取授权，同时请注明本文来自王国杰科学网博客。
链接地址：https://blog.sciencenet.cn/blog-569118-504684.html

上一篇：收到一封邮件...
下一篇：科学家，还是作家？