|
Last change
on this file since 3331 was
3331,
checked in by acozic, 7 years ago
|
|
Add modification for isotopes
|
-
Property svn:executable set to
*
|
|
File size:
877 bytes
|
| Line | |
|---|
| 1 | SUBROUTINE SURRTAB |
|---|
| 2 | |
|---|
| 3 | ! ----------------------------------------------------------------- |
|---|
| 4 | ! * E.C.M.W.F. PHYSICS PACKAGE ** AER'S RRTM LW RADIATION ** |
|---|
| 5 | |
|---|
| 6 | ! J.-J. MORCRETTE E.C.M.W.F. 98/07/14 |
|---|
| 7 | |
|---|
| 8 | ! ----------------------------------------------------------------- |
|---|
| 9 | |
|---|
| 10 | USE PARKIND1 ,ONLY : JPIM ,JPRB |
|---|
| 11 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
|---|
| 12 | |
|---|
| 13 | USE YOERRTAB , ONLY : TRANS, BPADE |
|---|
| 14 | |
|---|
| 15 | IMPLICIT NONE |
|---|
| 16 | |
|---|
| 17 | INTEGER(KIND=JPIM) :: ITR |
|---|
| 18 | |
|---|
| 19 | REAL(KIND=JPRB) :: ZTAU, ZTFN |
|---|
| 20 | REAL(KIND=JPRB) :: ZHOOK_HANDLE |
|---|
| 21 | |
|---|
| 22 | IF (LHOOK) CALL DR_HOOK('SURRTAB',0,ZHOOK_HANDLE) |
|---|
| 23 | BPADE=1.0_JPRB/0.278_JPRB |
|---|
| 24 | TRANS(0) =1.0_JPRB |
|---|
| 25 | TRANS(5000)=0.0_JPRB |
|---|
| 26 | DO ITR=1,4999 |
|---|
| 27 | ZTFN=REAL(ITR)/5000._JPRB |
|---|
| 28 | ZTAU=BPADE*ZTFN/(1.0_JPRB-ZTFN) |
|---|
| 29 | TRANS(ITR)=EXP(-ZTAU) |
|---|
| 30 | ENDDO |
|---|
| 31 | |
|---|
| 32 | ! ----------------------------------------------------------------- |
|---|
| 33 | |
|---|
| 34 | IF (LHOOK) CALL DR_HOOK('SURRTAB',1,ZHOOK_HANDLE) |
|---|
| 35 | END SUBROUTINE SURRTAB |
|---|
Note: See
TracBrowser
for help on using the repository browser.
