经纬度倒置

指甲钳

画500百帕的高度场，设置为兰勃脱投影后就变成这样了。是什么原因？

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;***************************************
begin
   grib_file="E:\20130727\fnl_20130728_00_00_c.grib2"
   fin = addfile(grib_file,"r")
   names = getfilevarnames(fin)        
;   print(names)

   atts = getfilevaratts(fin,names(0))
   dims = getfilevardims(fin,names(0))
;   print(atts)
;   print(dims)

; get vars

    hgt= fin->HGT_3_ISBL
    hgt=hgt/10
   lon= fin->lon_3
   lat= fin->lat_3
   lev=fin->lv_ISBL3

; create plot

    wks = gsn_open_wks("png","500")                 ; open wk station


  gsn_define_colormap(wks,"BlAqGrYeOrRe")        ; choose colormap

  res                       = True     ; plot mods desired

  res@cnFillOn              = False     ; turn on color fill
  res@cnLinesOn             = True    ; turn of contour lines
res@cnMonoLineColor  = False
;  res@cnLevelSpacingF       = 0.5      ; contour spacing
  res@gsnSpreadColors       = True     ; use full range of color map
  res@cnLevelSpacingF = 2
  res@lbLabelStride         = 4
res@cnMinLevelValF = 500
res@cnMaxLevelValF = 592
res@pmTickMarkDisplayMode = "Always"; use NCL default lat/lon labels

  res@gsnAddCyclic          = True    ; data already has cyclic point

  res@mpMinLatF            =  0      ; range to zoom in on
  res@mpMaxLatF            = 70
  res@mpMinLonF            = 40
  res@mpMaxLonF            = 160



;------------------------------map---------------------------


res@mpOutlineOn               = True              ; turn off outline
  res@mpFillOn                  = False             ; turn off map fill
res@mpDataSetName         = "Earth..4"   ; This database contains
                                           res@mpDataBaseVersion     = "MediumRes"  ; Medium resolution database
  res@mpGeophysicalLineThicknessF= 2.      ; double the thickness of geophysical boundaries
  res@mpNationalLineThicknessF= 2.         ; double the thickness of national boundaries
   res@mpOutlineSpecifiers = (/"China","India","Nepal","Bhutan","Bangladesh","Myanmar", \
                                    "Thailand","Cambodia","Vietnam","Taiwan"/)

res@mpLimitMode = "LatLon"
res@mpProjection = "LambertConformal"
res@mpLambertParallel1F = .001      ;Default: .001 ;可以自己改一改，看看投影有什么不同，挺有趣的
  res@mpLambertParallel2F = 89.999    ;Default: 89.999
;---------------------------------------------plot-------------------------------------------

   plot = gsn_csm_contour_map(wks,hgt(13,:,:), res)
;************************************************
end

kongfeng0824

感谢分享              

言深深

kongfeng0824 发表于 2013-9-10 16:35

                               
登录/注册后可看大图

感谢分享

亲，楼主这是求助帖哇···

尽头的尽头

兰伯特投影画图有两种形式，不知道楼主想画哪种
不知道是不是想要这样的效果

kongfeng0824

言深深 发表于 2013-9-10 17:02

                               
登录/注册后可看大图

亲，楼主这是求助帖哇···

是啊，看错帖子，以为是分享资源的呵呵不好意思

指甲钳

尽头的尽头 发表于 2013-9-10 17:41

                               
登录/注册后可看大图

兰伯特投影画图有两种形式，不知道楼主想画哪种
不知道是不是想要这样的效果

就是这种图，为什么我的这么丑？

尽头的尽头

指甲钳 发表于 2013-9-11 08:27

                               
登录/注册后可看大图

就是这种图，为什么我的这么丑？

加上这句 res@gsnMaskLambertConformal = True  

指甲钳

尽头的尽头 发表于 2013-9-11 09:09

                               
登录/注册后可看大图

加上这句 res@gsnMaskLambertConformal = True

十分感谢。 我刚用res@mpLambertMeridianF  = 115这个貌似也可以，官网上说是用来设置中轴线的意思？

Mr_Q

指甲钳 发表于 2013-9-11 09:29
十分感谢。 我刚用res@mpLambertMeridianF  = 115这个貌似也可以，官网上说是用来设置中轴线的意思？

我也遇到同样的问题，但是我设置了res@mpLambertMeridianF之后出来的图无法显示经纬度标注，请问您有遇到过这样的情况吗？

xyan88

学习了学习了