Frotran 删除字符串空格
Trim 删除字符串结尾处的空格
AdjustL 删除字符串开头的空格,左对齐
AdjustR 删除字符串结尾的空格,右对齐
write/read(UNIT=number,FMT=format,NML=namelist,REC=record,IOSTAT=stat,ERR=errlabel,END=endlabel)
(1)UNIT=number:指定read/write所使用的输入/出位置
(2)FMT=format,:指定输入输出格式
(3)NML=namelist：读写某个namelist
(4)REC=record：在直接读取文件中，设置索要读写的文件模块位置
(5)IOSTAT=stat:设置一个整数值给变量stat，用来说明文件的读写状态
stat>0,表示读取操作发生错误
stat=0,表示读取操作正常
stat<0,表示文件终了
(6)ERR=errlabel：在读写过程中发生错误时，转移到某个行代码来继续执行程序。
(7)END=endlabel：在读写到文件末尾时，转移到某个行代码来继续执行程序。
实例1
功能：挑出NLLJs的个例。
program main
implicit none
integer,parameter:: nline=100
integer i,hght(nline),relh(nline),drct(nline),sknt(nline)
integer max_lev,min_lev
real wind_max,wind_min
real pres(nline),temp(nline),dwpt(nline),mixr(nline),thta(nline),thte(nline),thtv(nline)
integer yr,imon,iday,ihr,days(12)
character(len=160) filename
character(len=4) cyr
character(len=2) cmon,cday,hours(2)
data hours/"00","12"/ !!!数组赋初值
data days/31,28,31,30,31,30,31,31,30,31,30,31/
open(13,file="lsd_select.log")
do yr = 2017, 2018
write(cyr,fmt='(i4.4)') yr
do imon = 1, 12
write(cmon,fmt='(i2)') imon
do iday = 1,days(imon)
write(cday,fmt='(i2)') iday
do ihr = 1,2
filename = "seasia_54511_sounding_"//cyr//"_"//trim(adjustl(cmon))//"_"//trim(adjustl(cday))//"_"//hours(ihr)//".txt"
!!!读取filename
open(12,file="/home/lsd/sound/"//filename)
do i=1,6
read(12,*,end=300) !!!!读写到文件末尾时，转移到300行代码继续执行
enddo
do i= 1,nline
read(12,200,err=100) pres(i),hght(i),temp(i),dwpt(i),relh(i),mixr(i),drct(i),sknt(i),thta(i),thte(i),thtv(i)
!write(*,200) pres(i),hght(i),temp(i),dwpt(i),relh(i),mixr(i),drct(i),sknt(i),thta(i),thte(i),thtv(i)
! error=100,读写过程中出现错误时，转移到100行代码继续执行.
enddo
200 FORMAT(1x,f6.1,1x,i6,2(1x,f6.1),1x,i6,1x,f6.1,2(1x,i6),3(1x,f6.1))
100 continue
close(12) !!!读取文件结束
!!! 判断3000m以下最大风速及最大风速的位置
wind_max=sknt(1)
max_lev=1
do i=1,nline
if(hght(i)>3000) exit
if(sknt(i)>wind_max) then
wind_max = sknt(i)
max_lev=i
endif
enddo
!!!在最大风速高度以上3000m以下找最小风速及高度
wind_min=wind_max
min_lev=max_lev
do i=max_lev,nline
if(hght(i)>3000) exit
if(sknt(i)<wind_min) then
wind_min = sknt(i)
min_lev=i
endif
enddo
print*,filename
print *,wind_max,hght(max_lev)
print *,wind_min,hght(min_lev)
print *,"----------------------------------"
if(wind_max * 0.514 >= 6 .and. wind_min <= wind_max*0.5) write(13,fmt='(A40,1x,2(f6.1,1x,i6))') filename,wind_max*0.514,hght(max_lev),wind_min*0.514,hght(min_lev)
300 continue
enddo
enddo
enddo
enddo
stop
end program main
实例2：
!数据：ASCII码
!目的：将单个站点的历史信息写入一个文件
program testio
implicit none
integer yr,imon,iday,ihr,iocode,num,ista,target,search_code
integer k,i,hr(4)
real days(12)
real lsd(25)
character(len=2) cmon,cdays,hours(4),chr
character(len=120) filename,dirname,ctarget
character(len=4) cyr
data days/31,28,31,30,31,30,31,31,30,31,30,31/
data hours/'02','08','14','20'/
data hr/2,8,14,20/
num = 20000
call getarg(1,ctarget)
read(ctarget,*) target
!target=40811
dirname = "."
open(13,file=trim(ctarget)//'.log') !!!trim将字符串末尾（右端）的空格删掉
open(22,file='./quece/'//trim(ctarget)//'_quece.log')
do yr = 2017, 2017
write(cyr,fmt='(i4.4)') yr
do imon = 1,12
write(cmon,fmt='(i2.2)') imon
do iday = 1,days(imon)
write (cdays,fmt='(i2.2)') iday
do ihr = 0,23
write (chr,fmt='(i2.2)') ihr
filename = cyr//cmon//cdays//chr//'00.000'
print*,filename
search_code=0
open(12,file=trim(dirname)//'/'//trim(filename), status='old',form='formatted',iostat=iocode)
if(iocode==0) then
do k=1,2
read(12,*,end=88) !如果文件为空，跳到88行语句执行
enddo
do k=1,num
read(12,*,end=88) ista, (lsd(i),i=1,25) !一共26列
if(ista==target) then
search_code=1
exit
endif
enddo
else
write(22,*) trim(filename)
endif
88 continue
if(search_code/=1) lsd=9999
write(13,fmt='(A,1x,25(f6.1,3x))') cyr//cmon//cdays//chr, (lsd(i),i=1,25)
close(12)
!call sleep( 1 )
enddo
enddo
enddo
enddo
end program testio
网友评论