@@ -138,17 +138,17 @@ contains
138138
139139 if (p > 0 ) then
140140 ! Compute 3D waves with phase shifts.
141- call s_instability_wave(2 * pi* 4.0 / Ldomain, 2 * pi* 4.0_wp / Ldomain, wave_tmp, 2 * pi* 11._wp / 31._wp )
141+ call s_instability_wave(2 * pi* 4.0_wp / Ldomain, 2 * pi* 4.0_wp / Ldomain, wave_tmp, 2 * pi* 11._wp / 31._wp )
142142 wave2 = wave2 + wave_tmp
143- call s_instability_wave(2 * pi* 2.0 / Ldomain, 2 * pi* 2.0_wp / Ldomain, wave_tmp, 2 * pi* 13._wp / 31._wp )
143+ call s_instability_wave(2 * pi* 2.0_wp / Ldomain, 2 * pi* 2.0_wp / Ldomain, wave_tmp, 2 * pi* 13._wp / 31._wp )
144144 wave2 = wave2 + wave_tmp
145- call s_instability_wave(2 * pi* 1.0 / Ldomain, 2 * pi* 1.0_wp / Ldomain, wave_tmp, 2 * pi* 17._wp / 31._wp )
145+ call s_instability_wave(2 * pi* 1.0_wp / Ldomain, 2 * pi* 1.0_wp / Ldomain, wave_tmp, 2 * pi* 17._wp / 31._wp )
146146 wave2 = wave2 + wave_tmp
147- call s_instability_wave(2 * pi* 4.0 / Ldomain, - 2 * pi* 4.0_wp / Ldomain, wave_tmp, 2 * pi* 19._wp / 31._wp )
147+ call s_instability_wave(2 * pi* 4.0_wp / Ldomain, - 2 * pi* 4.0_wp / Ldomain, wave_tmp, 2 * pi* 19._wp / 31._wp )
148148 wave2 = wave2 + wave_tmp
149- call s_instability_wave(2 * pi* 2.0 / Ldomain, - 2 * pi* 2.0_wp / Ldomain, wave_tmp, 2 * pi* 23._wp / 31._wp )
149+ call s_instability_wave(2 * pi* 2.0_wp / Ldomain, - 2 * pi* 2.0_wp / Ldomain, wave_tmp, 2 * pi* 23._wp / 31._wp )
150150 wave2 = wave2 + wave_tmp
151- call s_instability_wave(2 * pi* 1.0 / Ldomain, - 2 * pi* 1.0_wp / Ldomain, wave_tmp, 2 * pi* 29._wp / 31._wp )
151+ call s_instability_wave(2 * pi* 1.0_wp / Ldomain, - 2 * pi* 1.0_wp / Ldomain, wave_tmp, 2 * pi* 29._wp / 31._wp )
152152 wave2 = wave2 + wave_tmp
153153 wave = wave + 0.15_wp * wave2
154154 end if
0 commit comments