Hello,
I am writing a program that loops through a file with many instances of data taken at different times. I use a derived data type to create the linked list. See the code below.
module types
type
:: wavedat
real, dimension(:, :) :: espt
real :: winddir
real :: windspeed
real :: time
type (wavedat), pointer :: p
end type wavedat
end module
types
So I then have a module that reads in a data file and gets values for winddir, windspeed, and time. The value for wavedat%espt is given inReadOWI by the transpose of spec2d. Spec2d is a rank-2 array of size (nfreq,nang).How can I initialize wavedat%espt to dimension(nfreq,nang)?
module
ReadOWI
contains
subroutine
OWI(freqz,owidat)
use
types
integer :: OPEN_STATUS,yr,mo,da,hr,mm,nfreq,nang,count
real
:: TWD,U10,stime,dep
real
, allocatable, dimension(:,:) :: dummy,spec2d
real
, allocatable, dimension(:), intent(out) :: freqz
real
, allocatable, dimension(:) :: thetad
type
(wavedat), pointer :: tail
type
(wavedat), pointer :: head
type
(wave
dat), pointer :: ptr
type
(wavedat), allocatable, dimension(:), intent(out) :: owidat
!Initialize loop count
count = 0
!Open file for read only
open(UNIT=2,FILE=sample.asc,STATUS='OLD',&
ACTION='READ',IOSTAT=OPEN_STATUS)
if(OPEN_STATUS>0)STOP "----Error, File not opened properly-----"
! This section reads in the number of frequencies and directions (angles)
read(2,'(11X,I,12X,I)'), nfreq,nang !These values define size of owidat%espt
read(2,'(A)') !Blank Line
!Now allocate dimensions to arrays
allocate (freqz(nfreq))
allocate (thetad(nang))
allocate (spec2d(nang,nfreq))
allocate (dummy(nang,nfreq+1))
! Main loop
do
count = count+1
if (count == 1) then
read(2,'(64X,F9.4)'),dep !Depth
else
read(2,'(A)') !First Header Line
end if
read(2,'(A)') !Second Header Line
!Now read year, month, day, hour, minute, wind direction, and wind speed
read(2,'(1X,I4,I2,1X,3I2,7X,F6.1,F7.3)'),yr,mo,da,hr,mm,TWD,U10
read(2,'(A)') !Fourth Header Line
read(2,'(7X)',ADVANCE='NO') !Unwanted Characters
stime = datenum(yr,mo,da,hr,mm)
!Serial date number
!Read in the frequencies to vector freq
if (count==1) then
do p=1,nfreq
read(2,'(F7.4)',ADVANCE='NO') , freqz(p)
end do
end if
!Read directions and spec2d data to 'dummy' variable
do k=1,nang
read(2,'(A)')
do j=1,nfreq+1
if (j==1) then
read(2,'(F7.2)',ADVANCE='NO'),dummy(k,j)
else
read(2,'(F7.4)',ADVANCE='NO'),dummy(k,j)
end if
end do
end do
!Sort Directions and spectrum
dummy = Sort(dummy,nang,nfreq+1)
!Sorted Directions vector
thetad = dummy(:,1)
!Rotate spectrum angles to direction 'From'
spec2d(1:(nang/2),:) = dummy((nang/2)+1:nang,2:nfreq+1)
spec2d((nang/2)+1:nang,:) = dummy(1:nang/2,2:nfreq+1)
read(2,'(A)',IOSTAT=OPEN_STATUS) !Blank Line
!Update owidat structure with this observation
if (.NOT. ASSOCIATED(head)) then
allocate (head,STAT=istat)
! allocate (head%espt(nfreq,nang))
tail => head
nullify (tail%p)
tail%espt =
transpose(spec2d)
tail%winddir = TWD
tail%windspeed = U10
tail%time = stime
else
allocate(tail%p,STAT=istat)
tail => tail%p
nullify (tail%p)
tail%espt =
transpose(spec2d)
tail%winddir = TWD
tail%windspeed = U10
tail%time = stime
end if
!End loop at the end of file
end do
!Determine if loop terminated because of a READ error or the end of the input file
if (OPEN_STATUS>0) then
print *,'An Error Occurred While Reading File'
else
print *,'End of file Reached'
end if
allocate
(owidat(count))
ptr => head
do
i=1,count
owidat(i)%espt = ptr%espt
owidat(i)%winddir = ptr%winddir
owidat(i)%windspeed = ptr%windspeed
owidat(i)%time = ptr%time
ptr => ptr%p
end do
close
(2)
end subroutine
OWI
end module ReadOWI
bobfrat@gmail.com:
