/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/darwin_readparams.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/darwin/darwin_readparams.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.4 by jahn, Thu Aug 9 21:18:53 2012 UTC revision 1.8 by jahn, Tue Apr 16 20:21:57 2013 UTC
# Line 78  CEOP Line 78  CEOP
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
# Line 167  c default periodic forcing to same as fo Line 167  c default periodic forcing to same as fo
167    
168  #ifdef WAVEBANDS  #ifdef WAVEBANDS
169        DO ilam=1,tlam+1        DO ilam=1,tlam+1
170          darwin_wavebands(ilam) = 0          darwin_wavebands(ilam) = -1
171        ENDDO        ENDDO
172  #ifdef DAR_CALC_ACDOM  #ifdef DAR_CALC_ACDOM
173        darwin_Sdom = 0.014 _d 0        darwin_Sdom = 0.014 _d 0
174          darwin_lambda_aCDOM = 450 _d 0
175          darwin_aCDOM_fac = 0.2 _d 0
176  #endif  #endif
177  #ifdef DAR_DIAG_ACDOM  #ifdef DAR_DIAG_ACDOM
178  c value >= 100 will be converted to index in wavebands_init_fixed        c value >= 100 will be converted to index in wavebands_init_fixed      
# Line 189  c value >= 100 will be converted to inde Line 191  c value >= 100 will be converted to inde
191          darwin_bbphy(i) = 0 _d 0          darwin_bbphy(i) = 0 _d 0
192        enddo        enddo
193        darwin_radtrans_kmax = Nr        darwin_radtrans_kmax = Nr
194        darwin_radtrans_niter = 1        darwin_radtrans_niter = -2
195        darwin_part_size_P = 1 _d -15   ! mmol P per particle        darwin_part_size_P = 1 _d -15   ! mmol P per particle
196  #endif  #endif
197  #endif /* WAVEBANDS */  #endif /* WAVEBANDS */
# Line 368  C--   Print a summary of parameter value Line 370  C--   Print a summary of parameter value
370        CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid )        CALL PRINT_MESSAGE( msgBuf, iUnit, SQUEEZE_RIGHT, myThid )
371    
372  #if defined(WAVEBANDS) || defined(OASIM)  #if defined(WAVEBANDS) || defined(OASIM)
373        CALL WRITE_1D_I( darwin_waves, tlam, 0,        CALL WRITE_1D_RL( darwin_waves, tlam, 0,
374       &   'darwin_waves =',       &   'darwin_waves =',
375       &'   /* "central" wavelengths of wavebands */')       &'   /* "central" wavelengths of wavebands */')
376  #endif  #endif
# Line 416  C--   Print a summary of parameter value Line 418  C--   Print a summary of parameter value
418       &'   /* deepest level in which to compute irradiances */')       &'   /* deepest level in which to compute irradiances */')
419        CALL WRITE_0D_I( darwin_radtrans_niter, INDEX_NONE,        CALL WRITE_0D_I( darwin_radtrans_niter, INDEX_NONE,
420       &   'darwin_radtrans_niter =',       &   'darwin_radtrans_niter =',
421       &'   /* number of "radtrans improvement" iterations */')       &'   /* iterations/method for solveing 3-stream equations */')
422  #endif /* DAR_RADTRANS */  #endif /* DAR_RADTRANS */
423  #endif /* WAVEBANDS */  #endif /* WAVEBANDS */
424    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22