Porting S Stroteich`s attempt to add an explicit fourth order dissipa…

…tion operator. Initial investigations suggest that the feature in the current form is not adequate to stabilise the numerical stability of interest. In addition, the diffusion operator is likely called in the wrong place in the time advance (and should instead be called within advance_explicit).

Makefile.depend

@@ -210,7 +210,8 @@ fields_arrays.o: \
 	 common_types.o 
 file_utils.o: \
 	 command_line.o 
-finite_differences.o: 
+finite_differences.o: \
+	 mp.o 
 flow_shear.o: \
 	 constants.o \
 	 fields.o \
@@ -241,6 +242,17 @@ g_tofrom_h.o: \
 	 stella_layouts.o \
 	 vpamu_grids.o \
 	 zgrid.o 
+g_tofrom_h.old.o: \
+	 dist_fn_arrays.o \
+	 gyro_averages.o \
+	 kt_grids.o \
+	 physics_flags.o \
+	 run_parameters.o \
+	 species.o \
+	 stella_geometry.o \
+	 stella_layouts.o \
+	 vpamu_grids.o \
+	 zgrid.o 
 gauss_quad.o: \
 	 constants.o \
 	 file_utils.o \
@@ -261,30 +273,35 @@ gyro_averages.o: \
 	 zgrid.o 
 hyper.o: \
 	 dist_fn_arrays.o \
+	 dist_redistribute.o \
+	 extended_zgrid.o \
 	 file_utils.o \
+	 finite_differences.o \
 	 kt_grids.o \
 	 mp.o \
 	 physics_flags.o \
+	 redistribute.o \
 	 stella_geometry.o \
 	 stella_layouts.o \
 	 stella_time.o \
+	 vpamu_grids.o \
 	 zgrid.o 
 implicit_solve.o: \
+	 constants.o \
 	 dist_fn_arrays.o \
 	 extended_zgrid.o \
 	 fields.o \
 	 fields_arrays.o \
-	 g_tofrom_h.o \
 	 gyro_averages.o \
 	 job_manage.o \
 	 kt_grids.o \
 	 linear_solve.o \
 	 mp.o \
 	 neoclassical_terms.o \
 	 parallel_streaming.o \
+	 physics_flags.o \
 	 run_parameters.o \
 	 species.o \
-	 stella_geometry.o \
 	 stella_layouts.o \
 	 stella_time.o \
 	 vpamu_grids.o \
@@ -342,7 +359,11 @@ mirror_terms.o: \
 	 constants.o \
 	 dist_fn_arrays.o \
 	 dist_redistribute.o \
+	 fields.o \
+	 fields_arrays.o \
 	 finite_differences.o \
+	 g_tofrom_h.o \
+	 gyro_averages.o \
 	 job_manage.o \
 	 kt_grids.o \
 	 mp.o \
@@ -427,6 +448,29 @@ redistribute.o: \
 	 mp.o 
 response_matrix.o: \
 	 common_types.o \
+	 dist_fn_arrays.o \
+	 extended_zgrid.o \
+	 fields.o \
+	 fields_arrays.o \
+	 gyro_averages.o \
+	 implicit_solve.o \
+	 job_manage.o \
+	 kt_grids.o \
+	 linear_solve.o \
+	 mp.o \
+	 mp_lu_decomposition.o \
+	 physics_flags.o \
+	 physics_parameters.o \
+	 run_parameters.o \
+	 species.o \
+	 stella_geometry.o \
+	 stella_layouts.o \
+	 system_fortran.o \
+	 vpamu_grids.o \
+	 zgrid.o 
+response_matrix.writeout.o: \
+	 common_types.o \
+	 dist_fn_arrays.o \
 	 extended_zgrid.o \
 	 fields.o \
 	 fields_arrays.o \
@@ -437,20 +481,20 @@ response_matrix.o: \
 	 linear_solve.o \
 	 mp.o \
 	 mp_lu_decomposition.o \
-	 parallel_streaming.o \
 	 physics_flags.o \
+	 physics_parameters.o \
 	 run_parameters.o \
 	 species.o \
 	 stella_geometry.o \
 	 stella_layouts.o \
-	 stella_time.o \
 	 system_fortran.o \
 	 vpamu_grids.o \
 	 zgrid.o 
 run_parameters.o: \
 	 file_utils.o \
 	 mp.o \
 	 physics_flags.o \
+	 physics_parameters.o \
 	 text_options.o 
 sfincs_interface.o: \
 	 mp.o \
@@ -503,6 +547,7 @@ stella.o: \
 	 fields_arrays.o \
 	 file_utils.o \
 	 git_version.o \
+	 implicit_solve.o \
 	 init_g.o \
 	 job_manage.o \
 	 kt_grids.o \
@@ -568,6 +613,55 @@ stella_geometry.o: \
 	 vmec_geo.o \
 	 vmec_to_stella_geometry_interface.o \
 	 zgrid.o 
+stella_geometry.hack.o: \
+	 common_types.o \
+	 constants.o \
+	 file_utils.o \
+	 inputprofiles_interface.o \
+	 millerlocal.o \
+	 mp.o \
+	 physics_flags.o \
+	 physics_parameters.o \
+	 text_options.o \
+	 vmec_geo.o \
+	 vmec_to_stella_geometry_interface.o \
+	 zgrid.o 
+stella_hack.o: \
+	 dissipation.o \
+	 dist_fn.o \
+	 dist_fn_arrays.o \
+	 dist_redistribute.o \
+	 extended_zgrid.o \
+	 fields.o \
+	 fields_arrays.o \
+	 file_utils.o \
+	 git_version.o \
+	 implicit_solve.o \
+	 init_g.o \
+	 job_manage.o \
+	 kt_grids.o \
+	 mirror_terms.o \
+	 mp.o \
+	 multibox.o \
+	 parallel_streaming.o \
+	 physics_flags.o \
+	 physics_parameters.o \
+	 ran.o \
+	 redistribute.o \
+	 response_matrix.o \
+	 run_parameters.o \
+	 sources.o \
+	 species.o \
+	 stella_diagnostics.o \
+	 stella_geometry.o \
+	 stella_layouts.o \
+	 stella_save.o \
+	 stella_time.o \
+	 stella_transforms.o \
+	 time_advance.o \
+	 volume_averages.o \
+	 vpamu_grids.o \
+	 zgrid.o 
 stella_io.o: \
 	 convert.o \
 	 dist_fn_arrays.o \
@@ -600,7 +694,7 @@ stella_save.o: \
 	 vpamu_grids.o \
 	 zgrid.o 
 stella_time.o: \
-         mp.o 
+	 mp.o 
 stella_transforms.o: \
 	 fft_work.o \
 	 physics_flags.o \
@@ -621,6 +715,7 @@ time_advance.o: \
 	 file_utils.o \
 	 finite_differences.o \
 	 flow_shear.o \
+	 g_tofrom_h.o \
 	 gyro_averages.o \
 	 hyper.o \
 	 implicit_solve.o \