module theta
  implicit none
  integer,parameter,public :: ntheta=138
  real,parameter,public    :: theta=(/&
&270., 280., 290.,&
&300., 310., 320., 330., 340., 350., 360., 370., 380., 390.,&
&400., 410., 420., 430., 440., 450., 460., 470., 480., 490.,&
&500., 510., 520., 530., 540., 550., 560., 570., 580., 590.,&
&600., 625., 650., 675., 700., 725., 750., 775., 800., 825.,&
&850., 875., 900., 925., 950., 975.,&
&1000., 1050., 1100., 1150., 1200., 1250., 1300., 1350., 1400., 1450.,&
&1500., 1550., 1600., 1650., 1700., 1750., 1800., 1850., 1900., 1950.,&
&2000., 2100., 2200., 2300., 2400., 2500., 2600., 2700., 2800., 2900.,&
&3000., 3100., 3200., 3300., 3400., 3500., 3600., 3700., 3800., 3900.,&
&4000., 4250., 4500., 4750., 5000., 5250., 5500., 5750., 6000., 6250.,&
&6500., 6750., 7000., 7250., 7500., 7750., 8000., 8250., 8500., 8750.,&
&9000., 9250., 9500., 9750.,&
&10000., 10500., 11000., 11500., 12000., 12500., 13000., 13500., 14000., 14500.,&
&15000., 15500., 16000., 16500., 17000., 17500., 18000., 18500., 19000., 19500.,&
&20000., 21000., 22000., 23000., 25000./)
end module theta