Fortran Coder

标题: 如何将浮点型数组转化为字符串型数组 [打印本页]

作者: hpa    时间: 2020-3-22 09:34
标题: 如何将浮点型数组转化为字符串型数组
如题,需要将一个浮点型数组转化为字符串型数组
知道可以用read或者write直接读,但是具体的格式是啥不太清楚
  read(s(i),'(i)') a(i) 这个是将字符串型数组改为浮点型数组的语句,是不是只需要将'(i)'修改一下就好了

另外,还有一个问题,就是写了一段代码,将这个转化后的数组元素覆盖另外一个内容均为undefined的数组,但要判断转化完成后的字符型数组内的元素是否是74.40000 , 76.60000 , 70.00000 ,如果是的话就不予覆盖
下面是代码
[Fortran] 纯文本查看 复制代码
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)


作者: hpa    时间: 2020-3-22 14:23
我把所有的代码放在这里,希望路过的大佬帮我看一下,得不到正确的结果
[Fortran] 纯文本查看 复制代码
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
为什么将实型数组转为字符串,目的?
作者: hpa    时间: 2020-3-23 09:26
li913 发表于 2020-3-23 08:41
为什么将实型数组转为字符串,目的?

前辈,是这样的,我的老师让我先建一个内容均为undefined的字符串型数组,然后再让我将实型数据逐个填入将undefined替换掉,但是有一些数据是缺报,漏报的,在使用grads画图时时需要保留这些undefined的内容,所以我想到就是先将实型的数组转化为字符串型的,然后再实现数组内容的替换。
作者: liudy02    时间: 2020-3-23 11:45
hpa 发表于 2020-3-23 09:26
前辈,是这样的,我的老师让我先建一个内容均为undefined的字符串型数组,然后再让我将实型数据逐个填入 ...

感觉你们这个思路很奇怪
正常的想法应该是先定义一个正常不可能取到的值作为初值
哪怕做个一个同样形状的mask矩阵也比你们这种操作合理……
作者: li913    时间: 2020-3-24 11:55
比如统计一年365天的降雨数据,现在只有300天的实测数据。1、开辟数组A(365),2、赋初值为一个不可能的降雨值 a=1.0e10,3、读取300个数据到A中相应位置;4、判断数据是否缺失(值很大),如果缺失,输出XXX, 否则输出XXX。
作者: hpa    时间: 2020-3-25 09:51
li913 发表于 2020-3-24 11:55
比如统计一年365天的降雨数据,现在只有300天的实测数据。1、开辟数组A(365),2、赋初值为一个不可能的降雨 ...

前辈,确实是这种思路,我刚开始进行数据的处理,不知道可以赋给他不可能的数值,麻烦您了
作者: hpa    时间: 2020-3-25 09:53
liudy02 发表于 2020-3-23 11:45
感觉你们这个思路很奇怪
正常的想法应该是先定义一个正常不可能取到的值作为初值
哪怕做个一个同样形状的 ...

前辈,确实是这种思路,我刚开始进行数据的处理,不知道可以赋给他不可能的数值,麻烦您了




欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2