2 |
|
|
3 |
#include "CTRL_CPPOPTIONS.h" |
#include "CTRL_CPPOPTIONS.h" |
4 |
|
|
5 |
|
CBOP |
6 |
|
C !ROUTINE: ctrl_pack |
7 |
|
C !INTERFACE: |
8 |
|
subroutine ctrl_pack( myiter, mytime, mythid ) |
9 |
|
|
10 |
|
C !DESCRIPTION: \bv |
11 |
|
c *================================================================= |
12 |
|
c | SUBROUTINE ctrl_pack |
13 |
|
c | Pack the control vector |
14 |
|
c | * All control variable and adjoint variable fields are |
15 |
|
c | read from disk. |
16 |
|
c | * Wet points are extracted, and elements are |
17 |
|
c | normalized (optional) |
18 |
|
c | * A single control vector containing only (normalized |
19 |
|
c | wet points is written to file. |
20 |
|
c *================================================================= |
21 |
|
C \ev |
22 |
|
|
23 |
subroutine ctrl_pack( |
C !USES: |
|
I myiter, |
|
|
I mytime, |
|
|
I mythid |
|
|
& ) |
|
|
|
|
|
c ================================================================== |
|
|
c SUBROUTINE ctrl_pack |
|
|
c ================================================================== |
|
|
c |
|
|
c o Compress the control vector such that only ocean points are |
|
|
c written to file. |
|
|
c |
|
|
c started: Christian Eckert eckert@mit.edu 10-Mar=2000 |
|
|
c |
|
|
c changed: Patrick Heimbach heimbach@mit.edu 06-Jun-2000 |
|
|
c - Transferred some filename declarations |
|
|
c from here to namelist in ctrl_init |
|
|
c |
|
|
c Patrick Heimbach heimbach@mit.edu 16-Jun-2000 |
|
|
c - single file name convention with or without |
|
|
c ALLOW_ECCO_OPTIMIZATION |
|
|
c |
|
|
c |
|
|
c ================================================================== |
|
|
c SUBROUTINE ctrl_pack |
|
|
c ================================================================== |
|
|
|
|
24 |
implicit none |
implicit none |
25 |
|
|
26 |
c == global variables == |
c == global variables == |
|
|
|
27 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
28 |
#include "SIZE.h" |
#include "SIZE.h" |
29 |
#include "PARAMS.h" |
#include "PARAMS.h" |
30 |
#include "GRID.h" |
#include "GRID.h" |
|
|
|
31 |
#include "ctrl.h" |
#include "ctrl.h" |
32 |
#include "cost.h" |
#include "cost.h" |
33 |
#include "optim.h" |
#include "optim.h" |
34 |
|
|
35 |
|
C !INPUT/OUTPUT PARAMETERS: |
36 |
c == routine arguments == |
c == routine arguments == |
|
|
|
37 |
integer myiter |
integer myiter |
38 |
_RL mytime |
_RL mytime |
39 |
integer mythid |
integer mythid |
40 |
|
|
41 |
|
C !LOCAL VARIABLES: |
42 |
c == local variables == |
c == local variables == |
43 |
|
|
44 |
integer bi,bj |
integer bi,bj |
51 |
integer jtlo,jthi |
integer jtlo,jthi |
52 |
integer jmin,jmax |
integer jmin,jmax |
53 |
integer imin,imax |
integer imin,imax |
54 |
|
integer cbuffindex |
55 |
|
integer cunit |
56 |
|
integer prec |
57 |
|
|
58 |
logical doglobalread |
logical doglobalread |
59 |
logical ladinit |
logical ladinit |
|
integer cbuffindex |
|
60 |
|
|
|
integer cunit |
|
61 |
_RL cbuff( snx*nsx*npx*sny*nsy*npy ) |
_RL cbuff( snx*nsx*npx*sny*nsy*npy ) |
62 |
_RL globfld3d( snx,nsx,npx,sny,nsy,npy,nr ) |
_RL globfld3d( snx,nsx,npx,sny,nsy,npy,nr ) |
63 |
_RL globfld2d( snx,nsx,npx,sny,nsy,npy ) |
_RL globfld2d( snx,nsx,npx,sny,nsy,npy ) |
67 |
character*(128) cfile |
character*(128) cfile |
68 |
character*( 80) fname |
character*( 80) fname |
69 |
|
|
|
integer prec |
|
|
|
|
70 |
c == external == |
c == external == |
|
|
|
71 |
integer ilnblnk |
integer ilnblnk |
72 |
external ilnblnk |
external ilnblnk |
73 |
|
|
74 |
c == end of interface == |
c == end of interface == |
75 |
|
CEOP |
76 |
|
|
77 |
prec = precFloat64 |
prec = precFloat64 |
78 |
tmpvar = -9999. _d 0 |
tmpvar = -9999. _d 0 |