fortran
1 program readflux4 2 3 character (len=150) :: outfile, fname 4 real,allocatable, dimension (:,:,:,:,:,:,:) :: flux 5 real,allocatable, dimension (:,:,:) :: areanorth 6 real,allocatable, dimension (:,:,:) :: areaeast 7 8 open (9, file='fluxlist')910 read (9,'(a)',end=90) fname 10 11 open (10,file=trim(fname),form='unformatted',ACCESS='stream')12 print*, 'reading ',trim(fname)13 write (outfile,900) trim(fname)14 900 format (a,'.dat')15 16 open (20, file=trim(outfile) )17 read (10) 1.e12*flux(:,:,:,:,:,:,:)/areanorth(:,:,:)/outstep 18 write (20,901) 1.e12*flux(:,:,:,:,:,:,:)/areanorth(:,:,:)/outstep 19 901 format (f12.5)20 21go to 102290 close (9)23 24 end25
0 コメント