----赖格英-----

记忆不好了,记录工作中的点点滴滴....

  博客园  :: 首页  :: 新随笔  :: 联系 :: 订阅 订阅  :: 管理

一. Fortran 字符串与 C 字符串的区别
  Fortran的字符串处理能力其实很弱,关于字符串的语法还很落后。它与 C 字符串最大的区别就是:Fortran字符串是固定长度的,没有 \0 结束符。另外,Fortran 也不区分字符和字符串。即 'abc' 与 "abc" 是没有差别的。

二. Fortran 字符串的定义。
  Fortran字符串是固定长度的。因此,在声明时,就必须指定长度。(如不指定,大多数编译器认为是长度为 1 的字符串)
  声明时,可以按照这些格式进行:

1 integer , parameter :: L = 20
2 character :: s1*20 = "fcode.cn"
3 character*20 :: s2 = "fcode.cn"
4 character(20) :: s3 = "fcode.cn"
5 character( len = 20 ) :: s4 = "fcode.cn"
6 character( len = L ) :: s5 = "fcode.cn"

目前,笔者不建议使用第2行和第3行的方式。第4行的方式还算凑合。建议最好是使用第5行的方式,第6行演示了长度可以使用 parameter 常数,这样适合于需要大量定义相同长度的字符串。 

三. Fortran字符串的长度问题
    需要注意的是,字符串的长度一旦确定了,就无法改变。这一点笔者也感到语法应该改进,很不方便。

例如上例。character( Len = 20 ) :: s4 = "fcode.cn" 尽管后面的值 "fcode.cn" 只有 8 个字符,但是定义了变量是 20 长度,所以,实际 s4 的内容为 "fcode.cn            "(后面有12个空格)

如果我们要输出或使用字符串,往往需要去掉后面的空格,此时可使用 trim 函数。

write(*,*) trim(s4)

这样,就只会输出有内容的 8 个字符,后面的 12 个空格就不会输出。

特别需要注意的是,s4 = trim(s4) 这样的句子是没有任何意义的。因为虽然这个式子的右侧,trim(s4) 的结果是 8 个字符,但赋值给式子左侧的字符串,它依然是 20 长度。它等效于 s4 = "fcode.cn",所以 s4 的内容依然是 “fcode.cn            ”(读者可回去看上面的绿色文字)

这让我们很犯难。因此,trim 语句必须出现在每一次使用这个字符串的时候。

四. Fortran 字符串的连接(append)
  我们经常会挂靠字符串,比如 "fcode.cn" 在后面挂靠一个 "/bbs",很多初学者就会使用这样的代码:
  character( Len = 20 ) :: s4 = "fcode.cn"
  s4 = s4 // "/bbs"
  输出以后,发现 s4 的内容依然没有改变。并没有变成期望中的 "fcode.cn/bbs"。究其原因,其实就是第三个问题导致的。我们来观察上面的等式:s4 = s4 // "/bbs",等号左边 s4 有 20 长度,等号右边 s4 的 20 长度加上 "/bbs" 的 4 长度,一共 24 个长度。左边20个长度根本容纳不下右边的24个长度,于是,s4 依然是 s4,后面挂靠的 "/bbs" 由于存储不下而被丢弃。

  正确的挂靠方式是: s4 = trim(s4) // "/bbs"

  而反之,s4 = "http://" // s4 就不需要写 trim,请读者朋友自己思考为什么?

五. 字符串的左右对齐(AdjustL,AdjustR)
  我们经常会遇到这样的问题,除了后面的空格外,有时候字符串前面也有空格。例如 s4 = "    FortranCoder  ",单一的 trim 结果会是 "     FortranCoder"。
  此时我们就需要 AdjustL 函数,我们可以这样写: trim(AdjustL(s4)) , 或者 AdjustL(trim(s4)) 这两者的顺序无所谓,结果都是一样的。
  AdjustR 函数使用的情况非常少,通常不用。它可以把 "  abc  " 变为 "    abc"

