一个循环嵌套的设计问题.

再见大麋鹿

鄙人 Fortran 小白, 此问题查阅大量资料后无果, 且时间紧张, 故前来请教.
我的程序的思路是:
做两个自变量: time, voltage 的循环
就是先对 time 做一个循环，然后在每一个 time 数值做一个 voltage 的循环
最后得出2个自变量的函数 current=current(voltage, time)

而且我在给定 voltage 的时候是成功的, 也就是说, 原本只有一个 time 的循环, 现在我要在 time 的循环下加上一个 voltage 的循环.

原代码循环部分如下:

[Fortran] 纯文本查看 复制代码

查看代码

打印?

001	time = 0.0d0

002

do

003	if (time .ge. time_end) exit

004

005	if (time <= time_center) then

006	e_rabi_0 = e_rabi_i*exp(-(time-time_center)**2/time_width**2)       ! = exp(-4*(time-time_center)**2/time_center**2)

007

else

008	e_rabi_0 = e_rabi_i

009

end if

010	e_rabi_1 = e_rabi_0

011

012	call exci_rates(de,sq_fcf,gam,kop,e_rabi_0,e_rabi_1)

013	do mc = 0, 1

014	do me = 0, 1

015	do mv = 0, mv_max

016	do nc = 0, 1

017	do ne = 0, 1

018	do nv = 0, mv_max

019	k_ref(mc,me,mv,nc,ne,nv) = kivr(mc,me,mv,nc,ne,nv) + kop(mc,me,mv,nc,ne,nv) + klu(mc,me,mv,nc,ne,nv)

020

end do

021

end do

022

end do

023

end do

024

end do

025

end do

026

027	voltage = voltage_give

028

029

030	call mol_lead_rates(voltage,de,v_l,v_r,sq_fcf,klw,klw_l)

031	do mc = 0, 1

032	do me = 0, 1

033	do mv = 0, mv_max

034	do nc = 0, 1

035	do ne = 0, 1

036	do nv = 0, mv_max

037	k_tot(mc,me,mv,nc,ne,nv) = k_ref(mc,me,mv,nc,ne,nv) + klw(mc,me,mv,nc,ne,nv)

038

end do

039

end do

040

end do

041

end do

042

end do

043

end do

044

045	if (time == 0.0d0) then

046	do mc = 0, 1

047	do me = 0, 1

048	do mv = 0, mv_max

049	pop(mc,me,mv) = 0.d0

050

end do

051

end do

052

end do

053	pop(0,0,0) = 1.d0

054

055	stop_time = time_final

056	100 do mc = 0, 1

057	do me = 0, 1

058	do mv = 0, mv_max

059	pop_old(mc,me,mv) = pop(mc,me,mv)

060

end do

061

end do

062

end do

063	call multi_to_single(pop_old,fin)

064	call rk4(fin,time,time_step,fout,k_tot)

065	call single_to_multi(fout,pop)

066

067	time = time + time_step

068	if (time < stop_time) goto 100

069

070	error_pop = 0.d0

071	do mc = 0, 1

072	do me = 0, 1

073	do mv = 0, mv_max

074	error_pop = error_pop + (pop(mc,me,mv) - pop_old(mc,me,mv))**2

075

end do

076

end do

077

end do

078	error_pop = SQRT(error_pop)/length_coeff

079	if (error_pop > max_error) then

080	stop_time = stop_time + 0.5d0*time_final

081

end if

082	if (error_pop > max_error) goto 100

083	time = 0.0d0

084	time = time + time_step

085

else

086	do mc = 0, 1

087	do me = 0, 1

088	do mv = 0, mv_max

089	pop_old(mc,me,mv) = pop(mc,me,mv)

090

end do

091

end do

092

end do

093

094	call multi_to_single(pop_old,fin)

095	call rk4(fin,time,time_step,fout,k_tot)

096	call single_to_multi(fout,pop)

097	time = time + time_step

098

end if

099

100	total_pop = 0.d0

101	mean_vib = 0.d0

102	do mc = 0, 1

103	do me = 0, 1

104	pop_el(mc,me) = 0.d0

105

end do

106

end do

107

108	do mv = 0, mv_max

109	pop_vib(mv) = 0.d0

110

end do

111	do mc = 0, 1

112	do me = 0, 1

113	do mv = 0, mv_max

114	total_pop = total_pop + pop(mc,me,mv)

115	pop_el(mc,me) = pop_el(mc,me) + pop(mc,me,mv)

116

end do

117

end do

118

end do

119

120	mean_vib=0.0d0

121	do mc = 0, 1

122	do me = 0, 1

123	do mv = 0, mv_max

124	mean_vib = mean_vib + mv*pop(mc,me,mv)

