diff --git a/config-ex/config-mosaic.txt b/config-ex/config-mosaic.txt index 64fcdce330415092ab50a3f267ad8eb9192838f5..8750ecfc886f2e570a068364ac303ad3d4b1eb4f 100644 --- a/config-ex/config-mosaic.txt +++ b/config-ex/config-mosaic.txt @@ -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 diff --git a/srcF/sfx_main.f90 b/srcF/sfx_main.f90 index 4f22f5104d77f9a89b9367386bda24a79e4787f4..dcabad5f2ab4f35008ea870ca7b0f08ce1bdd13b 100644 --- a/srcF/sfx_main.f90 +++ b/srcF/sfx_main.f90 @@ -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