1) !---------------------------------------------------------------------------! !Restriction operators,f2c,从更多网格到少量网格,所以是限制 !还有多种限制器: !(1)简单的直接插值。 !(2)取前后左右4个点,加上中间点*4取平均值。5点格式 !(3)9点格式 !---------------------------------------------------------------------------! 2) 插值完了后,要处理一下f的边界 !update boundaries,还要进行边界更新,对的。这个好像忘记了 do i=0,nxh f(i,0) = r(2*i,0) f(i,nyh) = r(2*i,nyf) end do do j=0,nyh f(0,j) = r(0,2*j) f(nxh,j) = r(nxf,2*j) end do |
说点什么...