From c9a797f0839c6e46c69e0bc8b169e4615134d08a Mon Sep 17 00:00:00 2001 From: Sergey Yurchenko Date: Tue, 28 May 2024 16:38:17 +0100 Subject: [PATCH] Modifying the state line list name to include J1 and J2 into the file name --- tran.f90 | 11 ++++++----- 1 file changed, 6 insertions(+), 5 deletions(-) diff --git a/tran.f90 b/tran.f90 index a5c3a4a..1129d79 100644 --- a/tran.f90 +++ b/tran.f90 @@ -727,7 +727,10 @@ subroutine read_eigenval(njval, jval, error) write(linelistname, '(a, 2i5)') 'states for J=', intensity%J(1:2) call iostart(trim(linelistname), states_unit) ! - filename = trim(intensity%linelist_file)//'.states' + write(j1char,'(i3)') intensity%J(1) + write(j2char,'(i3)') intensity%J(2) + ! + filename = trim(intensity%linelist_file)//'_'//trim(adjustl(j1char))//'_'//trim(adjustl(j2char))//'.states' open(unit = states_unit, action = 'write',status='new',file = filename,err=21) ! if (.false.) then @@ -736,7 +739,8 @@ subroutine read_eigenval(njval, jval, error) ! use time as a random number to create a unique name timenow = get_real_time() write(char_timenow,'(f16.2)') timenow - filename = trim(intensity%linelist_file)//'.states.'//trim(adjustl(char_timenow)) + filename = trim(intensity%linelist_file)//'_'//trim(adjustl(j1char))//'_'//trim(adjustl(j2char))//& + '.states.'//trim(adjustl(char_timenow)) open(unit = states_unit, action = 'write',status='new',file = filename) ! endif @@ -745,9 +749,6 @@ subroutine read_eigenval(njval, jval, error) write(linelistname, '(a, 2i5)') 'trans for J=', intensity%J(1:2) call iostart(trim(linelistname), trans_unit) ! - write(j1char,'(i3)') intensity%J(1) - write(j2char,'(i3)') intensity%J(2) - ! filename = trim(intensity%linelist_file)//'_'//trim(adjustl(j1char))//'_'//trim(adjustl(j2char))//'.trans' open(unit = trans_unit, action = 'write',status='new',file = filename,err=22) !