问题1：函数的实参和虚参的类型不一致。这个问题你给出的代码上并看不出来具体问题所在，因为很多东西（RHO_T1,DELTA_X,COLLECTION等）没有给出定义，我无法知道它是数组，还是函数，如果是函数，每个参数代表什么意思。问题2：有一个子程序，你当做函数在用了。子程序subroutine没有返回值，而函数function有返回值。
