Running the model

1.) Set control parameters, i.e. edit dc.steer in the project directory at present the following control parameters are available

Important Notice!
The first line in dc.steer is ignored, but must be present.           

The variables must be given in the prescribed order.
Important Notice!

line1 will ignored
line2 integer, model time step in seconds
line3 real, length of integration time in days
line4 two integers separated by a blank
number 1: time interval in number of time steps between two output files
number 2: 0 if output shall be a snapshot of the state vector at the output time step 1 if output shall be the mean over time interval between two outputs
line5 scale of orography separated by a comma from the name of the the file containing the orography
scale = 0.0 no orography, scale can be used to switch on orography
slowly
line6

line17 read 3 numbers L1,F1,L12 as comma separated list containing information about the vertical model half levels in the first floor of the atmosphere. We have eta(kk-1) = F1*eta(kk).
Setting F1 < 0 equidistant levels will be produced
For more details see compost.f in dclib
L1 : number of levels
F1 : factor
L12 : number of levels to be used for a smooth transition
: the next floor of the atmosphere
line18 similar to line17 but for second floor of the atmosphere
L2 : number of levels
F2 : factor
L23 : number of levels to be used for a smooth transition
: the next floor of the atmosphere
line19 similar to line17 but for third floor of the atmosphere
L3 : number of levels
F3 : factor
L34 : number of levels to be used for a smooth transition
: the next floor of the atmosphere
line19 similar to line17 but for upper most floor of the atmosphere
F4 : factor
The number of levels L4 is given by L4 = Levels-(L1+L2+L3)
If L1+L2+L3 > Levels equidistant levels will be produced.

2.) If not available, produce the file dc.bini including an initial state vector, i.e. compile inidc.f by calling cfi, and execute inidc.exe. Reply to the requests. You may answer all requests by returns.

The resulting initial state vector is stored in dc.bini.
The initial state is a thermally balanced purely zonal flow
corresponding to the equilibrium temperature distribution Te,
where all asymmetries about the equator are omitted, i.e.
inidc.f assumes always feq=0, Tsum=0 and Twin=0.
The original intension of inidc.f was to provied initial conditions for adiabatic life cycle experiments. In order let the
baroclinic waves come early into play, wave excitations should be specified according to the “Standardwellenanregung".

3.) If the initial state vector shall be kept, copy dc.bini to a file with a different name, since the model will overwrite dc.bini after a completed model run with the state vectors needed for a restart.

4.) Run the model by calling kcmc.exe. After completion of the run two output files will be produced:

a.) dc.bini containing the state vectors for a model restart in double precision.
 Important notice!
If the restart conditions shall be kept, dc.bini should be copied to a file with a different name, since dc.bini will
be overwritten after completion of the next run.

b.) dc.new containing the model output (state vector in single precision, as specified in dc.steer)
 Important notice!
If model output shall be kept, dc.new should be copied to a file with a different name, since dc.new will
be overwritten after completion of the next run.

runl is an example of a shell script running kmcm several times starting from a given initial state using the same steering file.

 Seitenanfang