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