[ncl-talk] wind and emission flux on the same plot
Adam Phillips
asphilli at ucar.edu
Fri Aug 30 12:10:15 MDT 2019
Hi Setareh,
Please include ncl-talk on all replies. We cannot run the wind-dust2.ncl
script as you never included the files you use at the top (which are
different than the ones used in wind.ncl). I substituted the files used in
your wind.ncl script into wind-dust2.ncl, and I get vectors on the
resulting image. I've attached the script and output plot. As one set of
files works and another doesn't I would look closely at the differences
between the files.
Adam
On Fri, Aug 30, 2019 at 11:58 AM Setareh Rahimi <setareh.rahimi at gmail.com>
wrote:
> Dear Adam,
> I already changed the script name as "wind-dust2.ncl", you can find it
> through previous emails. I also attached the result I got. As you can see,
> there are emission flux plus boundary of the study area, but no wind
> vectors.
> Thank you.
>
> On Fri, Aug 30, 2019 at 10:16 PM Adam Phillips <asphilli at ucar.edu> wrote:
>
>> 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>
>>
>
>
> --
> 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/dcd4a6e8/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wind-dust2.mod.ncl
Type: text/x-ncl
Size: 5159 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190830/dcd4a6e8/attachment-0001.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: rcm_2018122812.png
Type: image/png
Size: 404502 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190830/dcd4a6e8/attachment-0001.png>
More information about the ncl-talk
mailing list