[ncl-talk] drawing satellite passes(14 per day)
Krishna C
chandrakrishna.90 at gmail.com
Thu Nov 19 04:38:44 MST 2015
Hi
I have also plotted the full swath, obviously it takes a lot of time
I have attached the full swath plot as well
Hope this helps
Regards
-Krishna-
On Thu, Nov 19, 2015 at 4:59 PM, Krishna C <chandrakrishna.90 at gmail.com>
wrote:
> 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/09c905d3/attachment-0001.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: s6.pdf
Type: application/pdf
Size: 2009568 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151119/09c905d3/attachment-0001.pdf
More information about the ncl-talk
mailing list