ncl绘图求教

叁玖壹

各位大佬帮帮忙，用FORTRAN将ENSO冷年的强度距平做合成分析和区域变化显著性检验后作图，在处理后的txt文件中按列为台站号、维度、经度、距平合成、显著性检验结果，利用NCL作图时显著性检验结果始终无法叠加在距平分布图上（斜线为通过显著性检验区域）。下面的txt文件为数据处理结果，冷年的显著性检验临界值为1.691
我的NCL绘图脚本：
begin
data = asciiread("g:/jg/htsj/Lqdjp.txt",(/663,5/),"float")
lon=data(:,2)
lat=data(:,1)
jp=data(:,3)
tp=data(:,4)
olon=new(127,"float")
olat=new(81,"float")
data_p=new((/81,127/),"float")
data_t=new((/81,127/),"float")
do i=0,63*2
olon(i)=72+i*0.5
end do
do l=0,40*2
olat(l)=17+l*0.5
end do
lon!0="lon"
lon@long_name="lonitude"
lon@units="degrees_east"
lon&lon=lon
lat!0="lat"
lat@long_name="latitude"
lat@units="degrees_north"
lat&lat=lat
olon!0="lon"
olon@long_name="lonitude"
olon@units="degrees_east"
olon&lon=olon
olat!0="lat"
olat@long_name="latitude"
olat@units="degrees_north"
olat&lat=olat
rscan=(/10,5,3/)
data_p=obj_anal_ic_deprecated(lon,lat,jp,olon,olat,rscan,False)
data_t=obj_anal_ic_deprecated(lon,lat,tp,olon,olat,rscan,False)
wks=gsn_open_wks("png","qdjpL")
res=True
res@gsnAddCyclic=False
res@gsnMaximize=True
res@mpDataBaseVersion="MediumRes"
res@mpDataSetName="Earth..4"
res@mpOutlineOn=True
res@mpFillOn = True
res@mpOutlineSpecifiers=(/"China:states","Taiwan"/)
res@mpGeophysicalLineThicknessF=2.0
res@mpNationalLineThicknessF=2.0
res@mpMinLatF=15.0
res@mpMaxLatF=55.0
res@mpMinLonF=105
res@mpMaxLonF=135.0
res@cnLinesOn=True   
res@mpFillDrawOrder   ="PreDraw"
res@cnFillDrawOrder   ="PreDraw"
res@cnFillOn                        = True
res@cnLineLabelsOn=True
res@gsnDraw = False
res@gsnFrame = False
plot=gsn_csm_contour_map(wks,data_p,res)
res2 = True                        
  res2@gsnDraw             = False     
  res2@gsnFrame            = False            
  res2@cnInfoLabelOn       = False      
  res2@cnLinesOn           = False      
  res2@cnLineLabelsOn      = False      
  plot2   = gsn_csm_contour(wks,data_t, res2)
opt = True
opt@gsnShadeFillType = "pattern"   
opt@gsnShadeMid =3           
plot2=gsn_contour_shade(plot2,1.691,2.6,opt)
overlay(plot,plot2)
draw(plot)
frame(wks)
end
绘图效果：

一大碗年糕

之前也是用过gsn_contour_shade这个函数不知道为啥检验出不来后来莫得办法只能自己手动contour用explicit模式了。。

叁玖壹

一大碗年糕 发表于 2020-4-11 15:55
之前也是用过gsn_contour_shade这个函数不知道为啥检验出不来后来莫得办法只能自己手动contour用explicit模 ...

之前可以的，我把数据重新整理了一下就弄不出来了

66kqkq

看着像是经纬度不对应吧