Download krc

Latest Version

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)

Previous Versions

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 03112009

• Test1: same results for constant K

loadt 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 ?? 28 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

Beta Version

Download a beta version