用linint_warp插值后，海温和地图不重合，求大神指教

哈特虎tiger

脚本如下

f=addfile("F:。。。.nc", "r")
p=f->SST
printVarSummary(p)

a=addfile("F:/。。.nc", "r")
lon_new=a->lon
lat_new=a->lat
printVarSummary(lat_new)
printVarSummary(lon_new)

lat_old=fspan(-90,90,384)
lon_old=fspan(0,358.75,320)

lat_lon = linint2_Wrap(lon_old, lat_old, p, True, lon_new, lat_new, 0)
lat_lon&lat@units = "degrees_north"
lat_lon&lon@units = "degrees_east"
printVarSummary(lat_lon)

wks = gsn_open_wks("png","F:/test")               
gsn_define_colormap(wks, "GMT_drywet")

res = True                                    
res@gsnDraw      =False
res@gsnFrame     =False
res@gsnAddCyclic =False
res@gsnMaximize  =True

res@cnFillOn             = True        
res@cnLinesOn            = True        
res@cnLineLabelsOn       = True        
res@cnLineThicknessF = 0.5
res@cnLineColor = "black"
res@cnLevelSelectionMode = "ManualLevels"   

res@lbLabelBarOn         = True      
res@gsnPanelLabelBar     = False
res@lbBoxLinesOn         = True

res@cnMinLevelValF       = 0               
res@cnMaxLevelValF       =  28               
res@cnLevelSpacingF      =   1              
res@cnLineLabelDensityF = 3
res@cnLineLabelInterval = 1
res@cnLineLabelFont = 0.15

res@gsnContourNegLineDashPattern =6
res@gsnContourPosLineDashPattern =0
res@gsnContourZeroLineThicknessF = 0

res@mpFillOn  = True
res@mpMinLatF =-90
res@mpMaxLatF =90
res@mpMinLonF =0
res@mpMaxLonF =360
res@mpCenterLonF = 180
plot = gsn_csm_contour_map(wks,lat_lon(0,0,:,:), res)  ; create plot

draw(plot)
frame(wks)

stepdance

我猜是因为linint_warp函数插值计算时，只要原数据有FillVallue或者missing data的点，新插值出来的结果就是FillVallue，所以你的图上，空白大陆的区域比实际的大陆区域要大，而且偏东和偏南，我猜是因为你的lon_new和lat_new排列顺序的问题

哈特虎tiger

现在解决了，是我用的函数不对，不过也是感谢回复