source: LMDZ6/trunk/libf/phylmd/ecrad/surrtab.F90 @ 4629

Last change on this file since 4629 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 
1SUBROUTINE 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
10USE PARKIND1  ,ONLY : JPIM     ,JPRB
11USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
12
13USE YOERRTAB , ONLY : TRANS, BPADE
14
15IMPLICIT NONE
16
17INTEGER(KIND=JPIM) :: ITR
18
19REAL(KIND=JPRB) :: ZTAU, ZTFN
20REAL(KIND=JPRB) :: ZHOOK_HANDLE
21
22IF (LHOOK) CALL DR_HOOK('SURRTAB',0,ZHOOK_HANDLE)
23BPADE=1.0_JPRB/0.278_JPRB
24TRANS(0)   =1.0_JPRB
25TRANS(5000)=0.0_JPRB
26DO ITR=1,4999
27  ZTFN=REAL(ITR)/5000._JPRB
28  ZTAU=BPADE*ZTFN/(1.0_JPRB-ZTFN)
29  TRANS(ITR)=EXP(-ZTAU)
30ENDDO
31
32!     -----------------------------------------------------------------
33
34IF (LHOOK) CALL DR_HOOK('SURRTAB',1,ZHOOK_HANDLE)
35END SUBROUTINE SURRTAB
Note: See TracBrowser for help on using the repository browser.