source: codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/config/ppsrc/dyn3d_common/inigrads.f @ 224

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


      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

!
! $Header$
!
      integer nfmx,imx,jmx,lmx,nvarmx
      parameter(nfmx=10,imx=200,jmx=150,lmx=200,nvarmx=1000)

      real xd(imx,nfmx),yd(jmx,nfmx),zd(lmx,nfmx),dtime(nfmx)

      integer imd(imx),jmd(jmx),lmd(lmx)
      integer iid(imx),jid(jmx)
      integer ifd(imx),jfd(jmx)
      integer unit(nfmx),irec(nfmx),itime(nfmx),nld(nvarmx,nfmx)

      integer nvar(nfmx),ivar(nfmx)
      logical firsttime(nfmx)

      character*10 var(nvarmx,nfmx),fichier(nfmx)
      character*40 title(nfmx),tvar(nvarmx,nfmx)

      common/gradsdef/xd,yd,zd,dtime,
     s   imd,jmd,lmd,iid,jid,ifd,jfd,
     s   unit,irec,nvar,ivar,itime,nld,firsttime,
     s   var,fichier,title,tvar

c     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'
     s   ,FORM='unformatted',
     s   ACCESS='direct'
     s  ,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