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

Setareh Rahimi setareh.rahimi at gmail.com
Sat Aug 31 03:16:07 MDT 2019 

Dear Adam,
 It works now. Many thank for your help.
Best wishes,

On Sat, Aug 31, 2019 at 12:38 AM Adam Phillips <asphilli at ucar.edu> wrote:

> Hi Setareh,
> You need to set the vfXArray/vfYArray resources in your resv resource
> list. gsn_vector is not throwing an error message in this case, and I
> believe it should, as the u/v arrays have no coordinate variables attached
> to them. Once you set those resources you get the vectors to show up.
> Script and plot are attached.
> Adam
>
> On Fri, Aug 30, 2019 at 1:09 PM Setareh Rahimi <setareh.rahimi at gmail.com>
> wrote:
>
>>  Yazd-2-2019_DUST04.2018120100.nc
>> <https://drive.google.com/file/d/1ozyderWAFvQO4admiktJsALa2BQvhRpZ/view?usp=drive_web>
>>  Yazd-2-2019_SRF.2018120100.nc
>> <https://drive.google.com/file/d/1Xr_Omj6LYkwiDRLRuiOnYQStHpNzKE8N/view?usp=drive_web>
>> Dear Adam,
>> I attached the files that I am using to plot.  As long as there is no
>> shapefile added to the script, wind vectors appear. However, the problem
>> starts from where I add shapefile to the script.
>> I need to fix this problem.
>> Thank you,
>>
>> On Fri, Aug 30, 2019 at 10:40 PM Adam Phillips <asphilli at ucar.edu> wrote:
>>
>>> 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>
>>>
>>
>>
>> --
>> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190831/51db9ec1/attachment.html>

More information about the ncl-talk mailing list