source: LMDZ.3.3/trunk/libf/dyn3d/inigrads.F @ 338

Last change on this file since 338 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.6 KB
RevLine 
[2]1      subroutine inigrads(if,im
2     s  ,x,fx,xmin,xmax,jm,y,ymin,ymax,fy,lm,z,fz
3     s  ,dt,file,titlel)
4
5
6      implicit none
7
8      integer if,im,jm,lm,i,j,l,lnblnk
9      real x(im),y(jm),z(lm),fx,fy,fz,dt
10      real xmin,xmax,ymin,ymax
11
12      character file*10,titlel*40
13
14#include "gradsdef.h"
15
16      data unit/66,32,34,36,38,40,42,44,46,48/
17      data nf/0/
18
19      if (if.le.nf) stop'verifier les appels a inigrads'
20
21      print*,'Entree dans inigrads'
22
23      nf=if
24      title(if)=titlel
25      ivar(if)=0
26
27      fichier(if)=file(1:lnblnk(file))
28
29      firsttime(if)=.true.
30      dtime(if)=dt
31
32      iid(if)=1
33      ifd(if)=im
34      imd(if)=im
35      do i=1,im
36         xd(i,if)=x(i)*fx
37         if(xd(i,if).lt.xmin) iid(if)=i+1
38         if(xd(i,if).le.xmax) ifd(if)=i
39      enddo
40      print*,'On stoke du point ',iid(if),'  a ',ifd(if),' en x'
41
42      jid(if)=1
43      jfd(if)=jm
44      jmd(if)=jm
45      do j=1,jm
46         yd(j,if)=y(j)*fy
47         if(yd(j,if).gt.ymax) jid(if)=j+1
48         if(yd(j,if).ge.ymin) jfd(if)=j
49      enddo
50      print*,'On stoke du point ',jid(if),'  a ',jfd(if),' en y'
51
52      print*,'Open de dat'
53      print*,'file=',file
54      print*,'fichier(if)=',fichier(if)
55
56      print*,4*(ifd(if)-iid(if))*(jfd(if)-jid(if))
57      print*,file(1:lnblnk(file))//'.dat'
58
59      OPEN (unit(if)+1,FILE=file(1:lnblnk(file))//'.dat',
60     s   FORM='UNFORMATTED',
61     s   ACCESS='DIRECT'
62     s  ,RECL=4*(ifd(if)-iid(if)+1)*(jfd(if)-jid(if)+1))
63
64      print*,'Open de dat ok'
65
66      lmd(if)=lm
67      do l=1,lm
68         zd(l,if)=z(l)*fz
69      enddo
70
71      irec(if)=0
72
73      print*,if,imd(if),jmd(if),lmd(if)
74      print*,'if,imd(if),jmd(if),lmd(if)'
75
76      return
77      end
Note: See TracBrowser for help on using the repository browser.