应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 求助Fortran编程
91/1返回列表
查看: 2758  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

liangyje

银虫 (著名写手)

[求助] 求助Fortran编程

长方形图形见附件，长方形尺寸已确定1mm*2mm，网格间距10微米，图形内的弧线上点的坐标值已知，如何判定图形内的节点是属于a区域还是属于b区域呢？
回复此楼

» 收录本帖的淘帖专辑推荐

Monte Carlo 模拟（fortran 编程） fortran 原子距离

» 本帖已获得的红花（最新10朵）

liangyje

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-09-15 18:32:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

【答案】应助回帖

★
jjdg(金币+1): 感谢参与 2011-09-17 00:58:53
ben_ladeng: 2011-09-21 12:42:29
liangyje(金币+3): 2011-09-28 21:54:26
我参想到就是比较节点坐标与 a 和 b 之间弧线上点的坐标值……

但这个问题的关键部分与 Fortran 编程没有任何关系……
赞一下(2人) 回复此楼
2楼2011-09-15 19:34:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liangyje

银虫 (著名写手)
送鲜花一朵
лл
回复此楼
3楼2011-09-20 22:12:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sd116202

木虫 (著名写手)

【答案】应助回帖

★
liangyje(金币+2): 2011-09-21 09:09:43
jjdg(金币+1): 感谢你的建议 2011-09-22 00:32:42
liangyje(金币+3): 2011-09-28 21:54:34
节点坐标与相同x,或者y坐标的弧线点进行比较, 确定一个,比较另一个.

这个是算法问题
赞一下(1人) 回复此楼
4楼2011-09-21 08:58:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superstat54

木虫 (小有名气)

【答案】应助回帖

★ ★ ★
余泽成(金币+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
赞一下(1人) 回复此楼
5楼2011-09-24 17:01:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liangyje

银虫 (著名写手)
非常感谢，你做晶粒长大模拟吗？想向你学习学习
赞一下 回复此楼
6楼2011-09-25 12:03:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superstat54

木虫 (小有名气)
引用回帖:
6楼: Originally posted by liangyje at 2011-09-25 12:03:43:
非常感谢，你做晶粒长大模拟吗？想向你学习学习

没有做过晶粒长大模拟，和这个问题有关吗？我做过原子模拟（MC、MD）方面的，以后应该做相图计算。
你在做晶粒长大模拟吗？
相互学习！
赞一下 回复此楼
7楼2011-09-25 19:28:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liangyje

银虫 (著名写手)
引用回帖:
7楼: Originally posted by superstat54 at 2011-09-25 19:28:04:
没有做过晶粒长大模拟，和这个问题有关吗？我做过原子模拟（MC、MD）方面的，以后应该做相图计算。
你在做晶粒长大模拟吗？
相互学习！

哇，很厉害啊，我现在用MC做晶粒长大模拟，你知道怎么处理多相和相变的问题?
赞一下 回复此楼
8楼2011-09-28 21:53:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

HXC211

新虫 (初入文坛)
引用回帖:
7楼: Originally posted by superstat54 at 2011-09-25 19:28:04
没有做过晶粒长大模拟，和这个问题有关吗？我做过原子模拟（MC、MD）方面的，以后应该做相图计算。
你在做晶粒长大模拟吗？
相互学习！...

您好！我是刚接触模拟计算方面，想向您请教一下
如何用Fortran编写出原子配位数的程序？可以帮帮我吗？非常感谢！
赞一下 回复此楼
9楼2017-03-27 15:13:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liangyje 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 情人节自我反思：在爱情中有过遗憾吗？ +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错，点赞！ +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题，请指教！ 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
关闭
请填处理意见
关闭