|
Last change
on this file since 5225 was
3908,
checked in by idelkadi, 3 years ago
|
|
Online implementation of the radiative transfer code ECRAD in the LMDZ model.
- Inclusion of the ecrad directory containing the sources of the ECRAD code
- interface routine : radiation_scheme.F90
- Adaptation of compilation scripts :
- compilation under CPP key CPP_ECRAD
- compilation with option "-rad ecard" or "-ecard true"
- The "-rad old/rtm/ecran" build option will need to replace the "-rrtm true" and "-ecrad true" options in the future.
- Runing LMDZ simulations with ecrad, you need :
- logical key iflag_rrtm = 2 in physiq.def
- namelist_ecrad (DefLists?)
- the directory "data" containing the configuration files is temporarily placed in ../libfphylmd/ecrad/
- Compilation and execution are tested in the 1D case. The repository under svn would allow to continue the implementation work: tests, verification of the results, ...
|
|
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.
