KRC for Eclipsing
From krc
(Difference between revisions)
Line 1: | Line 1: | ||
= Learn about Eclipsing = | = Learn about Eclipsing = | ||
− | [[Image:Bias_def.png | + | [[Image:Bias_def.png]] |
+ | |||
Revision as of 16:37, 14 January 2019
Contents |
Learn about Eclipsing
To learn about eclipsing: Eclipsing User Guide (PDF)
Or, visit: Eclipses on Wikipedia
Common bodies involving eclipsing
Earch/Moon
Taken from MoonA.inp
0 0 1 / KOLD: season to start with; KEEP: continue saving data in same disk file 0 0 0 0 0 0 debug codes Earth and Moon lunar eclipse 5 lats ALBEDO EMISS INERTIA COND2 DENS2 PERIOD SPEC_HEAT DENSITY .12 .95 55.0 0.070 1800.0 29.53 800. 1100. CABR AMW SatPrA PTOTAL FANON TATM TDEEP SpHeat2 0.11 43.5 27.9546 0.001 .055 200. 180.0 800. TAUD DUSTA TAURAT TWILI Heny-Gren ARC3=Safe SLOPE SLOAZI 0.3 .90 0.25 0.0 0.5 0.801 0.0 45. TFROST CFROST AFROST FEMIS AF1 AF2 FROEXT SatPrB 146.0 589944. .65 0.95 0.54 0.0009 50. 3182.48 RLAY FLAY CONVF DEPTH DRSET PhotoF GGT DTMAX 1.1500 .100 3.0 0.0 0.0 0.0 0.1 0.1 DJUL DELJUL SOLARDEC DAU LsubS SOLCON GRAV AtmCp 7305.00 35.511810 00.0 1.465 .0 1368. 3.727 735.9 ConUp0 ConUp1 ConUp2 ConUp3 ConLo0 ConLo1 ConLo2 ConLo3 0.0023471 0.0012847 5.8396e-4 8.8478e-5 0.0045689 0.0025008 0.0011367 0.0001722 SphUp0 SphUp1 SphUp2 SphUp3 SphLo0 SphLo1 SphLo2 SphLo3 790.701 578.508 208.610 66.0619 790.701 578.508 208.610 66.0619 N1 N2 N3 N4 N5 N24 IIB IC2 22 1536 15 5 20 96 0 999 NRSET NMHA NRUN JDISK IDOWN FlxP14 TUN/Flx15 KPREF 3 24 0 18 0 45 65 1 K4OUT JBARE Notif IDISK2 end -3 9999 5 0 0 LP1 LP2 LP3 LP4 LP5 LP6 LPGLOB LVFA LVFT LkofT F T F F F F F F F F LPORB LKEY LSC LZONE LOCAL Prt76 LPTAVE Prt78 Prt79 L_ONE T F F F T F F F F F Latitudes: in 10F7.2 _____7 _____7 _____7 _____7 _____7 _____7 _____7 -60. -30.0 0 30.0 60.0 _____7 _____7 _____7 Elevations: in 10F7.2 ____7 _____7 _____7 _____7 0.0 0.0 0.0 0.0 0.0 PORB:2014jun10 2018 Feb 21 21:13:11 IPLAN,TC= 103.0 0.20000 Earth:Moon 103.0000 0.2000000 0.000000 -0.4545965E-04 1.797730 0.1670245E-01 1.000004 0.4090926 0.000000 1.161374 4.712311 0.000000 0.000000 365.2584 6942.444 655.7199 0.000000 1.261200 0.3769714E-03 0.000000 0.000000 0.3046746 -0.9524564 -0.3590488E-03 0.9524565 0.3046746 0.1148536E-03 -0.000000 -0.3769713E-03 0.9999999 2 2 1536 'N2' / 8 5 0 './out/MoonA.t52' / Disk file name for Run 16 3 '00N' / Tsurf every time step 0 / 2 7 20. 'IIB' / heat flow. 0 / 15 .05157 .005 20 0.011 0.011 0 12 / Earth flux onto Moon at lunar midday 0 / 14 1 1. 6371.008 384.4e3 1737. 29.53 0.0 6000 12. 7 / Moon lunar 0 /^^^^ 4 14 1 1. 6371.008 384.4e3 1737. 29.53 0.5 6000 12. 7 / Moon lunar 0 / 14 3 1. 6371.008 384.4e3 1737. 29.53 0.5 6000 12. 7 / Moon lunar 0 / 1 3 200. 'Inertia' / 0 / 1 3 55. 'Inertia' / reset 14 3 1. 6371.008 384.4e3 1737. 29.53 0.5 6000 15. 7 / Moon lunar 15 .05157 .005 20 0.011 0.011 0 15 / Earth flux onto Moon at 15 H 0 / ^^^^ 8 14 1 1. 6371.008 384.4e3 1737. 29.53 0.5 6000 12. 7 / Moon lunar midday 15 .05157 .005 20 0.011 0.011 0 12 / Earth flux onto Moon at lunar midday 2 8 7 'IC2' / First layer of lower zone 1 8 1275. 'Dens' / upper zone 1 5 1625. 'Dens' / lower zone 3 10 1 'LKofT / use T-dep properties 0 / 0 / ======================= end of run 0 / 1 23 30. 'Slope' / I find Haynes c_p derived from his surface inertia, total conductivity (at 220K) and density to be 179.675, much too low. His polynomial for sp.heat at 220 is 791 Haynes density changes most rapidly about 0.05m Haynes K-phonon changes most rapidly about 0.05m above = .0015, below= .0029 ps Quat= 1100.00 1275.00 1450.00 1625.00 1800.00 kz Quat= 0.000740000 0.00140500 0.00207000 0.00273500 0.00340000 switch depth 1/2=exp(-zz/capH) 0.0415888 kk1 0.00234712 0.00128470 0.000583956 8.84783e-05 kk3 0.00456894 0.00250083 0.00113674 0.000172234 qt2=[0.0045689,0.0025008,0.0011367,0.0001722] ; rounded for KRC input 5.29950e-08 3.02750e-08 mean and stdev of resid after rounding
Mars/Phobos
Taken from PhoH.inp
0 0 1/ KOLD: season to start with; KEEP: continue saving data in same disk file 0 0 0 0 0 0 debug codes Version 355 Eclipse test. Mars and Phobos ALBEDO EMISS INERTIA COND2 DENS2 PERIOD SPEC_HEAT DENSITY .25 1.00 200.0 2.77 928.0 1.0275 647. 1600. CABR AMW SatPrA PTOTAL FANON TATM TDEEP SpHeat2 0.11 43.5 27.9546 546.0 .055 200. 180.0 1711. TAUD/PHT DUSTA TAURAT TWILI Hen-Gre ARC3/Safe SLOPE SLOAZI 0.3 .94 0.204 0.0 0.5 0.801 0.0 90. TFROST CFROST AFROST FEMIS AF1 AF2 FROEXT SatPrB 146.0 589944. .65 0.95 0.54 0.0009 50. 3182.48 RLAY FLAY CONVF DEPTH DRSET PhotoF GGT DTMAX 1.1500 .100 3.0 0.0 0.0 0.0 0.1 0.1 DJUL DELJUL SOLARDEC DAU LsubS SOLCON GRAV AtmCp -1222.69 17.174822 00.0 1.465 .0 1368. 3.727 735.9 ConUp0 ConUp1 ConUp2 ConUp3 ConLo0 ConLo1 ConLo2 ConLo3 0.038640 -0.002145 0.002347 -0.000750 2.766722 -1.298966 0.629224 -0.527291 SphUp0 SphUp1 SphUp2 SphUp3 SphLo0 SphLo1 SphLo2 SphLo3 646.6275 246.6678 -49.8216 7.9520 1710.648 721.8740 57.44873 24.37532 N1 N2 N3 N4 N5 N24 IIB IC2 28 1536 15 3 10 96 0 999 NRSET NMHA NRUN JDISK IDOWN FlxP14 TUN/Flx15 KPREF 3 24 0 1 0 45 65 1 K4OUT JBARE Notif IDISK2 end -1 9999 5 0 0 LP1 LP2 LP3 LP4 LP5 LP6 LPGLOB LVFA LVFT LkofT F T F F F F F F F F LPORB LKEY LSC LZONE LOCAL Prt76 LPTAVE Prt78 Prt79 L_ONE T F F F T F F F F F Latitudes: in 10F7.2 _____7 _____7 _____7 _____7 _____7 _____7 _____7 0.00 30.00 60.00 30.00 60.00 _____7 _____7 _____7 Elevations: in 10F7.2 ____7 _____7 _____7 _____7 1.22 1.17 0.17 -1.99 -4.08 2013 Jul 24 11:28:09=RUNTIME. IPLAN AND TC= 104.0 0.10000 Mars:Mars 104.0000 0.1000000 0.8644665 0.3226901E-01 -1.281586 0.9340198E-01 1.523712 0.4090926 0.000000 0.9229373 5.544402 0.000000 0.000000 686.9929 3397.977 24.62296 0.000000 -1.240317 0.000000 0.000000 0.000000 0.3244965 0.8559126 0.4026359 -0.9458869 0.2936298 0.1381285 0.000000 -0.4256703 0.9048783 8 5 0 './out/PhoH.t52' / Type 52 file to write 0 / base mars 14 3 1.52368 11.27 9376. 3396.2 0.3189 0.001 7000. 12. 17 0. / solar, odd 16 1 '00N' / Tsurf every time step 0 / 14 4 1.52368 11.27 9376. 3396.2 0.3189 0.001 7000. 12. 17 0. / solar, even 0 / 14 7 1.52368 11.27 9376. 3396.2 0.3189 0.001 7000. 12. 17 0. / solar, max 0 / 16 0 'none' / Tfine off 14 0 / no eclipse 1 12 0.1 'Ptotal' / no atm 1 6 0.3189 'Period' / sol 1 3 50 'inertia' / Phobos base 0 / 14 1 1.52368 3396.2 9376. 11.27 0.3189 0.001 7000. 12. 17 0. / Lunar 16 1 '00N' / Tsurf every time step 0/ 14 2 1.52368 3396.2 9376. 11.27 0.3189 0.001 7000. 12. 17 0. / rare 0 / 14 3 1.52368 3396.2 9376. 11.27 0.3189 0.95 7000. 12. 17 0. / rare big bias 0 / 0 / =============== end of run
Europa/Jupiter
Taken from EurH.inp
0 0 1 / KOLD: season to start with; KEEP: continue saving data in same disk file 0 0 1 0 7 0 debug codes Version 352 Europa test of eclipse and planetary load. 3 lats ALBEDO EMISS INERTIA COND2 DENS2 PERIOD SPEC_HEAT DENSITY .67 1.00 200.0 2.77 928.0 3.5511810 647. 1600. CABR AMW SatPrA PTOTAL FANON TATM TDEEP SpHeat2 0.11 43.5 27.9546 0.001 .055 200. 180.0 1711. TAUD DUSTA TAURAT TWILI ARC2/Pho ARC3=Safe SLOPE SLOAZI 0.3 .90 0.25 0.0 0.0 0.801 0.0 45. TFROST CFROST AFROST FEMIS AF1 AF2 FROEXT SatPrB 146.0 589944. .65 0.95 0.54 0.0009 50. 3182.48 RLAY FLAY CONVF DEPTH DRSET DDT GGT DTMAX 1.1500 .100 3.0 0.0 0.0 .0020 0.1 0.1 DJUL DELJUL SOLARDEC DAU LsubS SOLCON GRAV AtmCp 7305.00 35.511810 00.0 1.465 .0 1368. 3.727 735.9 ConUp0 ConUp1 ConUp2 ConUp3 ConLo0 ConLo1 ConLo2 ConLo3 0.038640 -0.002145 0.002347 -0.000750 2.766722 -1.298966 0.629224 -0.527291 SphUp0 SphUp1 SphUp2 SphUp3 SphLo0 SphLo1 SphLo2 SphLo3 646.6275 246.6678 -49.8216 7.9520 1710.648 721.8740 57.44873 24.37532 N1 N2 N3 N4 N5 N24 IIB IC2 22 1536 15 3 20 96 0 999 NRSET NMHA NRUN JDISK IDOWN FlxP14 TUN/Flx15 KPREF 3 24 0 18 0 45 65 1 K4OUT JBARE Notif IDISK2 end -3 9999 5 0 0 LP1 LP2 LP3 LP4 LP5 LP6 LPGLOB LVFA LVFT LkofT F T F F F F F F F F LPORB LKEY LSC LZONE LOCAL Prt76 LPTAVE Prt78 Prt79 L_ONE T F F F T F F F F F Latitudes: in 10F7.2 _____7 _____7 _____7 _____7 _____7 _____7 _____7 0.0 30.0 60.0 _____7 _____7 _____7 Elevations: in 10F7.2 ____7 _____7 _____7 _____7 0.0 0.0 0.0 0.0 0.0 PORB:2014jun10 2017 Mar 16 22:11:24 IPLAN,TC= 105.0 0.20000 Jupiter:Europa 105.0000 0.2000000 1.754315 0.2275961E-01 -1.496513 0.4835973E-01 5.202864 0.4090926 0.000000 1.125922 4.678847 0.000000 0.000000 4334.724 4095.201 85.22835 0.000000 2.136487 0.5414962E-01 0.000000 0.000000 -0.5359988 -0.8429813 -0.4569179E-01 0.8442187 -0.5352132 -0.2900995E-01 0.000000 -0.5412316E-01 0.9985343 2 2 6144 'N2' / 8 5 0 './out/EurH.t52' / Disk file name for Run 0 / 1 23 30. 'Slope' / 0 / 1 23 0. 'Slope' / flat 15 0.156 0. 0. 0.464 0.464 0. 12. / Jupiter heat load on Europa for 12 H 0 / 2 7 100. 'IIB' / heat flow. 0 /^^^^ 4 14 1 5.2026 71492. 0.6711D6 3121.6 3.551 0.00 7500. 12. 1 / Europa Daily 0 / 2 7 0 'IIB' / no heat flow. 0 / 14 1 5.2026 71492. 0.6711D6 3121.6 3.551 0.00 7500. 13. 0 / Europa Daily 0 / ^^^ 14 1 5.2026 66854. 0.6711D6 3121.6 3.551 0.63 7500. 12. 0 / Europa Daily 0 / --------------------------- 16 1 '00N' / Tsurf every time step 2 7 100. 'IIB' / heat flow. Need non-zero to test TFINE base options 14 3 5.2026 71492. 0.6711D6 3121.6 3.551 0.00 7500. 12. 1 / rare 0/ 2 7 100. 'IIB' / heat flow. Need non-zero to test TFINE base options 14 3 5.2026 71492. 0.6711D6 3121.6 3.551 0.00 -7500. 12. 0 / rare, base constant 0/ 0 / ======================= end of run 15 0.156 0. 0. 0.464 0.464 0. 12. / Jupiter heat load on Europa for 12 H 0 / 15 0 / planHeat off 2 7 100. 'IIB' / heat flow. 14 3 5.2026 71492. 0.6711D6 3121.6 3.551 0.01 7500. 13. 2 / Europa Rare 16 1 '00N' / output TOUT (input file stem will be added in front) 0 /^^^^7 14 3 5.2026 71492. 0.6711D6 3121.6 3.551 0.01 -7500. 12. 1 / Rare keep bot T 0 / 16 0 'off' / turn off output TOUT 14 1 5.2026 71492. 0.6711D6 3121.6 3.551 0.01 7500. 13. 1 / Europa Daily 15 0.156 0. 0. 0.464 0.464 0. 13. / Jupiter heat load on Europa for 14 H 1 23 30. 'Slope' / 0 / 14 0 / eclipse off 15 0 / PlanHeat off 2 7 0 'IIB' / no heat flow. Case should be identical to case 2 0 / DJUL 7305.00 is 2020 jan 01 DELJUL 35.511810 is 10 Europa days, or 1/122.004 of Jovian year Jupiter inclination 3.13 degree Europa inclination .466 degree, node precession 30.1 years