Skip to content
Snippets Groups Projects
Commit 018be472 authored by Evgeny Mortikov's avatar Evgeny Mortikov
Browse files

io implementation update

parent 493f2299
No related branches found
No related tags found
No related merge requests found
......@@ -7,7 +7,7 @@ dataset {
# --- redefine dataset
# filename = "data/mosaic.txt";
# surface_type = "land";
# h = 3.8; # meauserement height [m]
# h = 3.8; # measurement height [m]
# z0_m = 0.01; # aerodynamic roughness [m]
# z0_h = -1; # no prescribed value
# -> using scheme assigned by surface type
......
......@@ -50,7 +50,7 @@ subroutine run_dataset(filename_out, dataset, model)
! local variables
! --------------------------------------------------------------------------------
integer :: i
integer :: status
integer :: io, status
! --------------------------------------------------------------------------------
......@@ -66,7 +66,7 @@ subroutine run_dataset(filename_out, dataset, model)
!< @brief define number of cells
open(32, file = dataset%filename, iostat = status, status ='old')
open(newunit = io, file = dataset%filename, iostat = status, status ='old')
if (status /= 0) then
write(*, *) ' FAILURE! > unable to open file: ', trim(dataset%filename)
return
......@@ -75,12 +75,12 @@ subroutine run_dataset(filename_out, dataset, model)
num = 0
status = 0
do while (status.eq.0)
read (32, *, iostat = status) meteo_cell%U, meteo_cell%dT, meteo_cell%Tsemi, meteo_cell%dQ
read (io, *, iostat = status) meteo_cell%U, meteo_cell%dT, meteo_cell%Tsemi, meteo_cell%dQ
num = num + 1
enddo
num = num - 1
close(32)
close(io)
! --- print number of elements in dataset
write(*, '(a, g0)') ' size = ', num
......@@ -100,13 +100,13 @@ subroutine run_dataset(filename_out, dataset, model)
meteo_cell%z0_m = dataset%z0_m
!< @brief read input data
open(32, file = dataset%filename, iostat = status, status = 'old')
open(newunit = io, file = dataset%filename, iostat = status, status = 'old')
if (status /= 0) then
write(*, *) ' FAILURE! > unable to open file: ', trim(dataset%filename)
return
end if
do i = 1, num
read(32, *) meteo_cell%U, meteo_cell%dT, meteo_cell%Tsemi, meteo_cell%dQ
read(io, *) meteo_cell%U, meteo_cell%dT, meteo_cell%Tsemi, meteo_cell%dQ
meteo%h(i) = meteo_cell%h
meteo%U(i) = meteo_cell%U
......@@ -115,7 +115,7 @@ subroutine run_dataset(filename_out, dataset, model)
meteo%dQ(i) = meteo_cell%dQ
meteo%z0_m(i) = meteo_cell%z0_m
enddo
close(32)
close(io)
!< @brief calling flux module
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment