24小时热门版块排行榜    

查看: 1504  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 鲤鱼飞 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-18 01:04 by cITaGg3p5edr
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-18 00:28 by cITaGg3p5edr
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-13 9/450 2026-09-17 23:05 by cuPZTDXS3VOd
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 LwdutQ8HoqWP 2026-09-13 8/400 2026-09-17 22:53 by cuPZTDXS3VOd
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 7/350 2026-09-17 22:51 by cuPZTDXS3VOd
[有机交流] 各位大神,目前国内有哪些比较好用的逆合成软件? 20+3 zapen 2026-09-15 7/350 2026-09-17 16:40 by 北京莱茵编辑
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-17 14:27 by ObVzqQQh1rrE
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 40ms4Wlo4umh 2026-09-14 6/300 2026-09-17 13:53 by ObVzqQQh1rrE
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-17 13:39 by ObVzqQQh1rrE
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 23jxep3nCNZb 2026-09-14 6/300 2026-09-17 13:29 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 mibUvS8DDCwf 2026-09-15 4/200 2026-09-17 06:26 by k06BKNblrGRH
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-17 03:14 by Tql5LQhh5rLK
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 3n8v2C8RimXI 2026-09-13 8/400 2026-09-17 01:38 by Tql5LQhh5rLK
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-16 03:14 by DgNGHc3h5tPl
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
信息提示
请填处理意见