[dennou-ruby:000156] Re: reading binary files

[Numaguti Atusi <numa@xxxxxx>]

沼口です．

ほんとにとりあえずですが，

ftp://dennou-t.ms.u-tokyo.ac.jp/
GFD_Dennou_Club/ftp/arch/davis/ex/1999-11-02-numa
の 

mcdcl.991104.tar.gz   :   mcarray を用いた dcl インターフェース

をアップデートしました．

dcl の (oldsrcを除く) ルーチン全部のエントリを入れてあります．
ただし，値が返ってくるやつはまだ使えませんし．
他も全く無保証です．動かないの多々あるでしょう．
単に mksrc.rb で dcl 全体を自動変換しただけですので．

一応 grph2 も一部は動くことを確かめました．

-----
#!/usr/local/bin/ruby
require "mcarray"
require "mcdcl"
include Math

#     program u2df01

      Nx=19; Ny=19
      Xmin=0.0; Xmax=360.0; Ymin=-90.0; Ymax=90.0
      Drad=PI/180; Dz=0.05

      p = FloatArray.new(Nx,Ny)

      for j in 0...Ny
        for i in 0...Nx
          alon = ( Xmin + (Xmax-Xmin) * i / (Nx-1) ) * Drad
          alat = ( Ymin + (Ymax-Ymin) * j / (Ny-1) ) * Drad
          slat = sin(alat)
          p[i,j] = cos(alon) * (1-slat**2) * sin(2*PI*slat) + Dz
        end
      end

      Dcl.gropn(1)
      Dcl.grfrm

      Dcl.grswnd( Xmin, Xmax, Ymin, Ymax )
      Dcl.grsvpt( 0.2, 0.8, 0.2, 0.8 )
      Dcl.grstrn( 1 )
      Dcl.grstrf

      Dcl.usdaxs

      Dcl.udcntr( p, Nx, Nx, Ny )

      Dcl.grcls

-- 
沼口  敦＠北海道大学大学院地球環境科学研究科  大気海洋圏環境科学専攻
tel:011-706-2365 fax:011-706-4865 〒060-0810 札幌市北区北10条西5丁目  
mailto:numa@xxxxxx   http://wwwoa.ees.hokudai.ac.jp/~numa