pkg/radtrans/radtrans_init_fixed.F

b55e95f1ff Oliv*0001 #include "RADTRANS_OPTIONS.h"
                0002 #include "EXF_OPTIONS.h"
                0003
                0004 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0005 CBOP
                0006 C     !ROUTINE: RADTRANS_INIT_FIXED
                0007
                0008 C     !INTERFACE:
                0009       SUBROUTINE RADTRANS_INIT_FIXED( myThid )
                0010
                0011 C     !DESCRIPTION:
                0012 C     Initialize radtrans constants for exf field loading
                0013
                0014 C     !USES:
                0015       IMPLICIT NONE
                0016 #include "SIZE.h"
                0017 #include "EEPARAMS.h"
                0018 #include "PARAMS.h"
                0019 #include "EXF_PARAM.h"
                0020 #include "RADTRANS_SIZE.h"
                0021 #include "RADTRANS_PARAMS.h"
                0022
                0023 C     !INPUT PARAMETERS:
                0024       INTEGER myThid
                0025 CEOP
                0026
                0027 C     !LOCAL VARIABLES:
                0028 C     === Local variables ===
                0029 C     msgBuf    :: Informational/error message buffer
                0030 C     errCount  :: error counter
                0031       CHARACTER*(MAX_LEN_MBUF) msgBuf
                0032       INTEGER errCount, l, oUnit
                0033       _RL planck, c, hc, oavo, hcoavo, rlamm
                0034
                0035 #ifdef ALLOW_RADTRANS
                0036
                0037 #ifdef ALLOW_DIAGNOSTICS
                0038       IF ( useDiagnostics ) THEN
                0039         CALL RADTRANS_DIAGNOSTICS_INIT( myThid )
                0040       ENDIF
                0041 #endif
                0042
                0043       _BEGIN_MASTER( myThid )
                0044
                0045       planck = 6.6256 _d -34   !Plancks constant J sec
                0046       c = 2.998 _d 8                 !speed of light m/sec
                0047       hc = 1.0/(planck*c)
                0048       oavo = 1.0/6.023 _d 23   ! 1/Avogadros number
                0049       hcoavo = hc*oavo
                0050       DO l = 1,nlam
                0051         rlamm = RT_wbRefWLs(l)*1 _d -9  !lambda in m
                0052         RT_WtouEins(l) = 1 _d 6*rlamm*hcoavo      !Watts to uEin/s conversion
                0053       ENDDO
                0054
                0055 #ifdef ALLOW_EXF
                0056       errCount = 0
                0057
                0058       IF ( RT_Edfile(1) .NE. ' ' ) THEN
                0059 # ifdef ALLOW_DEBUG
                0060        IF (debugMode) CALL DEBUG_CALL('GETFFIELD_START RT_E',myThid)
                0061 # endif
a720384a8b Oliv*0062        CALL EXF_GETFFIELD_START( useEXFYearlyFields,
b55e95f1ff Oliv*0063      I                           'radtrans', 'RT_E', RT_E_period,
                0064      I                           RT_E_startdate1, RT_E_startdate2,
                0065      U                           RT_E_StartTime, errCount,
                0066      I                           myThid )
                0067       ENDIF
                0068
                0069
                0070       IF ( RT_icefile .NE. ' ' ) THEN
                0071 # ifdef ALLOW_DEBUG
                0072        IF (debugMode) CALL DEBUG_CALL('GETFFIELD_START RT_ice',myThid)
                0073 # endif
a720384a8b Oliv*0074        CALL EXF_GETFFIELD_START( useEXFYearlyFields,
b55e95f1ff Oliv*0075      I                           'RT', 'RT_ice', RT_iceperiod,
                0076      I                           RT_icestartdate1, RT_icestartdate2,
                0077      U                           RT_iceStartTime, errCount,
                0078      I                           myThid )
                0079       ENDIF
                0080
                0081       IF ( errCount.GE.1 ) THEN
                0082        WRITE(msgBuf,'(A,I3,A)')
                0083      &     'RADTRANS_INIT_FIXED: detected', errCount,' fatal error(s)'
                0084        CALL PRINT_ERROR( msgBuf, myThid )
                0085        CALL ALL_PROC_DIE( 0 )
                0086        STOP 'ABNORMAL END: S/R RADTRANS_INIT_FIXED'
                0087       ENDIF
                0088
                0089 #endif /* ALLOW_EXF */
                0090
                0091       _END_MASTER( myThid )
                0092       _BARRIER
                0093
                0094 #endif /* ALLOW_RADTRANS */
                0095
                0096       RETURN
                0097       END