文件操作

深度狙杀

fortran的文件操作里open（1，file=...）file后面可以是字符数组里的一个数组元素吗
试了一下貌似不可以
求指导呀

mofangbao

你的代码怎么写的？贴出来看看

深度狙杀

mofangbao 发表于 2013-3-16 11:29

                               
登录/注册后可看大图

你的代码怎么写的？贴出来看看

00041300.00U   00041300.00V   00041300.10H   00041300.10Q  文件名是这样的，有很多，220个
open(1,file='D:\pra\pro\name.txt')
     do i=1,nn
            read(1,*) name(i)
                write(*,*) name(i)
      end do
open(2,file='D:\pra\pro\name(1)')
      read(2,*)     
      do j=1,y
               read(2,100)(h(i,j),i=1,x)  
          end do
我把那些名字写到name.txt里。然后读到数组里 read(1,*) name(i)
然后我open  name（1） （就是00041300.00U ），第一个元素表示的数组名里的内容，但打不开

PS：character*12,dimension(220)::name
   

mofangbao

看这个帖子吧 你写的有问题 name(i)直接放到引号里面就是一个字符串而已，不会被解析为变量
http://bbs.06climate.com/forum.php?mod=viewthread&tid=3079

godenflame135

深度狙杀 发表于 2013-3-16 14:44

                               
登录/注册后可看大图

00041300.00U   00041300.00V   00041300.10H   00041300.10Q  文件名是这样的，有很多，220个
open(1,f ...

主要思想如下，其余请自我发挥。

!------------------------------------------------------------
!00041300.00U   00041300.00V   00041300.10H   00041300.10Q  文件名是这样的，有很多，220个
...
character(len=256)  ::   c_temp1     =  "  "
character(len=100)  ::   c_big_path  =  "  "
...
integer(kind=4)       ::   nn              = 0
integer(kind=4)       ::   n_files         = 0

nn      = ?
n_files = ?

open(1,file='D:\pra\pro\name.txt')
     do i=1,nn
            read(1,*) name(i)
                write(*,*) name(i)
     end do
close(1)

c_big_path="D:\pra\pro\"


do ii=1,n_files  

c_temp1=adjustl(trim(c_big_path))//adjustl(trim(name(i)))

open(2,file=adjustl(trim(c_temp1)))
      read(2,*)     
      do j=1,y
               read(2,100)(h(i,j),i=1,x)  
      end do !y
close(2)
end do !ii


...

!--------到此结束-------------------------

深度狙杀

godenflame135 发表于 2013-3-16 15:50

                               
登录/注册后可看大图

主要思想如下，其余请自我发挥。

!------------------------------------------------------------

谢谢了，编出来了
subroutine readdata
use  param
z=1
tt=1
a=1
b=1
open(1,file='D:\pra\pro\name.txt')
     do i=1,nn
            read(1,*) name(i)
        !        write(*,*) name(i)
      end do
          c_big_path="D:\pra\pro\"




    do n=1,nn
          
           c=mod(n,5)
          
          if(a==0) then
                     z=z+1
                         !print*,z
                    end if



                        if(b==0) then
                     tt=tt+1
                          z=1
                         !print*,tt
                    end if


             a=mod(n,5)
                 b=mod(n,55)
                 
                  




          c_temp1=adjustl(trim(c_big_path))//adjustl(trim(name(n)))
          select case(c)
            case(1)
                open(2,file=adjustl(trim(c_temp1)))
                          read(2,*)     
                         do j=1,y
                             read(2,100) (h(i,j,z,tt),i=1,x)
                                                   end do
                                    
                close(2)
   
             
                case(2)
                          open(2,file=adjustl(trim(c_temp1)))
                       read(2,*)     
                         do j=1,y
                           read(2,100) (q(i,j,z,tt),i=1,x)
                                                   end do
                                  
                close(2)
   
                case(3)

                 open(2,file=adjustl(trim(c_temp1)))
                       read(2,*)     
                         do j=1,y
                           read(2,100) (t(i,j,z,tt),i=1,x)
                                                   end do
                                    
                close(2)


                case(4)

                         open(2,file=adjustl(trim(c_temp1)))
                       read(2,*)     
                         do j=1,y
                           read(2,100) (u(i,j,z,tt),i=1,x)
                                                   end do
                                    
                close(2)
                case(0)


                         open(2,file=adjustl(trim(c_temp1)))
                       read(2,*)     
                         do j=1,y
                           read(2,100) (v(i,j,z,tt),i=1,x)
                                                   end do
                                  
                close(2)


               
          end select
                     100        format(10f8.2)               
                        end do



         
      do j=1,y
         lat(j)=10.0+float(j-1)*1.0   
      enddo
          
          
end subroutine readdata
  

深度狙杀

mofangbao 发表于 2013-3-16 15:04

                               
登录/注册后可看大图

看这个帖子吧 你写的有问题 name(i)直接放到引号里面就是一个字符串而已，不会被解析为变量
http://bbs.06 ...

恩恩，下载了，但是看不懂

ybc8012y

好文章，谢谢分享！