JGonzalez/vlaplex
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Fix fort.f40 mistake

Browse source
This commit is contained in:
JHendrikx 2025-02-24 13:57:18 +01:00
commit 03dc8671a1
1 changed files with 14 additions and 12 deletions
Download patch file Download diff file Expand all files Collapse all files

26
moduleOutput.f90
View file

@ -13,6 +13,7 @@ module output
integer, parameter:: dataF_id = 30 integer, parameter:: dataF_id = 30
integer, parameter:: dataPhi_id = 40 integer, parameter:: dataPhi_id = 40
integer, parameter:: dataCum_id = 50 integer, parameter:: dataCum_id = 50
integer, parameter:: dataTime_id = 60
character(len=*), parameter :: formatInt = 'I10' character(len=*), parameter :: formatInt = 'I10'
character(len=7), parameter:: formatFloat = 'ES0.6e3' character(len=7), parameter:: formatFloat = 'ES0.6e3'
character(len=3), parameter:: formatSep = '","' character(len=3), parameter:: formatSep = '","'
@ -71,8 +72,8 @@ module output
open(unit=dataF_id, file=pathOutput // filename) open(unit=dataF_id, file=pathOutput // filename)
write(dataF_id, '(A)') "t (s)" write(dataF_id, '(A)') "t (s)"
write(dataF_id, '('//formatFloat//')') t*dt*t_ref write(dataF_id, '('//formatFloat//')') t*dt*t_ref
write(dataPhi_id, '(A)') "Z" write(dataF_id, '(A)') "Z"
write(dataPhi_id, '('//formatFloat//')') Z_list(j) write(dataF_id, '('//formatFloat//')') Z_list(j)
write(myfmt, "(I0)") nr write(myfmt, "(I0)") nr
myfmt = '(A,' // trim(myfmt) // '(' // formatSep // ',' // formatFloat // '))' myfmt = '(A,' // trim(myfmt) // '(' // formatSep // ',' // formatFloat // '))'
write(dataF_id, myfmt) "v (m/s) / r (m)", r*L_ref write(dataF_id, myfmt) "v (m/s) / r (m)", r*L_ref
@ -199,26 +200,27 @@ module output
end subroutine writeOutputBoundary end subroutine writeOutputBoundary
subroutine writeOutputTime(t, time) subroutine writeOutputTime(t, time, bins)
integer, intent(in):: t integer, intent(in):: t
real(dp), intent(in):: time real(dp), intent(in):: time
real(dp), intent(in):: bins
character(len=8), parameter:: filename = 'time.csv' character(len=8), parameter:: filename = 'time.csv'
logical:: res logical:: res
inquire(file=pathOutput // filename, exist=res) inquire(file=pathOutput // filename, exist=res)
if (.not. res) then if (.not. res) then
write (*, '(A, A)') 'Writing: ', filename write (*, '(A, A)') 'Writing: ', filename
open(unit=dataBC_id, file=pathOutput // filename, action='write', position='append') open(unit=dataTime_id, file=pathOutput // filename, action='write', position='append')
write(dataBC_id, '(A,1(' // formatSep // ',A))') 'timestep', 'duration (s)' write(dataTime_id, '(A,2(' // formatSep // ',A))') 'timestep', 'duration (s)', '#bins'
close(dataBC_id) close(dataTime_id)
end if end if
open(unit=dataBC_id, file=pathOutput // filename, action='write', position='append') open(unit=dataTime_id, file=pathOutput // filename, action='write', position='append')
write(dataBC_id, '(' // formatInt // ',' // formatSep // ',' // formatFloat // ')') & write(dataTime_id, '(' // formatInt // ',2('// formatSep // ',' // formatFloat // '))') &
t, time t, time, bins
close(dataBC_id) close(dataTime_id)
end subroutine writeOutputTime end subroutine writeOutputTime
@ -288,8 +290,8 @@ module output
open(unit=dataCum_id, file=pathOutput // filename) open(unit=dataCum_id, file=pathOutput // filename)
write(dataCum_id, '(A)') "t (s)" write(dataCum_id, '(A)') "t (s)"
write(dataCum_id, '('//formatFloat//')') t*dt*t_ref write(dataCum_id, '('//formatFloat//')') t*dt*t_ref
write(dataPhi_id, '(A)') "Z" write(dataCum_id, '(A)') "Z"
write(dataPhi_id, '('//formatFloat//')') Z_list(j) write(dataCum_id, '('//formatFloat//')') Z_list(j)
write(myfmt, "(I0)") 1 write(myfmt, "(I0)") 1
myfmt = '(A,' // trim(myfmt) // '(' // formatSep // ',' // formatFloat // '))' myfmt = '(A,' // trim(myfmt) // '(' // formatSep // ',' // formatFloat // '))'
write(dataCum_id, myfmt) "v (m/s) / r (m)", r*L_ref write(dataCum_id, myfmt) "v (m/s) / r (m)", r*L_ref