爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5129|回复: 0

画图时选兰伯特投影后经纬度范围设置出错

[复制链接]

新浪微博达人勋

发表于 2015-1-17 21:25:14 | 显示全部楼层 |阅读模式
NCL
系统平台:
问题截图: -
问题概况: NCL画图时不设置投影方式经纬度范围没有问题,但设置为兰伯特投影后就有偏移
我看过提问的智慧: 看过
自己思考时长(天): 3

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

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

x
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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

begin
   data=fbindirread("t-mon.dat", -1, (/239,163,283/), "float")
   printVarSummary(data)
  lon = fspan(69.75,140.25,283)
  lat = fspan(14.75,55.25,163)
  time = fspan(1,239,239)
  ;time = tostring(time1)
  lat!0     = "lat"
  lon!0     = "lon"
  time!0    = "time"

  lat@units = "degrees_north"
  lon@units = "degrees_east"
  time@units= "days since 1998-04-01 00:00:00"
  time@long_name = "Time"
  lat&lat   = lat
  lon&lon   = lon
  time&time   = time


  data!1 = "lat"
  data!2 = "lon"
  data!0 = "time"
  data&time = time
  data&lat = lat
  data&lon = lon
  data@_FillValue = -9999
  data@description=""
  data@units = ""

  wks_type="eps"
  wks=gsn_open_wks(wks_type,"obs-t")

  sum_t=data(66,:,:)
  ;printVarSummary(sum_t)
  sum_t = sum_t / 30

  res= True

  res@gsnAddCyclic = False
  res@mpDataSetName         = "Earth..4"   ; This new database contains
                                           ; divisions for other countries.
  res@mpDataBaseVersion     = "MediumRes"  ; Medium resolution database
  res@mpOutlineOn           = True         ; Turn on map outlines
  res@mpOutlineSpecifiers   = (/"China:states","Taiwan"/)       ;China:states
  
  
  res@mpProjection = "LambertConformal"  ; Lambert
  res@mpLambertMeridianF = 104.06
  res@mpLimitMode = "LatLon"
  res@mpLambertParallel1F = 30      ;Default: .001
  res@mpLambertParallel2F = 60    ;Default: 89.999

  res@mpMinLatF             =  30          ; Asia limits
  res@mpMaxLatF             =  31.5
  res@mpMinLonF             =  102.45
  res@mpMaxLonF             =  104.45


res@cnFillOn=True
res@cnLinesOn=True
res@cnLineLabelsOn = False
res@cnFillDrawOrder= "PostDraw"     
res@cnInfoLabelOn=False

      res@cnLineColor = "Black"
      res@cnLineLabelBackgroundColor = -1
      res@gsnContourLineThicknessesScale = 2.0

  res@pmTickMarkDisplayMode = "Always"
  res@tmXTOn = False
  res@tmYROn = False

  res@cnLevelSelectionMode="ExplicitLevels"
  res@cnLevels = (/8,10,12,14,16,19,20,22,24,26,28/)  
  ;res@cnLevels = (/100,150,200,250/)

  res@lbLabelBarOn = True       ;LabelBar显示

    map = gsn_csm_contour_map(wks,sum_t,res)

end

上面是我的脚本,区域设置是30-31.5,102.45到104.45,不设置投影方式出来没问题,但是设为兰伯特投影后就有偏差了,请问各位,这该怎么解决呢?

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

本版积分规则

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

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

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