1 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
4 ! Copyright (C) 1995-2004, Scientific Computing Division,
5 ! University Corporation for Atmospheric Research
6 ! Licensed under the GNU General Public License (GPL)
8 ! Authors: Paul N. Swarztrauber and Richard A. Valent
10 ! $Id: mradb3.f,v 1.2 2004/06/15 21:29:19 rodney Exp $
12 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
14 SUBROUTINE MRADB3 (M,IDO,L1,CC,IM1,IN1,CH,IM2,IN2,WA1,WA2)
15 REAL CC(IN1,IDO,3,L1) ,CH(IN2,IDO,L1,3), &
20 ARG=2.*4.*ATAN(1.0)/3.
27 CH(M2,1,K,1) = CC(M1,1,1,K)+2.*CC(M1,IDO,2,K)
28 CH(M2,1,K,2) = CC(M1,1,1,K)+(2.*TAUR)*CC(M1,IDO,2,K) &
29 & -(2.*TAUI)*CC(M1,1,3,K)
30 CH(M2,1,K,3) = CC(M1,1,1,K)+(2.*TAUR)*CC(M1,IDO,2,K) &
31 & +2.*TAUI*CC(M1,1,3,K)
34 IF (IDO .EQ. 1) RETURN
42 CH(M2,I-1,K,1) = CC(M1,I-1,1,K)+(CC(M1,I-1,3,K)+CC(M1,IC-1,2,K))
43 CH(M2,I,K,1) = CC(M1,I,1,K)+(CC(M1,I,3,K)-CC(M1,IC,2,K))
44 CH(M2,I-1,K,2) = WA1(I-2)* &
45 & ((CC(M1,I-1,1,K)+TAUR*(CC(M1,I-1,3,K)+CC(M1,IC-1,2,K)))- &
46 & (TAUI*(CC(M1,I,3,K)+CC(M1,IC,2,K)))) &
48 & ((CC(M1,I,1,K)+TAUR*(CC(M1,I,3,K)-CC(M1,IC,2,K)))+ &
49 & (TAUI*(CC(M1,I-1,3,K)-CC(M1,IC-1,2,K))))
50 CH(M2,I,K,2) = WA1(I-2)* &
51 & ((CC(M1,I,1,K)+TAUR*(CC(M1,I,3,K)-CC(M1,IC,2,K)))+ &
52 & (TAUI*(CC(M1,I-1,3,K)-CC(M1,IC-1,2,K)))) &
54 & ((CC(M1,I-1,1,K)+TAUR*(CC(M1,I-1,3,K)+CC(M1,IC-1,2,K)))- &
55 & (TAUI*(CC(M1,I,3,K)+CC(M1,IC,2,K))))
56 CH(M2,I-1,K,3) = WA2(I-2)* &
57 & ((CC(M1,I-1,1,K)+TAUR*(CC(M1,I-1,3,K)+CC(M1,IC-1,2,K)))+ &
58 & (TAUI*(CC(M1,I,3,K)+CC(M1,IC,2,K)))) &
60 & ((CC(M1,I,1,K)+TAUR*(CC(M1,I,3,K)-CC(M1,IC,2,K)))- &
61 & (TAUI*(CC(M1,I-1,3,K)-CC(M1,IC-1,2,K))))
62 CH(M2,I,K,3) = WA2(I-2)* &
63 & ((CC(M1,I,1,K)+TAUR*(CC(M1,I,3,K)-CC(M1,IC,2,K)))- &
64 & (TAUI*(CC(M1,I-1,3,K)-CC(M1,IC-1,2,K)))) &
66 & ((CC(M1,I-1,1,K)+TAUR*(CC(M1,I-1,3,K)+CC(M1,IC-1,2,K)))+ &
67 & (TAUI*(CC(M1,I,3,K)+CC(M1,IC,2,K))))