module m_constituents
      contains
      SUBROUTINE constituents(ifreq,nconst,freq2,tname2)

cGE      implicit double precision (a-h,o-z)
 
      CHARACTER*4  tname2(nconst)
      dimension freq2(nconst)
      INTEGER ifreq(nconst)


      PARAMETER(nf=30)
      character*11 tname
      dimension tname(nf),freq(nf),per(nf),hsm(nf),chi(nf)
c
      data (tname(i),freq(i),per(i),hsm(i),chi(i),i=1,nf) /
c              name        f(cycle/d)       period      hs (m)    phase
c long period
     .     'LP   55.565',  0.000147094,    6798.382,    0.0279,   180.0, ! 1
     .     'Sa   56.554',  0.002737779,     365.260,   -0.0049,     0.0, ! 2
     .     'Ssa  57.555',  0.005475819,     182.621,   -0.0310,     0.0, ! 3
     .     'TERa 58.554',  0.008213597,     121.749,   -0.0018,     0.0, ! 4
     .     'Mm   65.455',  0.036291647,      27.555,   -0.0352,     0.0, ! 5
     .     'Mf   75.555',  0.073202203,      13.661,   -0.0666,     0.0, ! 6
     .     'TERm 85.455',  0.109493850,       9.133,   -0.0128,     0.0, ! 7
     .     '93a  93.555',  0.140928587,       7.096,   -0.0020,     0.0, ! 8
c diurnal
     .     '2Q1 125.755',  0.856952384,      -6.904,   -0.0066,   270.0, ! 9
     .     'Q1  135.655',  0.893244048,      -9.213,   -0.0502,   270.0, !10
     .     'O1  145.555',  0.929535695,     -13.841,   -0.2622,   270.0, !11
     .     'M1  155.655',  0.966446233,     -28.297,    0.0206,    90.0, !12
     .     'P1  163.555',  0.997262079,    -221.060,   -0.1220,   270.0, !13
     .     'S1  164.556',  1.000000000,    -559.897,    0.0029,    90.0, !14
     .     'K1  165.555',  1.002737898,    1050.246,    0.3688,    90.0, !15
     .     'PHI1167.555',  1.008213699,     155.580,    0.0053,    90.0, !16
     .     'J1  175.455',  1.039029527,      26.850,    0.0206,    90.0, !17
     .     'OO1 185.555',  1.075940083,      13.485,    0.0113,    90.0, !18
     .     'NU1 195.455',  1.112231730,       9.054,    0.0022,    90.0, !19
c semi-diurnal
     .     '227 227.655',  1.828255527,      -5.704,    0.0047,     0.0, !20
     .     '2N2 235.755',  1.859690264,      -6.950,    0.0160,     0.0, !21
     .     'MU2 237.555',  1.864547173,      -7.193,    0.0193,     0.0, !22
     .     'N2  245.655',  1.895981946,      -9.295,    0.1210,     0.0, !23
     .     'NU2 247.455',  1.900838820,      -9.734,    0.0230,     0.0, !24
     .     'M2  255.555',  1.932273593,     -14.026,    0.6319,     0.0, !25
     .     'L2  265.455',  1.968565204,     -28.566,   -0.0179,   180.0, !26
     .     'T2  272.556',  1.997262163,    -158.475,    0.0172,     0.0, !27
     .     'S2  273.555',  2.000000000,    -279.994,    0.2940,     0.0, !28
     .     'K2  275.555',  2.005475795,     525.123,    0.0800,     0.0, !29
     .     '285 285.455',  2.041767407,      26.181,    0.0045,     0.0/ !30

      
      write(*,'(''No Name    Freq       Period'')' )
      write(*,'(''          cyc/day      hours'')' )
      DO i=1,nconst
       ii=ifreq(i)
       tname2(i)=tname(ii)(1:4)
       freq2(i)= freq(ii)
       dum=24.D0/freq2(i)
       write(*,'(I2,1x,A4,1x,F11.9,F11.6)')
     .    ii,tname2(i),freq2(i),dum
      ENDDO
      write(*,*)

      end subroutine
      end module