/[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.2 by jahn, Thu Sep 22 19:18:37 2011 UTC revision 1.3 by jahn, Thu Sep 22 19:23:31 2011 UTC
# Line 65  CEOP Line 65  CEOP
65       &           darwin_relaxscale,       &           darwin_relaxscale,
66       &           darwin_ForcingPeriod, darwin_ForcingCycle,       &           darwin_ForcingPeriod, darwin_ForcingCycle,
67       &           darwin_PARunits, darwin_W_to_uEins,       &           darwin_PARunits, darwin_W_to_uEins,
68         &           darwin_ironUnits,
69         &           darwin_dustFrac,
70  #ifdef ALLOW_PAR_DAY  #ifdef ALLOW_PAR_DAY
71       &           darwin_PARavPeriod,       &           darwin_PARavPeriod,
72  #endif  #endif
# Line 134  C Set defaults values for parameters in Line 136  C Set defaults values for parameters in
136        darwin_oasim_esFile=' '        darwin_oasim_esFile=' '
137        darwin_PARunits='Ein/m2/d        '        darwin_PARunits='Ein/m2/d        '
138        darwin_W_to_uEins=1. _d 0/0.2174 _d 0        darwin_W_to_uEins=1. _d 0/0.2174 _d 0
139          darwin_ironUnits='mol Fe/m2/s     '
140          darwin_dustFrac=0.035 _d 0
141        darwin_relaxscale=0. _d 0        darwin_relaxscale=0. _d 0
142        darwin_seed=0        darwin_seed=0
143  c default periodic forcing to same as for GCHEM  c default periodic forcing to same as for GCHEM
# Line 228  C factor for conversion to uEin/m2/s Line 232  C factor for conversion to uEin/m2/s
232          STOP 'unknown darwin_PARunits'          STOP 'unknown darwin_PARunits'
233        ENDIF        ENDIF
234    
235    C factor for conversion to mmol Fe/m2/s (bioavailable fraction is done later)
236          IF ( darwin_ironUnits(1:16) .EQ. 'mol Fe/m2/s     ' ) THEN
237            darwin_ironFileConv = 1000. _d 0
238          ELSEIF ( darwin_ironUnits(1:16) .EQ. 'kg dust/m2/s    ' ) THEN
239            darwin_ironFileConv = 1000.*darwin_dustFrac*1000./55.845 _d 0
240          ELSE
241            WRITE(msgBuf,'(2A)') 'S/R DARWIN_READPARMS:',
242         &    'darwin_ironUnits must be one of mol Fe/m2/s, kg dust/m2/s'
243            CALL PRINT_ERROR( msgBuf , 1)
244            STOP 'unknown darwin_Feunits'
245          ENDIF
246    
247  #ifdef DAR_DIAG_CHL  #ifdef DAR_DIAG_CHL
248  C default values  C default values
249  C Geider:  chl:c = max(chl2cmin, chl2cmax/(1+(chl2cmax*alphachl*PARday)/(2*Pcm)))  C Geider:  chl:c = max(chl2cmin, chl2cmax/(1+(chl2cmax*alphachl*PARday)/(2*Pcm)))
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22