1 |
|
2 |
Structure until now: |
3 |
-------------------- |
4 |
|
5 |
grdchk_main: |
6 |
|
7 |
- save unperturbed cost function |
8 |
fcref = fc |
9 |
|
10 |
- choose an index |
11 |
-> done by equidistant loop |
12 |
do icomp = nbeg, nend, nstep |
13 |
ichknum = (icomp - nbeg)/nstep + 1 |
14 |
|
15 |
- grdchk_loc: |
16 |
determine entry within given field, i.e. 1-dim array, determined by |
17 |
* ctrlvarindex |
18 |
* ncvargrd (ctrlvarindex) |
19 |
* ncvarnrmax (ctrlvarindex) |
20 |
* ncvarrecs (ctrlvarindex) |
21 |
-> determine |
22 |
* icvrec |
23 |
* jtile |
24 |
* itile |
25 |
* layer |
26 |
|
27 |
- grdchk_getxx: |
28 |
get control variable component from file |
29 |
-> only relevant for optimcycle > 0, i.e. when |
30 |
xx_... has values .NE. 0 |
31 |
|
32 |
- grdchk_setxx: |
33 |
* get control variable component (read from file) |
34 |
* add perturbation |
35 |
* write perturbed field back to file |
36 |
|
37 |
|
38 |
- want to reset perturbation component after grad. check |