haley, one more time help me.

Dear haley.
I'm korean. I don't know English.I hope you understand if I commit mistakes.
I don't understand your help. (ㅠ_ㅠ)
I call  "draw(plot)" after "overlay(plot,contour)" call. still draw each domain in two frame.
Hear is my code: 
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
  FILES   = systemfunc("ls ../geo_em*nc")  numFILES = dimsizes(FILES)
  type = "x11"  wks = gsn_open_wks(type,"geo_test")  gsn_define_colormap(wks,"gui_default")
  res = True  res_at_InitTime = False  res_at_Footer = False
  res_at_cnFillOn = True
  pltres = True  mpres = True  mpres_at_mpGeophysicalLineColor = "black"  mpres_at_mpGridLineColor = "black"
  lnres = True  lnres_at_gsLineThicknessF = 1.8  lnres_at_gsLineColor = "white"  lnres_at_gsLineLabelFontColor = "red"  lnres_at_gsLineLabelFontHeightF = 0.0145
    in_file = addfile(FILES(0) , "r")    hgt_m = wrf_user_getvar(in_file,"HGT_M",0)
    pltres_at_FramePlot = False    contour=wrf_contour(in_file,wks,hgt_m,res)    plot = wrf_map_overlays(in_file,wks,(/contour/),pltres,mpres)
  do ifil = 1, numFILES-1;    in_file = addfile(FILES(0) , "r");    hgt_m = wrf_user_getvar(in_file,"HGT_M",0);;    pltres_at_FramePlot = False;    contour=wrf_contour(in_file,wks,hgt_m,res);    plot = wrf_map_overlays(in_file,wks,(/contour/),pltres,mpres)
    in_file = addfile(FILES(ifil) , "r")
    if ( ifil .gt. 0 ) then
;      overlay(plot,contour);      draw(plot)      hgt = wrf_user_getvar(in_file,"HGT_M",0)
      lat = in_file->XLAT_M(0,:,:)      lon = in_file->XLONG_M(0,:,:)
      dims = dimsizes(hgt)      r2 = dimsizes(dims)      r21 = dims(r2-1)      r22 = dims(r2-2)
      xbox = (/lon(0,0), lon(0,r21-1), lon(r22-1, r21-1), \               lon(r22-1,0), lon(0,0)/)      ybox = (/lat(0,0), lat(0,r21-1), lat(r22-1, r21-1), \      y_out = new(dimsizes(ybox), typeof(ybox))
      if ( ifil .eq. 1 )        n_d0 = ifil + 1        lnres_at_gsLineLabelString = "d0"+n_d0      else        n_d0 = ifil + 1        lnres_at_gsLineLabelString = "d0"+n_d0      end if
      datatondc(plot, xbox, ybox, x_out, y_out)      gsn_polyline_ndc(wks, x_out, y_out, lnres)
      delete(hgt)      delete(lat)      delete(lon)    frame(wks)    end if  end do
end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Sep 03 2009 - 19:20:55 MDT