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

Contents of /MITgcm/optim/optim_initmod.F

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


Revision 1.4 - (show 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
2 #include "CTRL_OPTIONS.h"
3
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 & scalname(1:9),'_',yctrlid,'.',optimcycle,'.data'
55
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