5 real,
parameter :: h1=1.e0,h1d5=1.5e0,h2d5=2.5e0,h3000=3000.e0, &
6 h10e5=100000.e0,d00=0.e0,d125=.125e0,d50=.5e0,d608=.608e0, &
7 g=9.8e0,cp=1004.6e0,capa=0.28589641e0,rog=287.04/9.8, &
8 elwv=2.50e6,elivw=2.72e6,row=1.e3,epsq=2.e-12, &
9 a2=17.2693882e0,a3=273.16e0,a4=35.86e0, &
10 t0=273.16e0,t1=274.16e0,pq0=379.90516e0,stresh=1.10e0, &
11 stabs=1.0e0,stabd=.90e0,stabfc=1.00e0,dttop=0.0e0, &
13 rhf=0.10,epsup=1.00,epsdn=1.05,epsth=0.0, &
14 pbm=13000.,pqm=20000.,pno=1000.,pone=2500.,zsh=2000., &
15 pfrz=15000.,pshu=45000., &
48 fss=.85e0,efimn=.20e0,efmnt=.70e0,fcc=.50,fcp=h1-fcc, &
58 dspbfl=-4843.75e0,dsp0fl=-7050.00e0,dsptfl=-2250.0e0,fsl=.850e0, &
62 dspbfs=-3875.e0,dsp0fs=-5875.e0,dsptfs=-1875.e0, &
63 dspbsl=dspbfl*fsl,dsp0sl=dsp0fl*fsl,dsptsl=dsptfl*fsl, &
64 dspbss=dspbfs*fss,dsp0ss=dsp0fs*fss,dsptss=dsptfs*fss, &
76 trel=3000.,epsntp=.0010e0,efifc=5.0e0, &
77 avgefi=(efimn+1.e0)*.5e0,dspc=-3000.e0,epsp=1.e-7, &
82 slopbl=(dspbfl-dspbsl)/(h1-efimn), &
83 slop0l=(dsp0fl-dsp0sl)/(h1-efimn), &
84 sloptl=(dsptfl-dsptsl)/(h1-efimn), &
85 slopbs=(dspbfs-dspbss)/(h1-efimn), &
86 slop0s=(dsp0fs-dsp0ss)/(h1-efimn), &
87 slopts=(dsptfs-dsptss)/(h1-efimn), &
88 slope=(h1 -efmnt)/(h1-efimn)
90 a23m4l=a2*(a3-a4)*elwv, &
91 elocp=elivw/cp,cprlg=cp/(row*g*elwv),rcp=h1/cp
92 logical,
parameter :: &
93 unis=.false.,unil=.false.,oct90=.false.