/[MITgcm]/MITgcm/optim/optim_initmod.F
ViewVC logotype

Annotation of /MITgcm/optim/optim_initmod.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.4 - (hide annotations) (download)
Tue May 26 22:54:08 2015 UTC (8 years, 11 months ago) by gforget
Branch: MAIN
Changes since 1.3: +1 -1 lines
- add CTRL_OPTIONS.h that is needed to set maxcvars correctly
  when using generic controls (contributed by D. Amrhein).

1 heimbach 1.2
2 gforget 1.4 #include "CTRL_OPTIONS.h"
3 heimbach 1.2
4     subroutine optim_initmod(
5     I nn,
6     O xx
7     & )
8    
9     c ==================================================================
10     c SUBROUTINE optim_initmod
11     c ==================================================================
12     c
13     c o Optim_Initmod is called before the optimization. It sets a first guess
14     c of the parameter vector. It may also contain the initialization
15     c of the model.
16     c
17     c started: Christian Eckert eckert@mit.edu 15-Feb-2000
18     c
19     c changed: Christian Eckert eckert@mit.edu 10-Mar-2000
20     c
21     c - Added ECCO layout.
22     c
23     c changed: Patrick Heimbach heimbach@mit.edu 19-Jun-2000
24     c - finished, revised and debugged
25     c
26     c ==================================================================
27     c SUBROUTINE optim_initmod
28     c ==================================================================
29    
30     IMPLICIT NONE
31    
32     c == global variables ==
33    
34     #include "EEPARAMS.h"
35     #include "SIZE.h"
36     #include "ctrl.h"
37     #include "optim.h"
38    
39     c == routine arguments ==
40    
41     integer nn
42     _RL xx(nn)
43    
44     c == local variables ==
45    
46     character*(80) fname
47    
48     c == end of interface ==
49    
50     c-- Generate file name for the scaling of the control vector
51     c-- contributions.
52     write(fname(1:80),'(80a)') ' '
53     write(fname(1:80),'(4a,i10.10,a)')
54 dfer 1.3 & scalname(1:9),'_',yctrlid,'.',optimcycle,'.data'
55 heimbach 1.2
56     c-- Read the scaling for the control vector.
57    
58     c --> Presently not implemented, since the control vector and the
59     c gradient of the cost function can already be non-dimensionalized
60     c by *ecco_ControlPack* and *ecco_ControlUnPack*.
61     c
62     c To do this set:
63     c
64     c #define USE_NONDIMENSIONAL_CONTROL_IO
65     c
66     c in CPP_OPTIONS.h. The subroutine *ecco_ControlPack* will then
67     c non-dimensionalize its output and *ecco_ControlUnPack* will
68     c dimensionalize its input.
69    
70     return
71     end

  ViewVC Help
Powered by ViewVC 1.1.22