call freesurface(npts_w,p_w(1:npts_w,1:3),
*av_new(npts+1:npts+npts_w,1:3),det_t,pot_w(1:npts_w),
*itype(1:npts_w), phi_x(npts+1:npts+npts_w)
*av_new(npts+1:npts+npts_w,1:3),det_t,pot_w(1:npts_w),
*itype(1:npts_w), phi_x(npts+1:npts+npts_w)
*,phi_y(npts+1:npts+npts_w),phi_z(npts+1:npts+npts_w))
在CVF中上面的函数能够正常运行,而在IVF中报错。因为phi_x被定义为phi_x(npts).