! ! $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