24小时热门版块排行榜    

查看: 1559  |  回复: 4

鲤鱼飞

铁虫 (正式写手)

[求助] 求助:用Fortran写一个关于晶粒个数统计和编号的程序。卡了我好久了。

小生近来自学Fortran,开始编写模拟程序实现晶粒长大过程的模拟。但周围没人做这个,而且没人写代码。因此,没有人可以讨论,我是一个头两个大。
目的:求一个逐个统计并赋值的小程序。也就是常见的晶粒长大模拟。
事件:划定一个N*N二维计算区域,如512*512个格点的区域。每个格点上可存储Q个不同的数值pi(0~1之间的实数),该数值的演化由一个PDE计算得到,并交由一个二维数组P(i,j)来保存pi。经过一段时间演化之后,形成如下图所示的形貌组织,即晶粒形成了。演化到最后,虽然每个格点上可以有Q个值,但只有一个为1,其余全部为0,因此属于一个晶粒。但如果格点处在晶界上,这可能存在好多个pi的值(0~1)。
    现在,我需要定义一个变量Q(一维或二维?)来统计生成了多少个晶粒,并给每个晶粒进行编号,如1,2,…等。显然,每个晶粒内含有的格点数目不同。
    这里晶粒的定义:在一个晶粒内部,每个格点上的某一个pi相同(均为1,比如p41为1,其余所有Pi都为0);两个晶粒之间的晶界上0     是否是从(1,1)格点开始,看它的 pi值,如果为1,则记为1#晶粒;再看(1,2)格点,如果小于1,则记为边界;再看(1,3)格点,如果为1,则记为2#晶粒….  还是要看某格点周围8个近邻的值来定?

111.jpg

[ Last edited by 鲤鱼飞 on 2013-1-16 at 20:55 ]
回复此楼

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

资源收集

» 猜你喜欢

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

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

chembetsey

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
从你的叙述看,晶粒的个数就是1的个数啊,按你说的做就是了
2楼2013-01-16 23:14:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

鲤鱼飞

铁虫 (正式写手)

No.
    512*512个格点计算区域,或许只有3个晶粒,或许30个晶粒。  一个晶粒或许包含100个格点,或许包含499个格点,在这个晶粒内部,所有格点的某一个p值为1,其余所有p值为0。晶界的宽度为过渡区,例如包含3-7个格点。
CAEP
3楼2013-01-18 08:31:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chembetsey

木虫 (小有名气)

【答案】应助回帖

★
jjdg: 金币+1, 感谢参与 2013-01-19 03:20:38
数组既然在那里,想怎么统计就怎么统计,没有太大困难
你若是觉得有困难,那是你对统计的条件还没想清楚
4楼2013-01-18 22:31:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

鲤鱼飞

铁虫 (正式写手)

最近想了下。可能要用三维数组P(i,j,k)来存储相场值和晶粒编号。 其中i,j记录二维点阵,k(1:Q)记录每个格点上可以存储的Q个相场值。   
   不太清楚,如何将晶粒和格点联系起来,而且晶粒的编号要固定不变,也就是每个格点归属的晶粒编号不变。当有新的晶粒出现后,还要增加一个晶粒编号,但格点个数不变...
CAEP
5楼2013-01-26 21:00:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 鲤鱼飞 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[硕博家园] 极限拉扯,如何抉择 +9 Isak0320 2026-10-12 13/650 2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见