KRC for Mars

From krc
(Difference between revisions)
Jump to: navigation, search
(→‎Command Line Example)
(→‎Running in the Command Line)
 
(26 intermediate revisions by 2 users not shown)
Line 1: Line 1:
−
= General advice about using KRC on Mars =
+
==Notes==
  
−
When using KRC on Mars the best practice is to:
+
Basalt is the default material for Mars (Mat1 = "basalt"), and T_user = 220 (temperature at which the inertia is defined).
−
::*Stay away from poles
+
−
::*Be aware of what season it is
+
−
::*Be aware of the opacity settings
+
  
 +
By default, PTOTAL = 545 Pa.
  
−
 
+
=Running in the Command Line=
−
= Command Line Example =
+
  
 
KRC is run within Davinci.
 
KRC is run within Davinci.
Line 16: Line 13:
 
For Mars, it can be run with as little input as a single latitude, e.g. 25°N:
 
For Mars, it can be run with as little input as a single latitude, e.g. 25°N:
  
−
  dv> OUT = krc(lat = 25.)
+
  OUT = krc(lat=25.)
  
 
Model output is saved in structure 'OUT'. The decimal is required for parameters to be initialized as a floating point number. All other model inputs are retrieved from lookup tables (e.g., longitude, surface albedo, thermal inertia etc.).
 
Model output is saved in structure 'OUT'. The decimal is required for parameters to be initialized as a floating point number. All other model inputs are retrieved from lookup tables (e.g., longitude, surface albedo, thermal inertia etc.).
Line 24: Line 21:
 
  tsurf: 96x1x360 array of double, bsq format [276,480 bytes]
 
  tsurf: 96x1x360 array of double, bsq format [276,480 bytes]
  
−
Additional fields can be prescribed within the parentheses when calling krc(). E.g., the same latitude but for a longitude of 120° and surface albedo of 0.3:
+
Additional fields can be prescribed within the parentheses when calling krc. E.g., the same latitude but for a longitude of 120° and surface albedo of 0.3:
  
−
  dv> OUT = krc(lat = 25., lon = 120., ALBEDO = 0.3)
+
  OUT = krc(lat=25.,lon=120.,ALBEDO=0.3)
  
 
If a particular season provided as solar longitude (in units of degrees) is desired, the annual dimension (e.g., 360) is removed. E.g., for Ls = 90° (northern summer solstice):
 
If a particular season provided as solar longitude (in units of degrees) is desired, the annual dimension (e.g., 360) is removed. E.g., for Ls = 90° (northern summer solstice):
  
−
  dv> OUT = krc(lat = 25., ls = 90.)
+
  OUT = krc(lat=25.,ls=90.)
 +
 
 +
