- 积分
- 11
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-4-16
- 最后登录
- 1970-1-1
|
GrADS
系统平台: |
fortran |
问题截图: |
- |
问题概况: |
我提取和文本文件资料是网格点形式的,输出时变成了只有一列表达的形式,如何修改? |
我看过提问的智慧: |
看过 |
自己思考时长(天): |
1 |
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
PROGRAM zzq
PARAMETER M = 63 !经度格点数
PARAMETER N = 51 !纬度格点数
PARAMETER X1 = 2001 !起始年份
PARAMETER X2 = 2001 !终止年份
PARAMETER Y1 = 1 !起始月份
PARAMETER Y2 = 1 !终止月份
real*8 :: A(N, M), C(N, M)
real*8 :: B(M)
character*30 track,outpath
character*4 year
character*2 month
integer :: i, j, k, t
track='F:\MINE\data\CDR\' !数据存放的位置
outpath='F:\MINE\result\' !结果保存的路径
A = 0
C = (X2 - X1 + 1)*(Y2 - Y1 + 1)*1.0
open(3,file=trim(outpath)//'cloud.txt',status='replace')
do k = X1, X2
write(year, "(i4.4)") k
do t = Y1, Y2
write(month, "(i2.2)") t
open(2,file=trim(track)//trim(year)//trim(month)//trim('.txt'),status='old')
do i = 1, N
read(2, *) (B(j), j = 1, M)
do j = 1, M
if(B(j) == -9999) then
C(i, j) = C(i, j) - 1.0
B(j) = 0.0
endif
A(i, j) = A(i, j) + B(j)
enddo
enddo
enddo
close(2)
enddo
do i = 1, N
do j = 1, M
write(3, *) A(i, j)/C(i, j)*0.01
enddo
enddo
end
|
|