/[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.9 by jahn, Wed Dec 4 21:21:32 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 115  CEOP Line 115  CEOP
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=' '
# Line 167  c default periodic forcing to same as fo Line 169  c default periodic forcing to same as fo
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      
# Line 189  c value >= 100 will be converted to inde Line 193  c value >= 100 will be converted to inde
193          darwin_bbphy(i) = 0 _d 0          darwin_bbphy(i) = 0 _d 0
194        enddo        enddo
195        darwin_radtrans_kmax = Nr        darwin_radtrans_kmax = Nr
196        darwin_radtrans_niter = 1        darwin_radtrans_niter = -2
197        darwin_part_size_P = 1 _d -15   ! mmol P per particle        darwin_part_size_P = 1 _d -15   ! mmol P per particle
198  #endif  #endif
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)
# Line 368  C--   Print a summary of parameter value Line 371  C--   Print a summary of parameter value
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
# Line 416  C--   Print a summary of parameter value Line 419  C--   Print a summary of parameter value
419       &'   /* deepest level in which to compute irradiances */')       &'   /* deepest level in which to compute irradiances */')
420        CALL WRITE_0D_I( darwin_radtrans_niter, INDEX_NONE,        CALL WRITE_0D_I( darwin_radtrans_niter, INDEX_NONE,
421       &   'darwin_radtrans_niter =',       &   'darwin_radtrans_niter =',
422       &'   /* number of "radtrans improvement" iterations */')       &'   /* iterations/method for solveing 3-stream equations */')
423  #endif /* DAR_RADTRANS */  #endif /* DAR_RADTRANS */
424  #endif /* WAVEBANDS */  #endif /* WAVEBANDS */
425    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22