125

end do

126

end do

127

end do

128

129	i_photon = 0.d0

130	do mc = 0, 1

131	do mv = 0, mv_max

132	do nv = 0, mv_max

133	i_photon = i_photon + klu(mc,1,mv,mc,0,nv)*pop(mc,1,mv)

134

end do

135

end do

136

end do

137

138	current = 0.0d0

139	do me = 0, 1

140	do mv = 0, mv_max

141	do ne = 0, 1

142	do nv = 0, mv_max

143	current = current + klw_l(0,me,mv,1,ne,nv)*pop(0,me,mv) - klw_l(1,ne,nv,0,me,mv)*pop(1,ne,nv)

144

end do

145

end do

146

end do

147

end do

148

149	do mc = 0, 1

150	do me = 0, 1

151	do nc = 0, 1

152	do ne = 0, 1

153	c_in(mc,me,nc,ne) = 0.0d0

154	c_out(mc,me,nc,ne) = 0.0d0

155

end do

156

end do

157

end do

158

end do

159

160	do mv = 0, mv_max

161	do nv = 0, mv_max

162	c_in(0,0,1,0) = c_in(0,0,1,0) + klw_l(0,0,mv,1,0,nv)*pop(0,0,mv)

163	c_in(0,0,1,1) = c_in(0,0,1,1) + klw_l(0,0,mv,1,1,nv)*pop(0,0,mv)

164	c_in(0,1,1,0) = c_in(0,1,1,0) + klw_l(0,1,mv,1,0,nv)*pop(0,1,mv)

165	c_in(0,1,1,1) = c_in(0,1,1,1) + klw_l(0,1,mv,1,1,nv)*pop(0,1,mv)

166

167	c_out(1,0,0,0) = c_out(1,0,0,0) - klw_l(1,0,nv,0,0,mv)*pop(1,0,nv)

168	c_out(1,0,0,1) = c_out(1,0,0,1) - klw_l(1,0,nv,0,1,mv)*pop(1,0,nv)

169	c_out(1,1,0,0) = c_out(1,1,0,0) - klw_l(1,1,nv,0,0,mv)*pop(1,1,nv)

170	c_out(1,1,0,1) = c_out(1,1,0,1) - klw_l(1,1,nv,0,1,mv)*pop(1,1,nv)

171

end do

172

end do

173

174	!!!!    Output

175	write (*,*) time, voltage, current*243.4d0

176	write (40,*) time, voltage, current*243.4d0

177

178

end do

这时候代码是可以正确运行的, 但是我想取若干个:

[Fortran] 纯文本查看 复制代码

查看代码

打印?

1	voltage = voltage_give

所以我把这行改成了:
voltage = voltage_min
然后 Output 的前一行加上:

[Fortran] 纯文本查看 复制代码

查看代码

打印?

1	if (num_voltage > 0) then      voltage = voltage + (voltage_max - voltage_min)/num_voltageend if

再然后, 末尾再加一个:

[Fortran] 纯文本查看 复制代码

查看代码

打印?

1

end do

编译没错, 但是, 但是并不是每个 time 数值都做一个voltage的循环, 而是同时作循环:

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	1.000000000000000E-002   40.0000000000000       1.988160357654218E-016

02	2.000000000000000E-002   80.0000000000000       1.364425784870382E-004

03	3.000000000000000E-002   120.000000000000       1.370837425057141E-004

04	4.000000000000000E-002   160.000000000000       2.083606552928942E-004

05	5.000000000000000E-002   200.000000000000       3.832251385343335E-004

06	6.000000000000000E-002   240.000000000000       3.846367630381155E-004

07	7.000000000000001E-002   280.000000000000       4.131201522155061E-004

08	8.000000000000000E-002   320.000000000000       5.825466767173463E-004

09	9.000000000000000E-002   360.000000000000       5.852940467602035E-004

10	0.100000000000000        400.000000000000       5.948948519956353E-004

11	0.110000000000000        440.000000000000       7.259287371840456E-004

12	0.120000000000000        480.000000000000       7.304897383455301E-004

13

……

14	38.1800000000010        2000.00000000000        74.8222446587716

15

9

按理说应该是第一列先都是 1.000000000000000E-002 第二列从 00.0000000000000 每隔 40 数到 2000.000000000000   
然后第一列变成 2.000000000000000E-002 第二列重新从 00.0000000000000 每隔 40 数到 2000.000000000000   
……


近期是论文dead_line, 心态简直好到爆炸! T.T

最后感谢大家的帮助！

Jackdaw

你可以吧这部分voltage 循环改成子程序，然后让这个子程序对time进行循环调用，输出结果，time作为参数每次传入一个数值