evalN1 krc()

Function and Definition

Evaluate N1: This function takes several input to determine the best choice for N1 (the number of layers to be run in the model) to create seasonal skin depth stability. It is activated when the parameter thick is used.

Input Arguments

Start davinci as normal.

>davinci
dv>

Type the function into davinci without any input parameters to bring up help menu for this function.

dv> evalN1_krc()
 Calculates the minimum number of layers required given the number of years the model runs
 Return an integer
 $1  = RLAY, Layer thickness ratio
 $2  = FLAY, First layer thickness (in skin depths)
 $3  = INERTIA, Surface material thermal inertia
 $4  = SPEC_HEAT, Surface material specific heat
 $5  = DENSITY, surface material density
 $6  = DELJUL, delta julian date to increment by (Default = 1.9083
 $7  = N5, number of "seasons" to run (Default = 1080)
 $8  = JDISK, start "season" to output (Default = 721)
 $9  = MAXN1, Maximum number of layers KRC can accept
 $10 = PERIOD, Length of solar day in days (of 86400 seconds)
 $11 = INERTIA2, bottom material thermal inertia
 $12 = SpHeat2, bottom material specific heat
 $13 = DENS2, bottom material density
 $14 = Top Material Thickness in m
s.piqueux 5/5/14

These are all the input parameters used for the function. You must include all 14 parameters, even the parameters that say they use a default.

Output

The purpose of the function is to output the number of layers needed to create model stability if a thickness of the upper material is declared.

This function is called automatically when thick is used, however, you can input all values in the function yourself and see

Examples

Start davinci as normal.

>davinci
dv>

Type the function into davinci with your desired input parameters. You can not include the parameter names, Davinci will only accept a list of values. But hey must be in the correct order.

dv> evalN1_krc(1.2,.18,50,595,1000,1.9083,1440,1080,50,1.0275,300,674,2500,.007)
struct, 8 elements
    req: struct, 4 elements
        FLAY: 0.1134460494   
        RLAY: 1.200000048   
        IC: 5.000000000   
        N1: 33.00000000   
    DSD1: 0.01412622258   
    DSD2: 0.02992914990   
    SSD1: 0.3657733798   
    SSD2: 0.7749620080   
    Depth: 3.091261148   
    top: 1x33x1 array of float, bsq format [132 bytes]
    thick: 1x33x1 array of float, bsq format [132 bytes]

It is important to note, that this process is not necessary, and you will save time by simply allowing davinci to perform evalN1_krc() on its own. Especially if you do not know all of the input parameters for the function.