!*************************************************************************** SUBROUTINE RRTM_CMBGB1 !*************************************************************************** ! The subroutines CMBGB1->CMBGB16 input the absorption coefficient ! data for each band, which are defined for 16 g-points and 16 spectral ! bands. The data are combined with appropriate weighting following the ! g-point mapping arrays specified in RRTMINIT. Plank fraction data ! in arrays FRACREFA and FRACREFB are combined without weighting. All ! g-point reduced data are put into new arrays for use in RRTM. ! BAND 1: 10-250 cm-1 (low - H2O; high - H2O) !*************************************************************************** ! Parameters USE PARKIND1 ,ONLY : JPIM ,JPRB USE YOMHOOK ,ONLY : LHOOK, DR_HOOK USE YOERRTO1 , ONLY : KAO, KBO, SELFREFO, FORREFO, FRACREFAO,FRACREFBO USE YOERRTA1 , ONLY : KA , KB , SELFREF , FORREF , FRACREFA ,FRACREFB USE YOERRTRWT, ONLY : FREFA ,FREFB ,RWGT USE YOERRTFTR, ONLY : NGC ,NGN IMPLICIT NONE INTEGER(KIND=JPIM) :: IGC, IPR, IPRSM, JP, JT REAL(KIND=JPRB) :: Z_SUMF1, Z_SUMF2, Z_SUMK REAL(KIND=JPRB) :: ZHOOK_HANDLE IF (LHOOK) CALL DR_HOOK('RRTM_CMBGB1',0,ZHOOK_HANDLE) DO JT = 1,5 DO JP = 1,13 IPRSM = 0 DO IGC = 1,NGC(1) Z_SUMK = 0.0_JPRB DO IPR = 1, NGN(IGC) IPRSM = IPRSM + 1 Z_SUMK = Z_SUMK + KAO(JT,JP,IPRSM)*RWGT(IPRSM) ENDDO KA(JT,JP,IGC) = Z_SUMK ENDDO ENDDO DO JP = 13,59 IPRSM = 0 DO IGC = 1,NGC(1) Z_SUMK = 0.0_JPRB DO IPR = 1, NGN(IGC) IPRSM = IPRSM + 1 Z_SUMK = Z_SUMK + KBO(JT,JP,IPRSM)*RWGT(IPRSM) ENDDO KB(JT,JP,IGC) = Z_SUMK ENDDO ENDDO ENDDO DO JT = 1,10 IPRSM = 0 DO IGC = 1,NGC(1) Z_SUMK = 0.0_JPRB DO IPR = 1, NGN(IGC) IPRSM = IPRSM + 1 Z_SUMK = Z_SUMK + SELFREFO(JT,IPRSM)*RWGT(IPRSM) ENDDO SELFREF(JT,IGC) = Z_SUMK ENDDO ENDDO IPRSM = 0 DO IGC = 1,NGC(1) Z_SUMK = 0.0_JPRB Z_SUMF1 = 0.0_JPRB Z_SUMF2 = 0.0_JPRB DO IPR = 1, NGN(IGC) IPRSM = IPRSM + 1 Z_SUMK = Z_SUMK + FORREFO(IPRSM)*RWGT(IPRSM) Z_SUMF1= Z_SUMF1+ FRACREFAO(IPRSM) Z_SUMF2= Z_SUMF2+ FRACREFBO(IPRSM) ENDDO FORREF(IGC) = Z_SUMK FRACREFA(IGC) = Z_SUMF1 FRACREFB(IGC) = Z_SUMF2 ENDDO DO IGC = 1,NGC(1) FREFA(IGC,1) = FRACREFA(IGC) FREFB(IGC,1) = FRACREFB(IGC) ENDDO IF (LHOOK) CALL DR_HOOK('RRTM_CMBGB1',1,ZHOOK_HANDLE) END SUBROUTINE RRTM_CMBGB1