/[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.1.2.1 - (hide 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 heimbach 1.1.2.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