cmip5 tos 插值

天玖Y

用ncl对cmip5 tos资料进行插值后，画出的图不对哇~~求大神帮帮忙呀~~

以下是ncl插值的脚本：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
diri= "/cygdrive/j/historicalDLN/tos/"
diri2="/cygdrive/d/OBS/nino34/"
in2  =addfile(diri2+"HadISST-sst.nc","r")  
lon2 = in2->longitude       ；（-179.5-179.5）1.0
lat2 = in2->latitude      ；（-89.589.5）  1.0
;printVarSummary(lon2)
;print(lat2)
ff=addfile(diri+"tos_Omon_CCSM4_historical_r1i1p1_185001-200512.nc","r")
names=getfilevarnames(ff)
print(names)
ts=lonFlip(ff->tos(1548:1871,:,:) )   ;197901-200512
lon1   = ts&i
lat1   = ts&jprintVarSummary(ts)
;print(lon1)
print(lat1)
fout = addfile("CCSM44.nc","c")
fol = linint2_Wrap (lon1,lat1,ts,True,lon2,lat2(::-1),0)
fout->tos=folend
插值后的图形为：[img=0,1]file:///C:\Users\hp\Documents\Tencent Files\364768116\Image\C2C\S{QZB~77RF_S35){[]Y0R$4.jpg[/img]

天玖Y

Hurricane_Hu 发表于 2015-3-9 14:28
ncdump -v看一下？我是弄过tas的，tos是只有海洋？如果是的话，看一下陆地的格点是不是缺测或者缺失值？
...

hp@hp-PC /cygdrive/j/historicalDLN/tos
$ ncdump -v tos_Omon_IPSL-CM5A-MR_historical_r1i1p1_185001-200512.nc
-bash: ncdump: command not found


怎么用不起来呢

tos就只是海洋，陆地没有数据的

Hurricane_Hu

lat1 lon1的范围是多少？
看上去插值只插到了北半球的样子。。。

天玖Y

Hurricane_Hu 发表于 2015-3-9 13:30
lat1 lon1的范围是多少？
看上去插值只插到了北半球的样子。。。

lat1 和lon1 只说了格点数320 和384 ；

还有就是北半球陆地多 海温资料分布也不对呀  是不是纬度反了~

Hurricane_Hu

天玖Y 发表于 2015-3-9 13:54
lat1 和lon1 只说了格点数320 和384 ；

还有就是北半球陆地多 海温资料分布也不对呀  是不是纬度反了~ ...

ncdump -v看一下？我是弄过tas的，tos是只有海洋？如果是的话，看一下陆地的格点是不是缺测或者缺失值？

好猫

以前用MATLAB插值过tos，对NCL的操作不是很清楚，但是需要注意一些问题：
1、tos资料多为非均匀格点资料，即在低纬精度较高。
2、要搞清原始数据的存放顺序，维度的变化，比如有些资料的起始经度为0，但有些为-270，不弄清楚的话，在插值时容易出现小失误。
希望能帮到你。

youko

进来学习一下！希望楼主早点解决问题~~

天玖Y

youko 发表于 2015-3-9 14:43
进来学习一下！希望楼主早点解决问题~~

这位小哥 心真善~

杨晓月

我也很期待有人解决楼主的问题

天玖Y

杨晓月 发表于 2015-3-9 14:45
我也很期待有人解决楼主的问题

这位小妹  心也善~