fortran可以给二进制数据赋值吗

兰北

二进制数据ps（M,N），想用fortran程序，将小于0的值统一赋值为0 ，可以实现吗？
我试了一下write或者read来读取内部文件，发现不提示错误，但实现不了。
语句如下：

REAL Lon0,Lon1,Lat0,Lat1,Ld,a1
integer DAYSTRT,DAYEND
Parameter (Lon0=100.0,Lon1=130.0,Lat0=25.0,Lat1=45.0)
Parameter (Ld=1.0,N=(Lat1-Lat0)/Ld+1,M=(Lon1-Lon0)/Ld+1)
!**********DAYSTRT为开始日期，DAYEND为结束日期数**********
parameter(DAYSTRT=14,DAYEND=16)
REAL  H(M,N),PS(M,N),WSH(M,N)
     open(5,file='H:\result\200904\1.dat',FORM='UNFORMATTED',ACCESS='DIRECT',RECL=4*M*N)
     open(6,file='H:\result\200904\2.dat',FORM='UNFORMATTED',ACCESS='DIRECT',RECL=4*M*N)
!*******从打开的文件中读取数据************
       READ(5,rec=K) ((PS(K1,K2),K1=1,M),K2=1,N)
       READ(6,rec=K) (( H(K1,K2),K1=1,M),K2=1,N)

       read(PS(K1,K2),*)a1
       if(a1.lt.0.0)   PS(K1,K2)=0.0
     write(*,*)PS(K1,K2)
       WSH(K1,K2)=PS(K1,K2)*H(K1,K2)
       close(5)
       close(6)
end

lqouc

为什么要这么做呢?直接if对PS(K1,K2)判断啊,如果小于0就直接对其赋值为0就可以了.

又是那隻貓

谢谢1楼版主分享方法

兰北

lqouc 发表于 2013-7-20 14:45

                               
登录/注册后可看大图

为什么要这么做呢?直接if对PS(K1,K2)判断啊,如果小于0就直接对其赋值为0就可以了.

谢谢回复!刚开始我也这样写的，但是不起作用啊！因为PS(K1，K2)是二进制数据啊，二进制数据可以直接跟0.0比大小么？

lqouc

兰北 发表于 2013-7-20 17:07

                               
登录/注册后可看大图

谢谢回复!刚开始我也这样写的，但是不起作用啊！因为PS(K1，K2)是二进制数据啊，二进制数据可以直接跟0.0 ...

你试过以后不可以么?我以前都是这么做啊,FORTRAN处理的再分析数据不都是二进制么?你只要打开数据的时候声明UNFORMATTED就好了.
那如果像你说的不能比的话你都没法做运算,比如数据乘以2,你就要写成*10,这样显然是不对的吧.
不行你就先PRINGT一下PS(K1,K2)看看是怎么显示的.

兰北

lqouc 发表于 2013-7-21 09:00

                               
登录/注册后可看大图

你试过以后不可以么?我以前都是这么做啊,FORTRAN处理的再分析数据不都是二进制么?你只要打开数据的时候声 ...

也对哦，那我再试一下，谢谢斑竹！

hptru64

看起来太复杂了

兰北

lqouc 发表于 2013-7-21 09:00

                               
登录/注册后可看大图

你试过以后不可以么?我以前都是这么做啊,FORTRAN处理的再分析数据不都是二进制么?你只要打开数据的时候声 ...

已解决了，谢谢！

海豚∩_∩浩

求楼主贴程序学习