Lines Matching refs:pt
104 PetscScalar temp, pt
108 pt = 0.1
197 PetscScalar temp, pt
206 pt = 0.1
300 part_p(ip) = (an_r(ip)/an_t)*pt
310 - (an_r(4)**i_cc)*(an_r(5)**i_h2o)*((pt/an_t)**idiff)
393 PetscScalar temp, pt
415 pt = 0.1
552 pt_cube = pt*pt*pt
553 pt_four = pt_cube*pt
554 pt_five = pt_cube*pt*pt
565 pt1 = pt**(ai_pwr1)
567 pt_val1 = (pt/an_t)**(ai_pwr1)
570 pt2 = pt**(i_pwr2)
572 pt_val2 = (pt/an_t)**(i_pwr2)
577 *((pt/an_t)**idiff)*(-idiff/an_t)
589 (an_r(5)**i_h2o)*((pt/an_t)**idiff)
592 *(an_r(4)**i_cc)*((pt/an_t)**idiff)
594 d_eq(3, 1) = -(an_r(4)**2)*(an_r(5)**3)*(pt/an_t)*(-1.0/an_t)
603 d_eq(3, 4) = d_eq(3, 4) - 2*an_r(4)*(an_r(5)**3)*(pt/an_t)
605 d_eq(3, 5) = d_eq(3, 5) - 3*(an_r(5)**2)*(an_r(4)**2)*(pt/an_t)
613 d_eq(7, 1) = pt*an_r(11)*(-1.0)/const2 &
614 - k_eq(1)*sqrt(pt)*sqrt(an_r(14) + 1d-50)*(-0.5/const3)
620 d_eq(7, 11) = d_eq(7, 11) + pt/an_t
622 - k_eq(1)*sqrt(pt)*(0.5/(sqrt((an_r(14) + 1d-50)*an_t)))
624 d_eq(8, 1) = pt*an_r(8)*(-1.0)/const2 &
625 - k_eq(2)*sqrt(pt)*sqrt(an_r(3) + 1.0d-50)*(-0.5/const3)
632 - k_eq(2)*sqrt(pt)*(0.5/(sqrt((an_r(3) + 1.0d-50)*an_t)))
633 d_eq(8, 8) = d_eq(8, 8) + pt/an_t
635 d_eq(9, 1) = pt*an_r(7)*(-1.0)/const2 &
636 - k_eq(3)*sqrt(pt)*sqrt(an_r(6))*(-0.5/const3)
642 d_eq(9, 7) = d_eq(9, 7) + pt/an_t
644 - k_eq(3)*sqrt(pt)*(0.5/(sqrt(an_r(6)*an_t)))
646 d_eq(10, 1) = pt*an_r(10)*(-1.0)/const2 &
647 - k_eq(4)*(pt)*sqrt((an_r(3) + 1.0d-50) &
654 - k_eq(4)*(pt)*sqrt(an_r(14)) &
656 d_eq(10, 10) = d_eq(10, 10) + pt/an_t
658 - k_eq(4)*(pt)*sqrt(an_r(3) + 1.0d-50) &
661 d_eq(11, 1) = pt*an_r(9)*(-1.0)/const2 &
662 - k_eq(5)*(pt)*sqrt((an_r(3) + 1.0d-50)*an_r(6)) &
670 - k_eq(5)*(pt)*sqrt(an_r(6))*(0.5/ &
673 - k_eq(5)*(pt)*sqrt(an_r(3) + 1.0d-50) &
675 d_eq(11, 9) = d_eq(11, 9) + pt/an_t
677 d_eq(12, 1) = pt*an_r(5)*(-1.0)/const2 &
678 - k_eq(6)*(pt**1.5)*sqrt(an_r(3) + 1.0d-50) &
686 - k_eq(6)*(pt**1.5)*((an_r(14) + 1.0d-50)/const3) &
689 d_eq(12, 5) = d_eq(12, 5) + pt/an_t
691 - k_eq(6)*(pt**1.5)*(sqrt(an_r(3) + 1.0d-50)/const3)
693 d_eq(13, 1) = pt*an_r(4)*(-1.0)/const2 &
694 - k_eq(7)*(pt**1.5)*sqrt(an_r(3) + 1.0d-50) &
702 - k_eq(7)*(pt**1.5)*(an_r(13)/const3) &
705 d_eq(13, 4) = d_eq(13, 4) + pt/an_t
707 - k_eq(7)*(pt**1.5)*(sqrt(an_r(3) + 1.0d-50)/const3)
709 d_eq(14, 1) = pt*an_r(15)*(-1.0)/const2 &
710 - k_eq(8)*(pt**1.5)*sqrt(an_r(3) + 1.0d-50) &
718 - k_eq(8)*(pt**1.5)*(an_r(9)/const3) &
721 - k_eq(8)*(pt**1.5)*(sqrt(an_r(3) + 1.0d-50)/const3)
722 d_eq(14, 15) = d_eq(14, 15) + pt/an_t
724 d_eq(15, 1) = pt*an_r(16)*(-1.0)/const2 &
725 - k_eq(9)*(pt**1.5)*sqrt(an_r(14) + 1.0d-50) &
733 - k_eq(9)*(pt**1.5)*(sqrt(an_r(14) + 1.0d-50)/const3)
735 - k_eq(9)*(pt**1.5)*(an_r(3)/const3) &
737 d_eq(15, 16) = d_eq(15, 16) + pt/an_t
739 d_eq(16, 1) = pt*an_r(12)*(-1.0)/const2 &
740 - k_eq(10)*(pt**1.5)*sqrt(an_r(3) + 1.0d-50) &
748 - k_eq(10)*(pt**1.5)*(an_r(6)/const3) &
752 - k_eq(10)*(pt**1.5)*(sqrt(an_r(3) + 1.0d-50)/const3)
753 d_eq(16, 12) = d_eq(16, 12) + pt/an_t
758 d_eq(17, 1) = an_r(14)*an_r(18)*an_r(18)*(pt**3)*(-3/const_four) &
759 - k_eq(15)*an_r(17)*pt*(-1/const2)
763 d_eq(17, 14) = d_eq(17, 14) + an_r(18)*an_r(18)*(pt**3)/const_cube
764 d_eq(17, 17) = d_eq(17, 17) - k_eq(15)*pt/an_t
766 & *(pt**3)/const_cube
768 d_eq(18, 1) = an_r(13)*an_r(13)*(pt**2)*(-2/const_cube) &
770 *(pt*pt*pt)*(-3/const_four)
775 - k_eq(16)*an_r(18)*an_r(18)*pt*pt*pt/const_cube
777 + 2*an_r(13)*pt*pt/const2
779 & *2*an_r(18)*pt*pt*pt/const_cube
783 d_eq(19, 1) = an_r(3)*an_r(19)*(pt**2)*(-2/const_cube) &
784 - k_eq(17)*an_r(13)*an_r(10)*pt*pt*(-2/const_cube)
789 - k_eq(17)*an_r(10)*pt*pt/const2
791 - k_eq(17)*an_r(13)*pt*pt/const2
792 d_eq(19, 3) = d_eq(19, 3) + an_r(19)*pt*pt/const2
793 d_eq(19, 19) = d_eq(19, 19) + an_r(3)*pt*pt/const2
796 d_eq(20, 1) = an_r(21)*an_r(20)*(pt**2)*(-2/const_cube) &
797 - k_eq(18)*an_r(19)*an_r(8)*pt*pt*(-2/const_cube)
802 - k_eq(18)*an_r(19)*pt*pt/const2
804 - k_eq(18)*an_r(8)*pt*pt/const2
805 d_eq(20, 20) = d_eq(20, 20) + an_r(21)*pt*pt/const2
806 d_eq(20, 21) = d_eq(20, 21) + an_r(20)*pt*pt/const2
811 d_eq(21, 1) = an_r(21)*an_r(23)*(pt**2)*(-2/const_cube) &
812 - k_eq(19)*an_r(7)*an_r(8)*pt*pt*(-2/const_cube)
817 - k_eq(19)*an_r(8)*pt*pt/const2
819 - k_eq(19)*an_r(7)*pt*pt/const2
820 d_eq(21, 21) = d_eq(21, 21) + an_r(23)*pt*pt/const2
821 d_eq(21, 23) = d_eq(21, 23) + an_r(21)*pt*pt/const2
825 d_eq(22, 1) = an_r(5)*an_r(11)*(pt**2)*(-2/const_cube) &
826 - k_eq(20)*an_r(21)*an_r(22)*pt*pt*(-2/const_cube)
831 - k_eq(20)*an_r(22)*pt*pt/const2
833 - k_eq(20)*an_r(21)*pt*pt/const2
834 d_eq(22, 11) = d_eq(22, 11) + an_r(5)*pt*pt/(const2)
835 d_eq(22, 5) = d_eq(22, 5) + an_r(11)*pt*pt/(const2)
840 d_eq(23, 1) = an_r(24)*(pt)*(-1/const2) &
841 - k_eq(21)*an_r(21)*an_r(3)*pt*pt*(-2/const_cube)
846 - k_eq(21)*an_r(21)*pt*pt/const2
848 - k_eq(21)*an_r(3)*pt*pt/const2
849 d_eq(23, 24) = d_eq(23, 24) + pt/(an_t)
853 d_eq(24, 1) = an_r(3)*an_r(25)*(pt**2)*(-2/const_cube) &
854 - k_eq(22)*an_r(24)*an_r(8)*pt*pt*(-2/const_cube)
859 - k_eq(22)*an_r(24)*pt*pt/const2
861 - k_eq(22)*an_r(8)*pt*pt/const2
862 d_eq(24, 3) = d_eq(24, 3) + an_r(25)*pt*pt/const2
863 d_eq(24, 25) = d_eq(24, 25) + an_r(3)*pt*pt/const2
868 d_eq(25, 1) = an_r(26)*(pt)*(-1/const2) &
869 - k_eq(23)*an_r(21)*an_r(10)*pt*pt*(-2/const_cube)
874 - k_eq(23)*an_r(21)*pt*pt/const2
876 - k_eq(23)*an_r(10)*pt*pt/const2
877 d_eq(25, 26) = d_eq(25, 26) + pt/(an_t)