var/convertor/kma_wave2grid/RADF2M.inc

   1 C***********************************************************************
   2       SUBROUTINE RADF2M (INC,LOT,IDO,L1,CC,CH,WA1)
   3       DIMENSION       CH(INC,IDO,2,L1)  ,CC(INC,IDO,L1,2)  ,WA1(*)
   4       DO 101 K=1,L1
   5         DO 101 L=1,LOT
   6           CH(L,  1,1,K) = CC(L,1,K,1)+CC(L,1,K,2)
   7           CH(L,IDO,2,K) = CC(L,1,K,1)-CC(L,1,K,2)
   8   101   CONTINUE
   9       IF (MOD(IDO,2) .EQ. 0) THEN
  10         DO 102 K=1,L1
  11           DO 102 L=1,LOT
  12             CH(L,  1,2,K) = -CC(L,IDO,K,2)
  13             CH(L,IDO,1,K) =  CC(L,IDO,K,1)
  14   102     CONTINUE
  15       END IF
  16       IF (IDO .GT. 2) THEN
  17           IDP2 = IDO+2
  18         DO 104 K=1,L1
  19 *VOPTION NOFVAL
  20         DO 104 I=3,IDO,2
  21             IC = IDP2-I
  22 *VOPTION NOFVAL
  23           DO 103 L=1,LOT
  24             TR2 = WA1(I-2)*CC(L,I-1,K,2)+WA1(I-1)*CC(L,I  ,K,2)
  25             TI2 = WA1(I-2)*CC(L,I  ,K,2)-WA1(I-1)*CC(L,I-1,K,2)
  26             CH(L,I -1,1,K) = CC(L,I-1,K,1)+TR2
  27             CH(L,I   ,1,K) = TI2          +CC(L,I,K,1)
  28             CH(L,IC-1,2,K) = CC(L,I-1,K,1)-TR2
  29             CH(L,IC  ,2,K) = TI2          -CC(L,I,K,1)
  30   103     CONTINUE
  31   104   CONTINUE
  32       END IF
  33       RETURN
  34       END SUBROUTINE RADF2M
  35