KRC for Eclipsing

From krc
(Difference between revisions)
Jump to: navigation, search
(Blanked the page)
 
(6 intermediate revisions by 2 users not shown)
Line 1: Line 1:
−
= Learn about Eclipsing =
 
  
−
To learn about eclipsing: [http://krc.mars.asu.edu/svn/filedetails.php?repname=KRC&path=%2Ftags%2Fkrc_3.5.6%2Fdoc%2Feclipse.pdf Eclipsing User Guide (PDF)]
 
−
 
−
= Common bodies involving eclipsing =
 
−
 
−
=== Earch/Moon ===
 
−
 
−
Taken from [http://krc.mars.asu.edu/svn/filedetails.php?repname=KRC&path=%2Ftags%2Fkrc_3.5.6%2Frun%2FMoonA.inp 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 [http://krc.mars.asu.edu/svn/filedetails.php?repname=KRC&path=%2Ftags%2Fkrc_3.5.6%2Frun%2FPhoH.inp 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 [http://krc.mars.asu.edu/svn/filedetails.php?repname=KRC&path=%2Ftags%2Fkrc_3.5.6%2Frun%2FEurH.inp 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
 

Latest revision as of 15:51, 15 January 2019

Personal tools