c_ --------------------------------------------------------------------- c_ RCS lines preceded by "c_ " c_ --------------------------------------------------------------------- c_ c_ $Source: /home/orr/WWW/Biotic/boundcond/RCS/readmap.f,v $ c_ $Revision: 1.3 $ c_ $Date: 1999/04/04 02:32:29 $ ; $State: Exp $ c_ $Author: orr $ ; $Locker: $ c_ c_ --------------------------------------------------------------------- c_ $Log: readmap.f,v $ c_ Revision 1.3 1999/04/04 02:32:29 orr c_ Cleaned up and reformatted for clarity c_ c_ Revision 1.2 1999/03/22 14:04:37 orr c_ changed name of input file (po4map.dat -> po4mapnew.dat) c_ c_ Revision 1.1 1999/03/22 14:03:47 orr c_ Initial revision c_ c_ --------------------------------------------------------------------- c_ program check dimension a(180,90) dimension alat(90) dimension alon(180) open(21,file ='po4mapnew.dat',form='unformatted') c loop through the months do i = 1,12 c loop through the first 6 level of Levitus 0 to 75 m do j = 1,6 write(*,*) 'read month',i,' level',j read(21) a do ii = 1,90 c First latitude 90 N alat(ii) = 91. - 2.*float(ii) do jj = 1,180 c First longitude Greenwich alon(jj) = -1. + 2.*float(jj) c Example c write longitude / latitude / PO4 for first month (janu) and c first level 0 m if(i .eq. 1 .and. j .eq. 1) then write(2,*) alon(jj),alat(ii),a(jj,ii) endif end do end do end do end do 3 format(2(f7.1,1x),f8.3) stop END