/[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.1 by jahn, Wed Apr 13 18:56:24 2011 UTC revision 1.2 by jahn, Thu Sep 22 19:18:37 2011 UTC
# Line 64  CEOP Line 64  CEOP
64  #endif  #endif
65       &           darwin_relaxscale,       &           darwin_relaxscale,
66       &           darwin_ForcingPeriod, darwin_ForcingCycle,       &           darwin_ForcingPeriod, darwin_ForcingCycle,
67         &           darwin_PARunits, darwin_W_to_uEins,
68  #ifdef ALLOW_PAR_DAY  #ifdef ALLOW_PAR_DAY
69       &           darwin_PARavPeriod,       &           darwin_PARavPeriod,
70  #endif  #endif
# Line 131  C Set defaults values for parameters in Line 132  C Set defaults values for parameters in
132        darwin_acdomFile=' '        darwin_acdomFile=' '
133        darwin_oasim_edFile=' '        darwin_oasim_edFile=' '
134        darwin_oasim_esFile=' '        darwin_oasim_esFile=' '
135          darwin_PARunits='Ein/m2/d        '
136          darwin_W_to_uEins=1. _d 0/0.2174 _d 0
137        darwin_relaxscale=0. _d 0        darwin_relaxscale=0. _d 0
138        darwin_seed=0        darwin_seed=0
139  c default periodic forcing to same as for GCHEM  c default periodic forcing to same as for GCHEM
# Line 211  C Open and read the data.darwin file Line 214  C Open and read the data.darwin file
214  #ifdef ALLOW_PAR_DAY  #ifdef ALLOW_PAR_DAY
215        darwin_PARnav = NINT(darwin_PARavPeriod*nsubtime/dTtracerLev(1))        darwin_PARnav = NINT(darwin_PARavPeriod*nsubtime/dTtracerLev(1))
216  #endif  #endif
217    C factor for conversion to uEin/m2/s
218          IF ( darwin_PARunits(1:16) .EQ. 'uEin/m2/s       ' ) THEN
219            darwin_PARFileConv = 1. _d 0
220          ELSEIF ( darwin_PARunits(1:16) .EQ. 'Ein/m2/d        ' ) THEN
221            darwin_PARFileConv = 1. _d 6/86400. _d 0
222          ELSEIF ( darwin_PARunits(1:16) .EQ. 'W/m2            ' ) THEN
223            darwin_PARFileConv = darwin_W_to_uEins
224          ELSE
225            WRITE(msgBuf,'(2A)') 'S/R DARWIN_READPARMS:',
226         &    'darwin_PARunits must be one of Ein/m2/d, uEin/m2/s, W/m2'
227            CALL PRINT_ERROR( msgBuf , 1)
228            STOP 'unknown darwin_PARunits'
229          ENDIF
230    
231  #ifdef DAR_DIAG_CHL  #ifdef DAR_DIAG_CHL
232  C default values  C default values
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22