sw_zfr 2014-08-23 04:55
浏览 1182

请教这个fortran脚本怎么用?据说是把txt文件转化成cst文件 linux下怎么操作 谢谢

下面是.f文件具体内容
parameter(npointp=900000)
real x(900000),y(900000)
integer num(900000),linenum(900000),npl(900000),lnum,nplt,
* lasti,i,j

    open(10, file='coastline.txt') 
    open(11, file='coastline.cst') 

  do i=1,npointp 
     read(10,*,end=100) num(i), y(i), x(i) 
  enddo 

100 continue
npoint=i-2
close(10)
lnum=0
nplt=0
lasti=0
do i=1,npoint

     if(num(i).eq.3) then 
    lnum=lnum+1 
    npl(lnum)=i-lasti 
    lasti=i 
     end if 
  enddo 

  write(11, *) 'COAST' 
write(11,'(I5,f4.1)') lnum,0.0 
  j=1 
  do i=1, npoint 
  print*,i 
      if(num(i).eq.1) then 
    write(11,'(I5,I4)')npl(j),0 
    j=j+1 
  endif 
  write(11,30) x(i),y(i),0.0 
enddo 

30 format(1x,2f20.4,f5.1)

   close(11) 
end
  • 写回答

0条回答 默认 最新

    报告相同问题?

    悬赏问题

    • ¥15 矩阵加法的规则是两个矩阵中对应位置的数的绝对值进行加和
    • ¥15 活动选择题。最多可以参加几个项目?
    • ¥15 飞机曲面部件如机翼,壁板等具体的孔位模型
    • ¥15 vs2019中数据导出问题
    • ¥20 云服务Linux系统TCP-MSS值修改?
    • ¥20 关于#单片机#的问题:项目:使用模拟iic与ov2640通讯环境:F407问题:读取的ID号总是0xff,自己调了调发现在读从机数据时,SDA线上并未有信号变化(语言-c语言)
    • ¥20 怎么在stm32门禁成品上增加查询记录功能
    • ¥15 Source insight编写代码后使用CCS5.2版本import之后,代码跳到注释行里面
    • ¥50 NT4.0系统 STOP:0X0000007B
    • ¥15 想问一下stata17中这段代码哪里有问题呀