安装ncl6.4.0遇到问题

小黄鹿

您好，求救一下，因为ncl更新了，所以也想换成6.4的版本，但是更新安装完后，原本6.3正常运行的程序在6.4下出现错误，给我的感觉是bashrc环境变量没有设置对，但检查后没有问题，不知道为什么，求大神解救


load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
minlon      = 96
maxlon      = 127
minlat      = 36
maxlat      = 54
;读取观测降水数据
diri = "/cygdrive/d/byteam/"
filepath = diri + "2016081808.024"
lines = asciiread(filepath,-1,"string")
; station=stringtofloat(str_get_field(lines(5:),1," "))
lon=stringtofloat(str_get_field(lines(5:),2," "))
lat=stringtofloat(str_get_field(lines(5:),3," "))
; height=stringtofloat(str_get_field(lines(5:),4," "))
obspr=stringtofloat(str_get_field(lines(5:),5," "))
;print(lines)
;print(station)
;print(lon)
;print(lat)
;print(height)
;print(obspr)
;========================================
;*****************************************
;graphic
;*****************************************
print("ok")
wks=gsn_open_wks("png","/cygdrive/d/byteam/1")
print("ok")
gsn_define_colormap(wks,"gui_default")
resl   = True
;resl@gsnMaximize=True
resl@gsnDraw            = False
resl@gsnFrame           = False         
resl@gsnCenterString  = ""
resl@gsnLeftString    = ""
resl@gsnRightString   = ""
resl@tiMainString     = ""
resl@tmYROn        = False        
resl@tmXTOn        = False      
resl@gsnAddCyclic         = False
;--------------------------------------
;=============================
;---------------------------------------
map_res                   = True
map_res@gsnFrame          = False
map_res@gsnDraw           = False
map_res@tmYROn        = False        
map_res@tmXTOn        = False      
;map_res@gsnMaximize=True
map_res@pmTickMarkDisplayMode = "Always"
map_res@mpCenterLonF      = 180.0
map_res@mpMinLonF         = minlon   
map_res@mpMaxLonF         = maxlon
map_res@mpMinLatF            = minlat
map_res@mpMaxLatF            = maxlat
map_res@gsnMajorLonSpacing =5
map_res@gsnMajorLatSpacing =4
map_res@gsnMinorLonSpacing =1
map_res@gsnMinorLatSpacing =1
;resl@tmXBLabels = (/"96~S~o~N~E","100~S~o~N~E","104~S~o~N~E","108~S~o~N~E","112~S~o~N~E","116~S~o~N~E","120~S~o~N~E","124~S~o~N~E"/)
;resl@tmYLLabels = (/"36~S~o~N~N","38~S~o~N~N","40~S~o~N~N","42~S~o~N~N","44~S~o~N~N","46~S~o~N~N","48~S~o~N~N","50~S~o~N~N","52~S~o~N~N","54~S~o~N~N"/)
;resl@pmTickMarkDisplayMode = "Always"
map_res@gsnLeftString      = ""
map_res@gsnRightString     = ""
map_res@gsnCenterString    = ""
;map_res@gsnAddCyclic         = False
map_res@tmXBLabelFontHeightF=0.015
map_res@tmYLLabelFontHeightF=0.015
map_res@gsnLeftStringFontHeightF=0.02
map_res@gsnRightStringFontHeightF=0.02
map_res@gsnCenterStringFontHeightF=0.02                  
map_res@gsnLeftStringOrthogonalPosF=-0.01
map_res@gsnRightStringOrthogonalPosF=-0.01
map_res@gsnCenterStringOrthogonalPosF=-0.01
map_res@mpLandFillColor      = "white"
map_res@mpDataSetName         = "Earth..4"   
map_res@mpDataBaseVersion     = "MediumRes"  
map_res@mpGeophysicalLineThicknessF= 1.      
map_res@mpNationalLineThicknessF= 1.         
map_res@mpMaskAreaSpecifiers    = (/"China","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/)
map_res@mpOutlineSpecifiers=(/"China:states","Nei Mongol","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/)  
;map_res@mpOutlineSpecifiers=(/"China"/)  
map_res@mpOutlineOn                 = True
;resl@mpOutlineBoundarySets       = "AllBoundaries"
map_res@mpGeophysicalLineColor = "black"
map_res@mpGeophysicalLineThicknessF = 0.5
map_res@mpPerimLineColor       = "black"
map_res@mpPerimLineThicknessF       = 0.5
map_res@mpNationalLineColor    = "black"
map_res@mpNationalLineThicknessF    = 0.5
map_res@mpProvincialLineColor="black"   
map_res@mpProvincialLineThicknessF =0.5  
;map_res@mpOutlineDrawOrder             = "PreDraw"
map_res@gsnCenterString = "2016_08_18_08_24h"
map = gsn_csm_map(wks,map_res)
;---------------------------------------------
;------------------------------------  
mapdir="/cygdrive/d/work/map/"
  ;---加上市级区域地图   
  shpfn3=mapdir+"nmg.shp"
  gres3        = True  
  gres3@minlat = minlat  
  gres3@maxlat = maxlat  
  gres3@minlon = minlon  
  gres3@maxlon = maxlon  
  gres3@gsLineColor   ="black"     
  gres3@gsLineThicknessF =1.0
  shp_plot3     = gsn_add_shapefile_polylines(wks,map,shpfn3,gres3)  
  ;---加上河流地图   
  shpfn4=mapdir+"river.shp"  
  gres4        = True  
  gres4@minlat = minlat  
  gres4@maxlat = maxlat  
  gres4@minlon = minlon  
  gres4@maxlon = maxlon   
  gres4@gsLineColor   ="blue"  
;  gres4@gsLineColor   ="gray"
  gres4@gsLineThicknessF =2.0
  shp_plot4    = gsn_add_shapefile_polylines(wks,map,shpfn4,gres4)  
;=======================================
;画站点信息
  colors=(/"white","DarkOliveGreen1","Chartreuse3","deepskyblue","blue","magenta","firebrick4"/)
  nstation=dimsizes(obspr)
  mkres = True
  dum_1 = new(nstation,graphic)
  dum_2 = new(nstation,graphic)
  dum_3 = new(nstation,graphic)
  dum_4 = new(nstation,graphic)
  dum_5 = new(nstation,graphic)
  dum_6 = new(nstation,graphic)
  mkres@gsMarkerThicknessF = 2.0      
  mkres@gsMarkerIndex = 16         
  mkres@gsMarkerSizeF = 0.007
  do i=0,nstation-1
    if(obspr(i).ge.0.1.and.obspr(i).lt.10.) then
    mkres@gsMarkerColor = colors(1)
    dum_1(i) = gsn_add_polymarker(wks,map,lon(i),lat(i),mkres)
    end if
    if(obspr(i).ge.10..and.obspr(i).lt.25.) then
    mkres@gsMarkerColor = colors(2)
    dum_2(i) = gsn_add_polymarker(wks,map,lon(i),lat(i),mkres)
    end if
    if(obspr(i).ge.25..and.obspr(i).lt.50.) then
    mkres@gsMarkerColor = colors(3)
    dum_3(i) = gsn_add_polymarker(wks,map,lon(i),lat(i),mkres)
    end if
    if(obspr(i).ge.50..and.obspr(i).lt.100.) then
    mkres@gsMarkerColor = colors(4)
    dum_4(i) = gsn_add_polymarker(wks,map,lon(i),lat(i),mkres)
    end if
    if(obspr(i).ge.100..and.obspr(i).lt.250.) then
    mkres@gsMarkerColor = colors(5)
    dum_5(i) = gsn_add_polymarker(wks,map,lon(i),lat(i),mkres)
    end if
    if(obspr(i).ge.250.) then
    mkres@gsMarkerColor = colors(6)
    dum_6(i) = gsn_add_polymarker(wks,map,lon(i),lat(i),mkres)
    end if
  end do
mkres@gsMarkerSizeF = 0.007
txres               = True                           
txres@txFontHeightF = 0.01
leglat=fspan(45.,51.,6)
leglon=(/100., 100.3, 100.3, 100.3, 100.6, 100./)
text=(/"0.1-9.9 mm","10.0-24.9 mm","25.0-49.9 mm","50.0-99.9 mm","100.0-249.9 mm","~F34~3 ~F21~250.0 mm"/)
  dum1 = new(6,graphic)
  dum2 = new(6,graphic)
  do i=0,5
  mkres@gsMarkerColor = colors(i+1)
  dum1(i) = gsn_add_polymarker(wks,map,97.5,leglat(5-i),mkres)
  dum2(i) = gsn_add_text(wks,map,text(i),leglon(i),leglat(5-i),txres)
  end do
draw(wks)
frame(wks)
delete(wks)

end

小黄鹿

尤其是我试运行xy04n.ncl的时候还是正确的出图了

Lancelot

应该是X11没配置好

blessing

NCL6.4和当前的cygwin不相配  使用systemfunc这个函数出的错   gsn_code.ncl 3482行用了这个函数   我也没有合适的cygwin版本

小黄鹿

blessing 发表于 2017-8-11 19:17
NCL6.4和当前的cygwin不相配  使用systemfunc这个函数出的错   gsn_code.ncl 3482行用了这个函数   我也没 ...

那就是需要更新版本的Cygwin了 好不容易有热情更新6.4 结果还不行

小黄鹿

Lancelot 发表于 2017-8-11 19:01
应该是X11没配置好

为啥呢？？？？？

西风带

你好，看到你发的帖子“安装ncl6.4.0遇到问题”，不知道你更新cygwin后是否解决了问题。我这两天安装最新的cygwin和ncl6.4.0，也遇到类似的问题，运行xy04n.ncl没有问题，运行下载的era40_1.ncl就出现问题，也是“ 0 [main] ncl 7412 child_info_fork::abort: C:\cygwin\bin\cygkrb5-3.dll: Loaded to different address: parent(0x210000) != child(0x190000)”，能否予以指导，谢谢！