WRITE(1,*)x(i), y(j), rho(i,j), p(i,j), u(i,j), v(i,j), e(i,j)得到数据后导入tecplot却只会显示前三个构成的坐标关系

这可能是因为您在导入数据时没有正确指定数据的变量类型和排列顺序。在TECPLOT中，您需要使用“Data”-“Data Load”菜单选项来导入数据，并在“Variables”选项卡中指定每个变量的类型和排列顺序。

对于您提供的数据输出语句，其中x(i)、y(j)、rho(i,j)、p(i,j)、u(i,j)、v(i,j)和e(i,j)分别表示数据的x坐标、y坐标、密度、压力、x方向速度、y方向速度和内能。在TECPLOT中，您需要将这些变量正确地映射到坐标轴和场变量中，并在“Zones”选项卡中指定数据排列方式（例如按行顺序）。

如果您仍然无法正确显示所有变量，请确保您的数据输出语句中包含了所有变量，并且在TECPLOT中正确指定了它们的映射关系。

相关问题

DO j = 2, ny-1 rho_new(i,j) = rho(i,j) - dt/dx*u(i,j)*(rho(i,j)-rho(i-1,j))- dt/dy*v(i,j)*(rho(i,j)-rho(i,j-1))& + 0.5*dt/dx*u(i,j)*(1.0-u(i,j)*dt/dx)*(rho(i+1,j)-rho(i,j))- 0.5*dt/dx*u(i,j)*(1.0+u(i,j)*dt/dx)*(rho(i,j)-rho(i-1,j))& + 0.5*dt/dy*v(i,j)*(1.0-v(i,j)*dt/dy)*(rho(i,j+1)-rho(i,j))- 0.5*dt/dy*v(i,j)*(1.0+v(i,j)*dt/dy)*(rho(i,j)-rho(i,j-1))

这是一个二维的离散化的非压缩流体连续性方程，用于求解流体的密度分布。其中 rho_new 是新时刻的密度分布，rho 是当前时刻的密度分布，u 和 v 分别是 x 和 y 方向的速度分量，dt 是时间步长，dx 和 dy 分别是 x 和 y 方向的网格大小。这个方程描述了流体在时间上的演化过程，通过对速度和密度的计算，可以求解出流体在下一个时刻的状态。

if ((i > nx * 3 / 20 + 3) .and. j > (ny * 2 / 5 + 3) .and. j < (ny * 3 / 5 + 3)) then! 钝头体区 else rho_sub = U_sub1(i,j,1); u_sub = U_sub1(i,j,2) / rho_sub; v_sub = U_sub1(i,j,3) / rho_sub; P_sub = (gamma - 1) * (U_sub1(i,j,4) - 0.5 * rho_sub * (u_sub * u_sub + v_sub * v_sub)); A(i,j) = sqrt(abs(gamma * P_sub / rho_sub));

这段代码是什么意思？

这段代码是一个if-else语句块，判断给定的i和j是否在指定的区域内。如果在，则A(i,j)的值为一个常数；如果不在，则根据给定的公式计算出rho_sub、u_sub、v_sub和P_sub的值，然后计算出A(i,j)的值。其中，gamma表示气体的比热比。这段代码可能用于计算钝头体的特性线在给定区域内的速度。