var/external/wavelet/print_filter.F90

   1 ! Author: Aime' Fournier
   2 ! E-mail: fournier@ucar.edu
   3
   4 #include "realt.h"
   5
   6  SUBROUTINE print_filter(nam,ran,kin,f)
   7  IMPLICIT NONE
   8  CHARACTER nam                          ! wavelet name
   9  INTEGER k,kin,ran
  10  realt, DIMENSION(ran) :: f             ! wavelet filter
  11  realt fs,fse
  12 !
  13  fs=0.
  14  IF( kin==0 ) THEN
  15     fse=SQRT(2.)
  16  ELSEIF( kin==1 ) THEN
  17     fse=0.
  18  ELSE
  19     STOP
  20  ENDIF
  21  WRITE(*,                                                              &
  22        '("Filter ",A1,I2.2,"[",I1,"] has alpha=0, omega=",I2,", f={",99(1X,F10.6,:))',&
  23        ADVANCE="NO")nam,ran,kin,ran-1,f
  24  DO k=1,ran
  25     fs=fs+f(k);                         ! f sum
  26  ENDDO
  27  PRINT '("} and sum-error ",ES9.2)',fs-fse
  28  END SUBROUTINE print_filter