求助大佬怎么样用fortran正确合成一个新变量dat

Jury

最近在合成dat文件上有些难受，有两个dat，u.dat和v.dat，现在急着把他们合成ws.dat，但是合出来的画图好像有问题，大佬们可以看下是不是哪里有问题。fortran代码

program main
implicit none      
integer, parameter :: nx=321, ny=221, nt=36
real, dimension(nx, ny, nt) :: u10,v10,ws
integer :: ii, jj, tt, irec
open(11,file='E:\shuju\ltrend_scripts\u_yr.dat',status='old',form='unformatted',access='direct',recl=nx*ny*4)
irec = 1
do tt = 1, nt
     read(11,rec=irec) ((u10(ii,jj,tt),ii=1,nx),jj=1,ny)   
end do
close(11)   
open(12,file='E:\shuju\ltrend_scripts\v_yr.dat',status='old',form='unformatted',access='direct', recl=nx*ny*4)
irec = 1
do tt = 1, nt
     read(12,rec=irec) ((v10(ii,jj,tt),ii=1,nx),jj=1,ny)
     irec = irec + 1   
end do
close(12)
open(13,file='E:\shuju\ltrend_scripts\wind_yr.dat',status='replace',form='unformatted',access='direct',recl=nx*ny*4)
ws = sqrt(u10**2+v10**2)
irec = 1
do tt = 1, nt
     write(13,rec=irec) ((ws(ii,jj,tt),ii=1,nx),jj=1,ny)
     irec = irec + 1
end do
close(13)
end program

ctl描述代码
dset E:\shuju\ltrend_scripts\wind_yr.dat
undef  -9.99e+08
title  ws
xdef  321  linear  35  0.25
ydef  221  linear  -30  0.25
zdef  1  levels 0
tdef  36  linear  01JAN1979   1yr
vars 1
ws 0 99
endvars


一年一年的画的画，会发现有几年画不出来，是数据缺失了吗？

苏拉苏拉

可能是ctl配的缺测值不对，换一个缺测值

Jury

苏拉苏拉 发表于 2020-5-14 16:14
可能是ctl配的缺测值不对，换一个缺测值

我今天直接sdfopen打开nc，逐年画，u的没问题，v的有些出不来。是不是v.nc原数据有问题

苏拉苏拉

Jury 发表于 2020-5-14 20:35
我今天直接sdfopen打开nc，逐年画，u的没问题，v的有些出不来。是不是v.nc原数据有问题

出不来是什么意思？直接没有数据（图）？

Jury

苏拉苏拉 发表于 2020-5-15 14:40
出不来是什么意思？直接没有数据（图）？

已经搞定了，是下载的原数据缺少了数据，谢谢大佬！！！！！

苏拉苏拉

Jury 发表于 2020-5-17 16:06
已经搞定了，是下载的原数据缺少了数据，谢谢大佬！！！！！

恭喜，不客气