| 查看: 3199 | 回复: 8 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
求助Fortran编程
|
||||
长方形图形见附件，长方形尺寸已确定1mm*2mm，网格间距10微米，图形内的弧线上点的坐标值已知，如何判定图形内的节点是属于a区域还是属于b区域呢？ |
回复此楼» 收录本帖的淘帖专辑推荐
 Monte Carlo 模拟（fortran 编程） | fortran 原子距离 |
» 本帖已获得的红花（最新10朵）
liangyje» 猜你喜欢
大佬帮着分析今年的评审意见？该课题连续申请了三年了，明年要不要换个方向。
已经有11人回复
-
两块石头
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】如何写FORTRAN程序实现求平均最近邻距离
已经有34人回复
- 【求助】初学fortran，
已经有15人回复
- 【求助】ANSYS 中命令流编程（fortran）
已经有11人回复
- 【求助】求指点，在Fortran里面，怎样实现“数组维度可调”的数组?
已经有5人回复
- 【求助】请问一下fortran的可视化编程
已经有7人回复
- 【求助】Intel Visual Fortran与Compaq Visual Fortran的区别？
已经有18人回复
- 【求助】有关fortran进行数值模拟
已经有10人回复
8楼2011-09-28 21:53:47
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
2楼2011-09-15 19:34:16
4楼2011-09-21 08:58:38
superstat54
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3809.1
- 散金: 136
- 红花: 10
- 帖子: 256
- 在线: 146.8小时
- 虫号: 1327326
- 注册: 2011-06-20
- 专业: 催化化学
【答案】应助回帖
★ ★ ★
余泽成(金币+3): 谢谢参与应助！ 2011-09-24 20:18:03
liangyje(金币+20): 2011-09-25 11:21:13
liangyje(金币+5): 2011-09-28 21:54:07
余泽成(金币+3): 谢谢参与应助！ 2011-09-24 20:18:03
liangyje(金币+20): 2011-09-25 11:21:13
liangyje(金币+5): 2011-09-28 21:54:07
|
!程序说明 !1 初始文件准备：把弧线上的点的坐标存在arc.txt文件里，文件的第一行为点的个数， ! 下面的每一行分别为x、y坐标值；如下： ! 3 ! 0.01 0.5 ! 0.02 0.49 ! 0.03 0.47 !2 结果文件为result.txt，输出nw*nl的矩阵，对应各个网格，如果该网格属于a则值为1，否则为0 ! ! ! PROGRAM TEST implicit none integer :: i, j, k integer :: nw, nl, n, log real :: with, length, d integer,allocatable :: flag(:,  , point(real,allocatable :: arc(:, !============================================= write(*, *) '程序说明' write(*, *) ' 1 初始文件准备：把弧线上的点的坐标存在arc.txt文件里，文件的第一行为点的个数，' write(*, *) ' 下面的每一行分别为x、y坐标值；如下：' write(*, *) ' 3' write(*, *) ' 0.01 0.5 ' write(*, *) ' 0.02 0.49 ' write(*, *) ' 0.03 0.47 ' write(*, *) ' 2 结果文件为result.txt，输出nw*nl的矩阵，对应各个网格 ' write(*, *) '如果该网格属于a则值为1，否则为0 ' write(*, *) '是否准备好初始文件？（是，输入1；否，输入2)' read(*, *) log if (log.eq.1) then continue else write(*, *) '没有进行判断！' goto 3000 endif !============================================= !初始化，矩形长(length)、宽(with)，网格间距(d) with = 1 length = 2 d = 0.01 !计算nw，nl值，即共有nw*nl个网格 nw = with/d nl = length/d allocate(flag(nw, nl), point(nl)) !=================================== !读取弧线上点的坐标 open(11, file = 'arc.txt') read(11, *) n allocate(arc(n, 2)) do i = 1, n read(11, *) (arc(i, j), j = 1, 2) enddo close(11) !============================================= !找到属于a的区域，point(i)表示，第i列的前point(i)行的网格属于a point = 0 flag = 0 do j = 1, n point(j) = arc(j, 2)/d !write(*, *) point(j) !read(*, *) enddo !属于a的网格赋值为1 do j = 1, n do i = 1, point(j) flag(i, j) = 1 enddo enddo !输出结果 open(12, file ='result.txt') do i = 1, nw do j = 1, nl write(12, '((I4))', advance='no') flag(i, j) enddo write(12, *) enddo close(12) write(*, *) '结果输出完毕' 3000 continue END PROGRAM |
5楼2011-09-24 17:01:45
100