[RegCNET] reading lake output

Michael Notaro mnotaro at wisc.edu
Fri Jul 13 19:42:45 CEST 2007


I ran regcm3 for 2 days, Jan 1-2, 1996.  It outputted LAK.010100.
I am trying to write a fortran code to read the output.
It is included below.

What did I do wrong here?  It seems there is more data on each
line than I am reading, and some temperature values are huge.

Is it outputting an initial set of values and then again every hour?
It seems to me to have 49 time steps if I am looking at it right.

Thanks, Mike


       integer day, ilake(527), jlake(527), depth(527), skip, i
       real evap(527), hi(527), hice(527), hsnow(527), t(527,500)
       integer ilake1(527), jlake1(527), depth1(527)

       open(1,file="../fort.40",form="formatted",status="old")
       read(1,10) skip
       do i=1,527
         read(1,20) ilake(i), jlake(i), depth(i)
       end do
    10 format(I2)
    20 format(3I3)
       close(1)

       open(2,file="LAK.1996010100",form="unformatted",status="old")
       do inum=1,49
         do i=1,527
           read(2) skip, skip,ilake1(i), jlake1(i), depth1(i), evap(i),
      #      hi(i), hice(i), hsnow(i), (t(i,j),j=1,depth(i))
         end do
       end do

       close(2)

       end





More information about the RegCNET mailing list