Simulated One Point Mode

Simulated One Point Mode

This mode is currently only supported for Mars.

It requires to provide: 1: a Temperature T 2: @ a season ls 3: @ a LTST hour

 OUT = krc(lat = 0., lon = 0., T = 170., ls = 152.4, hour = 3.12)

The output is the equivalent thermal inertia in J m-2 K-1 s-1/2, in this case OUT = XXX.

The One-Point Mode outputs a TI (SI) for an array of input temperatures with the given parameters

It can output various values that indicate why a TI was not derived for failure cases:

-600: The Temperature provided is lower than the lowest model output. Possible Frost Presence

-500: No value was derived or returned. KRC may have failed to run all together

-400: Negative TI value. This is a non-physical return as a function of the fitting routine

-300: TI value larger than 2200 SI. 2200 SI is the top of the available lookup table and extrapolation is not permitted

-200: Returned TI value doesn't meet proximity criterion (TI_Guess_PCT)

-100: Bad TI_Guess / TI_Guess_PCT combination. This is most likely a user error

SOMETHING ABOUT TI_Guess_PCT and TI_Guess

HERE enter examples of a bunch of cases

OUT = krc(T=190,ls=0.,hour=4.,lat=12.,lon=12.,TI_Guess=300.,TI_Guess_PCT=0.05)            #Validity not met => -200.
OUT = krc(T=190,ls=0.,hour=4.,lat=12.,lon=12.,TI_Guess=300.,TI_Guess_PCT=0.2)             #Validity met => 334.
OUT = krc(T=190,ls=0.,hour=4.,lat=12.,lon=12.,TI_Guess=300.)                              #No validity criterion specified, good guess => 334.
OUT = krc(T=190,ls=0.,hour=4.,lat=12.,lon=12.,TI_Guess=900.)                              #No validity criterion specified, poor guess => 641. You shoot yourself in the foot
OUT = krc(T=190,ls=0.,hour=4.,lat=12.,lon=12.)                                            #No guess, slow, => 334.  
OUT = krc(T=190,ls=0.,hour=4.,lat=12.,lon=12.,TI_Guess_PCT=0.2)                           #Pathological Case; KRC runs, sets TI_Guess to 0., and returns a flag => -100.
OUT = krc(T=130,ls=0.,hour=4.,lat=12.,lon=12.)                                            #Below Frost point => -600. (The solution was a negative TI) 
OUT = krc(T=130,ls=0.,hour=4.,lat=12.,lon=12.,TI_Guess=200.)                              #Below Frost point, horrible guess but -600 
OUT = krc(T=130,ls=0.,hour=4.,lat=12.,lon=12.,TI_Guess=200.,TI_Guess_PCT=0.5)             #Below Frost point, horrible guess -200                                                                    
OUT = krc(T=500,ls=0.,hour=4.,lat=12.,lon=12.)                                            #Much too warm => flags high TI => -300.
OUT = krc(T=500,ls=0.,hour=4.,lat=12.,lon=12.,TI_Guess=200.)                              #Much too warm, poor guess => flags high TI => -300.
OUT = krc(T=500,ls=0.,hour=4.,lat=12.,lon=12.,TI_Guess=200.,TI_Guess_PCT=0.1)             #Much too warm, poor guess => flags poor guess => -200.
OUT = krc(T=190,ls=000.00,hour=4.,lat=12.,lon=12.)                                        #No guess, slow, => 334.
OUT = krc(T=190,ls=359.99,hour=4.,lat=12.,lon=12.)                                        #No guess, slow, => 334.
OUT = krc(T=190,ls=359.99,hour=4.,lat=12.,lon=12.,DELLS=0.)                                 #sets DELLS=0. to avoid this issue but error message: Please do not set both DELLS and DELJUL XXX CORRECTION MADE: Not an issue anymore since Christopher extended lookup table to Ls<0 and Ls>360.