fortran90でgrd

fortran90でのGMTのgrdファイルを作るサブルーチン．

fortran用のnetcdfライブラリをインストールしてください．

ただし，GMTver4.0以下対応，GMT4.1.3は未対応

これまでのものよりサブルーチンの引数を減らしました．

----------------------------------

zz：配列「!#############################」の部分の２箇所を編集してください．

nx,ny：配列の大きさ，整数

xyrange(4)：xmin,xmax,ymin,ymax，実数配列，配列の大きさは４

dx：格子幅，実数，ここではdx=dy

filegrd：作成する配列名：文字列

c=====================================================

subroutine mk_grd_z(nx,ny,xyrange,dx,zz,filgrd)

IMPLICIT REAL*8 (A-H,O-Z)

c dimension xrange(*),yrange(*),zrange(*),space(*)

dimension xyrange(*)

dimension xrange(2),yrange(2),zrange(2),space(2)

real*8 zz(nx,ny,2) !###################################

integer ncid,icon,idside,idxr,ixy

integer ixrng(1),ixsize(1),node(1),istart(1),icount(1),idim(2)

character filgrd*18

real*8 scale(1)

real*8, allocatable:: z(:,:)

include 'netcdf.inc'

allocate(z(nx,ny),stat=ier)

xrange(1)=xyrange(1);xrange(2)=xyrange(2);

yrange(1)=xyrange(3);yrange(2)=xyrange(4);

space(1)=dx;space(2)=dx;

idim(1)=nx

idim(2)=ny

ixy=idim(1)*idim(2)

do i=1,nx

do j=1,ny

z(i,j)=zz(i,j,2) !###################################

enddo

zrange(1)=minval(z);zrange(2)=maxval(z);

icon=nf_create(filgrd,NF_CLOBBER,ncid)

if(icon.ne.0) write(6,*) "creat",icon

icon=nf_def_dim(ncid,"side",2,idside)

icon=nf_def_dim(ncid,"xysize",ixy,ixys)

ixrng(1)=idside

icon=nf_def_var(ncid,'x_range',NF_DOUBLE,1,ixrng,idxr)