[ncl-talk] drawing satellite passes(14 per day)
Ipsita Putatunda
ipsita.putatunda at gmail.com
Thu Nov 19 04:02:53 MST 2015
Hi,
I have already attached the overlay plots which I tried earlier with the
previous scripts. The modified script is giving individual 2 plots (if am
using 2 files); the first 2 plots attached herewith (Mt2..).
I have tried to put a secondary resource like "res= True" and "res1=True"
and take "res" for the first plot and rest of the plots with "res1", now am
getting 2 plots, where 1st plot is with the first file (plotA say) and
second plot is a overlay plot (plotA, plotB) but for plotB it is coming in
black color.Please find the plots attached (last 2 plots Mt1...).
Any help in this issue wil be appreciated.
Thanks,
Ipsita
On Thu, Nov 19, 2015 at 12:33 AM, Adam Phillips <asphilli at ucar.edu> wrote:
> 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/20151119/282aa933/attachment-0001.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Mt2.000001.png
Type: image/png
Size: 153576 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151119/282aa933/attachment-0004.png
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Mt2.000002.png
Type: image/png
Size: 163783 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151119/282aa933/attachment-0005.png
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Mt1.000001.png
Type: image/png
Size: 153576 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151119/282aa933/attachment-0006.png
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Mt1.000002.png
Type: image/png
Size: 205321 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151119/282aa933/attachment-0007.png
More information about the ncl-talk
mailing list