/[MITgcm]/MITgcm/pkg/ctrl/ctrl_pack.F
ViewVC logotype

Diff of /MITgcm/pkg/ctrl/ctrl_pack.F

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

revision 1.11 by heimbach, Thu Oct 30 19:09:05 2003 UTC revision 1.12 by heimbach, Thu Nov 6 22:05:08 2003 UTC
# Line 2  C Line 2  C
2  C $Header$  C $Header$
3  C $Name$  C $Name$
4    
5    #include "PACKAGES_CONFIG.h"
6  #include "CTRL_CPPOPTIONS.h"  #include "CTRL_CPPOPTIONS.h"
7    
8        subroutine ctrl_pack( first, mythid )        subroutine ctrl_pack( first, mythid )
# Line 26  c Line 27  c
27  c              G. Gebbie, added open boundary control packing,  c              G. Gebbie, added open boundary control packing,
28  c                  gebbie@mit.edu  18 -Mar- 2003  c                  gebbie@mit.edu  18 -Mar- 2003
29  c  c
30  c              heimbach@mit.edu totally restrucured 28-Oct-2003  c              heimbach@mit.edu totally restructured 28-Oct-2003
31  c  c
32  c     ==================================================================  c     ==================================================================
33  c     SUBROUTINE ctrl_pack  c     SUBROUTINE ctrl_pack
# Line 41  c     == global variables == Line 42  c     == global variables ==
42  #include "PARAMS.h"  #include "PARAMS.h"
43  #include "GRID.h"  #include "GRID.h"
44    
 #include "ecco.h"  
45  #include "ctrl.h"  #include "ctrl.h"
46  #include "cost.h"  #include "cost.h"
47    
48    #ifdef ALLOW_ECCO
49    # include "ecco_cost.h"
50    #else
51    # include "ctrl_weights.h"
52    #endif
53    
54  #ifdef ALLOW_ECCO_OPTIMIZATION  #ifdef ALLOW_ECCO_OPTIMIZATION
55  #include "optim.h"  # include "optim.h"
56  #endif  #endif
57    
58  c     == routine arguments ==  c     == routine arguments ==
# Line 131  c Line 137  c
137  c--     Only the master thread will do I/O.  c--     Only the master thread will do I/O.
138        _BEGIN_MASTER( mythid )        _BEGIN_MASTER( mythid )
139    
       print *, 'ph-pack in pack '  
140        if ( first .AND. optimcycle .EQ. 0 ) then        if ( first .AND. optimcycle .EQ. 0 ) then
141  c     >>> Initialise control vector for optimcycle=0 <<<  c     >>> Initialise control vector for optimcycle=0 <<<
       print *, 'ph-pack in ctrl '  
142            lxxadxx   = .TRUE.            lxxadxx   = .TRUE.
143            ictrlgrad = 1            ictrlgrad = 1
144            fcloc     = fmin            fcloc     = fmin
# Line 142  c     >>> Initialise control vector for Line 146  c     >>> Initialise control vector for
146       &      ctrlname(1:9),'_',yctrlid(1:10),'.opt', optimcycle       &      ctrlname(1:9),'_',yctrlid(1:10),'.opt', optimcycle
147        else        else
148  c     >>> Write gradient vector <<<  c     >>> Write gradient vector <<<
       print *, 'ph-pack in cost '  
149            lxxadxx   = .FALSE.            lxxadxx   = .FALSE.
150            ictrlgrad = 2            ictrlgrad = 2
151            fcloc     = fc            fcloc     = fc
# Line 150  c     >>> Write gradient vector <<< Line 153  c     >>> Write gradient vector <<<
153       &    costname(1:9),'_',yctrlid(1:10),'.opt', optimcycle       &    costname(1:9),'_',yctrlid(1:10),'.opt', optimcycle
154         endif         endif
155    
       print *, 'ph-pack vor open ', optimcycle, cfile  
156         call mdsfindunit( cunit, mythid )         call mdsfindunit( cunit, mythid )
157         open( cunit, file   = cfile,         open( cunit, file   = cfile,
158       &      status = 'unknown',       &      status = 'unknown',
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22