| 查看: 589 | 回复: 1 | ||
[求助]
求助二维fortran程序的输出问题
|
|
fortran菜鸟一枚,求大家帮助…… 程序是拉氏坐标下求解二维问题,初始时刻是二维矩形形网格,随着时间推进网格坐标会变化。 所以定义了坐标的二维数组 x ( i , j ),y ( i , j ),i 和 j 的范围分别是 im,jm 然后每次循环的时候都是做 do i=1, im do j=1, jm …… end do end do 所以想问: 1、结果输出时,直接在这个循环中加入 write (10 , *) 语句就可以么? 2、如果想输出当 j 为某一个特定值的时候,按照如下进行输出可以吗? 运行的结果不对啊,在 j 取50的时候,输出的结果中,第50行的数据为NaN。 do i=1,im do j=1,jm tao_in(i,j)=tao_in(i,j)+(……) u_in(i,j)=u_in(i,j)-(……) v_in(i,j)=v_in(i,j)-(……) Ene_in(i,j)=Ene_in(i,j)-(……) if (j==50) then write(10,*) x_in(i,j),p_in(i,j),u_in(i,j),v_in(i,j),den_in(i,j),tao_in(i,j),lamta_in(i,j) end if end do end do 3、如果我想输出二维网格,就是画出某一个时间步的网格,该怎么输出?我用下面的试了,用tecplot画出来的明显不是正常的二维网格啊。 do i=1,im do j=1,jm …… write(20,*) i,j , x(i ,j) , y(i,j) end do end do 真心求教,希望有大神指导 |
» 猜你喜欢
奇怪,两个人的filecode固定段从头到尾一模一样
已经有4人回复
8月时间戳变的,举个手。玩一下,释放压力
已经有13人回复
国自然结果
已经有6人回复
我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人
已经有9人回复
好奇怪的filecode
已经有5人回复
应该是93bebmhtak前后十一个字符比较关键
已经有21人回复
综述论文作为代表作会不会影响评审专家的印象分?
已经有10人回复
据悉今年马上要出结果了
已经有8人回复
关于代码变化问题,想知道的进来
已经有23人回复
FileCode能看出啥?
已经有6人回复


2楼2016-03-05 18:33:19










回复此楼
100