fortran程序调用子程序问题

一颗梅子在南京

我是一个fortran菜鸟，我有一个关于fortran调用子程序的问题想请教一下各位大神：
我的子程序为subroutine cffc(z,k,h,A,B,C)其中的A,B,C是要用于主程序计算的
在主程序中需要调用cffc五次，那就得到5组不同的A,B,C的值，因为不想设置那么多的变量，就想用二维数组来定义A(16，5),B(16，5)，C(16，5)，但是这样定义后调用子程序call cffc(z,ku,vcoeff,a(ny,1),b(ny,1),c(ny,1))的时候结果就只有A（16,1），B（16,1），C（16,1）一层的数据，其余层的都是0，难道这个只能定义成A1，A2，A3，A4，A5吗？

sunlight

第一次调用call cffc(z,ku,vcoeff,a(ny,1),b(ny,1),c(ny,1))
第二次是：call cffc(z,ku,vcoeff,a(ny，2),b(ny,2),c(ny,2))吗？

一颗梅子在南京

sunlight 发表于 2012-12-20 15:51

                               
登录/注册后可看大图

第一次调用call cffc(z,ku,vcoeff,a(ny,1),b(ny,1),c(ny,1))
第二次是：call cffc(z,ku,vcoeff,a(ny，2),b ...

不是
第一次调用call cffc(z,ku,vcoeff,a(ny,1),b(ny,1),c(ny,1))
第二次调用call cffc(z,kt,nteff,a(ny,2),b(ny,2),c(ny,2))
     

sunlight

我觉得，是不是第二次调用时
第二次是：call cffc(z,ku,vcoeff,a(ny，2),b(ny,2),c(ny,2))
如果一直用call cffc(z,ku,vcoeff,a(ny,1),b(ny,1),c(ny,1))
这样只能一直重复书写第一维数
而第二维数等一直没有计算写入

一颗梅子在南京

sunlight 发表于 2012-12-20 21:05

                               
登录/注册后可看大图

我觉得，是不是第二次调用时
第二次是：call cffc(z,ku,vcoeff,a(ny，2),b(ny,2),c(ny,2))
如果一直用cal ...

不是，我运行程序的感觉是调用子程序的出的值只是最后一项a(16,1)的值，不能得到整个数组的值，而如果写成call cffc(z,ku,vcoeff,a,b,c)这样的话就能得到整个a数组的值

言深深

调用的时候只需要用数组名传递就可以了
call s1(x,y,z)
子程序如下
subroutine s1(x,y,x)
其中子程序中的xyz可以使输入或者输出变量···

一颗梅子在南京

言深深 发表于 2012-12-21 14:44

                               
登录/注册后可看大图

调用的时候只需要用数组名传递就可以了
call s1(x,y,z)
子程序如下

这我知道呀，可是我要调用好几次，然后我不想设置那么多的变量，我想把调用子程序得到的A1，A2，A3...等用一个数组A（n）表示，这样好像不行

言深深

你call的时候这边call cffc(z,ku,vcoeff,a(ny,1),b(ny,1),c(ny,1))
这就用abc，不要写成数组的形式

一颗梅子在南京

言深深 发表于 2012-12-22 09:48

                               
登录/注册后可看大图

你call的时候这边call cffc(z,ku,vcoeff,a(ny,1),b(ny,1),c(ny,1))
这就用abc，不要写成数组的形式

我明白你的意思，但是我要调用好几次，那你的意思是只能写成a1,b1,c1,a2,b2,c2这样子的，我不想设那么多变量，想着对a1,a2,a3这样的能不能设一个数组a(ny,3)然后这样调用呢？不知道你有米有明白我的意思

言深深

一颗梅子在南京 发表于 2012-12-22 13:20

                               
登录/注册后可看大图

我明白你的意思，但是我要调用好几次，那你的意思是只能写成a1,b1,c1,a2,b2,c2这样子的，我不想设那么多变 ...

额····我明白你的意思了，但是你没有明白我的意思···调用数组（将所有需要调用的参数都赋值到这个数组里面）的时候，调用形式是这样的call sub(a,b,c)
其中abc可以是任何数组······比如这样一个程序
real a(100)
do i=1,100
a(i)=i
enddo
call sub(a)
end

subroutine sub(a)
do i=1,100
print*,x(i)
enddo
end

这个程序就是将主程序中的数组a，传递到子程序中，并print出来