24 |
c |
c |
25 |
c |
c |
26 |
INTEGER klevel |
INTEGER klevel |
27 |
REAL xcov( ip1jmp1,klevel ), ycov( ip1jm,klevel ) |
REAL, intent(in):: xcov( ip1jmp1,klevel ), ycov( ip1jm,klevel ) |
28 |
REAL grx( ip1jmp1,klevel ), gry( ip1jm,klevel ) |
REAL grx( ip1jmp1,klevel ), gry( ip1jm,klevel ) |
29 |
real, intent(in):: crot |
real, intent(in):: crot |
30 |
c |
c |
35 |
c |
c |
36 |
c |
c |
37 |
c |
c |
38 |
CALL SCOPY ( ip1jmp1*klevel, xcov, 1, grx, 1 ) |
grx = xcov |
39 |
CALL SCOPY ( ip1jm*klevel, ycov, 1, gry, 1 ) |
gry = ycov |
40 |
c |
c |
41 |
DO 10 iter = 1,lr |
DO 10 iter = 1,lr |
42 |
CALL rotat (klevel,grx, gry, rot ) |
CALL rotat (klevel,grx, gry, rot ) |