[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

<none>



small error, indeed
i found it useful to have a look
at the contributors of df:

        sit=1.e-5
        pt=1.e-6
        dic=2.1e-3



	df1 = ((k1*dic*b) - k1*x*dic*db)/b2 
      df2= -    	     2.0*dic*k12*db/b2 
      df3= -	     bt/kb/(1.0+x/kb)**2. 
      df4= -	     kw/x2 
      df5=	     (pt*k12p*(a - x*da))/a2 
      df6= -	     2.0*pt*k123p*da/a2 
      df7=  -	     sit/ksi/(1.0+x/ksi)**2. 
      df8= -	     1.0/c 
      df9=	     (st/(1.0 + ks/x/c)**2)*(ks/c/x2) 
      df10=	     (ft/(1.0 + kf/x)**2)*kf/x2 
      df11= -	     pt*x2*(3.0*a-x*da)/a2
      write(6,613)df1,df2,df3,df4,df5,df6,df7,df8,df9,df10,
     1 df11
613   format(5e12.3)


  call co2calc with temperature 1.
   0.580E+04  -0.167E+05  -0.414E+04  -0.552E+02   0.245E+01
  -0.788E+01  -0.138E+02  -0.906E+00   0.114E+00   0.164E-01
  -0.119E-05
  call co2calc with temperature 10.
   0.968E+04  -0.233E+05  -0.508E+04  -0.144E+03   0.512E+01
  -0.126E+02  -0.211E+02  -0.865E+00   0.180E+00   0.188E-01
  -0.928E-06
  call co2calc with temperature 20.
   0.147E+05  -0.325E+05  -0.616E+04  -0.384E+03   0.895E+01
  -0.198E+02  -0.317E+02  -0.810E+00   0.289E+00   0.214E-01
  -0.746E-06
  call co2calc with temperature 30.
   0.204E+05  -0.434E+05  -0.722E+04  -0.939E+03   0.135E+02
  -0.287E+02  -0.451E+02  -0.749E+00   0.447E+00   0.236E-01
  -0.629E-06