Results 1 to 2 of 2
relevant part of the code:
do i=1,nsites
do ico=1,nc
if(ivarp(ico).ne.0)then
iprd=ivarp(ico)
apr4sens(1,i,iprd)=apr1dsens(i,iprd)
apr4sens(2,i,iprd)=apr2dsens(i,iprd)
apr4sens(3,i,iprd)=apr3dsens(i,iprd)
apr4sens(4,i,iprd)=apr4dsens(i,iprd)
endif
enddo
enddo
do i=1,nsites
do ico=1,nc
if(ivarp(ico).ne.0)then
iprd=ivarp(ico)
do k=1,4
write(128,*)apr4sens(k,i,iprd)
enddo
endif
enddo
...
- 06-03-2011 #1Just Joined!
- Join Date
- May 2010
- Posts
- 41
fortran 90 rank 3array
relevant part of the code:
do i=1,nsites
do ico=1,nc
if(ivarp(ico).ne.0)then
iprd=ivarp(ico)
apr4sens(1,i,iprd)=apr1dsens(i,iprd)
apr4sens(2,i,iprd)=apr2dsens(i,iprd)
apr4sens(3,i,iprd)=apr3dsens(i,iprd)
apr4sens(4,i,iprd)=apr4dsens(i,iprd)
endif
enddo
enddo
do i=1,nsites
do ico=1,nc
if(ivarp(ico).ne.0)then
iprd=ivarp(ico)
do k=1,4
write(128,*)apr4sens(k,i,iprd)
enddo
endif
enddo
enddo
I manage to write apr4sens into file,works fine.But how to write it into array,which I could later reshape and so on..
I tryied like this
group=apr4sens
It doesn't work.
- 06-04-2011 #2Just Joined!
- Join Date
- Dec 2008
- Location
- Montreal
- Posts
- 5
write command in f90
Is the Output/result was something like
that...?
If i=1...prd=1
value of apr4sens(1,1,1)
value of apr4sens(2,1,1)
value of apr4sens(3,1,1)
value of apr4sens(4,1,1)
...
If i=1...prd=2
value of apr4sens(1,1,2)
value of apr4sens(2,1,2)
value of apr4sens(3,1,2)
value of apr4sens(4,1,2)
If you want this result,
values of
If i=i and prd=1
apr4sens(1,i,1) apr4sens(2,i,1) apr4sens(3,i,1) apr4sens(4,i,1)
...following by prd=2
apr4sens(1,i,2) apr4sens(2,i,2) apr4sens(3,i,2) apr4sens(4,i,2)
you shall write this line
write(128,*)apr4sens(1,i,prd),apr4sens(2,i,prd),ap r4sens(3,i,prd),apr4sens(4,i,prd)


Reply With Quote