[ncl-talk] wind and emission flux on the same plot
Setareh Rahimi
setareh.rahimi at gmail.com
Sat Sep 21 11:30:24 MDT 2019
Dear Dennis,
Thank you for your help. But it does not show color bar from zero.
Best wishes,
On Sun, Sep 15, 2019 at 6:20 PM Dennis Shea <shea at ucar.edu> wrote:
> The easiest way:
>
> xx = x(nt,:,:)
> xx = where(xx.le.0 , xx at _FillValue, xx)
> plot=gsn_csm_vector_scalar_map_ce(wks,u(ut,:,:),v(ut,:,:), xx,res)
>
> ===
>
> On Sat, Sep 14, 2019 at 11:08 AM Setareh Rahimi via ncl-talk <
> ncl-talk at ucar.edu> wrote:
>
>> Dear NCL users,
>>
>> I have a question regarding the color bar, which is:
>> once I set res at cnLevels to start from zero ( as bellow), NCL indicates
>> zero values as very light blue color on the map. However, I would like that
>> zero values will be shown as white color on the map. So how can I do so,
>> please?
>> Look forward to your helpful suggestions.
>> Best wishes,
>>
>>
>>
>> res at cnLevels = (/ 0, 100, 250, 500, 750, 1000, 2000\
>> , 3000, 4000, 5000, 6000, 7000, 8000\
>> , 9000,10000,11000,12000,13000,14000\
>> , 15000,17250,20000,25000,30000,40000\
>> , 50000,60000,75000,100000 /)
>>
>> 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
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
--
S.Rahimi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190921/41fd6a06/attachment.html>
More information about the ncl-talk
mailing list