24小时热门版块排行榜    

查看: 793  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

宝莲灯

木虫 (小有名气)

虫虫


[交流] 【讨论】编程求助

在Linux系统下编译fortran。想取sufer中geo2km-gzdw中的部分数据,把它加高20米,然后再放回到geo2km-gzdw文件中原来的位置。该怎么做?
       我编的程序如下:
                   dimension a(100,100)
                  sed -n '547,623p' geo2km-gzdw.dat > ab.dat   *这句话有问题,下面的都运行过,没有问题*
                  open (10,file='ab.adt')
                  do j=1,76
                  read (10,*) (a(i,j),i=1,66)
                  do i=1,66
                   if (i.ge.23.and.i.le.32.and.j.gz.32.and.j.le.42) then
      20        format(10f20.9)
                   a(i,j)=a(i,j)+20
                   endif
                  end do
                   end do
                   do j=1,76
                   write(*,30) a(i,j),i=1,66)
      30          format(66F20.9)
                  end do
                    close(10)
                     stop
                     end
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yscot

木虫 (正式写手)

★ ★
宝莲灯(金币+1):谢谢参与
余泽成(金币+1):谢谢参与应助! 2010-09-04 22:24:48
宝莲灯(金币+10): 2010-09-06 08:49:22
建议楼主把 sed -n '547,623p' geo2km-gzdw.dat > ab.dat
直接在程序中用fortran语句表达出来
或者程序中直接读入处理后的输入文件

把linux命令直接写到程序里好像执行不了的啊~
3楼2010-09-04 22:09:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 宝莲灯 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见