SMS模型格網(wǎng)轉(zhuǎn)換為MIKE21的格網(wǎng)源代碼 program main !sms網(wǎng)格轉(zhuǎn)換成mike21網(wǎng)格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

?

posted on 2014-01-20 11:44?向北方 閱讀(...) 評論(...) 編輯 收藏

轉(zhuǎn)載于:https://www.cnblogs.com/China3S/p/3526673.html

總結(jié)

以上是生活随笔為你收集整理的SMS模型格网转换为MIKE21的格网源代码的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。