aboutsummaryrefslogtreecommitdiff
path: root/orrs/src/mconv.f
blob: 6e05c08fef63b2d500e3e91477dfe87154adf5ed (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
      PROGRAM MAPGEN
      PARAMETER (NMAX=257,NRX=81,NWX=81)
      REAL ETA(NMAX), F(NMAX), U(NMAX), S(NMAX)
      LOGICAL*1 FNAME(32)
      REAL AR(NRX,NWX), AI(NRX,NWX), X(NRX,NWX), Y(NRX,NWX)
      REAL RT(NRX),RTL(NRX), WS(NWX),WSL(NWX)
      CHARACTER*1 ANS
C
      READ(19) N, H
      READ(19) (ETA(I),I=1, N)
      READ(19) (U(I)  ,I=1, N)
      READ(19) (S(I)  ,I=1, N)
      READ(19) NRP, NWP
      READ(19) (RTL(IR),IR=1,NRP)
      READ(19) (WSL(IW),IW=1,NWP)
C
      DO 10 IW=NWP, 1, -1
        DO 101 IR=1, NRP
          READ(19,END=11) AR(IR,IW), AI(IR,IW)
  101   CONTINUE
   10 CONTINUE
   11 CONTINUE
C
C
      GEO = (ETA(3)-ETA(2)) / (ETA(2)-ETA(1))
      ETAE = ETA(N)
C
CCC   CALL FS(3,2,BU,H,N,ETAE,GEO,ETA,F,U,S)
C
      WRITE(6,*) 'GEO =', GEO, '   ETAE =', ETAE
      WRITE(6,*) 'H   =', H  , '   N    =', N
C
      WRITE(29) N, H
      WRITE(29) (ETA(I),I=1, N)
      WRITE(29) (U(I)  ,I=1, N)
      WRITE(29) (S(I)  ,I=1, N)
      WRITE(29)  NRP, NWP
      WRITE(29) (RTL(IR),IR=1,NRP)
      WRITE(29) (WSL(IW),IW=1,NWP)
C
      DO 20 IW=1, NWP
        WRITE(29) (AR(IR,IW),IR=1,NRP)
        WRITE(29) (AI(IR,IW),IR=1,NRP)
   20 CONTINUE
C
      STOP
      END