source: LMDZ6/trunk/libf/phylmd/ecrad.v1.5.1/srtm_cmbgb26.F90 @ 5451

Last change on this file since 5451 was 3908, checked in by idelkadi, 4 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.1 KB
Line 
1SUBROUTINE SRTM_CMBGB26
2
3!     BAND 26:  22650-29000 cm-1 (low - nothing; high - nothing)
4!-----------------------------------------------------------------------
5
6USE PARKIND1  ,ONLY : JPIM , JPRB
7USE YOMHOOK   ,ONLY : LHOOK, DR_HOOK
8
9USE YOESRTM  , ONLY : NGN
10USE YOESRTWN , ONLY : NGC, NGS, RWGT
11!USE YOESRTWN , ONLY : NGC, NGS, NGN, RWGT
12USE YOESRTA26, ONLY : SFLUXREF, RAYL, &
13                    & SFLUXREFC, RAYLC
14
15IMPLICIT NONE
16
17! Local variables
18INTEGER(KIND=JPIM) :: IGC, IPR, IPRSM
19REAL(KIND=JPRB)    :: ZSUMF1, ZSUMF2
20
21REAL(KIND=JPRB) :: ZHOOK_HANDLE
22!     ------------------------------------------------------------------
23IF (LHOOK) CALL DR_HOOK('SRTM_CMBGB26',0,ZHOOK_HANDLE)
24
25IPRSM = 0
26DO IGC = 1,NGC(11)
27  ZSUMF1 = 0.
28  ZSUMF2 = 0.
29  DO IPR = 1, NGN(NGS(10)+IGC)
30    IPRSM = IPRSM + 1
31    ZSUMF1 = ZSUMF1 + RAYL(IPRSM)*RWGT(IPRSM+160)
32    ZSUMF2 = ZSUMF2 + SFLUXREF(IPRSM)
33  ENDDO
34  RAYLC(IGC) = ZSUMF1
35  SFLUXREFC(IGC) = ZSUMF2
36ENDDO
37
38!     -----------------------------------------------------------------
39IF (LHOOK) CALL DR_HOOK('SRTM_CMBGB26',1,ZHOOK_HANDLE)
40END SUBROUTINE SRTM_CMBGB26
41
Note: See TracBrowser for help on using the repository browser.