Download krc
From 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