/[MITgcm]/MITgcm/pkg/ecco/cost_generic.F
ViewVC logotype

Diff of /MITgcm/pkg/ecco/cost_generic.F

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

revision 1.9 by heimbach, Sat Feb 6 02:43:03 2010 UTC revision 1.10 by heimbach, Sat Feb 6 11:30:16 2010 UTC
# Line 8  C $Name$ Line 8  C $Name$
8       &     nnzbar, localbarfile, localbar, xx_localbar_mean_dummy,       &     nnzbar, localbarfile, localbar, xx_localbar_mean_dummy,
9       &     nnzobs, localobsfile, mult_local,       &     nnzobs, localobsfile, mult_local,
10       &     nrecloc, localstartdate, localperiod,       &     nrecloc, localstartdate, localperiod,
11       &     localmask, localweight,       &     ylocmask, localweight,
12       &     spminloc, spmaxloc, spzeroloc,       &     spminloc, spmaxloc, spzeroloc,
13       &     objf_local, num_local,       &     objf_local, num_local,
14       &     myiter, mytime, mythid )       &     myiter, mytime, mythid )
# Line 31  c     == global variables == Line 31  c     == global variables ==
31  #include "EEPARAMS.h"  #include "EEPARAMS.h"
32  #include "SIZE.h"  #include "SIZE.h"
33  #include "PARAMS.h"  #include "PARAMS.h"
34    #include "GRID.h"
35  #ifdef ALLOW_CAL  #ifdef ALLOW_CAL
36  # include "cal.h"  # include "cal.h"
37  #endif  #endif
# Line 53  c     == routine arguments == Line 54  c     == routine arguments ==
54    
55        _RL localbar   (1-olx:snx+olx,1-oly:sny+oly,nnzbar,nsx,nsy)        _RL localbar   (1-olx:snx+olx,1-oly:sny+oly,nnzbar,nsx,nsy)
56        _RL localweight(1-olx:snx+olx,1-oly:sny+oly,nnzobs,nsx,nsy)        _RL localweight(1-olx:snx+olx,1-oly:sny+oly,nnzobs,nsx,nsy)
       _RL localmask  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)  
57        _RL xx_localbar_mean_dummy        _RL xx_localbar_mean_dummy
58        _RL mult_local        _RL mult_local
59        _RL mytime        _RL mytime
# Line 64  c     == routine arguments == Line 64  c     == routine arguments ==
64        _RL objf_local(nsx,nsy)        _RL objf_local(nsx,nsy)
65        _RL num_local(nsx,nsy)        _RL num_local(nsx,nsy)
66    
67          character*(1) ylocmask
68        character*(MAX_LEN_FNAM) localbarfile        character*(MAX_LEN_FNAM) localbarfile
69        character*(MAX_LEN_FNAM) localobsfile        character*(MAX_LEN_FNAM) localobsfile
70    
# Line 90  c     == local variables == Line 91  c     == local variables ==
91        _RL localcost        _RL localcost
92        _RL junk        _RL junk
93    
94          _RL localmask  (1-olx:snx+olx,1-oly:sny+oly,nr,nsx,nsy)
95        _RL localobs   (1-olx:snx+olx,1-oly:sny+oly,nnzobs,nsx,nsy)        _RL localobs   (1-olx:snx+olx,1-oly:sny+oly,nnzobs,nsx,nsy)
96        _RL cmask (1-olx:snx+olx,1-oly:sny+oly,nnzobs)        _RL cmask (1-olx:snx+olx,1-oly:sny+oly,nnzobs)
97    
# Line 144  c--   Initialise local variables. Line 146  c--   Initialise local variables.
146          enddo          enddo
147        enddo        enddo
148    
149    c--   Assign mask
150          if ( ylocmask .EQ. 'C' .OR. ylocmask .EQ. 'c' ) then
151             localmask(:,:,:,:,:) = maskC(:,:,:,:,:)
152          elseif ( ylocmask .EQ. 'S' .OR. ylocmask .EQ. 's' ) then
153             localmask(:,:,:,:,:) = maskS(:,:,:,:,:)
154          elseif ( ylocmask .EQ. 'W' .OR. ylocmask .EQ. 'w' ) then
155             localmask(:,:,:,:,:) = maskW(:,:,:,:,:)
156          else
157             STOP 'cost_generic: wrong ylocmask'
158          endif
159    
160  c--   First, read tiled data.  c--   First, read tiled data.
161        doglobalread = .false.        doglobalread = .false.
162        ladinit      = .false.        ladinit      = .false.

Legend:
Removed from v.1.9  
changed lines
  Added in v.1.10

  ViewVC Help
Powered by ViewVC 1.1.22