KRC for Eclipsing

From krc
(Difference between revisions)
Jump to: navigation, search
Line 6: Line 6:
  
 
=== Earch/Moon ===
 
=== 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
  
  

Revision as of 23:59, 18 July 2018

Contents

Learn about Eclipsing

To learn about eclipsing: Eclipsing User Guide (PDF)

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
Personal tools