| 查看: 2775 | 回复: 8 | ||||
[求助]
求助Fortran编程
|
长方形图形见附件，长方形尺寸已确定1mm*2mm，网格间距10微米，图形内的弧线上点的坐标值已知，如何判定图形内的节点是属于a区域还是属于b区域呢？ |
回复此楼» 收录本帖的淘帖专辑推荐
 Monte Carlo 模拟（fortran 编程） | fortran 原子距离 |
» 本帖已获得的红花（最新10朵）
liangyje» 猜你喜欢
一志愿985，本科211，0817化学工程与技术319求调剂
已经有4人回复
-
311求调剂
已经有6人回复
-
0703化学调剂 290分有科研经历，论文在投
已经有4人回复
-
化学工程321分求调剂
已经有6人回复
-
本人考085602 化学工程 专硕
已经有9人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有7人回复
-
326求调剂
已经有3人回复
-
焦虑
已经有8人回复
-
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】如何写FORTRAN程序实现求平均最近邻距离
已经有34人回复
- 【求助】初学fortran，
已经有15人回复
- 【求助】ANSYS 中命令流编程（fortran）
已经有11人回复
- 【求助】求指点，在Fortran里面，怎样实现“数组维度可调”的数组?
已经有5人回复
- 【求助】请问一下fortran的可视化编程
已经有7人回复
- 【求助】Intel Visual Fortran与Compaq Visual Fortran的区别？
已经有18人回复
- 【求助】有关fortran进行数值模拟
已经有10人回复
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
2楼2011-09-15 19:34:16
3楼2011-09-20 22:12:08
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
非常感谢，你做晶粒长大模拟吗？想向你学习学习 6楼2011-09-25 12:03:43
superstat54
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3809.1
- 散金: 136
- 红花: 10
- 帖子: 256
- 在线: 146.8小时
- 虫号: 1327326
- 注册: 2011-06-20
- 专业: 催化化学
7楼2011-09-25 19:28:04
8楼2011-09-28 21:53:47
9楼2017-03-27 15:13:04
