主题:缺测值的处理
fn后面到的是风速数据,有31*8行,前面31*4是2分钟数据,等号后面31*4行是10分钟数据,我需要读取10分钟的数据
fn
....................................
166033 146034 155031 134031 147043 132039.
127033 118033 110035 117034 107038 126025
125027 076014 125026 137028 148025 //////
////// ////// PPC000 ////// ////// //////
////// PPC012 ////// ////// ////// //////.
...............................
我想把数据中的//////和PPC???替换成999999,用下面的语句
do j=1,n
if(b(j)=='//////'.or.b(j)==‘PPC???’) b(j)='999999'
read(b(j),"(i6)")a(i,j)
enddo
为什么结果是这样啊
166033 146034 155031 134031 147043 132039
127033 118033 110035 117034 107038 126025
125027 076014 125026 137028 148025 126025
125027 076014 125026 137028 148025 126025
125027 076014 125026 137028 148025 126025
....................................
fn
....................................
166033 146034 155031 134031 147043 132039.
127033 118033 110035 117034 107038 126025
125027 076014 125026 137028 148025 //////
////// ////// PPC000 ////// ////// //////
////// PPC012 ////// ////// ////// //////.
...............................
我想把数据中的//////和PPC???替换成999999,用下面的语句
do j=1,n
if(b(j)=='//////'.or.b(j)==‘PPC???’) b(j)='999999'
read(b(j),"(i6)")a(i,j)
enddo
为什么结果是这样啊
166033 146034 155031 134031 147043 132039
127033 118033 110035 117034 107038 126025
125027 076014 125026 137028 148025 126025
125027 076014 125026 137028 148025 126025
125027 076014 125026 137028 148025 126025
....................................