[ncl-talk] drawing satellite passes(14 per day)
Adam Phillips
asphilli at ucar.edu
Wed Nov 18 12:03:21 MST 2015
Hi Ipsita,
First, a developer noted that when I stated earlier that cnRasterSmoothingOn
= True would slow things down that I was incorrect. That resource should
not slow things down (that much).
I was thinking of cnSmoothingOn = True, which can slow things down.
So you are able to draw each swath individually (+correctly). I had you do
that to make sure that you could draw each individually before overlaying.
When you go to overlay them you state "it is not giving me the expected
plot with multiple overlying orbital swath." Can you please send us a
picture showing the overlay plot and tell us what is incorrect about it?
As always respond to ncl-talk and not to me personally.
Adam
On Wed, Nov 18, 2015 at 1:15 AM, Ipsita Putatunda <
ipsita.putatunda at gmail.com> wrote:
> Thanks for your suggestions, I have modified my script as per Adam's
> suggestions though the script has speed up a bit, but it is giving two
> separate plots instead of overlaying plot. If I plot each file separately
> it is showing single orbital swath, but while overlaying it is not giving
> me the expected plot with multiple overlying orbital swath.
> Any help in this issue will be appreciated.
>
> Thanks,
> Ipsita
>
>
>
> On Tue, Nov 17, 2015 at 10:30 PM, Adam Phillips <asphilli at ucar.edu> wrote:
>
>> Hi Ipsita,
>> I do not see anything amiss with your script. However, I would definitely
>> comment out this line:
>> res at cnRasterSmoothingOn = True
>> As I believe that will really slow down your script.
>>
>> If that does not fix the issue I would suggest the following:
>> 1) set res at gsnDraw = True and res at gsnFrame = True.
>> 2) Comment out the overlay(PLOT,plot), draw(PLOT) and frame(wks) lines.
>> 3) Change the output format from pdf to png
>>
>> Do you get each swath drawn on its own frame correctly?
>>
>> As always, please respond to the ncl-talk email list.
>> Adam
>>
>>
>>
>> On Mon, Nov 16, 2015 at 2:57 AM, Ipsita Putatunda <
>> ipsita.putatunda at gmail.com> wrote:
>>
>>> Thanks a lot for your suggestions. I have modified my script as per your
>>> suggestions, but instead of overlaying plot am getting a single plot, and
>>> it is taking a lot of time for only two files (for testing purpose I am
>>> doing with 2 files out of 14) not able to understand where am I making the
>>> mistake. The plot and the modified script is attached herewith.
>>>
>>> Thanks,
>>> Ipsita
>>>
>>> On Sat, Nov 14, 2015 at 12:50 AM, Dennis Shea <shea at ucar.edu> wrote:
>>>
>>>> diri = "./" ; input directory
>>>> fili = systemfunc("cd "+diri+" ; ls MT1SAPS*h5") ; all desired
>>>> files
>>>> print(fili)
>>>> nfili = dimsizes(fili)
>>>>
>>>> res = True ; define once
>>>> res at gsnDraw = False
>>>> res at gsnFrame = False
>>>> res at gsnAddCyclic = False
>>>>
>>>> res1 = True
>>>> :
>>>> wks = gsn_open_wks ("PDF","xx1")
>>>>
>>>> do nf=0,nfili-1 ;
>>>> loop over each file
>>>> rh =
>>>> if (nf.eq.0) then
>>>> PLOT = gsn_csm_contour_map (wks, rh, res)
>>>> else
>>>> plot = gsn_csm_contour(wks, rh, res)
>>>> overlay(PLOT, plot)
>>>> end if
>>>> end do ; nf loop
>>>>
>>>> draw(PLOT)
>>>> frame(wks)
>>>>
>>>> ===============
>>>>
>>>> Why are you doing this
>>>>
>>>> delete(fn)
>>>> delete(f)
>>>> delete(vnlon)
>>>> delete(ulon)
>>>> delete(vnlat)
>>>> delete(ulat)
>>>> delete(vnrr)
>>>> delete(urr)
>>>> delete(rr)
>>>> delete(dims)
>>>> delete(nlon)
>>>> delete(nlat)
>>>>
>>>> If the variables are the same size, they will overwrite the previous
>>>> values.
>>>>
>>>> If they are not the same size:
>>>>
>>>> delete( [/ fn, f, vnlon, ulon,......... /] ) delete multiple
>>>> variable in one command
>>>>
>>>> or, use NCL's := reassignment syntax
>>>>
>>>> urr := f->$vnrr$ ; if urr exists; delete and write with new
>>>> values
>>>>
>>>>
>>>> http://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclVariables.shtml#Reassignment
>>>>
>>>> ====
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> On Thu, Nov 12, 2015 at 10:25 PM, Ipsita Putatunda
>>>> <ipsita.putatunda at gmail.com> wrote:
>>>> > Dear NCL users,
>>>> > I want to make a single plot for all satellite passes (15 swath
>>>> per
>>>> > day) data. the script and the plot am getting from this is as
>>>> attached. for
>>>> > testing purpose I made it with only two files, but still the plot is
>>>> not
>>>> > coming exactly what I want. Could you please help me to fix this
>>>> issue. Is
>>>> > there any simpler way to do this plot? Any help in this issue will be
>>>> > appreciated.
>>>> >
>>>> > Thanks in advance,
>>>> > Ipsita
>>>> >
>>>> > _______________________________________________
>>>> > ncl-talk mailing list
>>>> > ncl-talk at ucar.edu
>>>> > List instructions, subscriber options, unsubscribe:
>>>> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>> >
>>>>
>>>
>>>
>>> _______________________________________________
>>> 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>
>>
>
>
--
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/20151118/a8c265bf/attachment.html
More information about the ncl-talk
mailing list