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

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

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

revision 1.4 by jahn, Tue Apr 16 20:21:21 2013 UTC revision 1.5 by jahn, Tue Apr 16 20:21:57 2013 UTC
# Line 49  c local variables Line 49  c local variables
49        _RL ssf        _RL ssf
50  c      _RL planck, c, hc, oavo, hcoavo, rlamm  c      _RL planck, c, hc, oavo, hcoavo, rlamm
51  #ifdef DAR_CALC_ACDOM  #ifdef DAR_CALC_ACDOM
52        _RL rlam450,rlam        _RL rlam
53  #else  #else
54        _RL sacdom        _RL sacdom
55  #endif  #endif
# Line 278  c if no file given then CDOM is zero Line 278  c if no file given then CDOM is zero
278  #else /* DAR_CALC_ACDOM */  #else /* DAR_CALC_ACDOM */
279  c for 3-D or for direct comparison to RADTRANS would need the same formulation for CDOM as in radtrans.  c for 3-D or for direct comparison to RADTRANS would need the same formulation for CDOM as in radtrans.
280  c   CDOM absorption exponent  c   CDOM absorption exponent
281        rlam450 = 450.0 _d 0        nlaCDOM = 0
282        do ilam = 1,tlam        do ilam = 1,tlam
283         if (pwaves(ilam) .eq. 450) nl450 = ilam         if (pwaves(ilam) .eq. darwin_lambda_aCDOM) nlaCDOM = ilam
284         rlam = float(pwaves(ilam))         rlam = float(pwaves(ilam))
285         excdom(ilam) = exp(-darwin_Sdom*(rlam-rlam450))         excdom(ilam) = exp(-darwin_Sdom*(rlam-darwin_lambda_aCDOM))
286        enddo        enddo
287          if (nlaCDOM.eq.0) then
288            WRITE(msgBuf,'(A,I3,A)')
289         &      'WAVEBANDS_INIT_FIXED: no waveband found at ',
290         &      darwin_lambda_aCDOM, ' nm (needed for DAR_CALC_ACDOM).'
291            CALL PRINT_ERROR( msgBuf, myThid )
292            STOP 'ABNORMAL END: S/R WAVEBANDS_INIT_FIXED'
293          endif
294    
295        WRITE(msgBuf,'(A,1P1E20.12)')        WRITE(msgBuf,'(A,1P1E20.12)')
296       &      'WAVEBANDS_INIT_FIXED: darwin_aCDOM_fac = ',darwin_aCDOM_fac       &      'WAVEBANDS_INIT_FIXED: darwin_aCDOM_fac = ',darwin_aCDOM_fac
# Line 293  c   CDOM absorption exponent Line 300  c   CDOM absorption exponent
300       &      'WAVEBANDS_INIT_FIXED: darwin_Sdom = ', darwin_Sdom       &      'WAVEBANDS_INIT_FIXED: darwin_Sdom = ', darwin_Sdom
301        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
302       &                    SQUEEZE_RIGHT, 1 )       &                    SQUEEZE_RIGHT, 1 )
303        WRITE(msgBuf,'(A,I3)')        WRITE(msgBuf,'(A,I3,A,I4)')
304       &      'WAVEBANDS_INIT_FIXED: nl450 = ', nl450       &      'WAVEBANDS_INIT_FIXED: nlaCDOM = ', nlaCDOM, ', lambda = ',
305         &      pwaves(nlaCDOM)
306        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
307       &                    SQUEEZE_RIGHT, 1 )       &                    SQUEEZE_RIGHT, 1 )
308  #endif /* DAR_CALC_ACDOM */  #endif /* DAR_CALC_ACDOM */
# Line 308  c     find waveband index for acdom diag Line 316  c     find waveband index for acdom diag
316            goto 60            goto 60
317          endif          endif
318         enddo         enddo
319           WRITE(msgBuf,'(2A,I3,A)') 'WAVEBANDS_INIT_FIXED: ',
320         & 'darwin_diag_acdom_ilam =',darwin_diag_acdom_ilam,
321         & ' not found in darwin_waves'
322           CALL PRINT_ERROR( msgBuf, myThid )
323           STOP 'ABNORMAL END: S/R WAVEBANDS_INIT_FIXED'
324  60     continue  60     continue
325        endif        endif
326    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22