c program ognca7 c c data conversion c parameter (nx=1338,ny= 446,nz= 894,nxp=223) c parameter (nx=1338,ny= 446,nz= 894,nxp=350) parameter (nx= 888,ny= 222,nz= 446,nxp=294) c parameter (mx= 600,my= 200,mz= 400,mxp=100) parameter (mx= 800,my= 200,mz= 400,mxp=300) parameter (nx2=nx+2,ny2=ny+2,nz2=nz+2) parameter (mx2=mx+2,my2=my+2,mz2=mz+2) parameter (n1=nx+2,n2=n1*(ny+2),n3=n2*(nz+2)) parameter (m1=mx+2,m2=m1*(my+2),m3=m2*(mz+2)) parameter (nb=8,nbb=11) parameter (last=2,itap=1) c dimension gf(nx2,ny2),gg(mx2,my2) open(11,file='/home/fpfs/usr6/a41456a/fpfs0/fa10900379.data', 1 access='sequential',form='unformatted') open(12,file='/home/fpfs/usr6/a41456a/fpfs0/fa1090037a.data', 1 access='sequential',form='unformatted') c open(21,file='/home/fpfs/usr6/a41456a/fpfs0/fa109013g.data', 1 access='sequential',form='unformatted') open(22,file='/home/fpfs/usr6/a41456a/fpfs0/fa109013h.data', 1 access='sequential',form='unformatted') c c nmz=(nz-mz)/2 nmx=(nx-mx)/2-(nxp-mxp) nmz1=1+nmz nmz2=mz2+nmz c do 100 jj=1,last itapi=jj+10 itapo=jj+20 c do 410 ii=1,itap do 174 m=1,nb do 174 k=1,nz2 read(itapi) gf if(k.lt.nmz1) go to 174 if(k.gt.nmz2) go to 174 c do 10 j=1,my2 do 10 i=1,mx2 i1=i+nmx gg(i,j)=gf(i1,j) 10 continue write(itapo) gg c 174 continue 410 continue c 100 continue c c stop end