Download krc

Previous Versions

Previous Versions of KRC may be downloaded, but are no longer supported. It is strongly recommended that you get the latest version of KRC from the Main Page.

Version 01222010

Download the latest version of krc (01222010)

Checking hoe one-point mode switch occurs:
TCARD(1 called only by KRC in single place
TCARD(2 called by KRC(2 places) and TSEAS (1 place)

When change card 10 encountered, tcard renames FINPUT and sets IR=4
Card 10 should only be used after the initial full parameter block.
then krc, when it gets IR=4, will switch input files (close old, open new) read
the first line of new file and write that to IOSP, and read the second line of
input file (ignored) and write a column header line to ISOP.

Debug direct fortran writes. All have comment !dbw
fort.41 tday 191 J,DIFFI(J),DTIMI, BLAY(J),SCONVG(J)
fort.42 tlats L101 n1,j3,j4,j5,n1pib
fort.43 tday 307 jj,jjj,j5,kn,fac7,(dtj(j),j=2,kn)
fort.44 tlats L287 j3,j4,j5,ncase,efrost,ave_a,taud,pres
& ,DTMJ(J3),DTMJ(J3P1),TMIN(2),TMAX(2)

Version 01092010

run KofT compare after resetting upper I to 7 decimal places

fix tday.f to recompute fac7
make krc

modify master.inp

run krc
DTIME 2.2677 3.4385

IDL .rnew lookrc
11
52

Version 12242009

Decide to elimiate conflicting use or PERIOD and PI by changing porbcm.inc and all that reference it. Then will
be able to reference it in tseas and tyearp.

PERIOD >> OPERIOD and PI >> PICON
Need to modify: porbcm.inc tseas.f tyearp.f
ephemr.f porb1.f porb.f porbin.f porbmn.f porbqq.f
Remove probdi and porbls as obsolete
porbin.f is not used??

Version 03112009

Test1: same results for constant K
load CONUP0:CONLO3 to be constant at same value as L10=0 case
dtimes
1.3018 0.0080 5.1952 0.0370
1.3028 0.0080 5.2032 0.0290
1.2588 0.0080 5.1922 0.0360
1.2598 0.0080 5.1982 0.0290
Output in work1/krc/mars/kot3.t52
IDL program qk to test differences.
Number of convergence days differs by 1 for ??

Do the following for top and lower layers
1) decide on input paramters that influence k of T
2) run either kotab.f or spspread.pro to generate k of T table
3) fit cubic function for scaled x=0.015*(T-220)
   koftfit.pro
4) transfer 4 coefficents into krc input file

Version 09092006

Design for allowing seasonally variable albedo and TauD Without changing sizes of common, there are few places available for control
values. Choose to use ID22 indicate variable soil surface (not frost) albedo and
23 to indicate variable soil Taud. 0 will mean they are not variable. Files of
text table of value versus season can be read at the start of a run. These will
apply to ALL latitudes.

Requires logical unit to read file, closes before return

Requires additional file names, read via fcard.

Requires a way to know how many season points are used.

Implimentation is complex. but can turn these variable features on and off:
To start variable albedo, input card
8 22 0 ’AlbedoFileName’ / Variable albedo text file name
Can revert to constant albedo by hokey technique of giving it bad name. E.g.,
8 22 0 ’badName’ / turn variable albedo off
Variable Tau done the same way, with 22 being replaceed with 23

Found code error in tcard.f REAL*4 LSUBS ! function should have been ALSUBS

Version 07052004

Desire to output the heatflow into the bottom of a frost layer
This would require output of DAYCOM for every latitude and season.
Decide to define a new common to contain simply this information.
define Array (# latitudes, # seasons, # models)
Then bin5 output only the defined part.

In TDAY, isolate the heat conduction term at the surface.
Requires a counter for model!
In TDISK, write new section for KODE=54
Actually: move all dimension-defining parameters from other commons into KRCCOM,
and create new hatcom.inc to hold new items

Version 01212003

Notes following tutorial for Robin Fergeson, Josh Bandfield and Tim Titus.
One-point model temperatures are influenced by elevation if the opacity is >0
Note, negative opacities can cause very long (infinite?) run times.
If N5=0, then the values entered for DAU ans SEDC are used for a single season.
If KPREF=2 and N4>8, then will do global condensed frost integration, in which case you better have input a
reasonable set of montonically increasing latitudes.
The stored times-of-day are such that the last one is ”midnight”. E.g., Hour 1 is 1/24 of a sol after midnight. (in
tday.f, see use of JJH)

Frost amount units are kg/m2

Minor code changes to insure non-negative TAUD and DRSET, and to avoid opening output binary file if it would
not be written to.

Version 07112004

There is a new version of KRC which can output the
daily average upward heatflow into the surface. This is intended for study of
the ice caps. The file output is an additional style of .bin5, selected by
K4OUT=54

In the process, I somewhat modernized the code by moving all the
dimension-defining parameters from other commons into KRCCOM and removing the
redundant definitions. I also added a new common, hatcom.inc, to transfer the new
items between routines; I am reluctant to change any sizes in the existing
commons because the C code used to read them into IDL can tolerate no mismatch
in dimensions. This common would be the place to add other additional items.

As far as I know, the new version is backward-compatible with the circa 2002
version that was distributed to Robin.

I would like to get some mileage on this new code before distributing it.

I checked, and the frost temperature does change if KPREF and LVFT are set
properly.

Beta Version

Download a beta version