implicit none              

integer::i,n            

parameter(n=86400)            

external hardspike        

type tl              

   character(21) date          

   integer rn          

   real(8) ux

   real(8)uy

   real(8)uz

   real(8)co2

   real(8)h2o

   real(8)ts

   real(8)p    

   integer dia            

end type tl              

type(tl),dimension(n)::tuanliu        

open(10,file=\'c:\\湍流\\20-new.dat\',status=\'old\')        

do i=1,n

read(10,*)tuanliu(i).date,tuanliu(i).rn,tuanliu(i).ux,tuanliu(i).uy,tuanliu(i).uz,tuanliu(i).co2,tuanliu(i).h2o,tuanliu(i).ts,tuanliu(i).p,tuanliu(i).dia

enddo

call hardspike(tuanliu)            

close(10)              

end program tuanliu_ex            

subroutine hardspike(tuanliu)      

implicit none

integer i,n            

parameter(n=864000)

type tl              

   character(21) date          

   integer  rn          

   real(8)ux

   real(8)uy

   real(8)uz

   real(8)co2

   real(8)h2o

   real(8)ts

   real(8)p    

   integer dia            

end type tl              

type(tl),dimension(n)::tuanliu                    

do i=2,n-1              

   if(tuanliu(i).dia==1)then

   tuanliu(i).dia=0

   tuanliu(i).date=tuanliu(i).date

   tuanliu(i).rn=tuanliu(i).rn            

   tuanliu(i).ux=(tuanliu(i-1).ux+tuanliu(i+1).ux)*1.0/2.0        

   tuanliu(i).uy=(tuanliu(i-1).uy+tuanliu(i+1).uy)*1.0/2.0        

   tuanliu(i).uz=(tuanliu(i-1).uz+tuanliu(i+1).uz)*1.0/2.0        

   tuanliu(i).co2=(tuanliu(i-1).co2+tuanliu(i+1).co2)*1.0/2.0        

   tuanliu(i).h2o=(tuanliu(i-1).h2o+tuanliu(i+1).h2o)*1.0/2.0        

   tuanliu(i).ts=(tuanliu(i-1).ts+tuanliu(i+1).ts)*1.0/2.0        

   tuanliu(i).p=(tuanliu(i-1).p+tuanliu(i+1).p)*1.0/2.0        

   endif              

end do              

do i=1,n

print*,tuanliu(i).date,tuanliu(i).rn,tuanliu(i).ux,tuanliu(i).uy,tuanliu(i).uz,tuanliu(i).co2,tuanliu(i).h2o,tuanliu(i).ts,tuanliu(i).p,tuanliu(i).dia

enddo      

end subroutine hardspike

--------------------Configuration: tl - Win32 Debug--------------------

Compiling Fortran...

C:\\MSDEV\\Projects\\tl.f90

C:\\MSDEV\\Projects\\tl.f90 : warning FOR4227: argument TUANLIU (number 1) in reference to procedure HARDSPIKE from procedure TUANLIU_EX incorrect: has the wrong data type

tl.obj - 0 error(s), 1 warning(s)

到底啥问题哈