Problem in modifying the subroutine


Click here for full thread
Clicked A Few Times
Dear Dr. Sean,
thank you very much for your kind reply.

I tried the repeat tests under no changing the source code (rt_tddft_excite.F) and *.nw input file. I could not get consistent results as shown in below.
I am troubled why such results came out.
Please help me.

Note:  
rt_tddft_excite.F = original
*.nw input file = H-GRH_SrcCode_170121A.nw. This file was attached in Jan. 25th
2017 this Post.

<< Results >>
[1]Step 1: original
 Peak E-Field = - 7.250470266E7 (V/cm)

[2]Step 2: different from Step 1
 Peak E-Field = - 6.146795755E7 (V/cm)

[3]Step 3: equal to Step 1
 Peak E-Field = - 7.250470266E7 (V/cm)

[4]Step 4: equal to Step 1
 Peak E-Field = - 7.250470266E7 (V/cm)

[5]Step 5: Field was changed from minus value to plus value.
 Peak E-Field = + 7.386817508E7 (V/cm)