| 查看: 527 | 回复: 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 真心求教,希望有大神指导 |
» 猜你喜欢
情人节自我反思:在爱情中有过遗憾吗?
已经有6人回复
基金正文30页指的是报告正文还是整个申请书
已经有4人回复
今年春晚有几个节目很不错,点赞!
已经有6人回复
球磨粉体时遇到了大的问题,请指教!
已经有15人回复
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
江汉大学解明教授课题组招博士研究生/博士后
已经有3人回复


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













回复此楼