program main !sms网格转换成mike21网格DIMENSION X(60000),Y(60000),H(60000),NDNN(110000,3),ncbd(10000)dimension NBS(50),NOBD(5000,50),NSED(5000,50),id0(80000),nbi(50)dimension x0(50000),y0(50000),h0(50000)character*2 a*2,b*2,gnn*3,demo*500open(3,file='./2013.8.30(after).grd',status='old')read(3,*)read(3,*)nel,ndedo i=1,nderead(3,*)num,x(i),y(i),h(i)h(i)=h(i)+2.2enddodo i=1,nelread(3,*)num,nd,(ndnn(i,j),j=1,3)enddodo i=1,4 !read(3,*)enddoread(3,*)kbread(3,*)do k=1,kbread(3,*)nbs(k)do i=1,nbs(k)read(3,*)nobd(i,k)enddoenddoclose(3)
c-----------------------------------------------------------open(4,file='滑道1.grd',status='old')read(4,*)read(4,*)nel0,nde0do i=1,nde0read(4,*)num,x0(i),y0(i),h0(i)enddoclose(4)do 25 i=1,ndedo j=1,nde0dis=sqrt((x(i)-x0(j))**2+(y(i)-y0(j))**2)if(dis.le.5)thenh(i)=2.2-h0(j)goto 25endifenddo
25    continue    open(4,file='滑道2.grd',status='old')read(4,*)read(4,*)nel0,nde0do i=1,nde0read(4,*)num,x0(i),y0(i),h0(i)enddoclose(4)do 35 i=1,ndedo j=1,nde0dis=sqrt((x(i)-x0(j))**2+(y(i)-y0(j))**2)if(dis.le.5)thenh(i)=2.2-h0(j)goto 35endifenddo
35    continue
c----------------------------------------------------------id0=0open(18,file='closebd.txt',status='old')do nn=1,100read(18,*,end=181)nkwrite(*,*)nkread(18,*)(ncbd(n),n=1,nk)do i=1,nkid0(ncbd(i))=1 !+nnenddoenddo
181    continueclose(18)do k=1,kbdo i=1,nbs(k)id0(nobd(i,k))=k+1enddoenddowrite(*,*)demo=' PROJCS["Beijing_1954_3_Degree_GK_CM_111E",GEOGCS*["GCS_Beijing_1954",DATUM["D_Beijing_1954",SPHEROID["Krasovsky*_1940",6378245.0,298.3]],PRIMEM["Greenwich",0.0],UNIT["Degree",*0.0174532925199433]],PROJECTION["Gauss_Kruger"],PARAMETER["False*_Easting",500000.0],PARAMETER["False_Northing",0.0],PARAMETER*["Central_Meridian",120.0],PARAMETER["Scale_Factor",1.0],*PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]'
  open(30,file='2013.8.30(after).mesh')write(30,'(I8,a)')nde,demod0=0do i=1,ndewrite(30,'(i8,2f12.2,f10.4,I4)')i,x(i),y(i),h(i)*-1,id0(i)  !enddowrite(30,*)nel,'  3  21'do i=1,nelwrite(30,'(5I6)')i,(ndnn(i,j),j=1,3),d0enddoclose(30)end