プラズマ核融合学会誌４月【８４－４】／講座５

Size: px

Start display at page:

Download "プラズマ核融合学会誌４月【８４－４】／講座５"

ちえこさくいし
5 years ago
Views:

1 Practical Data Analysis Using Open Source Software 5. Ruby for Geophysical Fluid Sciences HORINOUCHI Takeshi corresponding author s horinout rish.kyoto-u.ac.jp

3 air lon lat level time global_temp.rb require numru/ggraph include NumRu iws = ARGV[0] 1 DCL.gropn(iws) DCL.sgpset( lcntl, false) gphys = GPhys::IO.open( air.mon.ltm.nc, air ) GGraph.contour( gphys.cut( level =>925) ) DCL.grcls % ruby global_temp.rb *** MESSAGE (SWDOPN) *** GRPH1 : STARTED / IWS = 1. *** WARNING (STSWTR) *** WORKSTATION VIEWPORT WAS MODIFIED. *** MESSAGE (SWPCLS) *** GRPH1 : PAGE = 1 COMPLETED. air.mon.ltm. nc air

4 global_temp.rb cut cut (true,true,925,true) global_temp2.rb PostScript include iws GGraph ggraph_startup.rb % irb r ggraph_startup.rb *** MESSAGE (SWDOPN) *** GRPH1 : STARTED / IWS = 1. irb(main):001:0> gphys = GPhys::IO.open( air.mon.ltm.nc, air ) => <GPhys grid=<4d grid <axis pos=< lon in air.mon.ltm.nc sfloat[144]>> <axis pos=< lat in air.mon.ltm.nc sfloat[73]>> <axis pos=< level in air.mon.ltm.nc sfloat[17]>> <axis pos=< time in air.mon.ltm.nc float[12]>>> data=< air in air.mon.ltm.nc sfloat[144, 73, 17, 12]>> irb(main):002:0> GGraph.contour( gphys.cut( level =>925) ) *** WARNING (STSWTR) *** WORKSTATION VIEWPORT WAS MODIFIED. => nil irb(main):003:0>

5 => require numru/ggraph include NumRu iws = ARGV[0] 1 DCL.swpset( iwidth,550) # DCL.swpset( iheight,700) # DCL.sgscmn(5) # DCL.gropn(iws) DCL.sldiv( y,1,2) # DCL.sgpset( isub, 96) # DCL.sgpset( lfull, true) # viewport DCL.uzfact(0.8) # 0.8 gphys = GPhys::IO.open( air.mon.ltm.nc, air ) #<> GGraph.set_fig( itr =>12, viewport =>[0.07,0.93,0.1,0.55]) GGraph.set_map( coast_world =>true, grid =>true) GGraph.tone( gphys.cut( level =>500), true, interval =>5, clr_min =>50, # title => Temperature (Jan, 500hPa), annotate =>false ) GGraph.contour( gphys.cut( level =>500), false, interval =>5 ) #GGraph.color_bar # #<> GGraph.set_fig( itr =>2, viewport =>[0.25,0.75,0.12,0.52]) GGraph.line( gphys.cut( lat =>35).mean( lon ), true, exchange =>true, title => Temperature (35N Jan), annotate =>false ) DCL.grcls 2 require numru/gphys include NumRu def usage \n USAGE: % ruby #{$0} filename variablename n end filename = ARGV[0] or raise(usage) # varname = ARGV[1] or raise(usage) # FFT GPhys::fft_ignore_missing(true, 0.0) # 0.0 gphys = GPhys::IO.open(filename, varname) pre = gphys.detrend(3).cos_taper(3) # &cosine sp = pre.fft(false,0,3).abs ** 2

6 # false true 0,3==> fft power = sp.rawspect2powerspect(0,3).spect_one_sided(3).spect_zero_centering(0) # 1-sided etc outfile = NetCDF.create( pw2d_#{varname}.nc ) # GPhys::IO.write(outfile, power) # outfile.close or = n moment numru/gphys require ruby moment.rb require power moment.rb moment.rb require numru/gphys module NumRu class GPhys # n def moment(n, *dim) ( ( self - self.mean(*dim) )**n ).mean(*dim) end end end ###### test program ##### if $0 == FILE include NumRu raise( need 3 or more args ) if ARGV.length<3 fnm, vnm, n, = ARGV dims = ARGV[3..-1].collect{ s s.to_i} gphys = GPhys::IO.open(fnm,vnm) mmt = gphys.moment(n.to_i, *dims) p mmt end

7 dims * moment n dims * + moment gphys.moment(2, 0) self - self.mean(*dim) n [] % ruby moment.rb filename varname n [dims..] NArray require narray include NMath # NArray Math na = NArray.float(4,3,2).random! # nb = na[true,1..-1,0] # nc = sin(nb)**2 * cos(nb) # nd = na.mean(0,1) # na[0,false] = na[-1,false] # # false m = NMatrix[ [1.0,2.0], [3.0, 4.0] ] # M y = NVector[ 1.0, 1.0 ] # y x = m.lu.solve(y) # LU Mx y

8 GPhys cut [] GPhys::IO.open open open

9 % rdoc --charset sjis --inline-source global_ temp.rb global_temp2.rb moment.rb pw2d.rb RDoc html

10 GUI gave

11 Gfdnavi Object-oriented software construction

netcdf

netcdf 1. Anetcdf.rb netcdf C ruby open new create NetCDF C filename String NetCDF NetCDF_open mode r r, w share false true or false open open netcdf filename String NetCDF NetCDF_create noclobber false true