liangyje

銀蟲 (著名寫手)

應助: 50 (小學生)
金幣: 77.6
散金: 1285
紅花: 7
帖子: 1293
在線: 1093.7小時
蟲號: 798407
注冊: 2009-06-24
專業(yè): 無機非金屬類光電信息與功

[求助] 求助Fortran編程

長方形圖形見附件，長方形尺寸已確定1mm*2mm，網(wǎng)格間距10微米，圖形內(nèi)的弧線上點的坐標值已知，如何判定圖形內(nèi)的節(jié)點是屬于a區(qū)域還是屬于b區(qū)域呢？

回復此樓

» 收錄本帖的淘帖專輯推薦

Monte Carlo 模擬（fortran 編程）

fortran 原子距離

» 本帖已獲得的紅花（最新10朵）

liangyje

» 猜你喜歡

286求調(diào)劑已經(jīng)有9人回復
一志愿南昌大學324求調(diào)劑已經(jīng)有6人回復
0703化學/290求調(diào)劑/本科經(jīng)歷豐富/工科也可已經(jīng)有10人回復
291求調(diào)劑已經(jīng)有5人回復
308求調(diào)劑已經(jīng)有7人回復
085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富已經(jīng)有4人回復
壓國家一區(qū)線，求導師收留，有恩必謝！已經(jīng)有3人回復
一志愿中南大學化學0703總分337求調(diào)劑已經(jīng)有4人回復
083000學碩274求調(diào)劑已經(jīng)有6人回復
304求調(diào)劑已經(jīng)有4人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

【求助】如何寫FORTRAN程序?qū)崿F(xiàn)求平均最近鄰距離已經(jīng)有34人回復
【求助】初學fortran，已經(jīng)有15人回復
【求助】ANSYS 中命令流編程（fortran）已經(jīng)有11人回復
【求助】求指點，在Fortran里面，怎樣實現(xiàn)“數(shù)組維度可調(diào)”的數(shù)組? 已經(jīng)有5人回復
【求助】請問一下fortran的可視化編程已經(jīng)有7人回復
【求助】Intel Visual Fortran與Compaq Visual Fortran的區(qū)別？已經(jīng)有18人回復
【求助】有關fortran進行數(shù)值模擬已經(jīng)有10人回復

1樓 2011-09-15 18:32:42

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

snoopyzhao

至尊木蟲 (職業(yè)作家)

程序強帖: 16
應助: 157 (高中生)
貴賓: 0.02
金幣: 18844.7
紅花: 29
帖子: 3803
在線: 1422.4小時
蟲號: 183750
注冊: 2006-02-13
專業(yè): 污染生態(tài)化學

【答案】應助回帖

★
jjdg(金幣+1): 感謝參與 2011-09-17 00:58:53
ben_ladeng: 2011-09-21 12:42:29
liangyje(金幣+3): 2011-09-28 21:54:26

我參想到就是比較節(jié)點坐標與 a 和 b 之間弧線上點的坐標值……

但這個問題的關鍵部分與 Fortran 編程沒有任何關系……

贊一下(2人)

回復此樓

2樓2011-09-15 19:34:16

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

liangyje

銀蟲 (著名寫手)

應助: 50 (小學生)
金幣: 77.6
散金: 1285
紅花: 7
帖子: 1293
在線: 1093.7小時
蟲號: 798407
注冊: 2009-06-24
專業(yè): 無機非金屬類光電信息與功

送鮮花一朵

лл

回復此樓

3樓2011-09-20 22:12:08

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

sd116202

木蟲 (著名寫手)

應助: 11 (小學生)
金幣: 536.7
散金: 706
紅花: 8
帖子: 1747
在線: 235.9小時
蟲號: 1111478
注冊: 2010-09-30
性別: GG
專業(yè): 工程熱物理相關交叉領域

【答案】應助回帖

★
liangyje(金幣+2): 2011-09-21 09:09:43
jjdg(金幣+1): 感謝你的建議 2011-09-22 00:32:42
liangyje(金幣+3): 2011-09-28 21:54:34

