|
Last change
on this file since 5407 was
4773,
checked in by idelkadi, 12 months ago
|
- Update of Ecrad in LMDZ
The same organization of the Ecrad offline version is retained in order to facilitate the updating of Ecrad in LMDZ and the comparison between online and offline results.
version 1.6.1 of Ecrad (https://github.com/lguez/ecrad.git)
- Implementation of the double call of Ecrad in LMDZ
|
|
File size:
885 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, JPHOOK |
|---|
| 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=JPHOOK) :: 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.
