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

Adam Phillips asphilli at ucar.edu
Tue Nov 17 10:00:11 MST 2015


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/20151117/2cc0178d/attachment.html 


More information about the ncl-talk mailing list