select case(u0_use_file)
case(.True.)
open(29,file="vy.dat")
do j=pml_cell+1,Ny-pml_cell
read(29,*) (u0y(i,j),i=ypml_cell+1,Nx-ypml_cell)
end do
close(29)
case(.False.)
u0x=u0x_const
u0y=u0y_const
end select
end if
select case(n0_use_file)
case(.True.)
open(30,file="density.dat")
do j=pml_cell+1,Ny-pml_cell
read(30,*) (n0(i,j),i=ypml_cell+1,Nx-ypml_cell)
end do
close(30)
case(.False.)
n0=n0_constant
end select
kyra 发表于 2024-5-23 12:55
u0y 是否定义了?
另外
欢迎光临 Fortran Coder (http://bbs.fcode.cn/) | Powered by Discuz! X3.2 |