From 018be4729ecf3baff2c418709597198156c39fd4 Mon Sep 17 00:00:00 2001
From: Evgeny Mortikov <evgeny.mortikov@gmail.com>
Date: Sat, 21 Sep 2024 02:08:08 +0300
Subject: [PATCH] io implementation update
---
config-ex/config-mosaic.txt | 2 +-
srcF/sfx_main.f90 | 14 +++++++-------
2 files changed, 8 insertions(+), 8 deletions(-)
diff --git a/config-ex/config-mosaic.txt b/config-ex/config-mosaic.txt
index 64fcdce..8750ecf 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 4f22f51..dcabad5 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
--
GitLab