[ncl-talk] wind and emission flux on the same plot
Adam Phillips
asphilli at ucar.edu
Fri Aug 30 11:45:48 MDT 2019
Hi Setareh,
If your script does not return any error messages, then either:
1 - You are using a different script than I did. If my suggestions did not
work please attach the latest version of your script (and you might want to
rename it to say wind-2.ncl) and attach the resulting plot.
2 - You have turned off or redirected all NCL error messages.
Regardless of which it is, the suggestions I included in my two previous
emails will fix the issue with the vectors not showing up.
Adam
On Fri, Aug 30, 2019 at 11:08 AM Setareh Rahimi <setareh.rahimi at gmail.com>
wrote:
> 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
>
>
--
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>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190830/61d664b6/attachment.html>
More information about the ncl-talk
mailing list