var/convertor/wave2grid_kma/BSSLZ1.inc

   1       SUBROUTINE BSSLZ1(BES,N)
   2 C
   3       IMPLICIT REAL*8(A-H,O-Z)
   4       DIMENSION BES(N)
   5       DIMENSION BZ(50)
   6 C
   7       DATA BZ         / 2.40482555770, 5.52007811030,
   8      $  8.65372791290,11.79153443910,14.93091770860,18.07106396790,
   9      $ 21.21163662990,24.35247153080,27.49347913200,30.63460646840,
  10      $ 33.77582021360,36.91709835370,40.05842576460,43.19979171320,
  11      $ 46.34118837170,49.48260989740,52.62405184110,55.76551075500,
  12      $ 58.90698392610,62.04846919020,65.18996480020,68.33146932990,
  13      $ 71.47298160360,74.61450064370,77.75602563040,80.89755587110,
  14      $ 84.03909077690,87.18062984360,90.32217263720,93.46371878190,
  15      $ 96.60526795100,99.74681985870,102.8883742540,106.0299309160,
  16      $ 109.1714896490,112.3130502800,115.4546126530,118.5961766300,
  17      $ 121.7377420880,124.8793089130,128.0208770050,131.1624462750,
  18      $ 134.3040166380,137.4455880200,140.5871603520,143.7287335730,
  19      & 146.8703076250,150.0118824570,153.1534580190,156.2950342680/
  20       NN=N
  21       IF(N.LE.50) GO TO 12
  22       BES(50)=BZ(50)
  23       DO 5 J=51,N
  24     5 BES(J)=BES(J-1)+PI
  25       NN=49
  26    12 DO 15 J=1,NN
  27    15 BES(J)=BZ(J)
  28       RETURN
  29       END SUBROUTINE BSSLZ1