Changeset 3533 for trunk/LMDZ.COMMON/libf/evolution

Timestamp:

Dec 4, 2024, 6:29:10 PM (3 weeks ago)

Author:

jbclement

Message:

PEM:
Fixing the writing format for "info_PEM.txt".
JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• info_PEM_mod.F90 (modified) (4 diffs)

trunk/LMDZ.COMMON/libf/evolution/changelog.txt

@@ -504 +504 @@
 == 04/12/2024 == JBC
 Removing unecessary module/subroutine "interpol_TI_PEM2PCM.F90" + Few small corrections/cleanings.
+
+== 04/12/2024 == JBC
+Fixing the writing format for "info_PEM.txt".

trunk/LMDZ.COMMON/libf/evolution/info_PEM_mod.F90

@@ -31 +31 @@
 logical       :: ok
 integer       :: cstat
-character(50) :: info_frmt
-character(20) :: ich1, ich2, ich3, ich4, fch1, fch2, fch3
+character(20) :: fch1, fch2, fch3
 
 !----- Code
@@ -40 +39 @@
     write(fch2,'(f'//int2str(nb_digits(n_myear) + 5)//'.4)') n_myear
     write(fch3,'(f6.4)') convert_years ! 4 digits to the right of the decimal point to respect the precision of Martian year in "launch_pem.sh"
-    write(ich1,'(i0)') iPCM
-    write(ich2,'(i0)') iPEM + 1
-    write(ich3,'(i0)') nPCM
-    write(ich4,'(i0)') nPCM_ini
-    call execute_command_line('sed -i "1s/.*/'//trim(fch1)//' '//trim(fch2)//' '//trim(fch3)//' '//trim(ich1)//' '//trim(ich2)//' '//trim(ich3)//' '//trim(ich4)//'/" info_PEM.txt',cmdstat = cstat)
+    call execute_command_line('sed -i "1s/.*/'//trim(fch1)//' '//trim(fch2)//' '//trim(fch3)//' '//int2str(iPCM)//' '//int2str(iPEM + 1)//' '//int2str(nPCM)//' '//int2str(nPCM_ini)//'/" info_PEM.txt',cmdstat = cstat)
     if (cstat > 0) then
         error stop 'info_PEM: command execution failed!'
@@ -53 +48 @@
     ! Martian date, Number of Martians years done by the PEM run, Number of Martians years done by the chainded simulation, Code of the stopping criterion
     ! The conversion ratio from Planetary years to Earth years is given in the header of the file
-    info_frmt = '(f'//int2str(nb_digits(year_bp_ini + i_myear) + 5)//'.4,'//'f'//int2str(nb_digits(i_myear_leg) + 5)//'.4,'//'f'//int2str(nb_digits(i_myear) + 5)//'.4,'//'i0)'
-    write(1,trim(adjustl(info_frmt))) year_bp_ini + i_myear, i_myear_leg, i_myear, stopPEM
+    write(1,'(f20.4,f20.4,f20.4,i20)') year_bp_ini + i_myear, i_myear_leg, i_myear, stopPEM
     close(1)
 else
@@ -67 +61 @@
 ! Function to convert an integer into a string
 
-integer, intent(in) :: i
-character(20)       :: str
+integer, intent(in)       :: i
+character(20)             :: str_tmp
+character(:), allocatable :: str
 
-if (nb_digits(real(i)) > len(str)) error stop 'int2str [info_PEM_mod]: invalid integer for conversion!'
-write(str,'(i0)') i
-str = trim(adjustl(str))
+if (nb_digits(real(i)) > len(str_tmp)) error stop 'int2str [info_PEM_mod]: invalid integer for conversion!'
+write(str_tmp,'(i0)') i
+str = trim(adjustl(str_tmp))
 
 END FUNCTION int2str