** Warning **
Issuing rollback() due to DESTROY without explicit disconnect() of DBD::mysql::db handle dbname=MITgcm at /usr/local/share/lxr/lib/LXR/Common.pm line 1224.
MITgcm/MITgcm/pkg/ecco/ecco_readparms.F
File indexing completed on 2023-03-29 05:10:07 UTC
view on github raw file Latest commit cda1c18f on 2023-03-28 22:31:47 UTC
8f7d13d0c9 Jean* 0001 #include "ECCO_OPTIONS.h "
6b47d550f4 Mart* 0002 #ifdef ALLOW_CTRL
0003 # include "CTRL_OPTIONS.h "
0004 #endif
0c0a2531ee Patr* 0005
cda1c18f72 Jean* 0006 SUBROUTINE ECCO_READPARMS ( myThid )
0c0a2531ee Patr* 0007
0008
0009
0010
0011
5001c65f45 Patr* 0012
0013
0c0a2531ee Patr* 0014
0015
0016
0017
0018
44d3986245 Jean* 0019 IMPLICIT NONE
0c0a2531ee Patr* 0020
0021
0022 #include "EEPARAMS.h "
5001c65f45 Patr* 0023 #include "SIZE.h "
0024 #include "GRID.h "
0c0a2531ee Patr* 0025 #include "PARAMS.h "
5001c65f45 Patr* 0026
c509d7e04a Gael* 0027 #ifdef ALLOW_ECCO
cf705a6c8e Mart* 0028 # include "ECCO_SIZE.h "
0029 # include "ECCO.h "
0030 # include "ecco_local_params.h "
c509d7e04a Gael* 0031 #endif
0032 #ifdef ALLOW_CTRL
0033 # include "CTRL_OBCS.h "
0034 #endif
49484c0542 Gael* 0035 #ifdef ALLOW_CAL
0036 # include "cal.h "
0037 #endif
0c0a2531ee Patr* 0038
0039
9f5240b52a Jean* 0040 INTEGER myThid
0c0a2531ee Patr* 0041
44d3986245 Jean* 0042
9f5240b52a Jean* 0043 INTEGER ilnblnk
0044 EXTERNAL ilnblnk
44d3986245 Jean* 0045
0c0a2531ee Patr* 0046
f8e779c983 antn* 0047
9aaf43452b Patr* 0048
0049 CHARACTER *(MAX_LEN_MBUF ) msgBuf
9f5240b52a Jean* 0050 INTEGER k , k2 , iUnit
e7db56ba12 Gael* 0051 #ifdef ALLOW_GENCOST_CONTRIBUTION
c17b89ca05 Gael* 0052 INTEGER IL , kk , gencost_k3d , gencost_msk_k3d , ioUnit
447bdc4b79 Gael* 0053 CHARACTER *(128) tempfile
0054 LOGICAL exst
0055 _RS dummyRS (1)
e7db56ba12 Gael* 0056 #endif
0057
d39bcf1629 An T* 0058
9f5240b52a Jean* 0059 INTEGER nRetired
cda1c18f72 Jean* 0060 CHARACTER *(2) cost_yftype
0061 CHARACTER *(MAX_LEN_FNAM ) topexmeanfile
0062 LOGICAL using_cost_bp , using_cost_scat
d39bcf1629 An T* 0063
0c0a2531ee Patr* 0064
0065
5001c65f45 Patr* 0066
951926fb9b Jean* 0067 namelist /ecco_cost_nml /
5001c65f45 Patr* 0068 & data_errfile ,
9b12a028cc Patr* 0069 & tbarfile , sbarfile ,sstbarfile ,
d534937279 Matt* 0070 & psbarfile , bpbarfile ,iestaubarfile ,
5001c65f45 Patr* 0071 & ubarfile , vbarfile , wbarfile ,
0072 & tauxbarfile , tauybarfile ,
8c157ed454 Patr* 0073 & hfluxmeanbarfile , sfluxmeanbarfile ,
d39bcf1629 An T* 0074 & costTranspDataFile , mdtdatfile ,
0075 & mdtstartdate1 , mdtstartdate2 ,
0076 & mdtenddate1 , mdtenddate2 ,
9f3ef6a237 Gael* 0077 & using_cost_altim ,
f8e779c983 antn* 0078 & using_cost_bp , using_cost_scat ,
2ac072a19d Gael* 0079 & using_cost_sst , using_cost_seaice ,
f1e0eff8fe An T* 0080 & using_cost_transp ,
49484c0542 Gael* 0081 & using_mdt , using_tpj , using_ers , using_gfo ,
0082 & using_topex ,
5001c65f45 Patr* 0083 & topexstartdate1 , topexstartdate2 ,
c9ba60667e Patr* 0084 & topexperiod , topexfile , tpTimeMaskFile ,
0085 & topexintercept , topexslope , ersTimeMaskFile ,
0086 & ersstartdate1 , ersstartdate2 , gfoTimeMaskFile ,
5001c65f45 Patr* 0087 & ersperiod , ersfile ,
69531f39dd Patr* 0088 & ersintercept , ersslope ,
4fd9553bc8 Patr* 0089 & gfostartdate1 , gfostartdate2 ,
0090 & gfoperiod , gfofile ,
69531f39dd Patr* 0091 & gfointercept , gfoslope ,
951926fb9b Jean* 0092 & scatstartdate1 , scatstartdate2 ,
5001c65f45 Patr* 0093 & scatperiod , scatxdatfile , scatydatfile ,
0094 & scatx_errfile , scaty_errfile ,
8f7d13d0c9 Jean* 0095 & ssh_errfile , sshv4cost_scalefile ,
087d6770fc Gael* 0096 & sshv4cost_errfile , sshv4cost_errfactor ,
350884e02b Patr* 0097 & tp_errfile , ers_errfile , gfo_errfile ,
5001c65f45 Patr* 0098 & tmistartdate1 , tmistartdate2 , tmidatfile ,
951926fb9b Jean* 0099 & sststartdate1 , sststartdate2 , sstdatfile ,
0100 & sssstartdate1 , sssstartdate2 , sssdatfile ,
0101 & bpstartdate1 , bpstartdate2 , bpdatfile ,
d534937279 Matt* 0102 & iesstartdate1 , iesstartdate2 , iesdatfile ,
c9ba60667e Patr* 0103 & tdatfile , sdatfile , ctdtfile , ctdsfile ,
0104 & ctdtclimfile , ctdsclimfile , xbtfile ,
5001c65f45 Patr* 0105 & argotstartdate1 , argotstartdate2 , argotfile ,
0106 & argosstartdate1 , argosstartdate2 , argosfile ,
0107 & udriftfile , vdriftfile ,
0108 & udrifterrfile , vdrifterrfile ,
0109 & curmtrufile ,curmtrvfile ,
0110 & salterrfile , temperrfile , velerrfile ,
2b3c0560b2 Patr* 0111 & salt0errfile , temp0errfile , vel0errfile ,
30125f96e6 Patr* 0112 & etan0errfile , uvel0errfile , vvel0errfile ,
d534937279 Matt* 0113 & ssterrfile , ssserrfile , bperrfile ,ieserrfile ,
5001c65f45 Patr* 0114 & hflux_errfile , sflux_errfile ,
0115 & tauu_errfile , tauv_errfile ,
951926fb9b Jean* 0116 & atemp_errfile , aqh_errfile ,
02568c8633 Patr* 0117 & precip_errfile , swflux_errfile , swdown_errfile ,
fed442e87e Patr* 0118 & lwflux_errfile , lwdown_errfile , evap_errfile ,
0119 & snowprecip_errfile , apressure_errfile ,
5001c65f45 Patr* 0120 & uwind_errfile , vwind_errfile ,
0121 & geoid_errfile , geoid_covariancefile ,
bad79706a2 Gael* 0122 & bottomdrag_errfile ,edtau_errfile ,
7e2482cabc Gael* 0123 & kapgm_errfile ,kapredi_errfile ,diffkr_errfile ,
55bf8021cb Gael* 0124 & usercost_datafile , usercost_errfile ,
5001c65f45 Patr* 0125 & mult_hflux , mult_sflux , mult_tauu , mult_tauv ,
4bf6063982 Patr* 0126 & mult_hfluxmm , mult_sfluxmm ,
951926fb9b Jean* 0127 & mult_atemp , mult_aqh ,
02568c8633 Patr* 0128 & mult_precip , mult_swflux , mult_swdown ,
fed442e87e Patr* 0129 & mult_lwflux , mult_lwdown , mult_evap ,
0130 & mult_snowprecip , mult_apressure ,
d7ee8fe52e Patr* 0131 & mult_uwind , mult_vwind ,
087d6770fc Gael* 0132 & mult_hmean , mult_h , mult_sshv4cost ,
951926fb9b Jean* 0133 & mult_tp , mult_ers , mult_gfo ,
706976e7e4 Patr* 0134 & mult_temp , mult_salt ,
e284aa03cc Patr* 0135 & mult_tmi , mult_sst , mult_sss , mult_bp ,
d534937279 Matt* 0136 & mult_ctdt , mult_ctds , mult_ies ,
951926fb9b Jean* 0137 & mult_ctdtclim , mult_ctdsclim ,
5001c65f45 Patr* 0138 & mult_xbt , mult_drift ,
0139 & mult_argot , mult_argos ,
0140 & mult_scatx , mult_scaty ,
0141 & mult_sdrift , mult_tdrift , mult_wdrift ,
0142 & mult_temp0 , mult_salt0 ,
30125f96e6 Patr* 0143 & mult_etan0 , mult_uvel0 , mult_vvel0 ,
c509d7e04a Gael* 0144 #if (defined (ALLOW_CTRL ) && defined (ALLOW_OBCS ))
5001c65f45 Patr* 0145 & mult_obcsn , mult_obcss ,
0146 & mult_obcsw , mult_obcse , mult_obcsvol ,
c509d7e04a Gael* 0147 & mult_ageos ,
0148 #endif
0149 & mult_curmtr ,
7e2482cabc Gael* 0150 & mult_kapgm , mult_kapredi , mult_diffkr ,
bad79706a2 Gael* 0151 & mult_edtau , mult_bottomdrag ,
55bf8021cb Gael* 0152 & mult_usercost , mult_transp ,
f2b10b358f Patr* 0153 & mult_smooth_ic , mult_smooth_bc ,
d7ee8fe52e Patr* 0154 & whflux0 , wsflux0 , wtau0 ,
bad79706a2 Gael* 0155 & wbottomdrag0 ,
951926fb9b Jean* 0156 & watemp0 , waqh0 ,
02568c8633 Patr* 0157 & wprecip0 , wswflux0 , wswdown0 , wwind0 ,
fed442e87e Patr* 0158 & wsnowprecip0 , wlwflux0 , wlwdown0 , wevap0 ,
0159 & wapressure0 , wdiffkr0 , wkapgm0 , wedtau0 ,
7e2482cabc Gael* 0160 & wkapredi0 , wmean_hflux , wmean_sflux , wmean_tau ,
951926fb9b Jean* 0161 & wmean_atemp , wmean_aqh , wmean_precip ,
02568c8633 Patr* 0162 & wmean_swflux , wmean_swdown , wmean_wind ,
fed442e87e Patr* 0163 & wmean_lwflux , wmean_lwdown , wmean_evap ,
0164 & wmean_snowprecip , wmean_apressure ,
906a61c194 Ou W* 0165 & cost_iprec , ecco_output_sterGloH ,
0166 & ecco_keepTSeriesOutp_open ,
0167 & cost_yftype , topexmeanfile
0c0a2531ee Patr* 0168
8c157ed454 Patr* 0169 #ifdef ALLOW_GENCOST_CONTRIBUTION
0170 namelist /ecco_gencost_nml /
f8e779c983 antn* 0171 & using_gencost ,
69f8f4c84c Patr* 0172 & gencost_barfile ,
8f7d13d0c9 Jean* 0173 & gencost_datafile ,
d181f9168b Gael* 0174 & gencost_name ,
68315fab79 Gael* 0175 & gencost_scalefile ,
d3a3bee25e Patr* 0176 & gencost_errfile ,
81e05fa829 Gael* 0177 & gencost_itracer ,
a565f49279 Jean* 0178 & gencost_kLev_select ,
5cce2b5d76 Gael* 0179 & gencost_preproc ,
3bd7a4549f Gael* 0180 & gencost_preproc_c ,
0181 & gencost_preproc_i ,
0182 & gencost_preproc_r ,
92d10a5666 Gael* 0183 & gencost_posproc ,
3bd7a4549f Gael* 0184 & gencost_posproc_c ,
0185 & gencost_posproc_i ,
0186 & gencost_posproc_r ,
92d10a5666 Gael* 0187 & gencost_outputlevel ,
d3a3bee25e Patr* 0188 & gencost_mask ,
69f8f4c84c Patr* 0189 & gencost_spmin ,
0190 & gencost_spmax ,
0191 & gencost_spzero ,
6b2230d510 Ou W* 0192 & gencost_wei1d ,
69f8f4c84c Patr* 0193 & gencost_avgperiod ,
92d10a5666 Gael* 0194 & gencost_nrecperiod ,
49484c0542 Gael* 0195 & gencost_startdate1 ,
0196 & gencost_startdate2 ,
0197 & gencost_enddate1 ,
0198 & gencost_enddate2 ,
89d9f37d02 An T* 0199 & gencost_smooth2Ddiffnbt ,
6b2230d510 Ou W* 0200 & gencost_is1d ,
5cce2b5d76 Gael* 0201 & gencost_is3d ,
17944dd1e8 Gael* 0202 & gencost_msk_is3d ,
aa93ca8e85 Ciar* 0203 & gencost_useDensityMask ,
0204 & gencost_refPressure ,
0205 & gencost_sigmaLow ,
0206 & gencost_sigmaHigh ,
0207 & gencost_tanhScale ,
8227586f14 An T* 0208 & gencost_timevaryweight ,
8c157ed454 Patr* 0209 & mult_gencost
5cce2b5d76 Gael* 0210 #endif /* ALLOW_GENCOST_CONTRIBUTION */
8c157ed454 Patr* 0211
58023ceed9 Gael* 0212 IF ( .NOT. useECCO ) THEN
0213
0214 _BEGIN_MASTER (myThid )
0215
0216
0217 CALL PACKAGES_UNUSED_MSG ( 'useECCO' , ' ' , ' ' )
0218 _END_MASTER (myThid )
0219 RETURN
0220 ENDIF
0221
5001c65f45 Patr* 0222 _BEGIN_MASTER ( myThid )
cda1c18f72 Jean* 0223 nRetired = 0
0c0a2531ee Patr* 0224
5001c65f45 Patr* 0225
984d1519c6 Gael* 0226
0227
9f3ef6a237 Gael* 0228 #if (defined (ALLOW_SSH_COST_CONTRIBUTION ) && \
0229
984d1519c6 Gael* 0230 using_cost_altim =.TRUE.
0231 #else
0232 using_cost_altim =.FALSE.
0233 #endif
86576896a0 Gael* 0234 #ifdef ALLOW_SEAICE_COST_CONTRIBUTION
0235 using_cost_seaice =.TRUE.
0236 #else
0237 using_cost_seaice =.FALSE.
0238 #endif
c2dc7fc123 An T* 0239
f1e0eff8fe An T* 0240
0241 using_cost_transp =.FALSE.
c2dc7fc123 An T* 0242
984d1519c6 Gael* 0243 #ifdef ALLOW_SST_COST_CONTRIBUTION
0244 using_cost_sst =.TRUE.
0245 #else
0246 using_cost_sst =.FALSE.
0247 #endif
0248
0249
5001c65f45 Patr* 0250 data_errfile = ' '
0251 tbarfile = 'tbar'
0252 sbarfile = 'sbar'
61d31abfa0 Jean* 0253 sstbarfile = 'sstbar'
5001c65f45 Patr* 0254 psbarfile = 'psbar'
e284aa03cc Patr* 0255 bpbarfile = 'bpbar'
d534937279 Matt* 0256 iestaubarfile = 'iestaubar'
5001c65f45 Patr* 0257 ubarfile = 'ubar'
0258 vbarfile = 'vbar'
0259 wbarfile = 'wbar'
0260 tauxbarfile = 'tauxbar'
0261 tauybarfile = 'tauybar'
8c157ed454 Patr* 0262 hfluxmeanbarfile = 'hfluxmeanbar'
0263 sfluxmeanbarfile = 'sfluxmeanbar'
7484238bfd Patr* 0264 costTranspDataFile = ' '
9f3ef6a237 Gael* 0265 using_tpj = .false.
5001c65f45 Patr* 0266 topexstartdate1 = 0
0267 topexstartdate2 = 0
0268 topexperiod = 0. _d 0
69531f39dd Patr* 0269 topexintercept = 0. _d 0
0270 topexslope = 0. _d 0
5001c65f45 Patr* 0271 topexfile = ' '
d39bcf1629 An T* 0272
0273
9f3ef6a237 Gael* 0274 using_mdt =.FALSE.
d39bcf1629 An T* 0275 mdtdatfile = ' '
0276 mdtstartdate1 = 19930101
0277 mdtstartdate2 = 0
0278 mdtenddate1 = 20041231
0279 mdtenddate2 = 0
5001c65f45 Patr* 0280 using_ers = .false.
0281 ersstartdate1 = 0
0282 ersstartdate2 = 0
0283 ersperiod = 0. _d 0
69531f39dd Patr* 0284 ersintercept = 0. _d 0
0285 ersslope = 0. _d 0
5001c65f45 Patr* 0286 ersfile = ' '
4fd9553bc8 Patr* 0287 using_gfo = .false.
0288 gfostartdate1 = 0
0289 gfostartdate2 = 0
0290 gfoperiod = 0. _d 0
69531f39dd Patr* 0291 gfointercept = 0. _d 0
0292 gfoslope = 0. _d 0
4fd9553bc8 Patr* 0293 gfofile = ' '
5001c65f45 Patr* 0294 scatstartdate1 = 0
0295 scatstartdate2 = 0
0296 scatperiod = 0. _d 0
0297 scatxdatfile = ' '
0298 scatydatfile = ' '
0299 ssh_errfile = ' '
350884e02b Patr* 0300 tp_errfile = ' '
0301 ers_errfile = ' '
0302 gfo_errfile = ' '
5001c65f45 Patr* 0303 scatx_errfile = ' '
0304 scaty_errfile = ' '
0305 tmistartdate1 = 0
0306 tmistartdate2 = 0
0307 tmidatfile = ' '
0308 sststartdate1 = 0
0309 sststartdate2 = 0
0310 sstdatfile = ' '
0311 sssstartdate1 = 0
0312 sssstartdate2 = 0
0313 sssdatfile = ' '
e284aa03cc Patr* 0314 bpstartdate1 = 0
0315 bpstartdate2 = 0
0316 bpdatfile = ' '
d534937279 Matt* 0317 iesstartdate1 = 0
0318 iesstartdate2 = 0
0319 iesdatfile = ' '
5001c65f45 Patr* 0320 tdatfile = ' '
0321 sdatfile = ' '
0322 ctdtfile = ' '
0323 ctdsfile = ' '
0324 ctdtclimfile = ' '
0325 ctdsclimfile = ' '
0326 curmtrufile = ' '
0327 curmtrvfile = ' '
0328 xbtfile = ' '
0329 argotstartdate1 = 0
0330 argotstartdate2 = 0
0331 argotfile = ' '
0332 argosstartdate1 = 0
0333 argosstartdate2 = 0
0334 argosfile = ' '
0335 udriftfile = ' '
0336 vdriftfile = ' '
0337 udrifterrfile = ' '
0338 vdrifterrfile = ' '
0339 salterrfile = ' '
0340 temperrfile = ' '
0341 velerrfile = ' '
76ae088f98 Patr* 0342 salt0errfile = ' '
0343 temp0errfile = ' '
30125f96e6 Patr* 0344 etan0errfile = ' '
0345 uvel0errfile = ' '
0346 vvel0errfile = ' '
76ae088f98 Patr* 0347 vel0errfile = ' '
da86325d70 Gael* 0348 ssterrfile = ' '
0349 ssserrfile = ' '
49b79243e8 Gael* 0350 bperrfile = ' '
d534937279 Matt* 0351 ieserrfile = ' '
5001c65f45 Patr* 0352 hflux_errfile = ' '
0353 sflux_errfile = ' '
0354 tauu_errfile = ' '
0355 tauv_errfile = ' '
0356 atemp_errfile = ' '
0357 aqh_errfile = ' '
d7ee8fe52e Patr* 0358 precip_errfile = ' '
5ceb56546e Patr* 0359 swflux_errfile = ' '
02568c8633 Patr* 0360 swdown_errfile = ' '
fed442e87e Patr* 0361 lwflux_errfile = ' '
0362 lwdown_errfile = ' '
0363 evap_errfile = ' '
0364 snowprecip_errfile = ' '
0365 apressure_errfile = ' '
5001c65f45 Patr* 0366 uwind_errfile = ' '
0367 vwind_errfile = ' '
0368 geoid_errfile = ' '
0369 geoid_covariancefile = ' '
bad79706a2 Gael* 0370 bottomdrag_errfile = ' '
0371 edtau_errfile = ' '
0372 kapgm_errfile = ' '
c9ba60667e Patr* 0373 kapredi_errfile = ' '
bad79706a2 Gael* 0374 diffkr_errfile = ' '
55bf8021cb Gael* 0375 do k =1,NUSERCOST
cda1c18f72 Jean* 0376 usercost_datafile (k ) = ' '
0377 usercost_errfile (k ) = ' '
55bf8021cb Gael* 0378 enddo
087d6770fc Gael* 0379 do k =1,NSSHV4COST
2a529ed9f1 Gael* 0380 sshv4cost_scalefile (k ) = ' '
0381 sshv4cost_errfile (k ) = ' '
0382 sshv4cost_errfactor (k ) = 1. _d 0
087d6770fc Gael* 0383 enddo
c9ba60667e Patr* 0384 tpTimeMaskFile = ' '
0385 ersTimeMaskFile = ' '
0386 gfoTimeMaskFile = ' '
0387
5001c65f45 Patr* 0388 mult_hflux = 0. _d 0
0389 mult_sflux = 0. _d 0
4bf6063982 Patr* 0390 mult_hfluxmm = 0. _d 0
0391 mult_sfluxmm = 0. _d 0
5001c65f45 Patr* 0392 mult_tauu = 0. _d 0
0393 mult_tauv = 0. _d 0
0394 mult_atemp = 0. _d 0
0395 mult_aqh = 0. _d 0
d7ee8fe52e Patr* 0396 mult_precip = 0. _d 0
5ceb56546e Patr* 0397 mult_swflux = 0. _d 0
02568c8633 Patr* 0398 mult_swdown = 0. _d 0
fed442e87e Patr* 0399 mult_lwflux = 0. _d 0
0400 mult_lwdown = 0. _d 0
0401 mult_evap = 0. _d 0
0402 mult_snowprecip = 0. _d 0
0403 mult_apressure = 0. _d 0
5001c65f45 Patr* 0404 mult_uwind = 0. _d 0
0405 mult_vwind = 0. _d 0
0406 mult_hmean = 0. _d 0
706976e7e4 Patr* 0407 mult_h = 1. _d 0
0408 mult_tp = 0. _d 0
0409 mult_ers = 0. _d 0
0410 mult_gfo = 0. _d 0
5001c65f45 Patr* 0411 mult_temp = 0. _d 0
0412 mult_salt = 0. _d 0
0413 mult_temp0 = 0. _d 0
0414 mult_salt0 = 0. _d 0
30125f96e6 Patr* 0415 mult_etan0 = 0. _d 0
0416 mult_uvel0 = 0. _d 0
0417 mult_vvel0 = 0. _d 0
5001c65f45 Patr* 0418 mult_tmi = 0. _d 0
0419 mult_sst = 0. _d 0
0420 mult_sss = 0. _d 0
e284aa03cc Patr* 0421 mult_bp = 0. _d 0
d534937279 Matt* 0422 mult_ies = 0. _d 0
5001c65f45 Patr* 0423 mult_ctdt = 0. _d 0
0424 mult_ctds = 0. _d 0
0425 mult_ctdtclim = 0. _d 0
0426 mult_ctdsclim = 0. _d 0
0427 mult_xbt = 0. _d 0
0428 mult_argot = 0. _d 0
0429 mult_argos = 0. _d 0
0430 mult_drift = 0. _d 0
0431 mult_tdrift = 0. _d 0
0432 mult_sdrift = 0. _d 0
0433 mult_wdrift = 0. _d 0
0434 mult_scatx = 0. _d 0
0435 mult_scaty = 0. _d 0
c509d7e04a Gael* 0436 #if (defined (ALLOW_CTRL ) && defined (ALLOW_OBCS ))
5001c65f45 Patr* 0437 mult_obcsn = 0. _d 0
0438 mult_obcss = 0. _d 0
0439 mult_obcsw = 0. _d 0
0440 mult_obcse = 0. _d 0
0441 mult_obcsvol = 0. _d 0
0442 mult_ageos = 0. _d 0
c509d7e04a Gael* 0443 #endif
5001c65f45 Patr* 0444 mult_curmtr = 0. _d 0
0445 mult_kapgm = 0. _d 0
cda1c18f72 Jean* 0446 mult_kapredi = 0. _d 0
5001c65f45 Patr* 0447 mult_diffkr = 0. _d 0
bad79706a2 Gael* 0448 mult_edtau = 0. _d 0
0449 mult_bottomdrag = 0. _d 0
55bf8021cb Gael* 0450 do k =1,NUSERCOST
cda1c18f72 Jean* 0451 mult_usercost (k ) = 0. _d 0
55bf8021cb Gael* 0452 enddo
087d6770fc Gael* 0453 do k =1,NSSHV4COST
cda1c18f72 Jean* 0454 mult_sshv4cost (k ) = 0. _d 0
087d6770fc Gael* 0455 enddo
f2b10b358f Patr* 0456 mult_smooth_ic = 0. _d 0
0457 mult_smooth_bc = 0. _d 0
d851bc4637 Patr* 0458 mult_transp = 0. _d 0
d7ee8fe52e Patr* 0459 watemp0 = 0. _d 0
0460 waqh0 = 0. _d 0
0461 wprecip0 = 0. _d 0
5ceb56546e Patr* 0462 wswflux0 = 0. _d 0
02568c8633 Patr* 0463 wswdown0 = 0. _d 0
fed442e87e Patr* 0464 wlwflux0 = 0. _d 0
0465 wlwdown0 = 0. _d 0
0466 wevap0 = 0. _d 0
0467 wsnowprecip0 = 0. _d 0
0468 wapressure0 = 0. _d 0
d7ee8fe52e Patr* 0469 wwind0 = 0. _d 0
0e094cd207 Patr* 0470 wdiffkr0 = 1. _d 0
0471 wkapgm0 = 1. _d 0
cda1c18f72 Jean* 0472 wkapredi0 = 1. _d 0
0e094cd207 Patr* 0473 wedtau0 = 1. _d 0
d7ee8fe52e Patr* 0474 whflux0 = 0. _d 0
0475 wsflux0 = 0. _d 0
0476 wtau0 = 0. _d 0
bad79706a2 Gael* 0477 wbottomdrag0 = 1. _d 0
eeaf8785b4 Patr* 0478 wmean_atemp = 0. _d 0
0479 wmean_aqh = 0. _d 0
0480 wmean_precip = 0. _d 0
0481 wmean_swflux = 0. _d 0
02568c8633 Patr* 0482 wmean_swdown = 0. _d 0
fed442e87e Patr* 0483 wmean_lwflux = 0. _d 0
0484 wmean_lwdown = 0. _d 0
0485 wmean_evap = 0. _d 0
0486 wmean_snowprecip = 0. _d 0
0487 wmean_apressure = 0. _d 0
eeaf8785b4 Patr* 0488 wmean_wind = 0. _d 0
0489 wmean_hflux = 0. _d 0
0490 wmean_sflux = 0. _d 0
0491 wmean_tau = 0. _d 0
6bec9496d6 Jean* 0492 cost_iprec = precFloat32
906a61c194 Ou W* 0493 ecco_output_sterGloH = .FALSE.
0494 ecco_keepTSeriesOutp_open = .FALSE.
5001c65f45 Patr* 0495
cda1c18f72 Jean* 0496
0497 cost_yftype = 'RL'
0498 topexmeanfile = ' '
0499 using_cost_bp = .FALSE.
0500 using_cost_scat = .FALSE.
d39bcf1629 An T* 0501
0502
9aaf43452b Patr* 0503 WRITE (msgBuf ,'(A)' ) 'ECCO_READPARMS: opening data.ecco'
0504 CALL PRINT_MESSAGE ( msgBuf , standardMessageUnit ,
44d3986245 Jean* 0505 & SQUEEZE_RIGHT , myThid )
9aaf43452b Patr* 0506
0507 CALL OPEN_COPY_DATA_FILE (
0508 I 'data.ecco' , 'ECCO_READPARMS' ,
0509 O iUnit ,
0510 I myThid )
0511
0512 READ ( iUnit , nml = ecco_cost_nml )
0513
8f7d13d0c9 Jean* 0514 WRITE (msgBuf ,'(A)' )
8c157ed454 Patr* 0515 & 'ECCO_READPARMS: finished reading #1: ecco_cost_nml'
9aaf43452b Patr* 0516 CALL PRINT_MESSAGE ( msgBuf , standardMessageUnit ,
44d3986245 Jean* 0517 & SQUEEZE_RIGHT , myThid )
0c0a2531ee Patr* 0518
6bec9496d6 Jean* 0519
cda1c18f72 Jean* 0520
0521 IF ( cost_yftype .NE. 'RL' ) THEN
0522 WRITE (msgBuf ,'(2A)' ) 'S/R ECCO_READPARMS: ' ,
0523 & '"cost_yftype" is no longer allowed in file "data.ecco"'
0524 CALL PRINT_ERROR ( msgBuf , myThid )
d39bcf1629 An T* 0525 nRetired = nRetired + 1
cda1c18f72 Jean* 0526 ENDIF
0527 IF ( using_cost_bp ) THEN
0528 WRITE (msgBuf ,'(2A)' ) 'S/R ECCO_READPARMS: ' ,
0529 & '"using_cost_bp" is no longer allowed in file "data.ecco"'
0530 CALL PRINT_ERROR ( msgBuf , myThid )
0531 nRetired = nRetired + 1
0532 ENDIF
0533 IF ( using_cost_scat ) THEN
0534 WRITE (msgBuf ,'(2A)' ) 'S/R ECCO_READPARMS: ' ,
0535 & '"using_cost_scat" is no longer allowed in file "data.ecco"'
0536 CALL PRINT_ERROR ( msgBuf , myThid )
0537 nRetired = nRetired + 1
0538 ENDIF
0539 IF ( topexmeanfile .NE. ' ' ) THEN
d39bcf1629 An T* 0540 WRITE (msgBuf ,'(A,A)' )
0541 & 'S/R ECCO_READPARMS: "topexmeanfile" ' ,
0542 & 'is no longer allowed in file "data.ecco"'
0543 CALL PRINT_ERROR ( msgBuf , myThid )
0544 WRITE (msgBuf ,'(A)' )
0545 & 'S/R ECCO_READPARMS: set "mdtdatfile" instead'
0546 CALL PRINT_ERROR ( msgBuf , myThid )
cda1c18f72 Jean* 0547 nRetired = nRetired + 1
d39bcf1629 An T* 0548 ENDIF
0549 IF ( nRetired .GT. 0 ) THEN
0550 WRITE (msgBuf ,'(2A)' ) 'S/R ECCO_READPARMS: ' ,
0551 & 'Error reading parameter file "data.ecco"'
0552 CALL PRINT_ERROR ( msgBuf , myThid )
0553 WRITE (msgBuf ,'(A)' )
cda1c18f72 Jean* 0554 & 'some out-of-date parameters were found in the namelist(s)'
d39bcf1629 An T* 0555 CALL PRINT_ERROR ( msgBuf , myThid )
cda1c18f72 Jean* 0556 CALL ALL_PROC_DIE ( 0 )
d39bcf1629 An T* 0557 STOP 'ABNORMAL END: S/R ECCO_READPARMS'
0558 ENDIF
6bec9496d6 Jean* 0559
8c157ed454 Patr* 0560 _END_MASTER ( myThid )
0561
3bd7a4549f Gael* 0562
0563
0564
0565 _BEGIN_MASTER ( myThid )
0566 do k2 =1,NGENPPROC
447bdc4b79 Gael* 0567 no_preproc (k2 ) = ' '
0568 no_preproc_c (k2 ) = ' '
0569 no_preproc_i (k2 ) = 0
0570 no_preproc_r (k2 ) = 0. _d 0
0571 no_posproc (k2 ) = ' '
0572 no_posproc_c (k2 ) = ' '
0573 no_posproc_i (k2 ) = 0
0574 no_posproc_r (k2 ) = 0. _d 0
1b2ca1f822 Gael* 0575 clim_preproc (k2 ) = ' '
0576 clim_preproc_i (k2 ) = 0
dca8cfd541 Gael* 0577 anom_preproc (k2 ) = ' '
3bd7a4549f Gael* 0578 enddo
1b2ca1f822 Gael* 0579 clim_preproc (1) = 'clim'
0580 clim_preproc_i (1) = 12
dca8cfd541 Gael* 0581 anom_preproc (1) = 'anom'
3bd7a4549f Gael* 0582 _END_MASTER ( myThid )
0583
8c157ed454 Patr* 0584 #ifdef ALLOW_GENCOST_CONTRIBUTION
0585
0586 _BEGIN_MASTER ( myThid )
0587
0588
0589 do k =1,NGENCOST
130273d46b Gael* 0590 using_gencost (k ) = .FALSE.
0591 gencost_flag (k ) = 0
8c157ed454 Patr* 0592 gencost_avgperiod (k ) = ' '
49484c0542 Gael* 0593 gencost_startdate1 (k ) = 0
0594 gencost_startdate2 (k ) = 0
0595 gencost_enddate1 (k ) = 0
0596 gencost_enddate2 (k ) = 0
d3a3bee25e Patr* 0597 gencost_datafile (k ) = ' '
d181f9168b Gael* 0598 gencost_name (k ) = 'gencost'
3bd7a4549f Gael* 0599 do k2 =1,NGENPPROC
0600 gencost_preproc (k2 ,k ) = ' '
0601 gencost_preproc_c (k2 ,k ) = ' '
0602 gencost_preproc_i (k2 ,k ) = 0
0603 gencost_preproc_r (k2 ,k ) = 0. _d 0
0604 gencost_posproc (k2 ,k ) = ' '
0605 gencost_posproc_c (k2 ,k ) = ' '
0606 gencost_posproc_i (k2 ,k ) = 0
0607 gencost_posproc_r (k2 ,k ) = 0. _d 0
0608 enddo
92d10a5666 Gael* 0609 gencost_outputlevel (k )= 0
d3a3bee25e Patr* 0610 gencost_errfile (k ) = ' '
81e05fa829 Gael* 0611 gencost_itracer (k ) = 1
a565f49279 Jean* 0612 gencost_kLev_select (k )= 1
47d80787ea Gael* 0613 gencost_mask (k ) = ' '
d181f9168b Gael* 0614 gencost_barfile (k ) = ' '
89d9f37d02 An T* 0615 gencost_spmin (k ) = 0. _d 0
0616 gencost_spmax (k ) = 0. _d 0
0617 gencost_spzero (k ) = 9876. _d 0
6b2230d510 Ou W* 0618 gencost_wei1d (k ) = 0. _d 0
bf3a68e762 Gael* 0619 mult_gencost (k ) = 1. _d 0
6b2230d510 Ou W* 0620 gencost_is1d (k ) = .FALSE.
5cce2b5d76 Gael* 0621 gencost_is3d (k ) = .FALSE.
aa93ca8e85 Ciar* 0622 gencost_useDensityMask (k ) = .FALSE.
0623 gencost_refPressure (k ) = 0. _d 0
0624 gencost_sigmaLow (k ) = 0. _d 0
0625 gencost_sigmaHigh (k ) = 1. _d 3
0626 gencost_tanhScale (k ) = 1. _d 5
447bdc4b79 Gael* 0627 gencost_pointer3d (k ) = 0
44d3986245 Jean* 0628 gencost_msk_is3d (k ) = .FALSE.
17944dd1e8 Gael* 0629 gencost_msk_pointer3d (k ) = 0
447bdc4b79 Gael* 0630
f8e779c983 antn* 0631 call ecco_zero (gencost_mskCsurf (1-OLx ,1-OLy ,1,1,k ),
0632 & 1,zeroRL ,myThid )
0633 call ecco_zero (gencost_mskWsurf (1-OLx ,1-OLy ,1,1,k ),
0634 & 1,zeroRL ,myThid )
0635 call ecco_zero (gencost_mskSsurf (1-OLx ,1-OLy ,1,1,k ),
0636 & 1,zeroRL ,myThid )
9f5240b52a Jean* 0637 do k2 = 1,Nr
447bdc4b79 Gael* 0638 gencost_mskVertical (k2 ,k )= 1. _d 0
f8e779c983 antn* 0639 enddo
447bdc4b79 Gael* 0640
0641 gencost_nrecperiod (k ) = 0
0642 gencost_scalefile (k ) = ' '
b8fab26f74 Gael* 0643 gencost_smooth2Ddiffnbt (k ) = 0
447bdc4b79 Gael* 0644 gencost_timevaryweight (k ) = .FALSE.
8c157ed454 Patr* 0645 enddo
0646
17944dd1e8 Gael* 0647 #ifdef ALLOW_GENCOST3D
0648 do k =1,NGENCOST3D
0649
f8e779c983 antn* 0650 call ecco_zero (gencost_mskC (1-OLx ,1-OLy ,1,1,1,k ),
0651 & Nr ,zeroRL ,myThid )
0652 call ecco_zero (gencost_mskW (1-OLx ,1-OLy ,1,1,1,k ),
0653 & Nr ,zeroRL ,myThid )
0654 call ecco_zero (gencost_mskS (1-OLx ,1-OLy ,1,1,1,k ),
0655 & Nr ,zeroRL ,myThid )
17944dd1e8 Gael* 0656 enddo
0657 #endif
0658
8c157ed454 Patr* 0659 READ ( iUnit , nml = ecco_gencost_nml )
0660
8f7d13d0c9 Jean* 0661 WRITE (msgBuf ,'(A)' )
8c157ed454 Patr* 0662 & 'ECCO_READPARMS: finished reading #2: ecco_gencost_nml'
0663 CALL PRINT_MESSAGE ( msgBuf , standardMessageUnit ,
44d3986245 Jean* 0664 & SQUEEZE_RIGHT , myThid )
8c157ed454 Patr* 0665
5cce2b5d76 Gael* 0666 gencost_k3d =1
17944dd1e8 Gael* 0667 gencost_msk_k3d =1
5cce2b5d76 Gael* 0668
130273d46b Gael* 0669 do k =1,NGENCOST
5cce2b5d76 Gael* 0670
0671
e7db56ba12 Gael* 0672 if ( (gencost_name (k ).EQ. 'sshv4-mdt' ).OR.
0673 & (gencost_name (k ).EQ. 'sshv4-tp' ).OR.
0674 & (gencost_name (k ).EQ. 'sshv4-ers' ).OR.
0675 & (gencost_name (k ).EQ. 'sshv4-gfo' ).OR.
0676 & (gencost_name (k ).EQ. 'sshv4-lsc' ).OR.
4e6b53b435 Gael* 0677 & (gencost_name (k ).EQ. 'sshv4-gmsl' ).OR.
49484c0542 Gael* 0678 & (gencost_name (k ).EQ. 'bpv4-grace' ).OR.
e7db56ba12 Gael* 0679 & (gencost_name (k ).EQ. 'sstv4-amsre' ).OR.
130273d46b Gael* 0680 & (gencost_name (k ).EQ. 'sstv4-amsre-lsc' ) ) then
0681 gencost_flag (k )=-1
0682 using_gencost (k )=.TRUE.
877e0f385b An T* 0683
0684
0685
0686
447bdc4b79 Gael* 0687 elseif ( (gencost_barfile (k )(1:9).EQ. 'm_boxmean' ).OR.
0688 & (gencost_barfile (k )(1:9).EQ. 'm_horflux' ) ) then
877e0f385b An T* 0689 gencost_flag (k )=-3
0690 using_gencost (k )=.TRUE.
0691 il =ilnblnk (gencost_name (k ))
447bdc4b79 Gael* 0692 WRITE (msgBuf ,'(A,I3,3A,L5)' ) 'gencost k_boxmean/horflux=' ,
0693 & k ,' ' ,gencost_name (k )(1:il ),': ' ,using_gencost (k )
877e0f385b An T* 0694 CALL PRINT_MESSAGE ( msgBuf , standardMessageUnit ,
44d3986245 Jean* 0695 & SQUEEZE_RIGHT , myThid )
877e0f385b An T* 0696
17944dd1e8 Gael* 0697 if ( gencost_msk_is3d (k ) ) then
0698 gencost_msk_pointer3d (k )=gencost_msk_k3d
0699 gencost_msk_k3d =gencost_msk_k3d +1
0700 endif
0701
877e0f385b An T* 0702
0703 elseif ( gencost_name (k )(1:6).EQ. 'transp' ) then
0704 gencost_flag (k )=-4
0705 using_gencost (k )=.TRUE.
0706 il =ilnblnk (gencost_name (k ))
447bdc4b79 Gael* 0707 WRITE (msgBuf ,'(A,i3,3A,L5)' ) 'gencost k_transp=' ,
0708 & k ,' ' ,gencost_name (k )(1:il ),': ' ,using_gencost (k )
877e0f385b An T* 0709 CALL PRINT_MESSAGE ( msgBuf , standardMessageUnit ,
44d3986245 Jean* 0710 & SQUEEZE_RIGHT , myThid )
877e0f385b An T* 0711
df462307fb Timo* 0712
0713 elseif ( gencost_name (k )(1:3).EQ. 'moc' ) then
0714 gencost_flag (k )=-5
0715 using_gencost (k )=.TRUE.
0716 il =ilnblnk (gencost_name (k ))
0717 WRITE (msgBuf ,'(A,i2,2A,L5)' ) 'gencost k_moc=' ,
0718 & k ,gencost_name (k )(1:il ),': ' ,using_gencost (k )
0719 CALL PRINT_MESSAGE ( msgBuf , standardMessageUnit ,
44d3986245 Jean* 0720 & SQUEEZE_RIGHT , myThid )
df462307fb Timo* 0721
5cce2b5d76 Gael* 0722
f8e779c983 antn* 0723 elseif (
a443ec37be Gael* 0724 & (gencost_name (k ).EQ. 'siv4-conc' ).OR.
0a8c2c2ff2 An T* 0725 & (gencost_name (k ).EQ. 'siv4-deconc' ).OR.
0726 & (gencost_name (k ).EQ. 'siv4-exconc' ) ) then
a443ec37be Gael* 0727 gencost_flag (k )=2
0728 using_gencost (k )=.TRUE.
130273d46b Gael* 0729
5cce2b5d76 Gael* 0730 elseif ( gencost_datafile (k ) .NE. ' ' ) then
130273d46b Gael* 0731 gencost_flag (k )=1
0732 using_gencost (k )=.TRUE.
0733 endif
49484c0542 Gael* 0734
bf3a68e762 Gael* 0735
0736 if ( (gencost_barfile (k )(1:7).EQ. 'm_theta' ).OR.
0737 & (gencost_barfile (k )(1:6).EQ. 'm_salt' ).OR.
0738 & (gencost_barfile (k )(1:8).EQ. 'm_diffkr' ).OR.
0739 & (gencost_barfile (k )(1:7).EQ. 'm_kapgm' ).OR.
447bdc4b79 Gael* 0740 & (gencost_barfile (k )(1:9).EQ. 'm_kapredi' ).OR.
0741 & (gencost_barfile (k )(1:7).EQ. 'm_trVol' ).OR.
c13fbfe5a9 Nora* 0742 & (gencost_barfile (k )(1:8).EQ. 'm_trHeat' ).OR.
0743 & (gencost_barfile (k )(1:8).EQ. 'm_trSalt' )
877e0f385b An T* 0744 & )
bf3a68e762 Gael* 0745 & gencost_is3d (k )=.TRUE.
0746
0747
5cce2b5d76 Gael* 0748 if ( gencost_is3d (k ) ) then
0749 gencost_pointer3d (k )=gencost_k3d
0750 gencost_k3d =gencost_k3d +1
0751 endif
0752
130273d46b Gael* 0753 enddo
0754
447bdc4b79 Gael* 0755
0756 do k =1,NGENCOST
c17b89ca05 Gael* 0757 kk =gencost_msk_pointer3d (k )
47d80787ea Gael* 0758 if ( ( gencost_mask (k ) .NE. ' ' ).AND.
0d7023b5ce Jean* 0759 & (gencost_flag (k ).EQ. -3 .or.
df462307fb Timo* 0760 & gencost_flag (k ).eq. -4 .or.
0761 & gencost_flag (k ).eq. -5 ) ) then
447bdc4b79 Gael* 0762
47d80787ea Gael* 0763 il = ilnblnk (gencost_mask (k ))
0764 write (tempfile (1:128),'(2A)' ) gencost_mask (k )(1:il ),'C'
447bdc4b79 Gael* 0765 inquire( file=tempfile (1:il +1), exist=exst )
17944dd1e8 Gael* 0766 if (exst .AND. (.NOT. gencost_msk_is3d (k ))) then
f8e779c983 antn* 0767 CALL READ_REC_3D_RL ( tempfile , cost_iprec , 1,
0768 & gencost_mskCsurf (1-OLx ,1-OLy ,1,1,k ), 1, 1, myThid )
17944dd1e8 Gael* 0769 #ifdef ALLOW_GENCOST3D
0770 elseif (exst .AND. (gencost_msk_pointer3d (k ).LE.
0771 & NGENCOST3D )) then
f8e779c983 antn* 0772 CALL READ_REC_3D_RL ( tempfile , cost_iprec , Nr ,
0773 & gencost_mskC (1-OLx ,1-OLy ,1,1,1,kk ), 1, 1, myThid )
17944dd1e8 Gael* 0774 #endif
447bdc4b79 Gael* 0775 endif
0776
47d80787ea Gael* 0777 il = ilnblnk (gencost_mask (k ))
0778 write (tempfile (1:128),'(2A)' ) gencost_mask (k )(1:il ),'W'
447bdc4b79 Gael* 0779 inquire( file=tempfile (1:il +1), exist=exst )
17944dd1e8 Gael* 0780 if (exst .AND. (.NOT. gencost_msk_is3d (k ))) then
f8e779c983 antn* 0781 CALL READ_REC_3D_RL ( tempfile , cost_iprec , 1,
0782 & gencost_mskWsurf (1-OLx ,1-OLy ,1,1,k ), 1, 1, myThid )
17944dd1e8 Gael* 0783 #ifdef ALLOW_GENCOST3D
0784 elseif (exst .AND. (gencost_msk_pointer3d (k ).LE.
0785 & NGENCOST3D )) then
f8e779c983 antn* 0786 CALL READ_REC_3D_RL ( tempfile , cost_iprec , Nr ,
0787 & gencost_mskW (1-OLx ,1-OLy ,1,1,1,kk ), 1, 1, myThid )
17944dd1e8 Gael* 0788 #endif
447bdc4b79 Gael* 0789 endif
0790
47d80787ea Gael* 0791 il = ilnblnk (gencost_mask (k ))
0792 write (tempfile (1:128),'(2A)' ) gencost_mask (k )(1:il ),'S'
447bdc4b79 Gael* 0793 inquire( file=tempfile (1:il +1), exist=exst )
17944dd1e8 Gael* 0794 if (exst .AND. (.NOT. gencost_msk_is3d (k ))) then
f8e779c983 antn* 0795 CALL READ_REC_3D_RL ( tempfile , cost_iprec , 1,
0796 & gencost_mskSsurf (1-OLx ,1-OLy ,1,1,k ), 1, 1, myThid )
17944dd1e8 Gael* 0797 #ifdef ALLOW_GENCOST3D
0798 elseif (exst .AND. (gencost_msk_pointer3d (k ).LE.
0799 & NGENCOST3D )) then
f8e779c983 antn* 0800 CALL READ_REC_3D_RL ( tempfile , cost_iprec , Nr ,
0801 & gencost_mskS (1-OLx ,1-OLy ,1,1,1,kk ), 1, 1, myThid )
17944dd1e8 Gael* 0802 #endif
447bdc4b79 Gael* 0803 endif
0804
47d80787ea Gael* 0805 il = ilnblnk (gencost_mask (k ))
0806 write (tempfile (1:128),'(2A)' ) gencost_mask (k )(1:il ),'K'
447bdc4b79 Gael* 0807 inquire( file=tempfile (1:il +1), exist=exst )
0808 if (exst ) then
0809 ioUnit = 0
f8e779c983 antn* 0810 CALL MDS_READVEC_LOC ( tempfile , cost_iprec , ioUnit , 'RL' ,
0811 & Nr , gencost_mskVertical (1,k ), dummyRS , 0, 0, 1, myThid )
447bdc4b79 Gael* 0812 endif
0813
0814 endif
0815 enddo
0816
8c157ed454 Patr* 0817
0818 do k =1,NGENCOST
69f8f4c84c Patr* 0819
9aeddbf372 Gael* 0820
985662a3b3 Gael* 0821 if (gencost_barfile (k ).NE. ' ' ) then
9aeddbf372 Gael* 0822 if (gencost_barfile (k )(1:2).NE. 'm_' ) then
985662a3b3 Gael* 0823 il =ilnblnk (gencost_barfile (k ))
0824 WRITE (msgBuf ,'(A,A)' ) 'm_' ,gencost_barfile (k )(1:il )
0825 gencost_barfile (k )=msgBuf (1:il +2)
9aeddbf372 Gael* 0826 endif
0827 if (gencost_barfile (k )(1:8).EQ. 'm_tauZon' ) then
0828 il =ilnblnk (gencost_barfile (k ))
0829 WRITE (msgBuf ,'(A,A)' ) 'm_ustress' ,gencost_barfile (k )(9:il )
0830 gencost_barfile (k )=msgBuf (1:il +1)
0831 endif
0832 if (gencost_barfile (k )(1:8).EQ. 'm_tauMer' ) then
0833 il =ilnblnk (gencost_barfile (k ))
0834 WRITE (msgBuf ,'(A,A)' ) 'm_vstress' ,gencost_barfile (k )(9:il )
0835 gencost_barfile (k )=msgBuf (1:il +1)
0836 endif
985662a3b3 Gael* 0837 endif
0838
877e0f385b An T* 0839
0840 if ( using_gencost (k )) then
0841 if ( gencost_flag (k ).GE. 1 ) then
7a7bc10d23 Patr* 0842
0843 if ( gencost_name (k ) .EQ. 'gencost' )
1fa8e77fb9 Gael* 0844 & gencost_name (k ) = gencost_datafile (k )
0845
7a7bc10d23 Patr* 0846 if ( gencost_avgperiod (k ) .NE. 'day ' .AND.
69f8f4c84c Patr* 0847 & gencost_avgperiod (k ) .NE. 'DAY ' .AND.
0848 & gencost_avgperiod (k ) .NE. 'month' .AND.
0849 & gencost_avgperiod (k ) .NE. 'MONTH' .AND.
e7d9258ace Gael* 0850 & gencost_avgperiod (k ) .NE. 'step' .AND.
0851 & gencost_avgperiod (k ) .NE. 'STEP' .AND.
985662a3b3 Gael* 0852 & gencost_avgperiod (k ) .NE. 'const' .AND.
0853 & gencost_avgperiod (k ) .NE. 'CONST' .AND.
69f8f4c84c Patr* 0854 & gencost_avgperiod (k ) .NE. 'year ' .AND.
0855 & gencost_avgperiod (k ) .NE. 'YEAR ' ) then
e7db56ba12 Gael* 0856
0857 il =ilnblnk (gencost_name (k ))
0858 WRITE (msgBuf ,'(A,I2,2A)' )
0859 & 'ERROR in ECCO_READPARMS: for gencost' ,k ,
0860 & ' ' ,gencost_name (k )(1:il )
0861 CALL PRINT_MESSAGE ( msgBuf , errorMessageUnit ,
44d3986245 Jean* 0862 & SQUEEZE_RIGHT , myThid )
8f7d13d0c9 Jean* 0863 STOP
8c157ed454 Patr* 0864 & 'ECCO_READPARMS: gencost_avgperiod not properly set'
44d3986245 Jean* 0865 endif
61aa767fd2 Gael* 0866
7a7bc10d23 Patr* 0867 if ( gencost_spmin (k ) .EQ. 0. _d 0 .AND.
69f8f4c84c Patr* 0868 & gencost_spmax (k ) .EQ. 0. _d 0 ) then
e7db56ba12 Gael* 0869 il =ilnblnk (gencost_name (k ))
0870 WRITE (msgBuf ,'(A,I2,2A)' )
0871 & 'ERROR in ECCO_READPARMS: for gencost' ,k ,
0872 & ' ' ,gencost_name (k )(1:il )
0873 CALL PRINT_MESSAGE ( msgBuf , errorMessageUnit ,
44d3986245 Jean* 0874 & SQUEEZE_RIGHT , myThid )
69f8f4c84c Patr* 0875 STOP
0876 & 'ECCO_READPARMS: gencost_spmin, gencost_spmax not set'
7a7bc10d23 Patr* 0877 endif
69f8f4c84c Patr* 0878
7a7bc10d23 Patr* 0879 if ( gencost_spzero (k ) .EQ. 9876. _d 0 ) then
e7db56ba12 Gael* 0880 il =ilnblnk (gencost_name (k ))
0881 WRITE (msgBuf ,'(A,I2,2A)' )
0882 & 'ERROR in ECCO_READPARMS: for gencost' ,k ,
0883 & ' ' ,gencost_name (k )(1:il )
0884 CALL PRINT_MESSAGE ( msgBuf , errorMessageUnit ,
44d3986245 Jean* 0885 & SQUEEZE_RIGHT , myThid )
69f8f4c84c Patr* 0886 STOP
0887 & 'ECCO_READPARMS: gencost_spzero not set'
7a7bc10d23 Patr* 0888 endif
69f8f4c84c Patr* 0889
877e0f385b An T* 0890
0891 elseif ( gencost_flag (k ) .eq. -1 ) then
0892
0893 elseif ( gencost_flag (k ).lt. -1) then
0894 if ( gencost_avgperiod (k ) .NE. 'day ' .AND.
0895 & gencost_avgperiod (k ) .NE. 'DAY ' .AND.
0896 & gencost_avgperiod (k ) .NE. 'month' .AND.
0897 & gencost_avgperiod (k ) .NE. 'MONTH' .AND.
0898 & gencost_avgperiod (k ) .NE. 'step' .AND.
0899 & gencost_avgperiod (k ) .NE. 'STEP' .AND.
0900 & gencost_avgperiod (k ) .NE. 'const' .AND.
0901 & gencost_avgperiod (k ) .NE. 'CONST' .AND.
0902 & gencost_avgperiod (k ) .NE. 'year ' .AND.
0903 & gencost_avgperiod (k ) .NE. 'YEAR ' ) then
0904
0905 il =ilnblnk (gencost_name (k ))
0906 WRITE (msgBuf ,'(A,I2,2A)' )
0907 & 'ERROR in ECCO_READPARMS: for gencost' ,k ,
0908 & ' ' ,gencost_name (k )(1:il )
0909 CALL PRINT_MESSAGE ( msgBuf , errorMessageUnit ,
44d3986245 Jean* 0910 & SQUEEZE_RIGHT , myThid )
877e0f385b An T* 0911 STOP
0912 & 'ECCO_READPARMS: gencost_avgperiod not properly set'
0913 endif
0914 endif
0915 endif
61aa767fd2 Gael* 0916
8c157ed454 Patr* 0917 end do
0918
0919 _END_MASTER ( myThid )
0920
0921 #endif /* ALLOW_GENCOST_CONTRIBUTION */
0922
0923 _BEGIN_MASTER ( myThid )
0924
7a77863887 Mart* 0925 #ifdef SINGLE_DISK_IO
0926 CLOSE (iUnit )
0927 #else
0928 CLOSE (iUnit ,STATUS ='DELETE' )
0929 #endif /* SINGLE_DISK_IO */
8c157ed454 Patr* 0930
0931 _END_MASTER ( myThid )
0932
8f7d13d0c9 Jean* 0933 WRITE (msgBuf ,'(A)' )
8c157ed454 Patr* 0934 & 'ECCO_READPARMS: done'
0935 CALL PRINT_MESSAGE ( msgBuf , standardMessageUnit ,
44d3986245 Jean* 0936 & SQUEEZE_RIGHT , myThid )
0c0a2531ee Patr* 0937
0938 _BARRIER
0939
f8e779c983 antn* 0940 RETURN
0941 END
