[ncl-talk] drawing satellite passes(14 per day)

[Ipsita Putatunda]

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>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151118/0371d768/attachment.html