Alternatively, the model can run for a specific Gregorian Date, (currently ranging from 1990-Jan-01 to 2040-Jan-01, format: ????-Mmm-DD,
 +
with Mmm:Jan, Feb, Mar, Apr, May, Jun, Jul, Aug, Sep, Oct, Nov, Dec
 +
 
 +
OUT = krc(lat=12.,GD="2010-Jan-05")
 +
 
 +
or a specific Julian Date:
 +
 
 +
OUT = krc(lat=12.,JD=2455201)
 +
 
 +
Note: the possibility to specify the date with GD is only currently available for Mars, the Moon, Bennu, and Europa.
  
 
Now 'OUT.tsurf' has the dimensions of:
 
Now 'OUT.tsurf' has the dimensions of:
Line 36: Line 44:
 
  tsurf: 96x1x1 array of double, bsq format [768 bytes]
 
  tsurf: 96x1x1 array of double, bsq format [768 bytes]
  
−
One can plot the diurnal temperature series with:
+
One can plot the diurnal temperature series against local true solar time (LTST) with:
  
−
  dv> plot(OUT.tsurf , xaxis = OUT.time , "25N, Ls=90" , w = 2 , color = 1)
+
  plot(OUT.tsurf,xaxis=OUT.time,"25N,Ls=90",w=4 ,color=1)
−
  dv> labelxy("LTST" , "Temperature (K)")
+
  labelxy("LTST","Temperature (K)")
  
−
Which should look like:
+
[[Image:mars_temp_1.png|border|800px]]
−
[[Image:mars_temp_1.jpg|left|500px]]
+
  
−
= Table of Input Parameters =
+
Alternatively, one can prescribe local time and output will be provided for a full Mars year, e.g., for a local time of 3 AM:
 +
 
 +
OUT = krc(lat=25.,hour=3.)
 +
 
 +
'OUT.tsurf' will now be sampled roughly once per degree of solar longitude, e.g.,:
 +
 
 +
tsurf: 1x1x360 array of double, bsq format [2,880 bytes]
 +
 
 +
And can be plotted against time (with the addition of 3PM local time as 'OUT_2') with:
 +
 
 +
plot(OUT.tsurf[12,1,],w=4,color=4,"Dawn",OUT.tsurf[60,1,],w=4,color=1,"Noon")
 +
labelxy("Solar Longitude","Temperature (K)")
 +
 
 +
[[Image:mars_temp_2.png|border|800px]]
 +
 
 +
== Table of Input Parameters ==
  
 
Other common fields that can be prescribed are included in the table below (*NOTE fields are case sensitive*):
 
Other common fields that can be prescribed are included in the table below (*NOTE fields are case sensitive*):
 
[table of parameters (include example ranges?)]
 
[table of parameters (include example ranges?)]
  
−
= KRC Fortran Input File Example =
+
{| class="wikitable"
−
+
|+Input Parameters
−
  0 0 / KOLD: season to start with;  KEEP: continue saving data in same disk file
+
|-
−
+
|'''Parameter'''
−
Version 222 default values. 19 latitudes with mean Mars elevations
+
|'''KRC Syntax'''
−
+
|'''Range'''
−
    [[ALBEDO|ALBEDO]]    EMISS  INERTIA    COND2    DENS2    PERIOD SPEC_HEAT  DENSITY
+
|'''Units'''
−
+
|-
−
        .25      1.00    200.0      2.77    928.0    1.0275      647.    1600.
+
|Latitude
−
+
|lat
−
      CABR      AMW  [ABRPHA    PTOTAL    FANON      TATM    TDEEP  SpHeat2
+
| -90–90
−
+
|-
−
      0.11      43.5    -0.00    546.0      .055      200.    180.0    1711.
+
|Longitude (°E)
−
+
|lon
−
      TAUD    DUSTA    TAURAT    TWILI      ARC2    [ARC3    SLOPE    SLOAZI
+
|0–360
−
+
|-
−
        0.3      .90      0.5      0.0      0.5    -0.00      0.0      90.
+
|Albedo
−
+
|ALBEDO
−
    TFROST    CFROST    AFROST    FEMIS      AF1      AF2    FROEXT    [FD32
+
|0–1
−
+
|-
−
      146.0  589944.      .65      0.95      0.54    0.0009      50.      0.0
+
|Thermal Inertia
−
+
|INERTIA
−
      RLAY      FLAY    CONVF    DEPTH    DRSET      DDT      GGT    DTMAX
+
|20–2000
−
+
|-
−
    1.2000    .1800    2.0000      0.0      0.0    .0020      0.1      0.1
+
|Elevation
−
+
|ELEV
−
      DJUL    DELJUL  SOLARDEC      DAU    LsubS    SOLCON      GRAV    AtmCp
+
|-
−
+
|Local True Solar Time
−
  -1222.69 17.174822      00.0    1.465        .0    1368.    3.727    735.9
+
|hour
−
+
|0–24
−
    ConUp0    ConUp1    ConUp2    ConUp3    ConLo0    ConLo1    ConLo2    ConLo3
+
|-
−
+
|Solar Longitude
−
  0.038640 -0.002145  0.002347 -0.000750  2.766722 -1.298966  0.629224 -0.527291
+
|ls
−
+
|0–360
−
    SphUp0    SphUp1    SphUp2    SphUp3    SphLo0    SphLo1    SphLo2    SphLo3
+
|-
−
+
|Local True Solar Time
−
  646.6275  246.6678  -49.8216    7.9520  1710.648  721.8740  57.44873  24.37532
+
|hour
−
+
|}
−
        N1        N2        N3        N4        N5      N24        IB        IC
+
 
−
+
=Deriving Thermal Inertia and the Simulated One Point Mode=
−
        20      384        15        19      120        48        0        9
+
−
+
−
      NRSET      NMHA      NRUN    JDISK    IDOWN    FlxP14    FlxP15    KPREF
+
−
+
−
        3        24        0        81        0        45        65        1
+
−
+
−
      K4OUT    JBARE    Notif    IDISK2                                    end
+
−
+
−
        52        0        20        0                                      0
+
−
+
−
    LP1    LP2    LP3    LP4    LP5    LP6 LPGLOB  LVFA  LVFT  LkofT
+
−
+
−
      F      T      F      F      F      F      F      F      F      T
+
−
+
−
  LPORB  LKEY    LSC  spare  LOCAL  Prt76 LPTAVE  Prt78  Prt79  L_ONE
+
−
+
−
      T      F      F      F      F      T      F      T      F      F
+
−
+
−
Latitudes: in 10F7.2  _____7 _____7 _____7 _____7 _____7 _____7 _____7
+
−
+
−
  -87.50 -80.00 -70.00 -60.00 -50.00 -40.00 -30.00 -20.00 -10.00  0.00
+
−
+
−
  10.00  20.00  30.00  40.00  50.00  60.00  70.00  80.00  87.50  -0.00
+
−
+
−
  _____7 _____7 _____7 Elevations: in 10F7.2 ____7 _____7 _____7 _____7
+
−
+
−
    3.51  2.01  1.39  1.22  0.38  0.48  1.17  1.67  1.26  0.17
+
−
+
−
  -0.94  -1.28  -1.99  -2.51  -3.52  -4.08  -4.51  -4.38  -2.57  -0.00
+
−
+
−
  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 0 0 'master222.t52' / Disk file name for Run 1
+
−
+
−
0/
+
−
+
−
3 10 1 'LkofT' / Temperature-dependant conductivity
+
−
+
−
0/
+
−
+
−
0/  ======================= end of run
+
−
+
−
+
−
+
−
If LkofT set to T, then
+
−
+
−
Upper material: weakly cemented particulates:
+
−
+
−
Grain: k: BasicRocks_Zoth88    2:4
+
−
+
−
            Cp: Chlorite_Bert07_Fe=0.89  6:4.89
+
−
+
−
Cement: k: Limestone Zoth88 2:1
+
−
+
−
        Cp:  Sphene, which has relatively strong T dependence  5:0
+
−
+
−
Cement fraction 1.e-8
+
−
+
−
Yields c_0 of 0.050087, This adjusted to 0.038640 to agree with I=200 at 220 K
+
−
+
−
+
−
+
−
lower material: H2O Ice
+
−
+
−
k: koftop: 48
+
−
+
−
fit to A+B/T fit      2.766722 -1.298966  0.629224 -0.527291 <k H2O:ice3sources
+
−
+
−
Cp: koftop: @ 49,491,5,33 yields:
+
−
+
−
    3.95779 >  1710.648  721.8740  57.44873  24.37532 <SpH H2O:Ice_3sources
+
  
 +
In order to derive thermal inertia values for a single temperature measurement one can use the simulated one-point mode in KRC. *''available only for Mars at present''
  
−
= Common Problems =
+
Please refer to the [[Simulated One Point Mode]] page for more information.

Latest revision as of 16:16, 9 April 2020

Contents

[edit] Notes

Basalt is the default material for Mars (Mat1 = "basalt"), and T_user = 220 (temperature at which the inertia is defined).

By default, PTOTAL = 545 Pa.

[edit] Running in the Command Line

KRC is run within Davinci. 

> davinci

For Mars, it can be run with as little input as a single latitude, e.g. 25°N: 

OUT = krc(lat=25.)

Model output is saved in structure 'OUT'. The decimal is required for parameters to be initialized as a floating point number. All other model inputs are retrieved from lookup tables (e.g., longitude, surface albedo, thermal inertia etc.).

By default, the output is stored into multidimensional arrays sampled at 96 values per sol and 360 values per Mars year. Hence the structure element surface temperature ('OUT.tsurf') appears as: 

tsurf: 96x1x360 array of double, bsq format [276,480 bytes]

Additional fields can be prescribed within the parentheses when calling krc. E.g., the same latitude but for a longitude of 120° and surface albedo of 0.3: 

OUT = krc(lat=25.,lon=120.,ALBEDO=0.3)

If a particular season provided as solar longitude (in units of degrees) is desired, the annual dimension (e.g., 360) is removed. E.g., for Ls = 90° (northern summer solstice): 

OUT = krc(lat=25.,ls=90.)

Alternatively, the model can run for a specific Gregorian Date, (currently ranging from 1990-Jan-01 to 2040-Jan-01, format: ????-Mmm-DD, with Mmm:Jan, Feb, Mar, Apr, May, Jun, Jul, Aug, Sep, Oct, Nov, Dec 

OUT = krc(lat=12.,GD="2010-Jan-05")

or a specific Julian Date: 

OUT = krc(lat=12.,JD=2455201)

Note: the possibility to specify the date with GD is only currently available for Mars, the Moon, Bennu, and Europa.

Now 'OUT.tsurf' has the dimensions of: 

tsurf: 96x1x1 array of double, bsq format [768 bytes]

One can plot the diurnal temperature series against local true solar time (LTST) with: 

plot(OUT.tsurf,xaxis=OUT.time,"25N,Ls=90",w=4 ,color=1)
labelxy("LTST","Temperature (K)")

mars temp 1.png

Alternatively, one can prescribe local time and output will be provided for a full Mars year, e.g., for a local time of 3 AM: 

OUT = krc(lat=25.,hour=3.)

'OUT.tsurf' will now be sampled roughly once per degree of solar longitude, e.g.,: 

tsurf: 1x1x360 array of double, bsq format [2,880 bytes]

And can be plotted against time (with the addition of 3PM local time as 'OUT_2') with: 

plot(OUT.tsurf[12,1,],w=4,color=4,"Dawn",OUT.tsurf[60,1,],w=4,color=1,"Noon")
labelxy("Solar Longitude","Temperature (K)")

mars temp 2.png

[edit] Table of Input Parameters

Other common fields that can be prescribed are included in the table below (*NOTE fields are case sensitive*): [table of parameters (include example ranges?)]

Input Parameters
Parameter KRC Syntax Range Units
Latitude lat -90–90
Longitude (°E) lon 0–360
Albedo ALBEDO 0–1
Thermal Inertia INERTIA 20–2000
Elevation ELEV
Local True Solar Time hour 0–24
Solar Longitude ls 0–360
Local True Solar Time hour

[edit] Deriving Thermal Inertia and the Simulated One Point Mode

In order to derive thermal inertia values for a single temperature measurement one can use the simulated one-point mode in KRC. *available only for Mars at present

Please refer to the Simulated One Point Mode page for more information.

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