?COSANG 16-JUN-1997 16:57!Digital Fortran 77 V7.1-107 0,$DATA$$BSS$t$IODATA$ COSMOL$ANGINIWHICHARTHURFIXip+$CODE$($LINK$BNCSIZCNTRCTMERRORLOWSIZ>PATH(J INICOS(J QQ DFOR$WRITE_SEQ_LIS DFOR$HANDLER(J@@ ANGSIZ(J  LOWLEV DFOR$WRITE_SEQ_FMT(J  BNGSIZ(J@ ODEINT(J@P BSSTEP(J`' ASDRHS(J 0 MMM(J# POLEX(JH( ETAZ DFOR$STOP MATH$LOG10_G MATH$SQRT_G MATH$POW_GG MATH$ACOS_G MATH$COS_G(COSANG L=4,  0@` =|t6AHVH76H76H76AHVH76AHV76AHVH7 @=( 'ՙ''ՙ''ՙ''ՙ' D=,$'ՙ'R $=  $=  (='ՙ' = #~@> p^4Gx^~޴^~ޜcGGH~QGPp{8@4[Bb0 G _&X@^"`?"h>R"b"X"G@Zk!GH8"T#"TGBb XB"`]G8"_&X"hݶ?"R"b"X"G@Zkè೨@¤% $ `T% &AT CT%DU`4G@TGEtGG__`G@GŲ4GG__TG@tGEG@G@ _lU`G%x__-dTBbpmT@Zk TtGiáW8B@b"FAT\@Zkl"8Bc\"@T@bX@ZkX=B,T@ T T T`GR_ TG%IpfT8B"@b@Zk4G? %=G7_l€BbAfTpT@ZkAT cTBb\A TTAAUaAT0W@ZkfT0"HBPb0TTX=@ZkPTHBPb`X@Zk@T`=HB\Pb!W W Ta@]X=@Zk@TX݄TGT`&\`T T? %%fFGE"Q7 EUƅ„KJUJjUKUJU*KU NUFU(@0@@/UB"Gb"?"G @c_&R"``Xݶ(b"X"@ZktGe&FeT%$BcTCCT@TA#T  T TTcH= _8}E1eP¦f %eTu`G ]WcGp]x]}ݤ]}݌#k;aAAT@ TAAU@`U T Tk@#~x^^~^~cGGH^GPq?"P[ X"G_&Xp{"X{"PR"P"G@Zk ""Bb@}@DG"X="P}G`_&p"h=?"R"b"P"G("@Zk# $Gx"G!5@D #b"4GX}Ps_ $PB]`TXb) 8PjT@4GCS @"TG_&X"Pݶ?"R"pb"P"G@ZkFTG#h""AtG"X=P}Gh_&`"?"R"0b"P"("-@T$PBXb?"+8 AT G$ B"TGX]_&P"R"pb"P"G@ZkG#p""tG"X=P}Gh_&`"?"R"0b"P"("BbG@Zk@G""#`G#"" 4GX=P}_&"?"R"PBXbpb"P"G@Zk# $O8"@BHb@TqATJ$TAaUJU!AU-T@Zk8B lTP TtG0B"4G?"8pb" @$PB `]_&X"R"X"GXb@ZkP8"8]@T"*DU@@="8]"b"8 "` T0 "b#_+W@]u@0G#$ "H0 QWqcGx]]}]}#kG# ~^^~޴>^~cGGHG@[P>G$ X^" @} @4D Ģ"Gh`>_&pP{"R"x>?"`"GH{@Zk$X"@B G"h"`]_&p?"R"b"`""HbG@ZkGGD$!@ @]!@1D DGHb"@b"G#h}_&`?"pR"Pb"x]`"G@B,"@Zku@u0¤%T@d lU T8}1_("Bb.TQ1qOBTPUO/VPPV/VU@Zk$8b`tG-TB T""4G_&b8=" `#p=h}?"R"b"h"G@Zk`"@@T` Tx"@]@="8]"b"G( "b#=@$GB""G" _&h="`}?"R"b"`"Gb@ZkDb'xCT"T"!W@="8]"@]b"8G'( "b#@$GB x""TG_&h="`}?"R"b"`"Gb@ZkGx" ( AB\`_e("BbcUɅ LTMULUMU,UU@Zk%8b #Th TtGB"4G?"8=b" $B x]_&p"R"p"Gb@Zk =h8=x"@T!$T@="@8]"b"G``T( "b#_*W \G@$BbG ""TG_&p=h}"?"R"b"h"G@ZkGx? T&( @`] T U8-__T) IBbNTOUNVO/V.UU@Zk$8bh#T$tG T"4GBb8=?" ` xp]_&"R"b"p"G@Zkh@T``T"x@="@8]"b"G( "b#W@(G #@@"TGpBbh]_&"?"R"b"h"G@Zk1x"TT@="8]"b"G( "@b#82@GB"!G_&p="b?"hR"b"h"G@Zk x( (`}WH}WPXWCe$cG]]}ݤ=]}#k;#~`^h^p~xcGG#rPG@H@T@lU*`U _MU@ `Ua THBH _U@\~pb" "$ " "8=8""pb#8[@P@=pBx""G_&@!5@xbR"P"GA P}X="?"b"@Zk9áU@ tGpB"P]_&X"?"R"b"xbP"G@Zk(àu@pBGP" _&P=xb"X=?"R"b"P"G@ZkpBG _&P"xb?"XR"b"P"G@ZkH=#C0"1QcG`]h]p}x#kGGG1R@#^4G ^(~08@޴H_PX>*AU`^h~p+Wx޵^@]~GG~#Gk]p]GG}(=GG X6q`Bx #=` B 7WwB B8 " #Xx? B@70B B"8#_4GX!"I"Gx 1"GtG@Zk#4G": 2q @ AB` C "(`pTxtԅ4_B___ B "JjTKkTLlTMmT _ _ _ _ U U  V -V AU aU U UTt B GpTt0B"_ BRrT_rVAVT_m0"7q @4G1`4Tm_;1`C {B_0`Bx #RqBV$#`m 3`J BSsBSsB3B1" ؆8B Bі # 1"1_ BX0B B#Q1"\T<dW;eW[W`#T"GG5BEHGx " T" 4"^]">G}"GTG@Zk ``T 0B K0@BK%T0A=خ"tG$$T ` _DT!lUp`I @G] ](}08@ݤHPX=`]h}pxݥ]}#k¦ Vc"]4Gx 9#`8q`@ C`مB Bڕ 9# Z#? B@90B B9#:Z#*q`@4G*0b!0@x k!I@) `%HE@ CAE@A @GKk@ @Ee k! @`ˆ0@ @k!Ŗ ? XÈbTGL SGP#;(^0^h! 8~@HP޴X`h>p^x~޵^G~#GG0[GG'G3q@4GG=( `B```ԅ1TB "q B 1" ! A!a_ B G0B" B1"! _8B h,!L!l!!4GOA޴"~"GGG"b#@]_^Gݱ ="G}"H"Gb#_WaAUA!T =u@4GGG0!qA /@"@! ``1PqxWЇ@ " @ 1"z[W_Vw8\|W_Vw:\`W_Vw\]_Vw`] @ 0@lU @"1" _Vw]?}0`CoVU"ղ}A_G(=Q0VP0VPVV0B4GZG(]0]8}@HPݤX`h=p]x}ݥ]#k=@† @VVVזà @8"YV/A-@ p=_W_UVWuVGG4#^4G^ ~(08޴@HPX^`~G~GGՠpxۦ"޴4qBGWb T1 B# "`` = 3Ss@@TAATJJTU@ R"KKT҅U @px s" U֕ " Uh6Vv @ 0@R"LLT @ps"" U?BT"D GGGtGTGT@ZkY@G#4GGG5 1q @ !B&% pxӅ`2Rr@`TAaT3SJjTKkTsBUU Bp U Ux3 R"3 s"Ss`Shs@ B pӆ0BTtT BpӖR"s"VD0A"VGGGtG@Zk Ac4G<`0q`@ B CW) p'x煈G``GeŅ%gKKT*T UNNT @ UQQT 0@ SVTTT+T UUTVA TLUOUR@V%E #@ pdž0 @ WWT #@ VVUV?G]] }(08ݤ@HX]`}Pp#ks ;#4G~G^G^VqBG $6qB BG `8Xx7 @: 6@W " # Z#8Zwz7:? 0@ GX0 @ 0@"W#ZZ#5`B! 6qB B4G"Gr2Rw @r 6@7 R" "Ww 0@20 @ 0@R"7"G ;8qB4Gh! tB4G CRr2A@CR XCa R" ! !A PCr0@C PCR"a! 41BnsaB3B;4G1"Q2 K RRBRRB:qBC4Gh#GR"G+@Uk+TL TJKUT PC?@>@@@@r҅+U/VPU @ R"-U2؅KKUWs.UPVUU :@xQU #MMUKUr s"WVOOU "NNUWӆVKV3ؕWҖ3s2 _ 0@@Gs4V0 @R" 0@r#s""OUTTUV_0 C 4C1"_6J" VB"qBVB P@C4G"@@G@` @ "@֕V v "֖ 0@ 0 @ 0@ 6"?G]] #kGGG#$(^G0^8~@HP^X~`G~GG^P[dDą{@BTATL@U JUAaUKmUJBU!,T@`U ^.Tc\@Zk^ ~T`"\*KU!xb#aAU "> "NT ^#@cTÅd$@KlU "]G(]dDU0]8}KmU@a"TBbUBNTbCTX}$T`A/T"W!ATP]0Hp#k#~H^P^X`^cGG00["TjPTW_{p@ WqG4@"_!b Wd@ WeGBTB\b]@ KT^T@"!b T`@TTO@  TPJ]BUbI@  TJG_u@ JTbT?@  T@ð@@JUbT7@  T8á@ JJUb JTT.@  T/ð@lbBb+T @ZkTb!@  T"G_BTB\b@ KTG_5@!JTJ]A!T8=_0JTJJU8]8B@b@Zk8=b@TT@  T T T @GcH]P]`]Xp#kGHBG _&8"Pb?"@R"hb"8"G@ZkXB`bGGTG@Zk<GDFOR$WRITE_SEQ_LIS<b# DFOR$WRITE_SEQ_LIS<@DFOR$WRITE_SEQ_LIS<GDDFOR$WRITE_SEQ_LIS<b# DFOR$WRITE_SEQ_LIS<@@DDFOR$WRITE_SEQ_LIS4G MATH$SQRT_G4b#  MATH$SQRT_G4@ MATH$SQRT_G4 G MATH$POW_GG4 b#  MATH$POW_GG4 @ MATH$POW_GG4 G MATH$POW_GG4 b#  MATH$POW_GG4 @ MATH$POW_GG4 HGX MATH$POW_GG4 Pb#  MATH$POW_GG4 T@X MATH$POW_GG4G MATH$SQRT_G 4b#  MATH$SQRT_G4@ MATH$SQRT_G4G MATH$ACOS_G4b#  MATH$ACOS_G4@ MATH$ACOS_G4 G MATH$COS_G4b#  MATH$COS_G4 @ MATH$COS_G4 G  MATH$COS_G4b#  MATH$COS_G4 @  MATH$COS_G4 G8 MATH$COS_G4b#  MATH$COS_G4 4@8 MATH$COS_G<GDFOR$WRITE_SEQ_LIS<b# DFOR$WRITE_SEQ_LIS<@DFOR$WRITE_SEQ_LIS<GDFOR$WRITE_SEQ_LIS<b# @DFOR$WRITE_SEQ_LIS<@DFOR$WRITE_SEQ_LIS< GLDFOR$WRITE_SEQ_FMT< b# @DFOR$WRITE_SEQ_FMT< H@LDFOR$WRITE_SEQ_FMT<G DFOR$WRITE_SEQ_LIS<b# @DFOR$WRITE_SEQ_LIS<@ DFOR$WRITE_SEQ_LIS<lGDFOR$WRITE_SEQ_LIS<pb# @DFOR$WRITE_SEQ_LIS<@DFOR$WRITE_SEQ_LIS<  G( DFOR$WRITE_SEQ_FMT<  b# @DFOR$WRITE_SEQ_FMT< $ @( DFOR$WRITE_SEQ_FMT<x G DFOR$WRITE_SEQ_LIS<| b# @DFOR$WRITE_SEQ_LIS< @ DFOR$WRITE_SEQ_LIS4 G MATH$SQRT_G4 b# @ MATH$SQRT_G4 @ MATH$SQRT_G< GH DFOR$WRITE_SEQ_LIS<@ b# @DFOR$WRITE_SEQ_LIS<D @H DFOR$WRITE_SEQ_LIS< l G DFOR$WRITE_SEQ_FMT<  b# DFOR$WRITE_SEQ_FMT<  @ DFOR$WRITE_SEQ_FMT<  G< DFOR$WRITE_SEQ_FMT< 0 b# DFOR$WRITE_SEQ_FMT< 8 @< DFOR$WRITE_SEQ_FMT< p b# DFOR$WRITE_SEQ_FMT<  G DFOR$WRITE_SEQ_FMT<  @ DFOR$WRITE_SEQ_FMT4 GL MATH$SQRT_G4 b#  MATH$SQRT_G4H @L MATH$SQRT_G<h G DFOR$WRITE_SEQ_LIS<| b# DFOR$WRITE_SEQ_LIS< @ DFOR$WRITE_SEQ_LIS< GPDFOR$WRITE_SEQ_LIS<Hb# DFOR$WRITE_SEQ_LIS<L@PDFOR$WRITE_SEQ_LIS<GDFOR$WRITE_SEQ_LIS<b# DFOR$WRITE_SEQ_LIS<@DFOR$WRITE_SEQ_LIS48Gl MATH$SQRT_G4<b#  MATH$SQRT_G4h@l MATH$SQRT_G<GDFOR$WRITE_SEQ_LIS<b# DFOR$WRITE_SEQ_LIS<@DFOR$WRITE_SEQ_LIS<0GxDFOR$WRITE_SEQ_LIS<4b# DFOR$WRITE_SEQ_LIS<t@xDFOR$WRITE_SEQ_LIS 4G MATH$SQRT_G4b#  MATH$SQRT_G4@ MATH$SQRT_G<GXDFOR$WRITE_SEQ_LIS<b# DFOR$WRITE_SEQ_LIS<T@XDFOR$WRITE_SEQ_LIS<GDFOR$WRITE_SEQ_LIS<b# DFOR$WRITE_SEQ_LIS<@DFOR$WRITE_SEQ_LIS<4GxDFOR$WRITE_SEQ_LIS<Xb# DFOR$WRITE_SEQ_LIS<t@xDFOR$WRITE_SEQ_LIS<G,DFOR$WRITE_SEQ_LIS<b# DFOR$WRITE_SEQ_LIS<(@,DFOR$WRITE_SEQ_LIS<<GpDFOR$WRITE_SEQ_LIS<`b# DFOR$WRITE_SEQ_LIS<l@pDFOR$WRITE_SEQ_LIS<|GDFOR$WRITE_SEQ_LIS<b# DFOR$WRITE_SEQ_LIS<@DFOR$WRITE_SEQ_LIS<GDFOR$WRITE_SEQ_LIS<b# DFOR$WRITE_SEQ_LIS<@DFOR$WRITE_SEQ_LIS4'G' MATH$SQRT_G4'b# MATH$SQRT_G4'@' MATH$SQRT_G4@*GT* MATH$LOG10_G4D*b# H MATH$LOG10_G4P*@T* MATH$LOG10_G4*G* MATH$SQRT_G4*b# H MATH$SQRT_G4*@* MATH$SQRT_G< +GX+DFOR$WRITE_SEQ_LIS<8+b# HDFOR$WRITE_SEQ_LIS<T+@X+DFOR$WRITE_SEQ_LIS4X+Gp+ DFOR$STOP4\+b# H DFOR$STOP4l+@p+ DFOR$STOP   MATH$LOG10_G= SOMETHING'S WRONG!!!!!!;; ;=:?6-C=0H(;=p( 7 DFOR$HANDLER= @  MATH$SQRT_G DFOR$WRITE_SEQ_LIS DFOR$STOP=8B=0;;=0(`';= p0 A@A=0#;= 8 `;=00;=p< =0(P;=\ ;=? p ;= @333333@ffffff =0@;= x <9KH>;;= CAN'T GET HERE =0$ERROR IN ODEINT: STEPSIZE TOO SMALL = ,@  @;=0`;=( 7 DFOR$HANDLER= ; ;=y@o z>׼HERROR IN ODEINT: TOO MANY STEPS>ƵERROR IN BSSTEPERROR 12: ERROR CALCULATING D32, EXPANDING PHASE =@3ERROR 12: ERROR CALCULATING D32, CONTRACTING PHASE =@3ERROR 12: ERROR CALCULATING D14, CONTRACTING PHASE P=: INICOS ERROR `=ERROR  h0;=0"ERROR 11: OVERFLOW ERROR POSSIBLE =@=0 ;=L 7 DFOR$HANDLER= ;=+88x DFOR$WRITE_SEQ_FMT=8-8 8ERROR 13: BNGSIZ INVALID FOR WMTYPE  =8  (=@3ERROR 12: ERROR CALCULATING D34, CONTRACTING PHASE `=: ANGSIZ ERROR p=@1ERROR 12: ERROR CALCULATING D14, EXPANDING PHASE =Z2  =Z1  ;=@2ERROR 4: COMPUTATION NOT RECOMMENDED IN WMTYPE 17 =PA @NOT REALLY LARGER THAN ZMAX BUT CLOSE ENOUGH FOR TROUBLE @=0x@;=4 7 DFOR$HANDLER=88 x;= ERROR 9: Z2 > ZMAX =@6YOU MUST CALL INICOS FIRST!!!!!!!!!!!!!!!!!!!!!!!!!!! = ERROR 7: Z1 > ZMAX = ERROR 5: Z1 < 0.0 = ERROR 6: Z2 < 0.0 =8)INICOS ERROR 3: ROUNDOFF ERROR, BAD ZMAX H=0&INICOS ERROR 1: OMEGA < 0.0; OMEGA =  p=8. x=?h?UUUUUU=0p;=P |7 DFOR$HANDLER=@R-8es MATH$POW_GG MATH$COS_G=8)INICOS ERROR 2: ETA OUT OF RANGE; ETA =   MATH$ACOS_G=@UUUUUU I=910COSANGDigital Fortran 77 V7.1-107 = 4 4= INICOS4=C;lll ll   + - 3 3 6 6 8ȅ D D D  W \ f j  v               & & '  7 3 A A S O e g k k k i i ] _ ](     "     X s $  " $ 6 "&  I D $= p4 := 4=  r= 4 4=QQ4=  := <4= 0 .= @4 @4= ANGSIZ@4=      Z                !    !  ' )  9 9 9 5 =8 S O S S 4 $= 4 :=  4=  L=  4 4= BNGSIZ 4= e e˰ e 1 e 1 %  = ? A A A A L O  w w } } y                                             " $  4 0 4 4 8  T  ^ T  d `   r v r  z x r  x        S       D $=  4 := 4=  = 4 4= LOWLEV4=                ,8 8  , $= <4 := 4= P B= @4 4= ODEINT@4=         7             0  P                      D } & & }d                      M                 $= 4 := 4=  = P4 @4= BSSTEPP4=  t  9  9  9  9  9                     $ (    <            K M P '  )  +  )  '  +  /  $= 4 := P4=  /= 04 4= MMM04=  i f i b  i j                                      >P  8     8 $= 4 := "4=  b= #4 4= POLEX#4= b  : ' W U U d  d  d  f  f  d  d  d  d   p  p  p  p  p  p  p  p   |  z  |  |  |  |  |  |                                            h            $ $= <#4 := <'4= ` X= `'4 4= ASDRHS`'4= i   4                          j    &  &  '  &  '  &  '  '  &  &  $= '4 := (4= p = (4 H4= ETAZ(4= 1 Z 1 X F F     > <  <     /  +  '   #                    $= (4 := *4=  =