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 ) |
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 |
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 |
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 |
|
|
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 |
|
|
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. |