我有一系列的值,我希望在数据系列变为正数时找到第一个时间步的时间或日期。我知道我可以把它读成python,R或Fortran脚本来做它,但我想知道我是否可以从BASH的命令行中使用nco或cdo来做到这一点?
我想用
cdo gec,0.0 in.nc out.nc
制作面具,但这对我没什么帮助。我需要以某种方式根据数据的符号拆分文件,然后我可以使用
简单地选择日期cdo showdate
用管道输入awk。
此处有一个指向小示例文件的链接:
http://clima-dods.ictp.it/Users/tompkins/se/trmm_per10_pc0_year2000_nc2.nc
任何提示?
答案 0 :(得分:1)
我绞尽脑汁,想起NCO没有银弹。
我已经为ncap2提出了以下代码段。
(循环在ncap2中未经优化)
运行代码段使用命令
ncap2 -v -O -S sign.nco trmm_per10_pc0_year2000_nc2.nc foo.nc out.nc
/****************** sign.nco***********************************/
lon_sz=$lon.size;
lat_sz=$lat.size;
time_sz=$time.size;
*precip_prev=precip(0,0,0);
*precip_cur=2.0;
for(*idx=0;idx<time_sz;idx++)
{
for(*jdx=0;jdx<lat_sz;jdx++)
{
for(*kdx=0;kdx<lon_sz;kdx++)
{
precip_cur= precip(idx,jdx,kdx);
if( precip_cur > 0.0 && precip_prev<0.0)
print(time(idx));
precip_prev=precip_cur;
}
}
}
/***************************************************************/
答案 1 :(得分:0)
编辑回答:
使用已添加到CDO v1.8.0以后的新timcumsum函数(https://code.zmaw.de/projects/cdo/embedded/index.html#x1-3300002.8.1),现在可以为网格化字段完成此任务:
cdo gec,0 precip.nc mask.nc # positive entries are 1, negative are zero
cdo timcumsum mask.nc maskcum.nc # sum the mask in the time direction
# invert the mask so that you have 1 at start, but then zero from 1st +ve val
cdo lec,0.5 maskcum.nc maskinv.nc
# sum this new mask, and add 1, result is timestep when first positive value occurs:
cdo addc,1 -timcumsum maskinc.nc stepfirstpos.nc
我认为这个功能可以在一行中完成所有这些
cdo addc,1 -timcumsum -lec,0.5 -timcumsum -gec,0 precip.nc stepfirstpos.nc
原始回答:
花了我一年的时间,但我找到了一种方法,通过NCO和CDO的组合来实现单点文件但不是网格化的字段:
#!/bin/bash
# code to find date of first positive file entry
file=trmm_per10_pc0_year2000_nc2.nc
# set negative values to missing
cdo -s setrtomiss,-1.e36,0.0 $file set.nc
ntime=`cdo -s ntime set.nc` # number of steps in file
# loop over steps to find first positive value
for i in `seq 1 ${ntime}` ; do
# nco counts from step 0, cdo from 1
ncoi=`expr $i - 1`
# print out the value in step i
op=`ncks -d time,$ncoi -s "%16.10f\n" -H -C -v precip set.nc`
if [[ $op != "_" ]] ; then # not missing
# print the date of timestep i
cdo -s showdate -seltimestep,$i set.nc
rm -f set.nc
exit
fi
done
# no positive value found
rm -f set.nc
echo "All values in $file are negative"