fnl1×1格点资料画天气图

南宫蓝天

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
fnl_file=addfile("/gpfshome/fupj/wrf/data/nc/fnl_20130615_18_00.grib2.nc","r")
u=fnl_file->UGRD_P0_L100_GLL0(13,89:144,59:139)
v=fnl_file->VGRD_P0_L100_GLL0(13,89:144,59:139)
temp=fnl_file->TMP_P0_L100_GLL0(13,89:144,59:139)
temp=temp-273.15
ght=fnl_file->HGT_P0_L100_GLL0(13,89:144,59:139)
wks=gsn_open_wks("x11","PressureLevel")


res                             = True                 ; plot mods desired
res@mpFillOn                  = False             ; turn off map fill
res@mpOutlineSpecifiers      = (/"Continents"/)
  res@mpDataSetName         = "Earth..4"   ; This database contains
                                           ; divisions for other countries.
  res@mpDataBaseVersion     = "MediumRes"  ; Medium resolution database
  res@mpMinLatF = 0.
  res@mpMaxLatF = 55.
  res@mpMinLonF = 60.
  res@mpMaxLonF = 140.

  res@mpGeophysicalLineThicknessF= 1.      ; double the thickness of geophysical boundaries
  res@mpNationalLineThicknessF= 1.         ; double the thickness of national boundaries
  delete(res@mpOutlineSpecifiers)          ; about to change size of attribute, must delete
  res@mpOutlineSpecifiers = (/"China:states","Philippines","India","Nepal","Bhutan","Bangladesh","Myanmar", \
                                    "Thailand","Cambodia","Vietnam","Taiwan"/)
;map = gsn_csm_map_ce(wks,res)

;temp
opts=True
  opts@ContourParameters      = (/5.0/) ; set the contour levels
  opts@tiMainString = "T/U @500hPa"            ; set the main title
  opts@cnLineColor="Red"
plot_temp=gsn_contour(wks,temp,opts)


;ght
opts_z = res
opts_z@cnLineColor = "Blue"
opts_z@gsnContourLineThicknessesScale = 3.0
plot_ght= gsn_contour(wks,ght,opts_z)

;wind vector
;  sres                      = True              ; set up a second resource list
;  sres@gsnDraw              = False             ; do not draw the plot
;  sres@gsnFrame             = False             ; do not advance the frame
;  sres@cnLevelSelectionMode = "ExplicitLevels"  ; use explicit levels
;  sres@cnLevels             = ispan(-5,35,5 )   ; set the contour levels
  plot_v = gsn_vector(wks,u,v,False)

plot=wrf_map_overlays(wks,(/plot_temp,plot_ght,plot_v/),res)
end


想画天气图，不得解。跪求大神帮忙改下程序。

longlivehj

南宫蓝天 发表于 2015-4-26 09:36
求叠加函数，gsn_contour_map，gsn_csm_contour_map_overlay只找到这俩，用哪个？

看函数说明，根据需求选择。
两个contour，一个vector，用overlay。

cjjj

谢谢~谢谢~

兰北

提示 有什么问题？

南宫蓝天

兰北 发表于 2015-4-25 14:46
提示 有什么问题？

上面木有叠加到一张图上

南宫蓝天

cjjj 发表于 2015-4-25 14:11
谢谢~谢谢~

这是半成品，要改

cjjj

南宫蓝天 发表于 2015-4-25 14:48
这是半成品，要改

看错了。。。我还以为是贡献程序的~

longlivehj

fnl数据，却用的wrf的叠加函数！

南宫蓝天

longlivehj 发表于 2015-4-25 21:37
fnl数据，却用的wrf的叠加函数！

求叠加函数，gsn_contour_map，gsn_csm_contour_map_overlay只找到这俩，用哪个？

南宫蓝天

上面问题已解决