Re: Problem Returning graphic object from user defined function

Dear Adam,
Thank you very much. I will working now, however, the final figure still not working for different issues. Anyway, I really appreciate your help.

On Sep 14, 2013, at 12:55 AM, Adam Phillips <asphilli@ucar.edu> wrote:

> Hi Mohamed,
> (Everyone: this initially went offline.)
>
> Thank you for sending on a simple example (offline). I think the issue is that you cannot pass the plot back outside your plotting function when you create the workstation within that function. If you want to bring the plot back outside the function the workstation must be created outside and passed into the function. Here are modified versions of call_plot.ncl and make_plot.ncl that show this process.
> Adam
>
> call_plot.ncl
> ;=======================================================================================
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "make_plot.ncl"
> ;=======================================================================================
> begin
> x = fspan(0,1,10)
> y = x^2
> wks = gsn_open_wks("png", "tmp_figure" )
> figure = make_plot(wks,x,y)
> draw(figure)
> frame(wks)
>
> draw(figure)
> frame(wks)
> end
> ;============
> make_plot.ncl
> ;=======================================================================================
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl"
> ;=======================================================================================
> undef("make_plot")
> function make_plot(wks:graphic, data_x:numeric,data_y:numeric)
> begin
>
> res = True
> res@gsnDraw = False
> res@gsnFrame = False
> plot=gsn_csm_xy (wks,data_x,data_y,res)
>
> print(plot)
> return(plot)
> end
>
> On 9/13/13 1:50 AM, Mohamed Abdel Kader wrote:
>> Dear All,
>> I would like to make a figure of different scatter plots super-imposed a global map. I have a user defined function that generates the scatter plots and I return the scatter plot graphic object. The returned graphic object has missing value.
>> I have to do in that way since the scatter plot is highly customized and it's being called many times from different scripts. So, I would like to keep one scatter plot script. Could you please help me on solving this problem.
>>
>> I attached parts from the script.
>> P.S.: The scatter function works fine as standalone
>>
>> Thanks a lot and let me know if you need more information.
>>
>> ========================================================================================================================
>> Scatter plot function:
>> function Plot_Scatter_function(x_title:string,y_title:string,data_x:numeric,data_y:numeric,f_name:string,data_set_string:string, ps_op:logical, plot:graphic)
>> local plot1, plot2
>> begin
>> .
>> .
>> .
>> . "many commands"
>> .
>> .
>> wks = gsn_open_wks("eps", f_name )
>>
>>
>> plot1 = gsn_csm_xy (wks,ndtooned(data_x),ndtooned(data_y),res)
>>
>> plot2 = gsn_csm_xy (wks,(/res@trXMinF,res@trXMaxF/),(/res@trYMinF,res@trYMaxF/),res2)
>>
>> overlay(plot1,plot2)
>> draw(plot1)
>> frame(wks)
>> print(plot1)
>> plot=plot1
>> return (plot)
>> end
>>
>> the print out from scatter function
>> Variable: plot
>> Type: obj
>> Total Size: 4 bytes
>> 1 values
>> Number of Dimensions: 1
>> Dimensions and sizes: [1]
>> Coordinates:
>> Number Of Attributes: 2
>> dataspec : /EMAC_map_scatter_____tmp_data xyDataSpecClass 89
>> data : /EMAC_map_scatter_____tmp_data coordArraysClass 67
>> (0) /EMAC_map_scatter_____tmp_xy xyPlotClass 88
>> ========================================================================================================================
>>
>> Parent script:
>> .
>> .
>> .
>> .
>> dum_poly = new(NREG,graphic)
>>
>> do n=0,NREG-1
>> dum_poly(n)= Plot_Scatter_function(x_title,y_title,data_xx,data_yy,f_name+"_____tmp",data_set_string, ps_op,dum_poly(n))
>> ann = gsn_add_annotation (map,dum_poly(n),amres)
>> end do
>> .
>> .
>> .
>> .
>> .
>> .
>> end
>>
>> the message error
>>
>> warning:NhlAddAnnotation: bad HLU id passed in, ignoring it
>> fatal:_NclAddHLURef: internal error: ncl object 647412112 does not exist for vname anno_id
>> fatal:_NclAddHLURef: internal error: ncl object 647412112 does not exist for vname ann
>>
>> The print out of dum_poly(n): from the Parent script.
>> Variable: dum_poly
>> Type: obj
>> Total Size: 4 bytes
>> 1 values
>> Number of Dimensions: 1
>> Dimensions and sizes: [1]
>> Coordinates:
>> Number Of Attributes: 3
>> data : -1
>> dataspec : -1
>> _FillValue : -1
>> (0) -1
>> ========================================================================================================================
>>
>>
>>
>>
>> Best Regards
>> Mohamed M. Abdel Kader
>> Atmosphere and Climate Modelling Group
>> Research Assistant/PhD Candidate
>> Environmental & Atmospheric Sciences
>> The Cyprus Institute, Cyprus
>> 20 Constantinou Kavafi Street
>> P.O. Box: 27456
>> Nicosia 2121, Cyprus
>> Tel. +357 222 08679
>> Fax. +357 222 08625
>> Mobile: +357 992 62108
>> e-mails: m.abdelkader@cyi.ac.cy
>> mmkader@gmail.com
>> Skype: mmkader
>>
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

Best Regards
Mohamed M. Abdel Kader
Atmosphere and Climate Modelling Group
Research Assistant/PhD Candidate
Environmental & Atmospheric Sciences
The Cyprus Institute, Cyprus
20 Constantinou Kavafi Street
P.O. Box: 27456
Nicosia 2121, Cyprus
Tel. +357 222 08679
Fax. +357 222 08625
Mobile: +357 992 62108
e-mails: m.abdelkader@cyi.ac.cy
mmkader@gmail.com
Skype: mmkader

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Sep 14 01:43:29 2013