Re: haley, one more time help me.

Hi,

I don't understand what you are trying to plot.

You said you wanted to plot two contour domains on one map.

However, I'm not sure what two domains you want to plot. From the
script below it looks like you are trying to plot "HGT_M" twice
on the same plot, for different time steps.

Is this correct?

I tried to run your script below with some dummy data, but
you have commented some code out, and it doesn't run as expected.

Can you provide me with your script as you had it before,
and one or two of your geo*nc files?

You can email me offline about this.

--Mary

On Sep 3, 2009, at 7:20 PM, 감은아 wrote:

> 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

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