Re: Explicit tick marks and labels not plotting when using overlay

You are probably right. I didn't look at this too closely, but it sounds like you picked it up just find, and found "dim_sum_n". :-)

You don't you need this:

> area_car_dpvdyLPF = car_dpvdyLPF

unless you are trying to copy coordinates and attributes. And even then, it might be better to remove the above line, and then copy the metadata *after* the calculation with:

 copy_VarMeta(car_dpvdyLPF,area_car_dpvdyLPF)

--Mary

On Jun 23, 2010, at 6:46 PM, Kay Shelton wrote:

> Hi Mary,
>
> Thanks so much for your help, both with the plotting (works wonderfully!) and with cleaning up the script, it runs much faster now I have eliminated almost all the do loops.
>
> For completeness, the final suggestion you had using the where function didn't work. I could be wrong, but I think it is because I needed to sum over all lat/lons where the value was not missing and less than zero for each day/lev, and I was mixing the number of dimensions in parts of the where function. Anyway, I did still manage to still get around using a do loop here.
>
> I replaced this code:
>
> do days = 0,car_dims(0)-1
> do levs = 0,car_dims(1)-1
>
> do car_lat = 0,car_dims(2)-1
> do car_lon = 0,car_dims(3)-1
>
> if ((.not.(ismissing(car_dpvdyLPF(days,levs,car_lat,car_lon)))).and.(car_dpvdyLPF(days,levs,car_lat,car_lon).lt.0.)) then
> box_area(0,days,levs) = box_area(0,days,levs) + (car_dpvdyLPF(days,levs,car_lat,car_lon)/((car_dims(2)-1)*(car_dims(3)-1)))
> end if
>
> end do ; car_lon
> end do ; car_lat
>
> end do ; levs
> end do ;days
>
> with this code:
>
> area_car_dpvdyLPF = car_dpvdyLPF
>
> area_car_dpvdyLPF = where(((ismissing(car_dpvdyLPF)).or.(car_dpvdyLPF.ge.0.)),0.,area_car_dpvdyLPF)
>
> box_area(0,:,:) = (dim_sum_n(area_car_dpvdyLPF,(/2,3/)))/((car_dims(2)-1)*(car_dims(3)-1))
>
> This all ends up being significantly faster. Thanks again,
>
> Kay
>
>
> On 24 June 2010 07:26, Mary Haley <haley@ucar.edu> wrote:
> BTW, I think you can get rid of some of the do loops to make your script run faster. Do loops will cause any scripting language to slow down, so you want to keep them to a minimum.
>
> For example, the code:
>
> do nl = 0,(dimsizes(pvclim&lat)-1)
> dlon = (pvclim&lon(2)-pvclim&lon(1))*0.0174533
> dx = 6378388*cos(0.0174533*pvclim&lat(nl))*dlon
> dpvdx_T(nl,:,:,:)= dpvdx_T(nl,:,:,:)/dx
> end do
>
> doesn't need to have the "dlon" calculation inside the loop, since it's the same every time through.
>
> I believe you can accomplish the above without a do loop. Try this:
>
> dlon = (pvclim&lon(2)-pvclim&lon(1))*0.0174533
> dx = conform_dims(dimsizes(dpvdx_T),6378388*cos(0.0174533*pvclim&lat)*dlon,0)
> dpvdx_T= dpvdx_T/dx
>
> I'm not entirely sure about the above. It will only work if the leftmost dimension of dpvdx_T is the same size as pvclim&lat.
>
> This code:
>
> do days = 0,car_dims(0)-1
> do levs = 0,car_dims(1)-1
>
> box_min(0,days,levs) = min(car_dpvdyLPF(days,levs,:,:))
> box_min(1,days,levs) = min(epac_dpvdyLPF(days,levs,:,:))
>
> end do ; levs
> end do ; days
>
>
> can be replaced by using "dim_min_n":
>
> box_min(0,days,levs) = dim_min_n(car_dpvdyLPF(days,levs,:,:),(/2,3/))
> box_min(1,days,levs) = dim_min_n(epac_dpvdyLPF(days,levs,:,:),(/2,3/))
>
> This code:
>
> do car_lat = 0,car_dims(2)-1
> do car_lon = 0,car_dims(3)-1
>
> if ((.not.(ismissing(car_dpvdyLPF(days,levs,car_lat,car_lon)))).and.(ca\
> r_dpvdyLPF(days,levs,car_lat,car_lon).lt.0.)) then
> box_area(0,days,levs) = box_area(0,days,levs) + (car_dpvdyLPF(days,le\
> vs,car_lat,car_lon)/((car_dims(2)-1)*(car_dims(3)-1)))
> end if
>
> end do ; car_lon
> end do ; car_lat
>
>
> can be replaced by using the "where" function:
>
> box_area(0,days,levs) = where((.not.(ismissing(car_dpvdyLPF(days,levs,:,:)))).and.(car_dpvdyLPF(days,levs,:,:).lt.0.), \
> box_area(0,days,levs) + (car_dpvdyLPF(days,levs,:,:)/((car_dims(2)-1)*(car_dims(3)-1)))
> box_area(0,days,levs))
>
> I haven't tested all this, so there may be some syntax errors.
>
> --Mary
>
>
>
> On Jun 22, 2010, at 11:26 PM, Kay Shelton wrote:
>
> > Hello,
> >
> > I am trying to overlay one xy plot onto another using the overlay function. Both plots have explicit tick marks and labels on the x-axis. When I overlay the plots, the tick marks and labels from the base plot show up but the ones from the transform plot do not. The data from the transform plot is plotted though.
> >
> > I have included a plot of what I am trying to achieve (ncl-help_draw). This was plotted without using overlay, but since I later wish to panel several plots I cannot call draw. I have also included the result from using the overlay function (ncl-help_overlay) and also what the individual plots look like (ncl-help_separate). The script used to produce ncl-help_overlay is also attached.
> >
> > Thank you for your help,
> >
> > Kay Shelton
> > <ncl-help_draw.pdf><ncl-help_overlay.pdf><ncl-help_separate.pdf><ncl-help.ncl>_______________________________________________
> > 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 Jun 24 14:09:26 2010