[ncl-talk] drawing satellite passes(14 per day)
Adam Phillips
asphilli at ucar.edu
Thu Nov 19 14:31:04 MST 2015
Hi Ipsita,
It is very strange that you can draw each swath individually but when
overlaid one swath goes black. Can you upload the two data files to our ftp
site so one of us can take a look? Please follow the directions here:
http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP
Note that for security reasons, you cannot list the contents of this
directory, and neither can we. We need to know the *exact *names of the
files in order to retrieve them.
Let me know when the files have been uploaded.
Thanks,
Adam
On Thu, Nov 19, 2015 at 4:02 AM, Ipsita Putatunda <
ipsita.putatunda at gmail.com> wrote:
> 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>
>>
>
>
--
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/2008e236/attachment.html
More information about the ncl-talk
mailing list