changes ordering (in energy.dat output)

Author: danielpeter

src/specfem3D/compute_energy.f90

@@ -453,8 +453,8 @@ subroutine compute_energy_generic_slow()
   ! write the total to disk from the main
   if (myrank == 0) then
     total_energy_glob = kinetic_energy_glob + potential_energy_glob
-    ! format: #timestep #kinetic_energy #potential_energy #total_energy
-    write(IOUT_ENERGY,*) it,sngl(kinetic_energy_glob),sngl(potential_energy_glob),sngl(total_energy_glob)
+    ! format: #timestep #total_energy #kinetic_energy #potential_energy
+    write(IOUT_ENERGY,*) it,sngl(total_energy_glob),sngl(kinetic_energy_glob),sngl(potential_energy_glob)
     ! flush to file
     flush(IOUT_ENERGY)
   endif
@@ -876,8 +876,8 @@ subroutine compute_energy_fast_Deville()
   ! write the total to disk from the main
   if (myrank == 0) then
     total_energy_glob = kinetic_energy_glob + potential_energy_glob
-    ! format: #timestep #kinetic_energy #potential_energy #total_energy
-    write(IOUT_ENERGY,*) it,sngl(kinetic_energy_glob),sngl(potential_energy_glob),sngl(total_energy_glob)
+    ! format: #timestep #total_energy #kinetic_energy #potential_energy
+    write(IOUT_ENERGY,*) it,sngl(total_energy_glob),sngl(kinetic_energy_glob),sngl(potential_energy_glob)
     ! flush to file
     flush(IOUT_ENERGY)
   endif

src/specfem3D/iterate_time.F90

@@ -100,7 +100,7 @@ subroutine iterate_time()
   if (OUTPUT_ENERGY .and. myrank == 0) then
     open(unit=IOUT_ENERGY,file=trim(OUTPUT_FILES)//'energy.dat',status='unknown',action='write')
     ! format: #timestep #kinetic_energy #potential_energy #total_energy
-    write(IOUT_ENERGY,*) "#timestep  #kinetic_energy  #potential_energy  #total_energy"
+    write(IOUT_ENERGY,*) "#timestep  #total_energy  #kinetic_energy  #potential_energy"
     flush(IOUT_ENERGY)
   endif
 

src/specfem3D/iterate_time_undoatt.F90

@@ -240,7 +240,7 @@ subroutine iterate_time_undoatt()
   if (OUTPUT_ENERGY .and. myrank == 0) then
     open(unit=IOUT_ENERGY,file=trim(OUTPUT_FILES)//'energy.dat',status='unknown',action='write')
     ! format: #timestep #kinetic_energy #potential_energy #total_energy
-    write(IOUT_ENERGY,*) "#timestep  #kinetic_energy  #potential_energy  #total_energy"
+    write(IOUT_ENERGY,*) "#timestep  #total_energy  #kinetic_energy  #potential_energy"
     flush(IOUT_ENERGY)
   endif