      SUBROUTINE HESSI_ROTMS
      real*4 alpha

      call GSROTM (1,90.0,116.1,90.0,206.1,0.,0.)
      call GSROTM (2,90.,51.34,90.,141.34,0.,0.)
      call GSROTM (3,90.,128.66,90.,218.66,0.,0.)
      call GSROTM (4,90.,63.9,90.,153.9,0.,0.)
      call GSROTM (5,90.,60.9,90.,150.9,0.,0.)
      call GSROTM (6,90.,146.6,90.,236.6,0.,0.)
      call GSROTM (7,90.,119.1,90.,209.1,0.,0.)
      call GSROTM (8,90.,90.,90.,180.,0.,0.)
      call GSROTM (9,90.,33.4,90.,123.4,0.,0.)     
      call GSROTM (10,90.,270.,90.,0.,0.,0.)      
      call GSROTM (11,90.0,333.9,90.0,63.9,0.0,0.0)
      call GSROTM (12,90.0,26.1,90.0,116.1,0.0,0.0)
      call GSROTM (13,90.0,0.0,180.0,0.0,90.0,90.0)
      call GSROTM (14,90.0,0.0,0.0,0.0,90.0,270.0)
      call gsrotm (16, 90., 336.58, 90., 66.58, 0., 0.)
      call gsrotm (17, 90., 0., 180., 0., 90., 90.)
      call gsrotm (18, 90., 23.42, 90., 113.42, 0., 0.)
      call gsrotm (19, 90., 60., 90., 150., 0., 0.)
      call gsrotm (20, 90., 300., 90., 30., 0., 0.)
      call gsrotm (21, 90., 0., 0., 0., 90., 270.)
      call gsrotm (22, 90., 240., 0., 0., 90., 150.)
      call gsrotm (23, 90., 120., 0., 0., 90., 30.)
      CALL GSROTM(24,90.,59.,90.,149.,0.,0.)
      CALL GSROTM(25,90.,121.,90.,211.,0.,0.)
      CALL GSROTM(26,90.,145.365,90.,235.365,0.,0.)
      CALL GSROTM(27,90.,34.635,90.,124.635,0.,0.)
      CALL GSROTM(28,90.,55.366,90.,145.366,0.,0.)
      CALL GSROTM(29,90.,124.636,90.,214.636,0.,0.)
      CALL GSROTM(30,90.,116.118,90.,206.118,0.,0.)
      CALL GSROTM(31,90.,63.897,90.,153.897,0.,0.)
      CALL GSROTM(32,90.,0.0,0.0,0.0,90.0,270.0) 
      CALL GSROTM(33,180.,0.,90.,90.,90.,0.)
      CALL GSROTM(34,111.62,0.,90.,90.,21.62,0.)
      CALL GSROTM(35,68.38,0.0,90.,90.,21.62,180.0)  
      CALL GSROTM(36,90.,0.,0.,0.,90.,270.)
      CALL GSROTM(37,90.,0.,30.,90.,60.,270.)
      CALL GSROTM(38,90.,321.34,90.,51.34,0.,0.)
      CALL GSROTM(39,90.,38.66,90.,128.66,0.,0.)
      CALL GSROTM(40,90.,333.90,90.,63.90,0.,0.)
      CALL GSROTM(41,90.,26.1,90.,116.1,0.,0.)
      CALL GSROTM(42,69.5,0.,90.,90.,20.5,180.)
      CALL GSROTM(43,0.,0.,90.,270.,90.,180.)
      CALL GSROTM(44,90.,330.,90.,60.,0.,0.)
      CALL GSROTM(45,90.,30.,90.,120.,0.,0.)
      CALL GSROTM(46,90.,150.,90.,240.,0.,0.)
      CALL GSROTM(47,90.,210.,90.,300.,0.,0.)

      ALPHA=75.
      DO N=0,11
        CALL GSROTM(48+N,90.,ALPHA+90.,90.,ALPHA+180.,0.,0.)
        ALPHA=ALPHA-30.
      END DO

      CALL GSROTM(60,90.,324.56,90.,54.56,0.,0.)
      CALL GSROTM(61,90.,95.44,90.,5.44,0.,0.)
      CALL GSROTM(62,90.,24.56,90.,114.56,0.,0.)
      CALL GSROTM(63,90.,155.39,90.,245.39,0.,0.)
      CALL GSROTM(64,90.,84.56,90.,174.56,0.,0.)
      CALL GSROTM(65,90.,35.44,90.,125.44,0.,0.)
      CALL GSROTM(66,90.,320.,0.,0.,90.,230.)
      CALL GSROTM(67,90.,320.,35.,50.,55.,230.)
      CALL GSROTM(68,90.,31.,90.,121.,0.,0.)
      CALL GSROTM(69,90.,329.,90.,59.,0.,0.)
      CALL GSROTM(70,90.,12.,90.,102.,0.,0.)
      CALL GSROTM(71,90.,45.1,90.,135.1,0.,0.)
      CALL GSROTM(72,90.,135.,90.,225.,0.,0.)
      CALL GSROTM(73,90.,75.12,90.,165.12,0.,0.)
      CALL GSROTM(74,90.,90.,90.,180.,0.,0.)
      call gsrotm(82,90.,261.622,180.,0.,90.,351.622)
      call gsrotm(83,90.,0.,180.,0.,90.,90.)
      call gsrotm(84,90.,278.35768,90.,188.35768,180.,0.0)
      call gsrotm(85,90.,258.9079,90.,168.9079,180.,0.0)
      call gsrotm(86,90.,281.092,90.,11.092,0.0,0.0)
      call gsrotm(87,90.,274.7634,90.,4.7634,0.0,0.0)
      call gsrotm(88,90.,265.2366,90.,175.2366,180.,0.0)
      call gsrotm(91,90.,90.,0.0,0.0,90.,0.0)
      call gsrotm(92,90.,90.,180.,0.,90.,180.)          
      call gsrotm(93,90.,180.,0.0,0.0,90.,90.)
      call gsrotm(94,90.,108.5882,0.,0.,90.,18.5882)     
      call gsrotm(95,90.,71.4118,180.,0.,90.,161.4118) 
      call gsrotm(96,90.,180.,90.,90.,180.,0.)
      call gsrotm(97,90.,51.3341,180.,0.,90.,141.3341)
      call gsrotm(98,90.,128.6659,0.,0.,90.,38.6659)
      call gsrotm(100,90.,108.89,180.,0.,90.,198.89)
      call gsrotm(101,90.,278.3781,0.,0.,90.,188.3781)
      call gsrotm(102,90.,261.642319,90.,351.642319,0.,0.)
      call gsrotm(103,90.,12.01464752,180.,0.0,90.,102.01464752)
      call gsrotm(104,90.,108.396472,90.,198.396472,0.0,0.0)
      call gsrotm(105,90.,262.686977,90.,172.686977,180.,0.0)
      call gsrotm(106,90.,71.11,0.,0.,90.,341.11)
      call gsrotm(133,90.,45.,90.,135.,0.,0.)

      call gsrotm(150,90.,-5.985,  90.,84.015,0.,0.)
      call gsrotm(151,90.,-12.13,  90.,77.87,0.,0.)
      call gsrotm(152,90.,-30.56, 90.,59.44,0.,0.)
      call gsrotm(153,90.,-42.85, 90.,47.15,0.,0.)
      call gsrotm(154,90.,-48.99, 90.,41.01,0.,0.)

