| 78 |
#ifdef WAVEBANDS |
#ifdef WAVEBANDS |
| 79 |
& ,darwin_wavebands |
& ,darwin_wavebands |
| 80 |
#ifdef DAR_CALC_ACDOM |
#ifdef DAR_CALC_ACDOM |
| 81 |
& ,darwin_Sdom |
& ,darwin_Sdom, darwin_lambda_aCDOM, darwin_aCDOM_fac |
| 82 |
#endif |
#endif |
| 83 |
#ifdef DAR_DIAG_ACDOM |
#ifdef DAR_DIAG_ACDOM |
| 84 |
& ,darwin_diag_acdom_ilam |
& ,darwin_diag_acdom_ilam |
| 115 |
& dic_pCO2, dic_int1, dic_int2, dic_int3, dic_int4 |
& dic_pCO2, dic_int1, dic_int2, dic_int3, dic_int4 |
| 116 |
#endif |
#endif |
| 117 |
|
|
| 118 |
|
_BEGIN_MASTER(myThid) |
| 119 |
|
|
| 120 |
C Set defaults values for parameters in DARWIN_IO.h |
C Set defaults values for parameters in DARWIN_IO.h |
| 121 |
darwin_iceFile=' ' |
darwin_iceFile=' ' |
| 122 |
darwin_ironFile=' ' |
darwin_ironFile=' ' |
| 169 |
|
|
| 170 |
#ifdef WAVEBANDS |
#ifdef WAVEBANDS |
| 171 |
DO ilam=1,tlam+1 |
DO ilam=1,tlam+1 |
| 172 |
darwin_wavebands(ilam) = 0 |
darwin_wavebands(ilam) = -1 |
| 173 |
ENDDO |
ENDDO |
| 174 |
#ifdef DAR_CALC_ACDOM |
#ifdef DAR_CALC_ACDOM |
| 175 |
darwin_Sdom = 0.014 _d 0 |
darwin_Sdom = 0.014 _d 0 |
| 176 |
|
darwin_lambda_aCDOM = 450 _d 0 |
| 177 |
|
darwin_aCDOM_fac = 0.2 _d 0 |
| 178 |
#endif |
#endif |
| 179 |
#ifdef DAR_DIAG_ACDOM |
#ifdef DAR_DIAG_ACDOM |
| 180 |
c value >= 100 will be converted to index in wavebands_init_fixed |
c value >= 100 will be converted to index in wavebands_init_fixed |
| 199 |
#endif /* WAVEBANDS */ |
#endif /* WAVEBANDS */ |
| 200 |
|
|
| 201 |
C Open and read the data.darwin file |
C Open and read the data.darwin file |
|
_BEGIN_MASTER(myThid) |
|
| 202 |
WRITE(msgBuf,'(A)') ' DARWIN_READPARMS: opening data.darwin' |
WRITE(msgBuf,'(A)') ' DARWIN_READPARMS: opening data.darwin' |
| 203 |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |
CALL PRINT_MESSAGE(msgBuf, standardMessageUnit, |
| 204 |
& SQUEEZE_RIGHT , 1) |
& SQUEEZE_RIGHT , 1) |
| 371 |
CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid ) |
CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid ) |
| 372 |
|
|
| 373 |
#if defined(WAVEBANDS) || defined(OASIM) |
#if defined(WAVEBANDS) || defined(OASIM) |
| 374 |
CALL WRITE_1D_I( darwin_waves, tlam, 0, |
CALL WRITE_1D_RL( darwin_waves, tlam, 0, |
| 375 |
& 'darwin_waves =', |
& 'darwin_waves =', |
| 376 |
&' /* "central" wavelengths of wavebands */') |
&' /* "central" wavelengths of wavebands */') |
| 377 |
#endif |
#endif |