| 查看: 584 | 回复: 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 真心求教，希望有大神指导 |
回复此楼» 猜你喜欢
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
5个%2F，非固定段22个字母
已经有13人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有7人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
2楼2016-03-05 18:33:19
100