*** Some useful rotation matrices (PSI:)
      CALL GSROTM(216, 90., 90., 90., 180.,  0.,  0.) !   90 deg around z
      CALL GSROTM(217, 90.,180., 90., 270.,  0.,  0.) !  180 deg around z
      CALL GSROTM(218, 90.,270., 90.,   0.,  0.,  0.) !  270 deg around z
      CALL GSROTM(219, 90.,  0.,  0.,   0., 90.,270.) !   90 deg around x
      CALL GSROTM(220, 90.,120., 90., 210.,  0.,  0.) !  120 deg around z
      CALL GSROTM(221, 90.,240., 90., 330.,  0.,  0.) !  240 deg around z
      CALL GSROTM(222, 90., 45., 90., 135.,  0.,  0.) !   45 deg around z
      CALL GSROTM(224, 180., 0., 90.,  90., 90.,  0.) !   90 deg around y      
      CALL GSROTM(228, 90.,120.,  0.,   0., 90., 30.) !   90 deg around x 
                                                    ! & 120 deg around z
      CALL GSROTM(229, 90.,240.,  0.,   0., 90.,150.) !   90 deg around x
                                                    ! & 240 deg around z
      CALL GSROTM(230,120.,  0., 90.,  90., 30.,  0.) !   30 deg around y
      CALL GSROTM(231, 60.,  0., 90.,  90., 30.,180.) !  -30 deg around y
      CALL GSROTM(232, 90.,180., 90.,  90.,180.,  0.) !  180 deg around y
      CALL GSROTM(233, 90.,300., 90., 210.,180.,  0.) !  180 deg around y
                                                    ! & 120 deg around z
      CALL GSROTM(234, 90., 60., 90., 330.,180.,  0.) !  180 deg around y
                                                    ! & 240 deg around z
      CALL GSROTM(256, 90., 45.,180.,  0., 90., 135.) !  -90 deg around x
                                                    ! &  45 deg around z
      CALL GSROTM(257, 90.,135.,180.,  0., 90., 225.) !  -90 deg around x
                                                    ! & 135 deg around z      
      CALL GSROTM(258, 90.,225.,180.,  0., 90., 315.) !  -90 deg around x
                                                    ! & 225 deg around z      
      CALL GSROTM(259, 90.,315.,180.,  0., 90.,  45.) !  -90 deg around x
                                                    ! & 315 deg around z
      CALL GSROTM(260, 90.,315., 90., 45.,  0.,  0.) !   -45 deg around z
      
      end

