SUBROUTINE save_Hollmann_params(params)
TYPE(KORC_PARAMS), INTENT(IN) :: params
CHARACTER(MAX_STRING_LENGTH) :: filename
CHARACTER(MAX_STRING_LENGTH) :: gname
CHARACTER(MAX_STRING_LENGTH), DIMENSION(:), ALLOCATABLE :: attr_array
CHARACTER(MAX_STRING_LENGTH) :: dset
CHARACTER(MAX_STRING_LENGTH) :: attr
INTEGER(HID_T) :: h5file_id
INTEGER(HID_T) :: group_id
INTEGER :: h5error
REAL(rp) :: units
if (params%mpi_params%rank .EQ. 0) then
filename = TRIM(params%path_to_outputs) // "experimental_distribution_parameters.h5"
call h5fcreate_f(TRIM(filename), H5F_ACC_TRUNC_F, h5file_id, h5error)
gname = "pdf_params"
call h5gcreate_f(h5file_id, TRIM(gname), group_id, h5error)
dset = TRIM(gname) // "/max_pitch_angle"
attr = "Maximum pitch angle in PDF (degrees)"
call save_to_hdf5(h5file_id,dset,h_params%max_pitch_angle,attr)
dset = TRIM(gname) // "/min_pitch_angle"
attr = "Minimum pitch angle in PDF (degrees)"
call save_to_hdf5(h5file_id,dset,h_params%min_pitch_angle,attr)
dset = TRIM(gname) // "/min_energy"
attr = "Minimum energy in PDF (eV)"
units = 1.0_rp/C_E
call save_to_hdf5(h5file_id,dset,units*h_params%min_energy,attr)
dset = TRIM(gname) // "/max_energy"
attr = "Maximum energy in PDF (eV)"
units = 1.0_rp/C_E
call save_to_hdf5(h5file_id,dset,units*h_params%max_energy,attr)
dset = TRIM(gname) // "/max_g"
attr = "Maximum momentum in PDF (me*c)"
call save_to_hdf5(h5file_id,dset,h_params%max_g,attr)
dset = TRIM(gname) // "/min_g"
attr = "Maximum momentum in PDF (me*c)"
call save_to_hdf5(h5file_id,dset,h_params%min_g,attr)
dset = TRIM(gname) // "/max_sampling_g"
attr = "Maximum momentum in PDF (me*c)"
call save_to_hdf5(h5file_id,dset,h_params%max_sampling_g,attr)
dset = TRIM(gname) // "/min_sampling_g"
attr = "Maximum momentum in PDF (me*c)"
call save_to_hdf5(h5file_id,dset,h_params%min_sampling_g,attr)
dset = TRIM(gname) // "/Zeff"
attr = "Effective atomic number of ions."
call save_to_hdf5(h5file_id,dset,h_params%Zeff,attr)
dset = TRIM(gname) // "/sigmaZeff"
attr = "Effective atomic number of ions."
call save_to_hdf5(h5file_id,dset,h_params%sigma_Z,attr)
dset = TRIM(gname) // "/E"
attr = "Parallel electric field in (Ec)"
call save_to_hdf5(h5file_id,dset,h_params%E,attr)
dset = TRIM(gname) // "/sigmaE"
attr = "Parallel electric field in (Ec)"
call save_to_hdf5(h5file_id,dset,h_params%sigma_E,attr)
dset = TRIM(gname) // "/lambda"
attr = "Wavelength used when PDF is weighted with the distribution of synchrotron radiation."
call save_to_hdf5(h5file_id,dset,h_params%lambda,attr)
dset = TRIM(gname) // "/Bo"
attr = "Magnetic field used when PDF is weighted with the distribution of synchrotron radiation."
call save_to_hdf5(h5file_id,dset,h_params%Bo,attr)
call h5gclose_f(group_id, h5error)
call h5fclose_f(h5file_id, h5error)
end if
END SUBROUTINE save_Hollmann_params