Download krc
From krc
(Difference between revisions)
Line 3: | Line 3: | ||
[[Media:krc01222010.tar | Download]] the latest version of krc ([[Version 01222010|01222010]]) | [[Media:krc01222010.tar | Download]] the latest version of krc ([[Version 01222010|01222010]]) | ||
− | + | Checking hoe one-point mode switch occurs: | |
− | TCARD(1 called only by KRC in single place | + | TCARD(1 called only by KRC in single place |
− | TCARD(2 called by KRC(2 places) and TSEAS (1 place) | + | TCARD(2 called by KRC(2 places) and TSEAS (1 place) |
− | + | ||
− | Card 10 should only be used after the initial full parameter block. | + | When change card 10 encountered, tcard renames FINPUT and sets IR=4 |
− | then krc, when it gets IR=4, will switch input files (close old, open new) read | + | Card 10 should only be used after the initial full parameter block. |
− | the first line of new file and write that to IOSP, and read the second line of | + | then krc, when it gets IR=4, will switch input files (close old, open new) read |
− | input file (ignored) and write a column header line to ISOP. | + | 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. | |
− | fort.41 tday 191 J,DIFFI(J),DTIMI, BLAY(J),SCONVG(J) | + | |
− | fort.42 tlats L101 n1,j3,j4,j5,n1pib | + | Debug direct fortran writes. All have comment !dbw |
− | fort.43 tday 307 jj,jjj,j5,kn,fac7,(dtj(j),j=2,kn) | + | fort.41 tday 191 J,DIFFI(J),DTIMI, BLAY(J),SCONVG(J) |
− | fort.44 tlats L287 j3,j4,j5,ncase,efrost,ave_a,taud,pres | + | fort.42 tlats L101 n1,j3,j4,j5,n1pib |
− | & ,DTMJ(J3),DTMJ(J3P1),TMIN(2),TMAX(2) | + | 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== | ==Previous Versions== | ||
Line 22: | Line 24: | ||
Version [[Media:krc09092006.tar | 09092006]] | Version [[Media:krc09092006.tar | 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 | + | 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 | + | 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 | + | text table of value versus season can be read at the start of a run. These will |
− | apply to ALL latitudes. | + | apply to ALL latitudes. |
− | + | ||
− | + | Requires logical unit to read file, closes before return | |
− | + | ||
− | + | Requires additional file names, read via fcard. | |
− | + | ||
− | 8 22 0 ’AlbedoFileName’ / Variable albedo text file name | + | Requires a way to know how many season points are used. |
− | Can revert to constant albedo by hokey technique of giving it bad name. E.g., | + | |
− | 8 22 0 ’badName’ / turn variable albedo off | + | Implimentation is complex. but can turn these variable features on and off: |
− | Variable Tau done the same way, with 22 being replaceed with 23 | + | 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 [[Media:krc03112009.tar | 03112009]] | Version [[Media:krc03112009.tar | 03112009]] |
Revision as of 20:11, 24 February 2010
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