/[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.1.2.1 - (show annotations) (download)
Tue Feb 5 20:34:35 2002 UTC (22 years, 2 months ago) by heimbach
Branch: ecco-branch
CVS Tags: ecco_c50_e32, ecco_c50_e33, ecco_c50_e30, ecco_c50_e31, ecco_c51_e34d, ecco_c51_e34e, ecco_c51_e34f, ecco_c51_e34g, ecco_c51_e34a, ecco_c51_e34b, ecco_c51_e34c, icebear5, icebear4, icebear3, icebear2, ecco_c50_e29, ecco_c50_e28, ecco_c44_e19, ecco_c44_e18, ecco_c44_e17, ecco_c44_e16, ecco_c50_e33a, ecco_c51_e34, ecco_ice2, ecco_ice1, ecco_c44_e25, ecco_c44_e22, ecco_c44_e23, ecco_c44_e20, ecco_c44_e21, ecco_c44_e26, ecco_c44_e27, ecco_c44_e24, ecco-branch-mod1, ecco-branch-mod2, ecco-branch-mod3, ecco-branch-mod4, ecco-branch-mod5
Branch point for: c24_e25_ice, icebear
Changes since 1.1: +72 -0 lines
o Updating adjoint/makefile to ECCO code
o Adding optim and lsopt for line search optimization.
o Adding verif. experiments for ECCO
Code will be tagged ecco-branch-mod1.

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

  ViewVC Help
Powered by ViewVC 1.1.22