KRC for Moon

From krc
(Difference between revisions)
Jump to: navigation, search
Line 8: Line 8:
 
== Command Line Example ==
 
== Command Line Example ==
  
−
== KRC Fortran Input File Example ==
+
== Earth/Moon Eclipsing KRC Fortran Input File Example ==
 +
 
 +
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
  
 
== Online Tool ==
 
== Online Tool ==
  
 
== Common Problems ==
 
== Common Problems ==

Revision as of 23:11, 18 July 2018

Contents

General advice about using KRC on The Moon

When using KRC on The Moon the best practice is to:

Command Line Example

Earth/Moon Eclipsing KRC Fortran Input File Example

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

Online Tool

Common Problems

Retrieved from "https://krc.mars.asu.edu/index.php?title=KRC_for_Moon&oldid=1762"
Personal tools