如何将浮点型数组转化为字符串型数组
如题,需要将一个浮点型数组转化为字符串型数组知道可以用read或者write直接读,但是具体的格式是啥不太清楚
read(s(i),'(i)') a(i) 这个是将字符串型数组改为浮点型数组的语句,是不是只需要将'(i)'修改一下就好了
另外,还有一个问题,就是写了一段代码,将这个转化后的数组元素覆盖另外一个内容均为undefined的数组,但要判断转化完成后的字符型数组内的元素是否是74.40000 , 76.60000 , 70.00000 ,如果是的话就不予覆盖
下面是代码
do s = 1,14
do t = 1,18627
if(g(s,t) .eq. 74.40000 .or. g(s,t) .eq. 76.60000 .or. g(s,t) .eq. 70.00000) then
cycle
end if
e(s,t) = g(s,t)
end do
end do
注:e(s,t)即为另外一个均为undefined的字符串型数组,g(s,t)为由浮点型数组转化成为的字符型数组。
运行时出现了错误
e(s,t) = g(s,t)
1
Error: Unexpected STATEMENT FUNCTION statement at (1)
if(g(s,t) .eq. 74.40000 .or. g(s,t) .eq. 76.60000 .or. g(s,t) .eq. 70.00000) then
1
Error: Operands of comparison operator '.eq.' at (1) are CHARACTER(1)/REAL(4)
我把所有的代码放在这里,希望路过的大佬帮我看一下,得不到正确的结果
Program main
use timer_supports
use shuzu
implicit none
Integer,parameter :: length=1000
character(len=10)::g(18627,14)
integer i , f , idays , m , n , o , p , q , r , s , t
integer d
Character(LEN=length) A,B,Filename
Integer(kind=4)::Column, Row, K,L
real , allocatable :: c( : , : )
idays = 0
Column=0;Row=0;K=0;L=1;filename="ztzl-1.txt"
Open (55,file=Filename)
Read(55,'(a1000)') A
column=GetDataN(A)!//改成这个函数,具体实现在下面的代码
CLOSE(55)
Open (55,file=Filename)
Do
Read(55,*,iostat=k) B
If( k /= 0 ) Exit
Row = Row + 1
End Do
Write(*,"(A5,I6,5X,A7,I2)") "Row=",Row,"Column=",Column!读取行数
allocate(c(column-1,row))
rewind(55)
Do i = 1,Row
read(55,*) d,c(:,i)!//1改成 i
End Do
Close(55)!将数据转化为数组
do n = 7,10
do m = 1,18627
c(n,m) = c(n,m) * 0.1
end do
end do
p = 12
do o = 1,18627
if(c(p,o)>=32000) then
c(p,o) = (c(p,o)-32000) * 0.1
else if (c(p,o)>=31000) then
c(p,o) = (c(p,o)-31000) * 0.1
else if (c(p,o)>=30000) then
c(p,o) = (c(p,o)-30000) * 0.1
else
c(p,o) = c(p,o) * 0.1
end if
write(*,"(F10.5)") c(p,o)
end do
do q = 13 ,14
do r = 1 , 18627
c(q,r) = c(q,r) * 0.1
end do
end do
write(g,'(A)') c
do s = 1,14
do t = 1,18627
if(g(s,t) == '74.40000' .or. g(s,t) == '76.60000' .or. g(s,t) == '70.00000') then
cycle
end if
e(s,t) = g(s,t)
end do
end do
call daysfromdate(1957,1,1,2007,12,31,idays)
print *, idays
call indexit(1)
print *, e
contains
Integer Function GetDataN( cStr )
Character( Len = * ) , Intent( IN ) :: cStr
Integer :: i
Logical :: bIsSeparator , bIsQuote
GetDataN = 0
bIsSeparator = .TRUE.
bIsQuote = .FALSE.
Do i = 1 , Len_Trim( cStr )
Select Case( cStr(i:i) )
Case( '"' , "'" ) !// 如果遇到引号
If ( .Not.bIsQuote ) GetDataN = GetDataN + 1!//如果不在引号中,则增加一个数据
bIsQuote = .Not.bIsQuote !// 引号结束或开始
bIsSeparator = .FALSE.
Case( " " , "," , char(9) ) !// 如果遇到分隔符
If ( .Not.bIsQuote ) then!// 分隔符如果不在引号中
bIsSeparator = .TRUE.
End If
Case Default
If ( bIsSeparator ) then
GetDataN = GetDataN + 1
End If
bIsSeparator = .FALSE.
End Select
End Do
End Function GetDataN
End Program Main
数组e来自module shuzu,定义的是字符串型数组 为什么将实型数组转为字符串,目的? li913 发表于 2020-3-23 08:41
为什么将实型数组转为字符串,目的?
前辈,是这样的,我的老师让我先建一个内容均为undefined的字符串型数组,然后再让我将实型数据逐个填入将undefined替换掉,但是有一些数据是缺报,漏报的,在使用grads画图时时需要保留这些undefined的内容,所以我想到就是先将实型的数组转化为字符串型的,然后再实现数组内容的替换。 hpa 发表于 2020-3-23 09:26
前辈,是这样的,我的老师让我先建一个内容均为undefined的字符串型数组,然后再让我将实型数据逐个填入 ...
感觉你们这个思路很奇怪
正常的想法应该是先定义一个正常不可能取到的值作为初值
哪怕做个一个同样形状的mask矩阵也比你们这种操作合理…… 比如统计一年365天的降雨数据,现在只有300天的实测数据。1、开辟数组A(365),2、赋初值为一个不可能的降雨值 a=1.0e10,3、读取300个数据到A中相应位置;4、判断数据是否缺失(值很大),如果缺失,输出XXX, 否则输出XXX。 li913 发表于 2020-3-24 11:55
比如统计一年365天的降雨数据,现在只有300天的实测数据。1、开辟数组A(365),2、赋初值为一个不可能的降雨 ...
前辈,确实是这种思路,我刚开始进行数据的处理,不知道可以赋给他不可能的数值,麻烦您了 liudy02 发表于 2020-3-23 11:45
感觉你们这个思路很奇怪
正常的想法应该是先定义一个正常不可能取到的值作为初值
哪怕做个一个同样形状的 ...
前辈,确实是这种思路,我刚开始进行数据的处理,不知道可以赋给他不可能的数值,麻烦您了
页:
[1]