221 |
c set si half sat |
c set si half sat |
222 |
val_ksatsi=1. _d 0 |
val_ksatsi=1. _d 0 |
223 |
c set nutrient ratios for phyto |
c set nutrient ratios for phyto |
224 |
|
#ifdef ONLY_P_CYCLE |
225 |
|
val_R_SiP_diatom=0.0 _d 0 |
226 |
|
val_R_NP=0.0 _d 0 |
227 |
|
val_RFeP=0.0 _d -3 |
228 |
|
val_R_NP_diaz=0.0 _d 0 |
229 |
|
val_RFeP_diaz=0.0 _d 0 * val_RFeP |
230 |
|
val_R_PC=0.0 _d 0 |
231 |
|
val_R_PICPOC=0.0 _d 0 |
232 |
|
#else |
233 |
val_R_SiP_diatom=16.0 _d 0 ! 32 for Fanny's runs |
val_R_SiP_diatom=16.0 _d 0 ! 32 for Fanny's runs |
234 |
val_R_NP=16.0 _d 0 |
val_R_NP=16.0 _d 0 |
235 |
val_RFeP=1.0 _d -3 |
val_RFeP=1.0 _d -3 |
237 |
val_RFeP_diaz=30.0 _d 0 * val_RFeP |
val_RFeP_diaz=30.0 _d 0 * val_RFeP |
238 |
val_R_PC=120.0 _d 0 |
val_R_PC=120.0 _d 0 |
239 |
val_R_PICPOC=0.8 _d 0 |
val_R_PICPOC=0.8 _d 0 |
240 |
|
#endif |
241 |
|
c |
242 |
#ifdef OLD_GRAZE |
#ifdef OLD_GRAZE |
243 |
c grazing hlaf saturation |
c grazing hlaf saturation |
244 |
kgrazesat = 0.1 _d 0 |
kgrazesat = 0.1 _d 0 |
334 |
|
|
335 |
C make sure we have reserved enough space in Ptracers |
C make sure we have reserved enough space in Ptracers |
336 |
IF ( nCompZooMax .LT. 4 ) THEN |
IF ( nCompZooMax .LT. 4 ) THEN |
337 |
|
#ifndef ONLY_P_CYCLE |
338 |
WRITE(msgBuf,'(A,A,I3)') |
WRITE(msgBuf,'(A,A,I3)') |
339 |
& 'MONOD_INIT_FIXED: ERROR: 4 zooplankton components, but ', |
& 'MONOD_INIT_FIXED: ERROR: 4 zooplankton components, but ', |
340 |
& 'nCompZooMax = ', nCompZooMax |
& 'nCompZooMax = ', nCompZooMax |
341 |
CALL PRINT_ERROR( msgBuf , 1) |
CALL PRINT_ERROR( msgBuf , 1) |
342 |
STOP 'ABNORMAL END: S/R MONOD_INIT_FIXED' |
STOP 'ABNORMAL END: S/R MONOD_INIT_FIXED' |
343 |
|
#endif |
344 |
ENDIF |
ENDIF |
345 |
DO nz = 1,nzmax |
DO nz = 1,nzmax |
346 |
iZooP (nz) = iZoo + (nz-1)*strideTypeZoo |
iZooP (nz) = iZoo + (nz-1)*strideTypeZoo |