節(jié)點坐標與相同x,或者y坐標的弧線點進行比較, 確定一個,比較另一個.

這個是算法問題

贊一下(1人)

回復此樓

4樓2011-09-21 08:58:38

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

superstat54

木蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 3809.1
散金: 136
紅花: 10
帖子: 256
在線: 146.8小時
蟲號: 1327326
注冊: 2011-06-20
專業(yè): 催化化學

【答案】應助回帖

★ ★ ★
余澤成(金幣+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文件里，文件的第一行為點的個數(shù)，
!  下面的每一行分別為x、y坐標值；如下：
!  3
!  0.01 0.5
!  0.02 0.49
!  0.03 0.47
!2 結(jié)果文件為result.txt，輸出nw*nl的矩陣，對應各個網(wǎng)格，如果該網(wǎng)格屬于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文件里，文件的第一行為點的個數(shù)，'
  write(*, *) ' 下面的每一行分別為x、y坐標值；如下：'
  write(*, *) ' 3'
  write(*, *) ' 0.01 0.5 '
  write(*, *) ' 0.02 0.49 '
  write(*, *) ' 0.03 0.47 '
  write(*, *) ' 2 結(jié)果文件為result.txt，輸出nw*nl的矩陣，對應各個網(wǎng)格 '
  write(*, *) '如果該網(wǎng)格屬于a則值為1，否則為0 '
  write(*, *) '是否準備好初始文件？（是，輸入1；否，輸入2)'
  read(*, *) log
  if (log.eq.1) then
continue
  else
write(*, *) '沒有進行判斷！'
goto 3000
  endif
  !=============================================
  !初始化，矩形長(length)、寬(with)，網(wǎng)格間距(d)
  with = 1
  length = 2
  d = 0.01
  !計算nw，nl值，即共有nw*nl個網(wǎng)格
  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的區(qū)域，point(i)表示，第i列的前point(i)行的網(wǎng)格屬于a
  point = 0
  flag = 0
  do j = 1, n
point(j) = arc(j, 2)/d
!write(*, *) point(j)
!read(*, *)
  enddo
  !屬于a的網(wǎng)格賦值為1
  do j = 1, n
do i = 1, point(j)
   flag(i, j) = 1
enddo
  enddo
  !輸出結(jié)果
  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(*, *) '結(jié)果輸出完畢'
3000 continue
END PROGRAM

贊一下(1人)

回復此樓

5樓2011-09-24 17:01:45

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

liangyje

銀蟲 (著名寫手)

應助: 50 (小學生)
金幣: 77.6
散金: 1285
紅花: 7
帖子: 1293
在線: 1093.7小時
蟲號: 798407
注冊: 2009-06-24
專業(yè): 無機非金屬類光電信息與功

非常感謝，你做晶粒長大模擬嗎？想向你學習學習

贊一下

回復此樓

6樓2011-09-25 12:03:43

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

superstat54

木蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 3809.1
散金: 136
紅花: 10
帖子: 256
在線: 146.8小時
蟲號: 1327326
注冊: 2011-06-20
專業(yè): 催化化學

引用回帖:

6樓: Originally posted by liangyje at 2011-09-25 12:03:43:

非常感謝，你做晶粒長大模擬嗎？想向你學習學習

沒有做過晶粒長大模擬，和這個問題有關嗎？我做過原子模擬（MC、MD）方面的，以后應該做相圖計算。
你在做晶粒長大模擬嗎？
相互學習！

贊一下

回復此樓

7樓2011-09-25 19:28:04

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

liangyje

銀蟲 (著名寫手)

應助: 50 (小學生)
金幣: 77.6
散金: 1285
紅花: 7
帖子: 1293
在線: 1093.7小時
蟲號: 798407
注冊: 2009-06-24
專業(yè): 無機非金屬類光電信息與功

引用回帖:

7樓: Originally posted by superstat54 at 2011-09-25 19:28:04:
沒有做過晶粒長大模擬，和這個問題有關嗎？我做過原子模擬（MC、MD）方面的，以后應該做相圖計算。
你在做晶粒長大模擬嗎？
相互學習！

哇，很厲害啊，我現(xiàn)在用MC做晶粒長大模擬，你知道怎么處理多相和相變的問題?

贊一下

回復此樓

8樓2011-09-28 21:53:47

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

HXC211

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 3
帖子: 14
在線: 6.5小時
蟲號: 6083230
注冊: 2017-03-21
專業(yè): 金屬基復合材料

引用回帖:

7樓: Originally posted by superstat54 at 2011-09-25 19:28:04
沒有做過晶粒長大模擬，和這個問題有關嗎？我做過原子模擬（MC、MD）方面的，以后應該做相圖計算。
你在做晶粒長大模擬嗎？
相互學習！...

您好！我是剛接觸模擬計算方面，想向您請教一下
如何用Fortran編寫出原子配位數(shù)的程序？可以幫幫我嗎？非常感謝！

贊一下

回復此樓

9樓2017-03-27 15:13:04

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 291求調(diào)劑 +5	HanBeiNingZC 2026-03-24	5/250	2026-03-28 07:44 by barnett0632
[考研] 083000學碩274求調(diào)劑 +6	Li李魚 2026-03-26	6/300	2026-03-28 07:25 by barnett0632
[考研] 352分化工與材料 +5	海納百川Ly 2026-03-27	5/250	2026-03-28 03:39 by fmesaito
[考研] 265求調(diào)劑 +8	小木蟲085600 2026-03-27	8/400	2026-03-27 22:16 by 無際的草原
[考研] 277跪求調(diào)劑 +5	1915668 2026-03-27	8/400	2026-03-27 19:53 by WYUMater
[考研] 調(diào)劑 +3	李嘉圖·S·路 2026-03-27	3/150	2026-03-27 11:19 by wangjy2002
[考研] 調(diào)劑推薦 +5	清酒714 2026-03-26	6/300	2026-03-27 11:12 by 不吃魚的貓
[考研] 一志愿陜師大生物學071000，298分，求調(diào)劑 +5	SYA！ 2026-03-23	5/250	2026-03-27 09:29 by 不吃魚的貓
[考研] 085602 289分求調(diào)劑 +8	WWW西西弗斯 2026-03-24	8/400	2026-03-26 16:33 by 不吃魚的貓
[考研] 一志愿南京郵電大學 288分材料考研求調(diào)劑 +3	jl0720 2026-03-26	3/150	2026-03-26 13:39 by zzll406
[考研] 303求調(diào)劑 +6	藍山月 2026-03-25	6/300	2026-03-25 22:47 by 418490947
[考研] 求調(diào)劑 +3	QiMing7 2026-03-25	3/150	2026-03-25 21:13 by 給你你注意休息
[考研] 26考研-291分-廈門大學（085601）-柔性電子學院材料工程專業(yè)求調(diào)劑 +3	min3 2026-03-24	4/200	2026-03-25 18:22 by xcjcqu
[考研] 0854AI CV方向招收調(diào)劑 +4	章小魚567 2026-03-23	4/200	2026-03-25 17:04 by CoderLoser
[考研] 0703化學求調(diào)劑 +6	奶油草莓. 2026-03-22	7/350	2026-03-25 10:00 by shangxh
[考研] 化學調(diào)劑 +6	yzysaa 2026-03-21	6/300	2026-03-25 09:27 by aa331100
[考研] 311求調(diào)劑 +3	冬十三 2026-03-24	3/150	2026-03-24 21:31 by peike
[考博] 26申博自薦 +3	whh869393 2026-03-24	3/150	2026-03-24 09:55 by 21018060
[論文投稿] 急發(fā)核心期刊論文 +3	賢達問津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹
[考研] 石河子大學（211、雙一流）碩博研究生長期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋懷