爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5786|回复: 0

求助关于linint_Wrap 的问题,插值画后不出图,

[复制链接]

新浪微博达人勋

发表于 2012-7-22 11:47:24 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
求助关于linint_Wrap 的问题,插值画后不出图,我用 linint_Wrap 这个函数想把ERA-interim 1.5*1.5 的插值成 NCEP格式2.5*2.5,结果插值好了,运行也没提示有错误或警告,但就是画不出等值线,只有底图,而色标,单位都有的,print插值后的数也是有值的,就是画不出来,请各位大虾帮忙解决!
下边是脚本:
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"
  begin
    latS=0.
    latN=60.
    lonL=60.
    lonR=160.
     diri="/home/clear/ncl/mon/"
     h=addfile(diri+"slp.mon.mean.nc","r")
     latout=h->lat
     lonout=h->lon
     
     
     h1=addfile(diri+"slp7-8.nc","r")
     TIME=h1->time
     YYYY =ut_calendar(TIME,0)
     iYYYY =ind(YYYY(:,1).eq.7)
     iYYYYY =ind((YYYY(:,0).ge.1979.and.YYYY(:,0).le.2011).and.(YYYY(:,1).eq.8))
     h7=short2flt(h1->msl(iYYYY,:,)
     h8=short2flt(h1->msl(iYYYYY,:,)
     
         
     
   
     h7ave=dim_avg(h7(latitude|:,longitude|:,time|)
     h8ave=dim_avg(h8(latitude|:,longitude|:,time|)
   
     h2=(h7ave+h8ave)/200
     copy_VarCoords(h7(0,:,,h2)
     latin=h2&latitude
     lonin=h2&longitude
     var=h2
     h3=linint2_Wrap(lonin,latin,var,True,lonout,latout,0)
     
  
    wks=gsn_open_wks("x11","slpera")
    gsn_define_colormap(wks,"gsdtol")

   
   
    resCN                      =True
    resCN@cnFillOn             =False
    resCN@cnLinesOn            =True
    resCN@gsnAddCyclic         =False
    resCN@mpOutlineBoundarySets       = "National"
    resCN@mpGeophysicalLineThicknessF =1.
    resCN@mpNationalLineThicknessF    =1.
    resCN@gsnSpreadColors      =False
    resCN@mpMinLatF            = latS
    resCN@mpMaxLatF            = latN
    resCN@mpMinLonF            = lonL
    resCN@mpMaxLonF            = lonR
    resCN@mpFillOn             =False
    resCN@pmTickMarkDisplayMode ="Conditional"
    resCN@cnLineLabelInterval             = 1
   
    resCN@cnLevelSpacingF        = 2
    resCN@gsnContourLineThicknessesScale=2
    resCN@cnLineLabelPlacementMode        = "Constant"
    resCN@cnInfoLabelOn        = False         
   
    plot=gsn_csm_contour_map(wks,h3(:,,resCN)
   
   
   

    end
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表