!*************************************************************************** SUBROUTINE RRTM_CMBGB9 !*************************************************************************** ! BAND 9: 1180-1390 cm-1 (low - H2O,CH4; high - CH4) !*************************************************************************** ! Parameters #include "tsmbkind.h" USE PARRRTM , ONLY : JPBAND ,JPG ,JPXSEC ,JPGPT USE YOERRTO9 , ONLY : KAO ,KBO ,SELFREFO ,FRACREFAO ,& &FRACREFBO, ABSN2OO USE YOERRTA9 , ONLY : KA ,KB ,SELFREF ,FRACREFA ,& & FRACREFB,ABSN2O ,& & ABSA ,ABSB ,NG9 USE YOERRTRWT, ONLY : FREFA ,FREFB ,FREFADF ,FREFBDF ,RWGT USE YOERRTFTR, ONLY : NGC ,NGS ,NGN ,NGB ,NGM , WT IMPLICIT NONE ! LOCAL INTEGER SCALARS INTEGER_M :: IGC, IPR, IPRSM, JN, JND, JNDC, JP, JT INTEGER_M :: MEQ, NEQ ! To force equivalence, HG, 13-DEC-2003 ! LOCAL REAL SCALARS REAL_B :: SUMF, SUMK DO JN = 1,11 DO JT = 1,5 DO JP = 1,13 IPRSM = 0 DO IGC = 1,NGC(9) SUMK = _ZERO_ DO IPR = 1, NGN(NGS(8)+IGC) IPRSM = IPRSM + 1 SUMK = SUMK + KAO(JN,JT,JP,IPRSM)*RWGT(IPRSM+128) ENDDO KA(JN,JT,JP,IGC) = SUMK ENDDO ENDDO ENDDO ENDDO DO JT = 1,5 DO JP = 13,59 IPRSM = 0 DO IGC = 1,NGC(9) SUMK = _ZERO_ DO IPR = 1, NGN(NGS(8)+IGC) IPRSM = IPRSM + 1 SUMK = SUMK + KBO(JT,JP,IPRSM)*RWGT(IPRSM+128) ENDDO KB(JT,JP,IGC) = SUMK ENDDO ENDDO ENDDO DO JT = 1,10 IPRSM = 0 DO IGC = 1,NGC(9) SUMK = _ZERO_ DO IPR = 1, NGN(NGS(8)+IGC) IPRSM = IPRSM + 1 SUMK = SUMK + SELFREFO(JT,IPRSM)*RWGT(IPRSM+128) ENDDO SELFREF(JT,IGC) = SUMK ENDDO ENDDO DO JN = 1,3 IPRSM = 0 DO IGC = 1,NGC(9) SUMK = _ZERO_ DO IPR = 1, NGN(NGS(8)+IGC) IPRSM = IPRSM + 1 JND = (JN-1)*16 SUMK = SUMK + ABSN2OO(JND+IPRSM)*RWGT(IPRSM+128) ENDDO JNDC = (JN-1)*NGC(9) ABSN2O(JNDC+IGC) = SUMK ENDDO ENDDO DO JP = 1,9 IPRSM = 0 DO IGC = 1,NGC(9) SUMF = _ZERO_ DO IPR = 1, NGN(NGS(8)+IGC) IPRSM = IPRSM + 1 SUMF = SUMF + FRACREFAO(IPRSM,JP) ENDDO FRACREFA(IGC,JP) = SUMF ENDDO ENDDO IPRSM = 0 DO IGC = 1,NGC(9) SUMF = _ZERO_ DO IPR = 1, NGN(NGS(8)+IGC) IPRSM = IPRSM + 1 SUMF = SUMF + FRACREFBO(IPRSM) ENDDO FRACREFB(IGC) = SUMF ENDDO DO JP = 1,9 DO IGC = 1,NGC(9) FREFA(NGS(8)+IGC,JP) = FRACREFA(IGC,JP) ENDDO ENDDO DO JP = 1,8 DO IGC = 1,NGC(9) FREFADF(NGS(8)+IGC,JP) = FRACREFA(IGC,JP+1) -FRACREFA(IGC,JP) ENDDO ENDDO DO IGC = 1,NGC(9) FREFB(NGS(8)+IGC,1) = FRACREFB(IGC) ENDDO ! +--Force the equivalence: BEGIN (HG, 13-DEC-2003) ! + ============================ ! +--ABSA ! + ^^^^ JN = 0 JT = 1 JP = 1 IGC = 1 DO NEQ=1,NG9 DO MEQ=1,715 JN = JN + 1 IF ( JN == 11 + 1) THEN JN = 1 JT = JT + 1 IF ( JT == 5 + 1 ) THEN JT = 1 JP = JP + 1 IF ( JP == 13 + 1 ) THEN JP = 1 IGC= IGC + 1 END IF END IF END IF ABSA(MEQ,NEQ) = KA(JN,JT,JP,IGC) ENDDO ENDDO ! +--ABSB ! + ^^^^ JT = 0 JP = 13 IGC = 1 DO NEQ=1,NG9 DO MEQ=1,235 JT = JT + 1 IF ( JT == 5 + 1 ) THEN JT = 1 JP = JP + 1 IF ( JP == 59 + 1 ) THEN JP = 13 IGC= IGC + 1 END IF END IF ABSB(MEQ,NEQ) = KB(JT,JP,IGC) ENDDO ENDDO ! +--Force the equivalence: END (HG, 13-DEC-2003) ! + ========================== RETURN END SUBROUTINE RRTM_CMBGB9