NCL新手遇到的问题，救命

黄小仙儿

刚刚清风大哥说我提问区不对，帖子被删了。。。改了版好不习惯
Huanglei@Huanglei-PC:/home/Huanglei>ncl /home/Huanglei/wrf_CrossSection1.ncl
/home/Huanglei/bin/ncl.exe: error while loading shared libraries: ?: cannot open shared object file: No such file or directory

我在运行官网那些简单的例子gsun..等是还可以运行并出图，现在参考官网上WRF-NCL自己画就出现了这个问题我的NCL是在cygwin下安装的
我查了ncl。exe。还有加载的两个数据库都有，但是就是出现这个问题。
曲折的学习之路，求路过的大神看看是什么原因


脚本：
就是问号。我刚刚查了，大家遇到的问题都是libc.so文件找不到我的就是问号
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

begin
;
; The WRF ARW input file.  
; This needs to have a ".nc" appended, so just do it.
  a = addfile("/home/Huanglei/data/d03"+".nc","r")


; We generate plots, but what kind do we prefer?
  type = "pdf"
; type = "pdf"
; type = "ps"
; type = "ncgm"
  wks = gsn_open_wks(type,"plt_CrossSection1")


; Set some basic resources
  res = True
  res@MainTitle = "REAL-TIME WRF"
  res@Footer = False
  
  pltres = True


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

  FirstTime = True
  times  = wrf_user_getvar(a,"times",-1) ; get times in the file
  ntimes = dimsizes(times)          ; number of times in the file

  mdims = getfilevardimsizes(a,"P") ; get some dimension sizes for the file
  nd = dimsizes(mdims)

;---------------------------------------------------------------

  do it = 0,ntimes-1,2             ; TIME LOOP

    print("Working on time: " + times(it) )
    res@TimeLabel = times(it)   ; Set Valid time to use on plots

    tc  = wrf_user_getvar(a,"tc",it)     ; T in C
    rh = wrf_user_getvar(a,"QICE",it)      ; relative humidity
    z   = wrf_user_getvar(a, "z",it)     ; grid point height

    if ( FirstTime ) then                ; get height info for labels
      zmin = 0.
      zmax = max(z)/1000.
      nz   = floattoint(zmax/2 + 1)
      FirstTime = False
    end if

;---------------------------------------------------------------

          angle = 0.
   
        plane = new(2,float)
        plane = (/ mdims(nd-1)/2, mdims(nd-2)/2 /)    ; pivot point is center of domain (x,y)
        opts = False  

        rh_plane = wrf_user_intrp3d(rh,z,"v",plane,angle,opts)
        tc_plane = wrf_user_intrp3d(tc,z,"v",plane,angle,opts)

        dim = dimsizes(rh_plane)                      ; Find the data span - for use in labels
        zspan = dim(0)

        
      ; Options for XY Plots
        opts_xy                         = res
        opts_xy@tiYAxisString           = "Height (km)"
        opts_xy@cnMissingValPerimOn     = True
        opts_xy@cnMissingValFillColor   = 0
        opts_xy@cnMissingValFillPattern = 11
        opts_xy@tmYLMode                = "Explicit"
        opts_xy@tmYLValues              = fspan(0,zspan,nz)                    ; Create tick marks
        opts_xy@tmYLLabels              = sprintf("%.1f",fspan(zmin,zmax,nz))  ; Create labels   
        opts_xy@tiXAxisFontHeightF      = 0.020
        opts_xy@tiYAxisFontHeightF      = 0.020
        opts_xy@tmXBMajorLengthF        = 0.02
        opts_xy@tmYLMajorLengthF        = 0.02
        opts_xy@tmYLLabelFontHeightF    = 0.015
        opts_xy@PlotOrientation         = tc_plane@Orientation


      ; Plotting options for RH
        opts_rh = opts_xy
        opts_rh@pmLabelBarOrthogonalPosF = -0.07
        opts_rh@cnFillOn                = True
        
      ; Plotting options for Temperature
        opts_tc = opts_xy
        opts_tc@cnInfoLabelOrthogonalPosF = 0.00
        opts_tc@ContourParameters  = (/ 5. /)


      ; Get the contour info for the rh and temp
        contour_tc = wrf_contour(a,wks,tc_plane,opts_tc)
        contour_rh = wrf_contour(a,wks,rh_plane,opts_rh)


      ; MAKE PLOTS         
        plot = wrf_overlays(a,wks,(/contour_rh,contour_tc/),pltres)

      ; Delete options and fields, so we don't have carry over
        delete(opts_tc)
        delete(opts_rh)
        delete(tc_plane)
        delete(rh_plane)

   
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

  end do        ; END OF TIME LOOP

end

andrewsoong

你这是画什么图？变量~~~

黄小仙儿

andrewsoong 发表于 2014-4-8 16:40
你这是画什么图？变量~~~

画剖面图，变量是Qice，再加上温度垂直分布

longlivehj

o(∩∩)o...哈哈，我说咋帖子瞬间就木有了。提问有专区了，也好！不过，还要适应一下。

你cd进入ncl可执行程序所在的bin目录，然后运行ldd ncl，看看输出里面有没有类似“not found”之类的，会提示你缺少什么库。

黄小仙儿

andrewsoong 发表于 2014-4-8 16:40
你这是画什么图？变量~~~

类似于这种

黄小仙儿

longlivehj 发表于 2014-4-8 16:42
o(∩∩)o...哈哈，我说咋帖子瞬间就木有了。提问有专区了，也好！不过，还要适应一下。

你cd进入ncl可执 ...

貌似没有啊

longlivehj

黄小仙儿 发表于 2014-4-8 16:46
貌似没有啊

给你发个链接，你用里面提到的一个cygcheck命令试试先。我没有cygwin，所以……
cygcheck

黄小仙儿

longlivehj 发表于 2014-4-8 16:56
给你发个链接，你用里面提到的一个cygcheck命令试试先。我没有cygwin，所以……
cygcheck

找到了
cygcheck: track_down:could not find cygheimntlm-0.dll
那接下来我该怎么办？在网上找到这个文件下载后，放到bin文件夹里吗？

andrewsoong

黄小仙儿 发表于 2014-4-8 16:43
类似于这种

没问题啊~~~你的有问题，如果系统一样，我可以拷贝一个你需要的库给你。我的是
Centos 6.4 64位的。

longlivehj

黄小仙儿 发表于 2014-4-8 17:03
找到了
cygcheck: track_down:could not find cygheimntlm-0.dll
那接下来我该怎么办？在网上找到这个 ...

Good！对，下载了放在bin里面！