[ncl-talk] drawing satellite passes(14 per day)
Krishna C
chandrakrishna.90 at gmail.com
Thu Nov 19 04:29:29 MST 2015
Hi
I also had the same problem as you are facing. I found a binning example in
ncl page
I have plotted BTs for some 3 overpass of MT for my domain of concern.
Please take a look if it is of any help
Cheers
Regards
-Krishna-
On Thu, Nov 19, 2015 at 4:32 PM, 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>
>>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151119/862310dd/attachment-0001.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: s6.pdf
Type: application/pdf
Size: 1012990 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151119/862310dd/attachment-0001.pdf
-------------- next part --------------
A non-text attachment was scrubbed...
Name: h5new.ncl
Type: application/octet-stream
Size: 2388 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151119/862310dd/attachment-0001.obj
More information about the ncl-talk
mailing list