(已解决)用grads打开两个nc文件做差值出错

东风急流

从NOAA上下载了nc格式的AVHRR的海表温度资料，想做两个时次的差值场，资料是全球的，分辨率为0.25*0.25，但在grads里面运行的时候报错：

直接做差值提示我时间和格点不匹配，于是想用define来存储第二个文件里面的海温，再用第一个sst去减，但依然提示我超出；我直接存储第一个文件的海温，用第二文件海温去减，再次提示超出：

这里是什么原因呢？两个文件时次不一样无法相减吗？如何用define实现呢？

东风急流

问题已经解决，赋上脚本：

'sdfopen D:\tu3\sst\amsr-avhrr-v2.20090716.nc'
'sdfopen D:\tu3\sst\amsr-avhrr-v2.20090718.nc'
'set dfile 1'
'set lon 105 130'
'set lat 10 35'
'set t 1'
'set lev 0'
'define a=sst'

'set dfile 2'
'set lon 105 130'
'set lat 10 35'
'set t 1'
'set lev 0'
'define b=sst.2'

'set gxout shaded'
'd b-a'
'cbarn'
'printim D:\tu3\sst\avhrr18-16.png white'

东风急流

两个文件经纬度一样的，为什么说X格点不匹配呢？

吴阿好大人

请问楼主有没有找到台风最大持续风速以及最大风圈半径资料？谢谢！！

东风急流

JTWC资料里面应该有

听海落雪

东风急流 发表于 2016-10-17 19:10
JTWC资料里面应该有

这样直接算好像不行，应该是用函数循环什么的。

东风急流

听海落雪 发表于 2016-10-17 20:55
这样直接算好像不行，应该是用函数循环什么的。

官网上有这两项：
VMAX - Maximum sustained wind speed in knots: 0 through 300.
MRD - radius of max winds, 0 - 999 nm.
我的理解是可以直接用的哈

liuliwa

楼主最后解决了吗？