[Thread Prev][Thread Next][Index]

Re: [ferret_users] limit of numbers of viewports



Dear Mick

Thank you very much for your suggestion.

I had tried but the result is the same as

XgksDuplicatePrimi()  300 Storage overflow has occurred in
GKS
XgksDuplicatePrimi()  300 Storage overflow has occurred in
GKS

Anyhow, it is good to avoid the problem accordint to your
suggestion.

Best Regards,

Takaya namba



--- Mick Spillane <Mick.Spillane@noaa.gov> からのメッセ
ージ:
> Hello Takaya,
> You might try adding " ; can memory/all " as the
> last command in the
> REPEAT loop.
> In my simple/small demo the same image is produced
> but without running
> it with a full scale dataset like yours, I can't be
> sure it will work in
> your situation.
> If this suggestion does not work, perhaps you could
> smooth and regrid
> the data to the largest size that will work (you say
> 100x100 does work).
> The full resolution image would need to be so big
> anyway in terms of
> pixels that it may be difficult to actually see the
> detail. If someone
> needed to see the full resolution 570x684 image for
> a particular
> year/month you could generate the monthly images as
> individual plots
> then treat the lower resolution "quilt" image with a
> little HTML as a
> set of low-res "thumbnails" which, with a click,
> would lead to the full
> image. Or perhaps an animation or "slideshow"
> stepping through the
> high-res images would be suitable.
> Mick
> 
> Namba Takaya wrote:
> 
> >Dear Mick
> >
> >When I try to do the grid of 570 
> >x 684, it has a problem.If I make the size of grid
> >smaller, ,eg, 100 x 100, it works very well. So the
> >problem looks like too much memory consumption.Are
> there
> >any  ways to refresh the memory in this script.Or
> any ways
> >to run this script with large grid.
> >
> >Best Regards,
> >
> >Takaya namba
> >
> >
> >--- Mick Spillane <mick.spillane@noaa.gov>
> からのメッセ
> >ージ:
> >  
> >
> >>Hi Takaya & Ansley
> >>This problem might be better solved with
> regridding
> >>than with multiple
> >>viewports. Below is a demo, for an admittedly
> small
> >>x-y grid, that does
> >>the job; a GIF image of the result. is attached.
> >>Good luck,
> >>Mick
> >>!-------------------
> >>! 6-Dec-05 : drawing 20 years of monthly data (on
> an
> >>x-y grid) without
> >>using multiple viewports
> >>! NOTE: the lines in the repeat loop end in
> >>"backslash" which the
> >>mailtool seems to alter
> >>
> >>! preliminary step - create some fake data 240
> >>months on a 9x7 grid
> >>def axis/x=1:9:1 xax ; def axis/y=1:7:1 yax ; def
> >>axis/t=15-dec-1981:15-dec-2000/npoints=240 tax
> >>def grid/x=xax/y=yax/t=tax grd
> >>let/title="240 months of data on 9x7 grid"
> >>
> >>    
> >>
>
>v=(1+randu(l[g=grd]))*(10-((x[g=grd]-5)^2+(y[g=grd]-4)^2)^0.5)
> >  
> >
> >>! ... and save it to a file
> >>save/file=testdata.nc v
> >>
> >>! remove the variable used to define the "data"
> >>can var v
> >>
> >>! ... then read back the "data" ... now the demo
> >>begins
> >>use testdata.nc
> >>
> >>! define a conversion from time index L to month
> and
> >>year
> >>let mod12=mod(l,12) ; set var/bad=0 mod12 ; let
> >>month=missing(mod12,12)
> >>let year=1981+(l-month)/12
> >>
> >>! create a base plot with appropriate axes to span
> >>the entire year and
> >>month range
> >>plot/vs/nolab/nokey/set 1,1
> >>ppl xaxis,1980,2001,1 ; ppl yaxis,0,13,1 ; ppl
> >>xfor,(i4) ; ppl yfor,(i2)
> >>; ppl plot
> >>
> >>! then, one by one, overlay the monthly gridded
> >>images
> >>repeat/l=1:240 (def
> >>axis/x=`year-0.45`:`year+0.45`/npoints=9 xnew ;\
> >>def axis/y=`month-0.4`:`month+0.4`/npoints=7 ynew
> ;\
> >>def grid/like=v/x=xnew/y=ynew gnew ; let
> >>vnew=v[g=gnew,gx=@asn,gy=@asn] ;\
> >>
> >>    
> >>
>
>shade/o/nolab/nokey/lev=(0,20,0.5)/hlim=1980:2001/vlim=0:13
> >  
> >
> >>vnew ; ppl
> >>shaset reset )
> >>
> >>frame/file=quilt.gif
> >>
> >>q
> >>
> >>
> >>    
> >>
> >
> >  
> >
> >>ATTACHMENT part 2 image/gif x-mac-type=0;
> >>    
> >>
> >x-mac-creator=0; name=quilt.gif
> >  
> >
> >> !  6-Dec-05 : drawing 20 years of monthly data
> (on
> >>an x-y grid) without using multiple viewports 
> >>
> >>! preliminary step - create some fake data 240
> >>months on a 9x7 grid
> >>def axis/x=1:9:1 xax ; def axis/y=1:7:1 yax ; def
> >>axis/t=15-dec-1981:15-dec-2000/npoints=240 tax
> >>def grid/x=xax/y=yax/t=tax grd
> >>let/title="240 months of data on 9x7 grid"
> >>
> >>    
> >>
>
>v=(1+randu(l[g=grd]))*(10-((x[g=grd]-5)^2+(y[g=grd]-4)^2)^0.5)
> >  
> >
> >>! ... and save it to a file
> >>save/file=testdata.nc v
> >>
> >>! remove the variable used to define the "data"
> >>can var v
> >>
> >>! ... then read back the "data" ... now the demo
> >>begins 
> >>use testdata.nc
> >>
> >>! define a conversion from time index L to month
> and
> >>year
> >>let mod12=mod(l,12) ; set var/bad=0 mod12 ; let
> >>month=missing(mod12,12) ; let
> year=1981+(l-month)/12
> >>
> >>! create a base plot with appropriate axes to span
> >>the entire year and month range
> >>plot/vs/nolab/nokey/set 1,1 ; ppl
> xaxis,1980,2001,1
> >>; ppl yaxis,0,13,1 ; ppl xfor,(i4) ; ppl yfor,(i2)
> ;
> >>ppl plot
> >>
> >>! then, one by one, overlay the monthly gridded
> >>images
> >>repeat/l=1:240 (def
> >>axis/x=`year-0.45`:`year+0.45`/npoints=9 xnew ;
> def
> >>axis/y=`month-0.4`:`month+0.4`/npoints=7 ynew ;\
> >>def grid/like=v/x=xnew/y=ynew gnew ; let
> >>vnew=v[g=gnew,gx=@asn,gy=@asn] ;\
> >>
> >>    
> >>
>
>shade/o/nolab/nokey/lev=(0,20,0.5)/hlim=1980:2001/vlim=0:13
> >  
> >
> >>vnew ; ppl shaset reset )
> >>frame/file=quilt.gif
> >>
> >>q
> 
=== message truncated ===

--------------------------------------
STOP HIV/AIDS.
Yahoo! JAPAN Redribbon Campaign 2005
http://pr.mail.yahoo.co.jp/redribbon/

[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement