[ncl-talk] wind and emission flux on the same plot

Setareh Rahimi setareh.rahimi at gmail.com
Fri Aug 30 11:07:54 MDT 2019 

Dear Adam,
I did not get any of those errors you have mentioned. The script returns an
image, showing dust emission and boundary of the study area, but not
showing wind vectors, which means the issue is something else!
Thank you,

On Fri, Aug 30, 2019 at 9:28 PM Adam Phillips <asphilli at ucar.edu> wrote:

> Hi Setareh,
> I'm not sure what script you are referring to anymore, but I just ran your
> wind.ncl script (I do not have the data files to run the wind-dust2.ncl
> script), and got the following error messages:
>  (0) check_for_y_lat_coord: Warning: Data either does not contain
> (0) a valid latitude coordinate array or doesn't contain one at all.
> (0) A valid latitude coordinate array should have a 'units'
> (0) attribute equal to one of the following values:
> (0)    'degrees_north' 'degrees-north' 'degree_north' 'degrees north'
> 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
> (0) check_for_lon_coord: Warning: Data either does not contain
> (0) a valid longitude coordinate array or doesn't contain one at all.
> (0) A valid longitude coordinate array should have a 'units'
> (0) attribute equal to one of the following values:
> (0)    'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
> 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'
> warning:cnLevels is not a valid resource in rcm_vector at this time
> warning:cnFillOn is not a valid resource in rcm_vector at this time
> warning:cnFillMode is not a valid resource in rcm_vector at this time
> warning:cnLineLabelsOn is not a valid resource in rcm_vector at this time
> warning:cnLevelSelectionMode is not a valid resource in rcm_vector at this
> time
> warning:cnLinesOn is not a valid resource in rcm_vector at this time
> warning:cnFillPalette is not a valid resource in rcm_vector at this time
> (0) gsn_csm_map_ce: Warning: you are calling one of the CE map functions,
> (0) but setting the map projection to something other than
> (0) 'CylindricalEquidistant'. You may get errors or unexpected results.
> warning:ContourPlotDraw: out of range coordinates encountered; standard
> line rendering method may be unreliable;
>  consider setting the resource trGridType to "TriangularMesh" if
> coordinates contain missing values
>
> Whenever you write to ncl-talk you should always include the error
> messages you are getting. You should also try to fix the errors before
> contacting ncl-talk. The "check_for_y_lat_coord message is a common one
> whose solution is easily found by searching for the error message:
> https://www.ncl.ucar.edu/Document/Graphics/error_msg.shtml
>
> The units are correctly set in your xlat and xlon arrays, so that is not
> the issue. However, while you are setting the sfXArray/sfYArray resources,
> you are not setting the vfXArray/vfYArray resources to set the coordinates
> for the vector array that you are passing to gsn_csm_vector_scalar_map. So,
> set this to fix the check_for_* errors:
>   res at vfXArray = xlon
>   res at vfYArray = xlat
>
> The set of warning:cn* messages can be ignored as you are calling
> gsn_csm_vector_scalar_map with a common set of resources for both contour
> and vectors.
>
> The gsn_csm_map_ce error message was referenced by Dennis in a previous
> email, change gsn_csm_vector_scalar_map_ce to gsn_csm_vector_scalar_map.
>
> Finally, I noted in a previous email to remove the gsn_panel call. When
> overlaying one plot on the other, the generic format of a script looks like
> this:
>
> res = True
> res at gsnDraw = False
> res at gsnFrame = False
> ....
> plot = gsn_csm_contour_map(...
>
> res2 = True
> res2 at gsnDraw = False
> res2 at gsnFrame = False
> ....
> oplot = gsn_csm_contour(...
> overlay(plot,oplot)
> draw(plot)
> frame(plot)
>
> If you throw a gsn_panel call in between the plot = and the oplot = lines,
> you will not be able to overlay the oplot on the plot.
>
> Please apply the suggestions above before responding to ncl-talk, and as
> always, only respond to the ncl-talk email list.
> Adam
>
> On Fri, Aug 30, 2019 at 7:47 AM Setareh Rahimi via ncl-talk <
> ncl-talk at ucar.edu> wrote:
>
>> Dear Dennis,
>> Thanks for your help. I added shapefile to the script, but it does not
>> show wind vectors again. just showing emission flux and border of desire. I
>> attached the script. What should I do please? I have tried:
>> overlay(plotc,shape1), but faced an error.
>> Thanks in advance,
>> Best wishes,
>>
>>
>> On Wed, Aug 28, 2019 at 8:17 PM Dennis Shea <shea at ucar.edu> wrote:
>>
>>>  You used:
>>>   plot=*gsn_csm_vector_scalar_map_ce*(wks,u(ut,:,:),v(ut,:,:),
>>> x(nt,:,:),res)
>>>
>>> The *_ce *means 'cylindrical equidistant' ....
>>> You *must* have seen a Warning message! Did you try to address the
>>> issue?
>>> ---
>>> The attached draws a contour over a map; then overlays a vector plot.
>>>
>>>
>>>
>>>
>>>
>>> *https://www.ncl.ucar.edu/Applications/overlay.shtml*
>>> <https://www.ncl.ucar.edu/Applications/overlay.shtml>
>>>
>>> You used
>>>
>>> On Wed, Aug 28, 2019 at 12:28 AM Setareh Rahimi via ncl-talk <
>>> ncl-talk at ucar.edu> wrote:
>>>
>>>> Dear Adam,
>>>> You and Rick are right, but the problem was not fixed with removing
>>>> gsn_panel. On the other hand, I tried to add res at gsnDraw = False and
>>>> res at gsnFrame =  False to my first res resource list, but I saw no
>>>> change. I need to plot wind vectors over surface dust emission flux.
>>>> Thank you for your email.
>>>> Best wishes
>>>>
>>>> On Mon, Aug 26, 2019 at 9:54 PM Adam Phillips <asphilli at ucar.edu>
>>>> wrote:
>>>>
>>>>> Hi Setareh,
>>>>> Rick is correct that you should remove the gsn_panel call, as
>>>>> gsn_panel is not needed if you are creating a single plot. You will also
>>>>> need to add res at gsnDraw = False and res at gsnFrame =  False to your
>>>>> first res resource list, as you do not want to draw the plot nor advance
>>>>> the frame until you have added the wind vector overlay to the original
>>>>> plot. Those changes should fix the issue you are having.
>>>>> Adam
>>>>>
>>>>> On Fri, Aug 23, 2019 at 11:17 AM Setareh Rahimi via ncl-talk <
>>>>> ncl-talk at ucar.edu> wrote:
>>>>>
>>>>>> Dear Rick,
>>>>>> I did your suggestion, but nothing changed.
>>>>>> Thank you
>>>>>>
>>>>>> On Fri, Aug 23, 2019 at 8:14 PM Rick Brownrigg <brownrig at ucar.edu>
>>>>>> wrote:
>>>>>>
>>>>>>> Hi,
>>>>>>>
>>>>>>> I'm a little confused about the call the gsn_panel() to draw a
>>>>>>> "panel" plot of the emflx variable in a 1x1 panel, and then the subsequent
>>>>>>> call to overlay() overlaying the wind speed on top of emflx.  You might try
>>>>>>> removing the call the gsn_panel().
>>>>>>>
>>>>>>> Rick
>>>>>>>
>>>>>>> On Fri, Aug 23, 2019 at 9:21 AM Setareh Rahimi via ncl-talk <
>>>>>>> ncl-talk at ucar.edu> wrote:
>>>>>>>
>>>>>>>> Dear all,
>>>>>>>> I need to plot my desirable varible called "emflx", and wind
>>>>>>>> direction ( wind at 10m)  on one plot. Once I run the script, it just plot
>>>>>>>> the wind direction, not 'emflx' . Would you please kindy advise me in this
>>>>>>>> regards.
>>>>>>>> data and  the script, attached.
>>>>>>>> Many thanks in advanc,
>>>>>>>> Best wishes,
>>>>>>>>  Haboob-May-28-first_ATM.2018052400.nc
>>>>>>>> <https://drive.google.com/file/d/11LHXIs8w-nI--3cLwgGs9ulREVuBwX7h/view?usp=drive_web>
>>>>>>>>  Haboob-May-28-first_DUST04.2018052400.nc
>>>>>>>> <https://drive.google.com/file/d/1ICE51X0dN-tM4ffBwk6vo-uBUDR_ASLd/view?usp=drive_web>
>>>>>>>>  wind.ncl
>>>>>>>> <https://drive.google.com/file/d/1grYpZhKNcr1Flt0wF2xNJVbUK9b1dK2b/view?usp=drive_web>
>>>>>>>>
>>>>>>>>
>>>>>>>> --
>>>>>>>> S.Rahimi
>>>>>>>>
>>>>>>>> _______________________________________________
>>>>>>>> ncl-talk mailing list
>>>>>>>> ncl-talk at ucar.edu
>>>>>>>> List instructions, subscriber options, unsubscribe:
>>>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>>>
>>>>>>>
>>>>>>
>>>>>> --
>>>>>> S.Rahimi
>>>>>>
>>>>>> _______________________________________________
>>>>>> ncl-talk mailing list
>>>>>> ncl-talk at ucar.edu
>>>>>> List instructions, subscriber options, unsubscribe:
>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>
>>>>>
>>>>>
>>>>> --
>>>>> Adam Phillips
>>>>> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
>>>>> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>>>>>
>>>>> <http://www.cgd.ucar.edu/staff/asphilli>
>>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>> --
>> S.Rahimi
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
> --
> Adam Phillips
> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>
> <http://www.cgd.ucar.edu/staff/asphilli>
>


-- 
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190830/3c622b65/attachment.html>

More information about the ncl-talk mailing list