module params_impy_coronal_m !! Parameters related with cornal impurity use precision_grillix_m, only : GP use error_handling_grillix_m, only: handle_error use status_codes_grillix_m, only : GRILLIX_ERR_NAMELIST use screen_io_m, only : get_stdout use technical_constants_m, only : PATHLEN_MAX implicit none real(GP), protected :: impy_concentration = - 1.0_GP !! Impurity concentration (fraction of electron density) !! <= 0 disables functionality character(len=PATHLEN_MAX), protected :: impy_path = 'coronal_nitrogen.dat' !! Filepath to radiation rate coefficient table integer, protected :: impy_tab_length = 121 !! Number of data points in the table public :: read_params_impy_coronal public :: write_params_impy_coronal namelist / impy_coronal / & impy_concentration, & impy_path, & impy_tab_length private impy_coronal contains subroutine read_params_impy_coronal(filename) !! Reads parameters related with impy_coronal character(len=*), intent(in) :: filename !! Filename, to read from integer :: io_unit, io_error character(len=256) :: io_errmsg open(newunit=io_unit, file=filename, status='old', action='read', & iostat=io_error, iomsg=io_errmsg) if (io_error /= 0) then call handle_error(io_errmsg, GRILLIX_ERR_NAMELIST, & __LINE__, __FILE__) endif read(io_unit, nml=impy_coronal, iostat=io_error, iomsg=io_errmsg) if (io_error /= 0) then call handle_error(io_errmsg, GRILLIX_ERR_NAMELIST, & __LINE__, __FILE__) endif close(io_unit) end subroutine subroutine write_params_impy_coronal(filename) !! Reads parameters related with impy_coronal character(len=*), intent(in), optional :: filename !! If present, filename where params are written to, !! if not present, writes to screen integer :: io_unit, io_error character(len=256) :: io_errmsg if (present(filename)) then open(newunit=io_unit, file=filename, status='unknown', & access='append', action='write', & iostat=io_error, iomsg=io_errmsg ) if (io_error /= 0) then call handle_error(io_errmsg, GRILLIX_ERR_NAMELIST, & __LINE__, __FILE__) endif else io_unit = get_stdout() endif write(io_unit, nml=impy_coronal, iostat=io_error, iomsg=io_errmsg) if (io_error /= 0) then call handle_error(io_errmsg, GRILLIX_ERR_NAMELIST, & __LINE__, __FILE__) endif if (present(filename)) then close(io_unit) endif end subroutine end module