Context Navigation

lect-eis.F @ 23

Last change on this file since 23 was 4, checked in by dumas, 10 years ago
initial import GRISLI trunk
File size: 1.5 KB

Line
1	!> \file lect-eis.F
2	!! Module pour la lecture des donnees
3	!<
4
5	!> SUBROUTINE: lect_eis
6	!! Lecture des coupes
7	!! \author ...
8	!! \date ...
9	!! @param nx dimension along x
10	!! @param ny dimension along y
11	!! @param data array for reading
12	!! @param filin
13	!! @param dirnameinp
14	!! @return data the value of Tab(nxx,nyy): Warning it is a real
15
16
17	subroutine lect_eis(nx,ny,data,filin,DIRNAMEINP)
18
19	implicit none
20	character(len=17) :: DIRNAMEINP
21	integer nx,ny,i,j,kdum
22	real data(nx,ny)
23	character(len=80) :: title,fmt
24	character(len=80) :: filin
25
26	integer :: num_coupe = 11
27
28	filin = TRIM(DIRNAMEINP)//TRIM(filin)
29
30
31	open(num_coupe,file=filin)
32
33	read(num_coupe,'(A80)') title
34	read(num_coupe,'(A80)') fmt
35
36	do J=1,NY
37	read(num_coupe,fmt) kdum,(data(i,j),i=1,NX)
38	end do
39
40	close(num_coupe)
41	return
42	end
43
44	subroutine lect_ieis(nx,ny,data,filin,DIRNAMEINP)
45
46	implicit none
47	character(len=17) :: DIRNAMEINP
48	integer nx,ny,i,j,kdum
49	integer data(nx,ny)
50	character(len=80) :: title,fmt
51	character(len=80) :: filin
52
53	integer :: num_coupe = 11
54	filin = TRIM(DIRNAMEINP)//TRIM(filin)
55
56	open(num_coupe,file=filin)
57
58	read(num_coupe,'(A80)') title
59	read(num_coupe,'(A80)') fmt
60
61	do J=1,NY
62	read(num_coupe,fmt) kdum,(data(i,j),i=1,NX)
63	end do
64
65	close(num_coupe)
66	return
67	end

Note: See TracBrowser for help on using the repository browser.