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

Diff of /MITgcm/optim/optim_initmod.F

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

revision 1.1 by heimbach, Tue Feb 5 20:34:35 2002 UTC revision 1.1.2.1 by heimbach, Tue Feb 5 20:34:35 2002 UTC
# Line 0  Line 1 
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
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.1.2.1
  ViewVC Help
Powered by ViewVC 1.1.22