source: LMDZ6/trunk/libf/phylmd/ecrad/srtm_kgb26.F90 @ 4640

Last change on this file since 4640 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: 1.5 KB
RevLine 
[3908]1SUBROUTINE SRTM_KGB26
2
3!     Originally by J.Delamere, Atmospheric & Environmental Research.
4!     Revision: 2.4
5!     BAND 26:  22650-29000 cm-1 (low - nothing; high - nothing)
6!     Reformatted for F90 by JJMorcrette, ECMWF
7!     G.Mozdzynski March 2011 read constants from files
8
9!     ------------------------------------------------------------------
10
11USE PARKIND1  , ONLY : JPRB
12USE YOMHOOK   , ONLY : LHOOK, DR_HOOK
13USE YOESRTA26 , ONLY : SFLUXREF, RAYL
14
15!     ------------------------------------------------------------------
16
17IMPLICIT NONE
18
19! KURUCZ
20REAL(KIND=JPRB) :: ZHOOK_HANDLE
21IF (LHOOK) CALL DR_HOOK('SRTM_KGB26',0,ZHOOK_HANDLE)
22
23SFLUXREF = (/ &
24 !  &     129.462_JPRB, 15*_ZERO_ /)
25 & 29.0079_JPRB,  28.4088_JPRB,     20.3099_JPRB,  13.0283_JPRB &
26 & ,  11.8619_JPRB,  9.95840_JPRB,     6.68696_JPRB,  5.38987_JPRB &
27 & ,  3.49829_JPRB, 0.407693_JPRB,    0.299027_JPRB, 0.236827_JPRB &
28 & , 0.188502_JPRB, 0.163489_JPRB, 4.64335E-02_JPRB, 2.72662E-03_JPRB /) 
29
30!     Rayleigh extinction coefficient at all v
31RAYL = (/ &
32 & 1.21263E-06_JPRB,1.43428E-06_JPRB,1.67677E-06_JPRB,1.93255E-06_JPRB &
33 & , 2.19177E-06_JPRB,2.44195E-06_JPRB,2.66926E-06_JPRB,2.85990E-06_JPRB &
34 & , 3.00380E-06_JPRB,3.06996E-06_JPRB,3.08184E-06_JPRB,3.09172E-06_JPRB &
35 & , 3.09938E-06_JPRB,3.10456E-06_JPRB,3.10727E-06_JPRB,3.10818E-06_JPRB /) 
36
37!     ------------------------------------------------------------------
38IF (LHOOK) CALL DR_HOOK('SRTM_KGB26',1,ZHOOK_HANDLE)
39END SUBROUTINE SRTM_KGB26
Note: See TracBrowser for help on using the repository browser.