| 查看: 570 | 回复: 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 真心求教,希望有大神指导 |
» 猜你喜欢
面上项目没有好文章就没希望了吗?
已经有14人回复
困死了
已经有8人回复
不知道还有没有招博士的学校了
已经有5人回复
售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有8人回复
材料博士申请
已经有5人回复
售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有9人回复
还有课题组有博士名额吗
已经有6人回复
关于水星近日点进动成因的质疑 与实证分析
已经有10人回复
博士申请
已经有3人回复


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











回复此楼