六. 字符串与整型,实型的相互转换
  Fortran 的字符串与整型实现转换,是个很有意思的事情。它不像其他语言那样,提供一个函数来进行。而使用 read 和 write 读写来实现。
  看下面的例子:

 

 1 Program www_fcode_cn
 2   Implicit None
 3   Real :: r
 4   Integer :: i
 5   Character( Len = 20 ) :: c
 6   c = "3.1415926"
 7   read( c , * ) r !// 将 字符串c 转换为实数 r
 8   write( * , * ) r + 1.0
 9   i = int( r - 2.0 ) !// 给 i 一个值
10   write( c , '(i0)' ) i !// 把整型i 转换为字符串 c
11   c = AdjustL(trim(c)) // "st"
12   write( * , * ) Trim( c )
13 End Program www_fcode_cn

这里的第 7 行,看起来是一个 read 语句,它其实是转换。意思是:从字符串 c 中读取 r 的值。读者可以把此时的字符串 c 想成是一个虚拟的文件。
这里的第 10 行,看起来是一个 write 语句,其实它还是转换。意思是:把 i 的值,写入字符串 c 中,读者依然可以想象成 c 是一个虚拟的文件。
简单的说,就是 read 和 write 语句,可以直接对字符串进行操作。字符串可以被认为是虚拟的文件,从字符串中获得数值,就是read;把值写入到字符串中,就是write。

七. 一个实例
  这里,我们列举一个实际字符串使用的例子。
  假设现在有一个文件,其内容如下:
// This is a sample file
// There may be some comment text
Time:13:23:42 RecordID:18 Weather:Sunny
Fcode.cn Date:2014-02-38

前两行是注释,而且有可能不是两行,不确定有多少行。
然后是一些数据组合在一起。我们需要找到 RecordID: 后面的数字,即,18。我们事先不能确定 RecordID 在第几行的什么位置里。
找到之后,我们需要给它加上 2014,即,2014+18=2032。
最后我们打开 File2031.dat 文件,开始读取其中的数据。 

上面的一系列操作,如果不使用字符串,恐怕会很难。下面我们来看看如何书写代码来完成。利用本文所提到的各种内容。

代码中有一部分注释,相信读者朋友能够理解

 1 Program www_fcode_cn
 2   Implicit None
 3   Integer , parameter :: MAX_PATH = 512
 4   Character( Len = MAX_PATH ) :: c
 5   Character( Len = * ) , parameter :: STR_FIND = "RecordID:"
 6   integer :: ID , i , iErr
 7   Open( 12 , File = "fcode.txt" )
 8   Do
 9     Read( 12 , '(a512)' , ioStat = iErr ) c !必须a512格式,否则遇到空格会终止
10     if ( iErr /= 0 ) Exit !// 如果读错了,比如遇到文件结束,则退出循环
11     if ( c(1:2) == "//" ) cycle !// 如果读取的是注释行,则直接读下一行
12     i = index( c , STR_FIND ) !// 在 c 中搜索 RecordID: 的位置
13     if ( i > 0 ) then !// 如果找到
14       i = i + Len(STR_FIND) !// 将位置移到RecordID:后面
15       c = c(i:) !// 使得 c 的内容为 i 位置之后的内容,即放弃i位置之前的内容
16       read( c , * ) ID !// 读取 ID 编号
17       write( c , * ) ID + 2014 !// ID 编号加上2014,再写回c中
18       c = "File" // Trim(AdjustL(c)) // ".dat" !组合成File2032.dat,注意trim
19       write( * , * ) Trim( c ) !// 输出看看对不对
20       !// Open( 13 , File = Trim( c ) )
21       Exit
22     end if
23   End Do
24   Close( 12 )
25 End Program www_fcode_cn

转自:http://fcode.cn/guide-48-1.html

 

posted on 2015-04-06 11:05  向北方  阅读(3508)  评论(0编辑  收藏  举报