C $Header: /u/gcmpack/MITgcm/pkg/ecco/ecco_readparms.F,v 1.25 2006/07/07 22:53:36 heimbach Exp $
#include "COST_CPPOPTIONS.h"
subroutine ECCO_READPARMS( mythid )
c ==================================================================
c SUBROUTINE ecco_readparms
c ==================================================================
c
c o This routine initialises the package cost.
c started: Ralf Giering 18-Jan-2001
c
c ==================================================================
c SUBROUTINE ecco_readparms
c ==================================================================
implicit none
c == global variables ==
#include "EEPARAMS.h"
#include "SIZE.h"
#include "GRID.h"
#include "PARAMS.h"
#include "ecco_cost.h"
c == routine arguments ==
integer mythid
c == local variables ==
C msgBuf - Informational/error meesage buffer
C iUnit - Work variable for IO unit number
CHARACTER*(MAX_LEN_MBUF) msgBuf
INTEGER iUnit, num_file, num_var
c == end of interface ==
c-- Read the namelist input.
namelist //ecco_cost_nml
& data_errfile,
& tbarfile, sbarfile, psbarfile,
& ubarfile, vbarfile, wbarfile,
& tauxbarfile, tauybarfile,
& hfluxbarfile, sfluxbarfile,
& using_topex,
& topexstartdate1, topexstartdate2,
& topexperiod, topexfile, topexmeanfile,
& topexintercept, topexslope,
& using_ers,
& ersstartdate1, ersstartdate2,
& ersperiod, ersfile,
& ersintercept, ersslope,
& using_gfo,
& gfostartdate1, gfostartdate2,
& gfoperiod, gfofile,
& gfointercept, gfoslope,
& scatstartdate1, scatstartdate2,
& scatperiod, scatxdatfile, scatydatfile,
& scatx_errfile, scaty_errfile,
& ssh_errfile,
& tmistartdate1, tmistartdate2, tmidatfile,
& sststartdate1, sststartdate2, sstdatfile,
& sssstartdate1, sssstartdate2, sssdatfile,
& tdatfile, sdatfile,
& ctdtfile, ctdsfile,
& ctdtclimfile, ctdsclimfile,
& xbtfile,
& argotstartdate1, argotstartdate2, argotfile,
& argosstartdate1, argosstartdate2, argosfile,
& udriftfile, vdriftfile,
& udrifterrfile, vdrifterrfile,
& curmtrufile,curmtrvfile,
& salterrfile, temperrfile, velerrfile,
& salt0errfile, temp0errfile, vel0errfile,
& hflux_errfile, sflux_errfile,
& tauu_errfile, tauv_errfile,
& atemp_errfile, aqh_errfile,
& precip_errfile, swflux_errfile, swdown_errfile,
& uwind_errfile, vwind_errfile,
& geoid_errfile, geoid_covariancefile,
& mult_hflux, mult_sflux, mult_tauu, mult_tauv,
& mult_hfluxmm, mult_sfluxmm,
& mult_atemp, mult_aqh,
& mult_precip, mult_swflux, mult_swdown,
& mult_uwind, mult_vwind,
& mult_hmean, mult_h,
& mult_tp, mult_ers, mult_gfo,
& mult_temp, mult_salt,
& mult_tmi, mult_sst, mult_sss,
& mult_ctdt, mult_ctds,
& mult_ctdtclim, mult_ctdsclim,
& mult_xbt, mult_drift,
& mult_argot, mult_argos,
& mult_scatx, mult_scaty,
& mult_sdrift, mult_tdrift, mult_wdrift,
& mult_temp0, mult_salt0,
& mult_obcsn, mult_obcss,
& mult_obcsw, mult_obcse, mult_obcsvol,
& mult_curmtr, mult_ageos,
& mult_kapgm, mult_diffkr,
& mult_smooth_ic, mult_smooth_bc,
& whflux0, wsflux0, wtau0,
& watemp0, waqh0,
& wprecip0, wswflux0, wswdown0, wwind0,
& wdiffkr0, wkapgm0, wedtau0,
& wmean_hflux, wmean_sflux, wmean_tau,
& wmean_atemp, wmean_aqh, wmean_precip,
& wmean_swflux, wmean_swdown, wmean_wind,
& cost_iprec, cost_yftype
_BEGIN_MASTER( myThid )
c-- Set default values.
data_errfile = ' '
tbarfile = 'tbar'
sbarfile = 'sbar'
psbarfile = 'psbar'
ubarfile = 'ubar'
vbarfile = 'vbar'
wbarfile = 'wbar'
tauxbarfile = 'tauxbar'
tauybarfile = 'tauybar'
hfluxbarfile = 'hfluxbar'
sfluxbarfile = 'sfluxbar'
using_topex = .false.
topexstartdate1 = 0
topexstartdate2 = 0
topexperiod = 0. _d 0
topexintercept = 0. _d 0
topexslope = 0. _d 0
topexfile = ' '
topexmeanfile = ' '
using_ers = .false.
ersstartdate1 = 0
ersstartdate2 = 0
ersperiod = 0. _d 0
ersintercept = 0. _d 0
ersslope = 0. _d 0
ersfile = ' '
using_gfo = .false.
gfostartdate1 = 0
gfostartdate2 = 0
gfoperiod = 0. _d 0
gfointercept = 0. _d 0
gfoslope = 0. _d 0
gfofile = ' '
scatstartdate1 = 0
scatstartdate2 = 0
scatperiod = 0. _d 0
scatxdatfile = ' '
scatydatfile = ' '
ssh_errfile = ' '
scatx_errfile = ' '
scaty_errfile = ' '
tmistartdate1 = 0
tmistartdate2 = 0
tmidatfile = ' '
sststartdate1 = 0
sststartdate2 = 0
sstdatfile = ' '
sssstartdate1 = 0
sssstartdate2 = 0
sssdatfile = ' '
tdatfile = ' '
sdatfile = ' '
ctdtfile = ' '
ctdsfile = ' '
ctdtclimfile = ' '
ctdsclimfile = ' '
curmtrufile = ' '
curmtrvfile = ' '
xbtfile = ' '
argotstartdate1 = 0
argotstartdate2 = 0
argotfile = ' '
argosstartdate1 = 0
argosstartdate2 = 0
argosfile = ' '
udriftfile = ' '
vdriftfile = ' '
udrifterrfile = ' '
vdrifterrfile = ' '
salterrfile = ' '
temperrfile = ' '
velerrfile = ' '
salt0errfile = ' '
temp0errfile = ' '
vel0errfile = ' '
hflux_errfile = ' '
sflux_errfile = ' '
tauu_errfile = ' '
tauv_errfile = ' '
atemp_errfile = ' '
aqh_errfile = ' '
precip_errfile = ' '
swflux_errfile = ' '
swdown_errfile = ' '
uwind_errfile = ' '
vwind_errfile = ' '
geoid_errfile = ' '
geoid_covariancefile = ' '
mult_hflux = 0. _d 0
mult_sflux = 0. _d 0
mult_hfluxmm = 0. _d 0
mult_sfluxmm = 0. _d 0
mult_tauu = 0. _d 0
mult_tauv = 0. _d 0
mult_atemp = 0. _d 0
mult_aqh = 0. _d 0
mult_precip = 0. _d 0
mult_swflux = 0. _d 0
mult_swdown = 0. _d 0
mult_uwind = 0. _d 0
mult_vwind = 0. _d 0
mult_hmean = 0. _d 0
mult_h = 1. _d 0
mult_tp = 0. _d 0
mult_ers = 0. _d 0
mult_gfo = 0. _d 0
mult_temp = 0. _d 0
mult_salt = 0. _d 0
mult_temp0 = 0. _d 0
mult_salt0 = 0. _d 0
mult_tmi = 0. _d 0
mult_sst = 0. _d 0
mult_sss = 0. _d 0
mult_ctdt = 0. _d 0
mult_ctds = 0. _d 0
mult_ctdtclim = 0. _d 0
mult_ctdsclim = 0. _d 0
mult_xbt = 0. _d 0
mult_argot = 0. _d 0
mult_argos = 0. _d 0
mult_drift = 0. _d 0
mult_tdrift = 0. _d 0
mult_sdrift = 0. _d 0
mult_wdrift = 0. _d 0
mult_scatx = 0. _d 0
mult_scaty = 0. _d 0
mult_obcsn = 0. _d 0
mult_obcss = 0. _d 0
mult_obcsw = 0. _d 0
mult_obcse = 0. _d 0
mult_obcsvol = 0. _d 0
mult_ageos = 0. _d 0
mult_curmtr = 0. _d 0
mult_kapgm = 0. _d 0
mult_diffkr = 0. _d 0
mult_smooth_ic = 0. _d 0
mult_smooth_bc = 0. _d 0
watemp0 = 0. _d 0
waqh0 = 0. _d 0
wprecip0 = 0. _d 0
wswflux0 = 0. _d 0
wswdown0 = 0. _d 0
wwind0 = 0. _d 0
wdiffkr0 = 1. _d 0
wkapgm0 = 1. _d 0
wedtau0 = 1. _d 0
whflux0 = 0. _d 0
wsflux0 = 0. _d 0
wtau0 = 0. _d 0
wmean_atemp = 0. _d 0
wmean_aqh = 0. _d 0
wmean_precip = 0. _d 0
wmean_swflux = 0. _d 0
wmean_swdown = 0. _d 0
wmean_wind = 0. _d 0
wmean_hflux = 0. _d 0
wmean_sflux = 0. _d 0
wmean_tau = 0. _d 0
c-- Next, read the cost data file.
WRITE(msgBuf,'(A)') 'ECCO_READPARMS: opening data.ecco'
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT , 1)
CALL OPEN_COPY_DATA_FILE(
I 'data.ecco', 'ECCO_READPARMS',
O iUnit,
I myThid )
READ( iUnit, nml = ecco_cost_nml )
WRITE(msgBuf,'(A)') 'ECCO_READPARMS: finished reading data.ecco'
CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
& SQUEEZE_RIGHT , 1)
CLOSE( iUnit )
#ifdef ALLOW_CAL
c-- Get the complete dates of the ...
c-- ... SST data.
call CAL_FULLDATE( sststartdate1, sststartdate2,
& sststartdate, mythid )
c-- ... TMI data.
call CAL_FULLDATE( tmistartdate1, tmistartdate2,
& tmistartdate, mythid )
c-- ... SSS data.
call CAL_FULLDATE( sssstartdate1, sssstartdate2,
& sssstartdate, mythid )
c-- ... T/P data.
call CAL_FULLDATE( topexstartdate1, topexstartdate2,
& topexstartdate, mythid )
c-- ... ERS data.
call CAL_FULLDATE( ersstartdate1, ersstartdate2,
& ersstartdate, mythid )
c-- ... GFO data.
call CAL_FULLDATE( gfostartdate1, gfostartdate2,
& gfostartdate, mythid )
c-- ... SCAT data.
call CAL_FULLDATE( scatstartdate1, scatstartdate2,
& scatxstartdate, mythid )
call CAL_FULLDATE( scatstartdate1, scatstartdate2,
& scatystartdate, mythid )
#endif
_END_MASTER( mythid )
_BARRIER
end