gradsからデータを読み込み、pminを探すf90

gradsからデータを読み込み、pminを探すf90

だたし、gxoutを使ってgradsにおいてp1.grdを作ってから走らせること。
@yoshi/analysis/f_find_pmin.f90


program find_pmin
 implicit none
 integer, parameter :: xgrid=400, ygrid=300
 integer :: n, i, j, kx, ky
 real(kind=4), dimension(xgrid,ygrid) :: p_d
 real(kind=4) :: p
open (unit=17,file='p1.grd',form='unformatted',access='direct',recl=xgrid*ygrid*4)
read (unit=17,rec=1) p_d

open (unit=7, file='psea.txt', form='FORMATTED')
write(7,*) 'write mini psea xgrid ygrid lon lat'


p=2000
do i = 1, xgrid
   do j = 1, ygrid
      if (p_d(i,j) < p) then
         kx=i
         ky=j
        p=p_d(i,j)
         write(*,*) i, j, n
         write(7,*) n, 'min psea', p, 'kx', kx, 'ky', ky 
      endif
   n=n+1  
   enddo
enddo

write(7,*) n, ' out put   min psea', p, 'kx', kx, 'ky', ky 
write(7,*) 'yoshiyoshi'

close(7)
close(17)

end program  find_pmin