Changeset 3690
- Timestamp:
- May 30, 2020, 4:14:30 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/StratAer/nucleation_tstep_mod.F90
r3689 r3690 529 529 jnuc_n= 2.1361182605986115E-1 + 3.3827029855551838*tln -3.2423555796175563E-2*tln**2 + & 530 530 & 7.0120069477221989E-5*tln**3 +8.0286874752695141/x_n + & 531 & -2.6939840579762231E-1*LOG(satratln) +1.6079879299099518*tln*LOG(satratln) +&532 & -1.9667486968141933E-2*tln**2*LOG(satratln) +&531 & (-2.6939840579762231E-1)*LOG(satratln) +1.6079879299099518*tln*LOG(satratln) + & 532 & (-1.9667486968141933E-2)*tln**2*LOG(satratln) + & 533 533 & 5.5244755979770844E-5*tln**3*LOG(satratln) + (7.8884704837892468*LOG(satratln))/x_n + & 534 534 & 4.6374659198909596*LOG(satratln)**2 - 8.2002809894792153E-2*tln*LOG(satratln)**2 + & 535 535 & 8.5077424451172196E-4*tln**2*LOG(satratln)**2 + & 536 & -2.6518510168987462E-6*tln**3*LOG(satratln)**2 +&536 & (-2.6518510168987462E-6)*tln**3*LOG(satratln)**2 + & 537 537 & (-1.4625482500575278*LOG(satratln)**2)/x_n - 5.2413002989192037E-1*LOG(satratln)**3 + & 538 538 & 5.2755117653715865E-3*tln*LOG(satratln)**3 + & 539 & -2.9491061332113830E-6*tln**2*LOG(satratln)**3 +&540 & -2.4815454194486752E-8*tln**3*LOG(satratln)**3 +&539 & (-2.9491061332113830E-6)*tln**2*LOG(satratln)**3 + & 540 & (-2.4815454194486752E-8)*tln**3*LOG(satratln)**3 + & 541 541 & (-5.2663760117394626E-2*LOG(satratln)**3)/x_n + & 542 542 & 1.6496664658266762*LOG(rhoaln) + & 543 & -8.0809397859218401E-1*tln*LOG(rhoaln) +&543 & (-8.0809397859218401E-1)*tln*LOG(rhoaln) + & 544 544 & 8.9302927091946642E-3*tln**2*LOG(rhoaln) + & 545 & -1.9583649496497497E-5*tln**3*LOG(rhoaln) +&545 & (-1.9583649496497497E-5)*tln**3*LOG(rhoaln) + & 546 546 & (-8.9505572676891685*LOG(rhoaln))/x_n + & 547 & -3.0025283601622881E+1*LOG(satratln)*LOG(rhoaln) +&547 & (-3.0025283601622881E+1)*LOG(satratln)*LOG(rhoaln) + & 548 548 & 3.0783365644763633E-1*tln*LOG(satratln)*LOG(rhoaln) + & 549 & -7.4521756337984706E-4*tln**2*LOG(satratln)*LOG(rhoaln) +&550 & -5.7651433870681853E-7*tln**3*LOG(satratln)*LOG(rhoaln) +&549 & (-7.4521756337984706E-4)*tln**2*LOG(satratln)*LOG(rhoaln) + & 550 & (-5.7651433870681853E-7)*tln**3*LOG(satratln)*LOG(rhoaln) + & 551 551 & (1.2872868529673207*LOG(satratln)*LOG(rhoaln))/x_n + & 552 & -6.1739867501526535E-1*LOG(satratln)**2*LOG(rhoaln) +&552 & (-6.1739867501526535E-1)*LOG(satratln)**2*LOG(rhoaln) + & 553 553 & 7.2347385705333975E-3*tln*LOG(satratln)**2*LOG(rhoaln) + & 554 & -3.0640494530822439E-5*tln**2*LOG(satratln)**2*LOG(rhoaln) +&554 & (-3.0640494530822439E-5)*tln**2*LOG(satratln)**2*LOG(rhoaln) + & 555 555 & 6.5944609194346214E-8*tln**3*LOG(satratln)**2*LOG(rhoaln) + & 556 556 & (-2.8681650332461055E-2*LOG(satratln)**2*LOG(rhoaln))/x_n + & 557 557 & 6.5213802375160306*LOG(rhoaln)**2 + & 558 & -4.7907162004793016E-2*tln*LOG(rhoaln)**2 +&559 & -1.0727890114215117E-4*tln**2*LOG(rhoaln)**2 +&558 & (-4.7907162004793016E-2)*tln*LOG(rhoaln)**2 + & 559 & (-1.0727890114215117E-4)*tln**2*LOG(rhoaln)**2 + & 560 560 & 5.6401818280534507E-7*tln**3*LOG(rhoaln)**2 + & 561 561 & (5.4113070888923009E-1*LOG(rhoaln)**2)/x_n + & 562 562 & 5.2062808476476330E-1*LOG(satratln)*LOG(rhoaln)**2 + & 563 & -6.0696882500824584E-3*tln*LOG(satratln)*LOG(rhoaln)**2 +&563 & (-6.0696882500824584E-3)*tln*LOG(satratln)*LOG(rhoaln)**2 + & 564 564 & 2.3851383302608477E-5*tln**2*LOG(satratln)*LOG(rhoaln)**2 + & 565 & -1.5243837103067096E-8*tln**3*LOG(satratln)*LOG(rhoaln)**2 +&565 & (-1.5243837103067096E-8)*tln**3*LOG(satratln)*LOG(rhoaln)**2 + & 566 566 & (-5.6543192378015687E-2*LOG(satratln)*LOG(rhoaln)**2)/x_n + & 567 & -1.1630806410696815E-1*LOG(rhoaln)**3 +&567 & (-1.1630806410696815E-1)*LOG(rhoaln)**3 + & 568 568 & 1.3806404273119610E-3*tln*LOG(rhoaln)**3 + & 569 & -2.0199865087650833E-6*tln**2*LOG(rhoaln)**3 +&570 & -3.0200284885763192E-9*tln**3*LOG(rhoaln)**3 +&569 & (-2.0199865087650833E-6)*tln**2*LOG(rhoaln)**3 + & 570 & (-3.0200284885763192E-9)*tln**3*LOG(rhoaln)**3 + & 571 571 & (-6.9425267104126316E-3*LOG(rhoaln)**3)/x_n 572 572 jnuc_n=EXP(jnuc_n) 573 573 574 ntot_n =-3.5863435141979573E-3 - 1.0098670235841110E-1*tln + 8.9741268319259721E-4*tln**2 - 1.4855098605195757E-6*tln**3 &575 & - 1.2080330016937095E-1/x_n + 1.1902674923928015E-3*LOG(satratln) - 1.9211358507172177E-2*tln*LOG(satratln) + &576 & 2.4648094311204255E-4*tln**2*LOG(satratln) - 7.5641448594711666E-7*tln**3*LOG(satratln) + &577 & (-2.0668639384228818E-02*LOG(satratln))/x_n - 3.7593072011595188E-2*LOG(satratln)**2 + &578 & 9.0993182774415718E-4 *tln*LOG(satratln)**2 + &579 & -9.5698412164297149E-6*tln**2*LOG(satratln)**2 + 3.7163166416110421E-8*tln**3*LOG(satratln)**2 + &580 & (1.1026579525210847E-2*LOG(satratln)**2)/x_n + 1.1530844115561925E-2 *LOG(satratln)**3 + &581 & - 1.8083253906466668E-4 *tln*LOG(satratln)**3 + 8.0213604053330654E-7*tln**2*LOG(satratln)**3 +&582 & -8.5797885383051337E-10*tln**3*LOG(satratln)**3 + (1.0243693899717402E-3*LOG(satratln)**3)/x_n + &583 & -1.7248695296299649E-2*LOG(rhoaln) + 1.1294004162437157E-2*tln*LOG(rhoaln) + &584 & -1.2283640163189278E-4*tln**2*LOG(rhoaln) + 2.7391732258259009E-7*tln**3*LOG(rhoaln) + &585 & (6.8505583974029602E-2*LOG(rhoaln))/x_n +2.9750968179523635E-1*LOG(satratln)*LOG(rhoaln) + &586 & -3.6681154503992296E-3*tln*LOG(satratln)*LOG(rhoaln) + 1.0636473034653114E-5*tln**2*LOG(satratln)*LOG(rhoaln)+ &587 & 5.8687098466515866E-9*tln**3*LOG(satratln)*LOG(rhoaln) + (-5.2028866094191509E-3*LOG(satratln)*LOG(rhoaln))/x_n+ &588 & 7.6971988880587231E-4*LOG(satratln)**2*LOG(rhoaln) - 2.4605575820433763E-5*tln*LOG(satratln)**2*LOG(rhoaln) + &589 & 2.3818484400893008E-7*tln**2*LOG(satratln)**2*LOG(rhoaln) + &590 & -8.8474102392445200E-10*tln**3*LOG(satratln)**2*LOG(rhoaln) + &591 & (-1.6640566678168968E-4*LOG(satratln)**2*LOG(rhoaln))/x_n - 7.7390093776705471E-2*LOG(rhoaln)**2 + &592 & 5.8220163188828482E-4*tln*LOG(rhoaln)**2 + 1.2291679321523287E-6*tln**2*LOG(rhoaln)**2 + &593 & -7.4690997508075749E-9*tln**3*LOG(rhoaln)**2 + (-5.6357941220497648E-3*LOG(rhoaln)**2)/x_n + &594 & -4.7170109625089768E-3*LOG(satratln)*LOG(rhoaln)**2 + 6.9828868534370193E-5*tln*LOG(satratln)*LOG(rhoaln)**2 + &595 & -3.1738912157036403E-7*tln**2*LOG(satratln)*LOG(rhoaln)**2 + &596 & 2.3975538706787416E-10*tln**3*LOG(satratln)*LOG(rhoaln)**2 + &597 & (4.2304213386288567E-4*LOG(satratln)*LOG(rhoaln)**2)/x_n + 1.3696520973423231E-3*LOG(rhoaln)**3 + &598 & -1.6863387574788199E-5*tln*LOG(rhoaln)**3 + 2.7959499278844516E-8*tln**2*LOG(rhoaln)**3 + &574 ntot_n =-3.5863435141979573E-3 - 1.0098670235841110E-1*tln + 8.9741268319259721E-4*tln**2 - 1.4855098605195757E-6*tln**3 & 575 & - 1.2080330016937095E-1/x_n + 1.1902674923928015E-3*LOG(satratln) - 1.9211358507172177E-2*tln*LOG(satratln) + & 576 & 2.4648094311204255E-4*tln**2*LOG(satratln) - 7.5641448594711666E-7*tln**3*LOG(satratln) + & 577 & (-2.0668639384228818E-02*LOG(satratln))/x_n - 3.7593072011595188E-2*LOG(satratln)**2 + & 578 & 9.0993182774415718E-4 *tln*LOG(satratln)**2 + & 579 & (-9.5698412164297149E-6)*tln**2*LOG(satratln)**2 + 3.7163166416110421E-8*tln**3*LOG(satratln)**2 + & 580 & (1.1026579525210847E-2*LOG(satratln)**2)/x_n + 1.1530844115561925E-2 *LOG(satratln)**3 + & 581 & (-1.8083253906466668E-4)*tln*LOG(satratln)**3 + 8.0213604053330654E-7*tln**2*LOG(satratln)**3 + & 582 & (-8.5797885383051337E-10)*tln**3*LOG(satratln)**3 + (1.0243693899717402E-3*LOG(satratln)**3)/x_n + & 583 & (-1.7248695296299649E-2)*LOG(rhoaln) + 1.1294004162437157E-2*tln*LOG(rhoaln) + & 584 & (-1.2283640163189278E-4)*tln**2*LOG(rhoaln) + 2.7391732258259009E-7*tln**3*LOG(rhoaln) + & 585 & (6.8505583974029602E-2*LOG(rhoaln))/x_n +2.9750968179523635E-1*LOG(satratln)*LOG(rhoaln) + & 586 & (-3.6681154503992296E-3)*tln*LOG(satratln)*LOG(rhoaln) + 1.0636473034653114E-5*tln**2*LOG(satratln)*LOG(rhoaln)+ & 587 & 5.8687098466515866E-9*tln**3*LOG(satratln)*LOG(rhoaln) + (-5.2028866094191509E-3*LOG(satratln)*LOG(rhoaln))/x_n+ & 588 & 7.6971988880587231E-4*LOG(satratln)**2*LOG(rhoaln) - 2.4605575820433763E-5*tln*LOG(satratln)**2*LOG(rhoaln) + & 589 & 2.3818484400893008E-7*tln**2*LOG(satratln)**2*LOG(rhoaln) + & 590 & (-8.8474102392445200E-10)*tln**3*LOG(satratln)**2*LOG(rhoaln) + & 591 & (-1.6640566678168968E-4*LOG(satratln)**2*LOG(rhoaln))/x_n - 7.7390093776705471E-2*LOG(rhoaln)**2 + & 592 & 5.8220163188828482E-4*tln*LOG(rhoaln)**2 + 1.2291679321523287E-6*tln**2*LOG(rhoaln)**2 + & 593 & (-7.4690997508075749E-9)*tln**3*LOG(rhoaln)**2 + (-5.6357941220497648E-3*LOG(rhoaln)**2)/x_n + & 594 & (-4.7170109625089768E-3)*LOG(satratln)*LOG(rhoaln)**2 + 6.9828868534370193E-5*tln*LOG(satratln)*LOG(rhoaln)**2 + & 595 & (-3.1738912157036403E-7)*tln**2*LOG(satratln)*LOG(rhoaln)**2 + & 596 & 2.3975538706787416E-10*tln**3*LOG(satratln)*LOG(rhoaln)**2 + & 597 & (4.2304213386288567E-4*LOG(satratln)*LOG(rhoaln)**2)/x_n + 1.3696520973423231E-3*LOG(rhoaln)**3 + & 598 & (-1.6863387574788199E-5)*tln*LOG(rhoaln)**3 + 2.7959499278844516E-8*tln**2*LOG(rhoaln)**3 + & 599 599 & 3.9423927013227455E-11*tln**3*LOG(rhoaln)**3 + (8.6136359966337272E-5*LOG(rhoaln)**3)/x_n 600 600 ntot_n=EXP(ntot_n) … … 673 673 ELSE 674 674 jnuc_i1 = 3.0108954259038608E+01+tli*6.1176722090512577E+01+(tli**2)*8.7240333618891663E-01+(tli**3)* & 675 & -4.6191788649375719E-03+(tli**(-1))*8.3537059107024481E-01 +&675 & (-4.6191788649375719E-03)+(tli**(-1))*8.3537059107024481E-01 + & 676 676 & (1.5028549216690628E+01+tli*(-1.9310989753720623E-01)+(tli**2)*8.0155514634860480E-04+(tli**3)* & 677 & -1.0832730707799128E-06+(tli**(-1))*1.7577660457989019)*(LOG(satratli)**(-2)) +&677 & (-1.0832730707799128E-06)+(tli**(-1))*1.7577660457989019)*(LOG(satratli)**(-2)) + & 678 678 & (-2.0487870170216488E-01 + tli * 1.3263949252910405E-03 + (tli**2) *(-8.4195688402450274E-06) + & 679 679 & (tli**3)*1.6154895940993287E-08 + (tli**(-1))*3.8734212545203874E+01) * (LOG(satratli)**(-2)*LOG(rhoali)) + & … … 690 690 & 3.8059489037584171E-10 + (tli**(-1)) * 1.5000492788553410E+02) * (LOG(satratli)**(-1) * LOG(rhoali)**3) + & 691 691 & (1.2250990965305315 + tli * 3.0495946490079444E+01 + (tli**2) * 2.1051563135187106E+01 + (tli**3) * & 692 & -8.2200682916580878E-02 + (tli**(-1)) * 2.9965871386685029E-02) * (LOG(rhoali)**(-2)) +&692 & (-8.2200682916580878E-02) + (tli**(-1)) * 2.9965871386685029E-02) * (LOG(rhoali)**(-2)) + & 693 693 & (4.8281605955680433 + tli * 1.7346551710836445E+02 + (tli**2) *(-1.0113602140796010E+01) + (tli**3) * & 694 694 & 3.7482518458685089E-02 + (tli**(-1)) *(-1.4449998158558205E-01)) * (LOG(rhoali)**(-1)) + & … … 702 702 & 1.3397152304977591E-02 + (tli**(-1)) *(-2.4155187776460030E-02)) * (LOG(satratli)* LOG(rhoali)**(-2)) + & 703 703 & (3.7592282990713963 + tli *(-1.5257988769009816E+02) + (tli**2) * 2.6113805420558802 + (tli**3) * & 704 & -9.0380721653694363E-03 + (tli**(-1)) *(-1.3974197138171082E-01)) * (LOG(satratli)* LOG(rhoali)**(-1)) +&704 & (-9.0380721653694363E-03) + (tli**(-1)) *(-1.3974197138171082E-01)) * (LOG(satratli)* LOG(rhoali)**(-1)) + & 705 705 & (1.8293600730573988E+01 + tli * 1.8344728606002992E+01 + (tli**2) *(-4.0063363221106751E-01) + (tli**3) & 706 706 & * 1.4842749371258522E-03 + (tli**(-1)) * 1.1848846003282287) * (LOG(satratli)) + & 707 707 & (-1.7634531623032314E+02 + tli * 4.9011762441271278 + (tli**2) *(-1.3195821562746339E-02) + (tli**3) * & 708 & -2.8668619526430859E-05 + (tli**(-1)) *(-2.9823396976393551E-01)) * (LOG(satratli)* LOG(rhoali)) +&708 & (-2.8668619526430859E-05) + (tli**(-1)) *(-2.9823396976393551E-01)) * (LOG(satratli)* LOG(rhoali)) + & 709 709 & (-3.2944043694275727E+01 + tli * 1.2517571921051887E-01 + (tli**2) * 8.3239769771186714E-05 + (tli**3) * & 710 710 & 2.8191859341519507E-07 + (tli**(-1)) *(-2.7352880736682319E+01)) * (LOG(satratli)* LOG(rhoali)**2) + & … … 714 714 & 3.7899330796456790E-06 + (tli**(-1)) * 7.1106427501756542E-01) * (LOG(satratli)**2 * LOG(rhoali)**(-1)) + & 715 715 & (-2.8901906781697811E+01 + tli *(-1.5356398793054860) + (tli**2) * 1.9267271774384788E-02 + (tli**3) * & 716 & -5.3886270475516162E-05 + (tli**(-1)) * 5.0490415975693426E-01) * (LOG(satratli)**2) +&716 & (-5.3886270475516162E-05) + (tli**(-1)) * 5.0490415975693426E-01) * (LOG(satratli)**2) + & 717 717 & (3.3365683645733924E+01 + tli *(-3.6114561564894537E-01) + (tli**2) * 9.2977354471929262E-04 + (tli**3) * & 718 718 & 1.9549769069511355E-07 + (tli**(-1)) *(-8.8865930095112855)) * (LOG(satratli)**2 * LOG(rhoali)) + & 719 719 & (2.4592563042806375 + tli *(-8.3227071743101084E-03) + (tli**2) * 8.2563338043447783E-06 + (tli**3) * & 720 & -8.4374976698593496E-09 + (tli**(-1)) *(-2.0938173949893473E+02)) * (LOG(satratli)**2 * LOG(rhoali)**2) +&720 & (-8.4374976698593496E-09) + (tli**(-1)) *(-2.0938173949893473E+02)) * (LOG(satratli)**2 * LOG(rhoali)**2) + & 721 721 & (4.4099823444352317E+01 + tli * 2.5915665826835252 + (tli**2) *(-1.6449091819482634E-02) + (tli**3) * & 722 722 & 2.6797249816144721E-05 + (tli**(-1)) * 5.5045672663909995E-01)* satratli … … 731 731 & (tli**(-1)) *(-1.9213956820114927E+03) + (tli**(3)) *(-9.3839114856129453E-06)) * LOG(satratli)**(-2) + & 732 732 & (3.9652478944137344E+00 + tli * 1.2469375098256536E-02 + (tli**2) *(-9.9837754694045633E-05) + (tli**(-1)) * & 733 & -5.1919499210175138E+02 + (tli**(3)) * 1.6489001324583862E-07) * LOG(satratli)**(-2) * LOG(rhoali) +&733 & (-5.1919499210175138E+02) + (tli**(3)) * 1.6489001324583862E-07) * LOG(satratli)**(-2) * LOG(rhoali) + & 734 734 & (2.4975714429096206E+02 + tli * 1.7107594562445172E+02 + (tli**2) *(-7.8988711365135289E-01) + (tli**(-1)) * & 735 & -2.2243599782483177E+01 + (tli**(3)) *(-1.6291523004095427E-04)) * LOG(satratli)**(-1) * LOG(rhoali)**(-2) +&735 & (-2.2243599782483177E+01) + (tli**(3)) *(-1.6291523004095427E-04)) * LOG(satratli)**(-1) * LOG(rhoali)**(-2) +& 736 736 & (-8.9270715592533611E+02 + tli * 1.2053538883338946E+02 + (tli**2) *(-1.5490408828541018E+00) + (tli**(-1)) * & 737 & -1.1243275579419826E+01 + (tli**(3)) * 4.8053105606904655E-03) * LOG(satratli)**(-1) * LOG(rhoali)**(-1) +&737 & (-1.1243275579419826E+01) + (tli**(3)) * 4.8053105606904655E-03) * LOG(satratli)**(-1) * LOG(rhoali)**(-1) + & 738 738 & (7.6426441642091631E+03 + tli *(-7.1785462414656578E+01) + (tli**2) * 2.3851864923199523E-01 + (tli**(-1)) * & 739 739 & 8.5591775688708395E+01 + (tli**(3)) *(-3.7000473243342858E-04)) * LOG(satratli)**(-1) + & … … 743 743 & 1.2625883141097162E+01 + (tli**(3)) *(-1.2728497822219101E-01)) * LOG(rhoali)**(-2) + & 744 744 & (-3.3763494256461472E+03 + tli * 3.1916579136391006E+03 + (tli**2) *(-2.7234339474441143E+01) + (tli**(-1)) * & 745 & -2.1897653262707397E+01 + (tli**(3)) * 5.1788505812259071E-02) * LOG(rhoali)**(-1) +&745 & (-2.1897653262707397E+01) + (tli**(3)) * 5.1788505812259071E-02) * LOG(rhoali)**(-1) + & 746 746 & (-1.8817843873687068E+03 + tli * 4.3038072285882070E+00 + (tli**2) * 6.6244087689671860E-03 + (tli**(-1)) * & 747 & -2.7133073605696295E+03 + (tli**(3)) *(-1.7951557394285043E-05)) * LOG(rhoali) +&747 & (-2.7133073605696295E+03) + (tli**(3)) *(-1.7951557394285043E-05)) * LOG(rhoali) + & 748 748 & (-1.7668827539244447E+02 + tli * 4.8160932330629913E-01 + (tli**2) *(-6.3133007671100293E-04) + (tli**(-1)) * & 749 749 & 2.5631774669873157E+04 + (tli**(3)) * 4.1534484127873519E-07) * LOG(rhoali)**(2) + & 750 750 & (-1.6661835889222382E+03 + tli * 1.3708900504682877E+03 + (tli**2) *(-1.7919060052198969E+01) + (tli**(-1)) * & 751 & -3.5145029804436405E+01 + (tli**(3)) * 5.1047240947371224E-02) * LOG(satratli)* LOG(rhoali)**(-2) +&751 & (-3.5145029804436405E+01) + (tli**(3)) * 5.1047240947371224E-02) * LOG(satratli)* LOG(rhoali)**(-2) + & 752 752 & (1.0843549363030939E+04 + tli *(-7.3557073636139577E+01) + (tli**2) * 1.2054625131778862E+00 + (tli**(-1)) * & 753 753 & 1.9358737917864391E+02 + (tli**(3)) *(-4.2871620775911338E-03)) * LOG(satratli)* LOG(rhoali)**(-1) + & … … 755 755 & 2.3709874548950634E+03 + (tli**(3)) * 1.4091851828620244E-04) * LOG(satratli) + & 756 756 & (5.2745372575251588E+02 + tli *(-2.6080675912627314E+00) + (tli**2) * 5.6902218056670145E-03 + (tli**(-1)) * & 757 & -3.2149319482897838E+04 + (tli**(3)) *(-5.4121996056745853E-06)) * LOG(satratli)* LOG(rhoali) +&757 & (-3.2149319482897838E+04) + (tli**(3)) *(-5.4121996056745853E-06)) * LOG(satratli)* LOG(rhoali) + & 758 758 & (-1.6401959518360403E+01 + tli * 2.4322962162439640E-01 + (tli**2) * 1.1744366627725344E-03 + (tli**(-1)) * & 759 & -8.2694427518413195E+03 + (tli**(3)) *(-5.0028379203873102E-06))* LOG(satratli)**(2) +&759 & (-8.2694427518413195E+03) + (tli**(3)) *(-5.0028379203873102E-06))* LOG(satratli)**(2) + & 760 760 & (-2.7556572017167782E+03 + tli * 4.9293344495058264E+01 + (tli**2) *(-2.6503456520676050E-01) + (tli**(-1)) * & 761 761 & 1.2130698030982167E+03 + (tli**(3)) * 4.3530610668042957E-04)* LOG(satratli)**2 * LOG(rhoali)**(-1) + & … … 767 767 & 1.8735938211539585E+02 + (tli**(3)) * 5.0974564680442852E-04)* (satratli) + & 768 768 & (3.6409071302482083E+01 + tli * 1.7919859306449623E-01 + (tli**2) *(-1.0020116255895206E-03) + (tli**(-1)) * & 769 & -8.3521083354432303E+03+ (tli**(3)) * 1.5879900546795635E-06)* satratli * LOG(rhoali)**(2))769 & (-8.3521083354432303E+03) + (tli**(3)) * 1.5879900546795635E-06)* satratli * LOG(rhoali)**(2)) 770 770 771 771 rc_i = (-3.6318550637865524E-08 + tli * 2.1740704135789128E-09 + (tli**2) * & 772 & -8.5521429066506161E-12 + (tli**3) *(-9.3538647454573390E-15)) +&772 & (-8.5521429066506161E-12) + (tli**3) *(-9.3538647454573390E-15)) + & 773 773 & (2.1366936839394922E-08 + tli *(-2.4087168827395623E-10) + (tli**2) * 8.7969869277074319E-13 + & 774 774 & (tli**3) *(-1.0294466881303291E-15))* LOG(satratli)**(-2) * LOG(rhoali)**(-1) + &
Note: See TracChangeset
for help on using the changeset viewer.
