source: LMDZ6/branches/contrails/libf/dyn3d_common/inigrads.f90 @ 5451

!
! $Header$
!
subroutine inigrads(if,im &
        ,x,fx,xmin,xmax,jm,y,ymin,ymax,fy,lm,z,fz &
        ,dt,file,titlel)
  USE gradsdef_mod_h
  implicit none

  integer :: if,im,jm,lm,i,j,l
  real :: x(im),y(jm),z(lm),fx,fy,fz,dt
  real :: xmin,xmax,ymin,ymax

  character(len=*),intent(in) :: file
  character(len=*),intent(in) :: titlel

  ! data unit/66,32,34,36,38,40,42,44,46,48/
  integer :: nf
  save nf
  data nf/0/

  unit(1)=66
  unit(2)=32
  unit(3)=34
  unit(4)=36
  unit(5)=38
  unit(6)=40
  unit(7)=42
  unit(8)=44
  unit(9)=46

  if (if.le.nf) stop 'verifier les appels a inigrads'

  print*,'Entree dans inigrads'

  nf=if
  title(if)=titlel
  ivar(if)=0

  fichier(if)=trim(file)

  firsttime(if)=.true.
  dtime(if)=dt

  iid(if)=1
  ifd(if)=im
  imd(if)=im
  do i=1,im
     xd(i,if)=x(i)*fx
     if(xd(i,if).lt.xmin) iid(if)=i+1
     if(xd(i,if).le.xmax) ifd(if)=i
  enddo
  print*,'On stoke du point ',iid(if),'  a ',ifd(if),' en x'

  jid(if)=1
  jfd(if)=jm
  jmd(if)=jm
  do j=1,jm
     yd(j,if)=y(j)*fy
     if(yd(j,if).gt.ymax) jid(if)=j+1
     if(yd(j,if).ge.ymin) jfd(if)=j
  enddo
  print*,'On stoke du point ',jid(if),'  a ',jfd(if),' en y'

  print*,'Open de dat'
  print*,'file=',file
  print*,'fichier(if)=',fichier(if)

  print*,4*(ifd(if)-iid(if))*(jfd(if)-jid(if))
  print*,trim(file)//'.dat'

  OPEN (unit(if)+1,FILE=trim(file)//'.dat' &
        ,FORM='unformatted', &
        ACCESS='direct' &
        ,RECL=4*(ifd(if)-iid(if)+1)*(jfd(if)-jid(if)+1))

  print*,'Open de dat ok'

  lmd(if)=lm
  do l=1,lm
     zd(l,if)=z(l)*fz
  enddo

  irec(if)=0

  print*,if,imd(if),jmd(if),lmd(if)
  print*,'if,imd(if),jmd(if),lmd(if)'

  return
end subroutine inigrads