[ncl-talk] wind and emission flux on the same plot
Adam Phillips
asphilli at ucar.edu
Fri Aug 30 10:58:12 MDT 2019
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>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190830/77339602/attachment.html>
More information about the ncl-talk
mailing list