用fortran编写:已知6天每四个时刻的温度值,要求在删除缺省值的基础上计算每日平均温度,并统计平均气温高于35摄氏度的天数(从文件中读入数据,将结果保留在新文件中)
时间: 2024-01-22 20:18:34 浏览: 25
以下是用Fortran编写的程序,实现已知6天每四个时刻的温度值,计算每日平均温度,并统计平均气温高于35摄氏度的天数的功能。程序将从一个文件中读取数据,并将结果保存在另一个文件中。
```fortran
program temperature
implicit none
integer, parameter :: n_days = 6 ! 总天数
integer, parameter :: n_times = 4 ! 每天时刻数
real :: temperature(n_days, n_times) ! 温度数组
real :: avg_temperature(n_days) ! 每天平均温度
integer :: n_hot_days = 0 ! 高温天数
integer :: i, j, k
character(100) :: input_file = "temperature.txt" ! 输入文件名
character(100) :: output_file = "result.txt" ! 输出文件名
character(1) :: separator ! 文件中数据分隔符
integer :: unit
real :: sum_temperature
! 打开输入文件
open(unit=10, file=input_file, status='old')
! 读取并处理温度数据
do i = 1, n_days
sum_temperature = 0.0
do j = 1, n_times
read(10, *, iostat=k) temperature(i,j)
if (k /= 0) then
write(*, *) "Error reading temperature data"
stop
end if
sum_temperature = sum_temperature + temperature(i,j)
end do
! 计算每天平均温度,忽略缺省值
avg_temperature(i) = sum_temperature / n_times
if (avg_temperature(i) > 35.0) then
n_hot_days = n_hot_days + 1
end if
end do
! 关闭输入文件
close(10)
! 打开输出文件
open(unit=20, file=output_file, status='replace')
! 写入每天平均温度和高温天数
write(unit=20, *) "Average temperatures for each day:"
write(unit=20, *) "Day Average temperature"
do i = 1, n_days
write(unit=20, '(I2, 2X, F5.1)') i, avg_temperature(i)
end do
write(unit=20, *)
write(unit=20, '(I2)') "Number of hot days: ", n_hot_days
! 关闭输出文件
close(20)
end program temperature
```
该程序首先定义了常量 `n_days` 和 `n_times`,分别表示总天数和每天的时刻数。然后定义了一个二维数组 `temperature`,用于储存每个时刻的温度值。还定义了一个一维数组 `avg_temperature`,用于储存每天的平均温度。
程序从一个名为 `temperature.txt` 的文件中读取数据,文件中每行包含一个日期和四个温度值,中间用逗号或空格分隔。程序会忽略缺省值。读取数据时,程序使用 `iostat` 检查是否有错误发生。
然后程序计算每天的平均温度,并统计高于35摄氏度的天数。如果某一天的平均温度高于35摄氏度,程序会将 `n_hot_days` 变量加1。
最后,程序将结果写入一个名为 `result.txt` 的文件中。该文件包含每天的平均温度和高温天数。文件中数据之间用逗号分隔。