source: LMDZ6/branches/blowing_snow/libf/phylmd/rrtm/val923.F90

Last change on this file was 1990, checked in by Laurent Fairhead, 11 years ago

Corrections à la version r1989 pour permettre la compilation avec RRTM
Inclusion de la licence CeCILL_V2 pour RRTM


Changes to revision r1989 to enable RRTM code compilation
RRTM part put under CeCILL_V2 licence

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 3.3 KB
Line 
1SUBROUTINE VAL923(LDNEW)
2
3USE PARKIND1  ,ONLY : JPIM     ,JPRB
4USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
5
6USE YOMCLI   , ONLY : NTPMER   ,NTPGLA   ,NTPDES   ,NTPLAC   ,&
7 & SMASK    ,SMANQ    ,STHER    ,SALBN    ,SALBX    ,&
8 & SALBM    ,SALBG    ,SALBB    ,SALBD    ,SEMIN    ,&
9 & SEMIX    ,SEMIM    ,SEMIG    ,SEMIB    ,SEMID    ,&
10 & SDEPN    ,SDEPX    ,SDEPD    ,SARGN    ,SARGX    ,&
11 & SARGD    ,SSABN    ,SSABX    ,SSABD    ,SRSMN    ,&
12 & SRSMX    ,SRSMD    ,SZZ0N    ,SZZ0M    ,SZZ0B    ,&
13 & SZZ0U    ,SZZ0D 
14USE YOMLUN   , ONLY : NULOUT
15
16#ifdef DOC
17!**** *GEO923*
18
19!     PURPOSE.
20!     --------
21!      Compute the constants (YOMCLI) which are used by configuration 923.
22
23!     INTERFACE.
24!     ----------
25!      CALL VAL923(LDNEW)
26!        LDNEW = .FALSE. if old fields required
27!        Results in YOMCLI.
28
29!     AUTHORS.
30!     --------
31!      D. Giard 97-05-06
32
33!     MODIFICATIONS.
34!     --------------
35!      D. Giard   99-07-02 add values for desert and prints
36!      S. Kertesz 99-12-16 new type for lakes
37!        M.Hamrud      01-Oct-2003 CY28 Cleaning
38
39#endif
40
41IMPLICIT NONE
42
43LOGICAL           ,INTENT(IN)    :: LDNEW
44REAL(KIND=JPRB) :: ZHOOK_HANDLE
45
46!  Threshold defining the mask
47IF (LHOOK) CALL DR_HOOK('VAL923',0,ZHOOK_HANDLE)
48SMASK= 0.5_JPRB
49!  Value for missing data + 1
50SMANQ=-9998._JPRB
51!  Land-use types for sea, ice-cap, desert, lakes
52NTPMER= 1
53NTPGLA= 2
54NTPDES= 3
55NTPLAC= 5
56!  Roughness length : minimum, sea, sea-ice, urban areas, desert
57SZZ0N= 0.001_JPRB
58SZZ0M= 0.001_JPRB
59SZZ0B= 0.001_JPRB
60SZZ0U= 2.500_JPRB
61SZZ0D= 0.001_JPRB
62!  Ration of thermal to kinetic roughness length
63STHER= 0.10_JPRB
64!  Albedo : minimum, maximum, sea, ice-cap, sea-ice, desert
65IF (LDNEW) THEN
66  SALBN= 0.05_JPRB
67  SALBX= 0.80_JPRB
68ELSE
69  SALBN= 0.07_JPRB
70  SALBX= 0.70_JPRB
71ENDIF
72SALBM= 0.07_JPRB
73SALBG= 0.75_JPRB
74SALBB= 0.65_JPRB
75SALBD= 0.10_JPRB
76!  Emissivity : minimum, maximum, sea, ice-cap, sea-ice, desert
77SEMIN= 0.90_JPRB
78SEMIX= 1.00_JPRB
79SEMIM= 0.96_JPRB
80SEMIG= 0.98_JPRB
81SEMIB= 0.97_JPRB
82SEMID= 0.943_JPRB
83!  Soil depth : minimum, maximum, desert
84SDEPN= 0.10_JPRB
85SDEPX= 8.00_JPRB
86SDEPD= 0.10_JPRB
87!  Percentage of clay : minimum, maximum, desert
88SARGN=  3._JPRB
89SARGX= 58._JPRB
90SARGD=  3._JPRB
91!  Percentage of sand : minimum, maximum, desert
92SSABN=  6._JPRB
93SSABX= 92._JPRB
94SSABD= 92._JPRB
95!  Minimum surface resistance : minimum, maximum, desert
96SRSMX=5000._JPRB
97SRSMN=   1.0_JPRB
98SRSMD=5000._JPRB
99
100WRITE(UNIT=NULOUT,FMT=111) SMASK,SMANQ,STHER,&
101 & NTPMER,NTPGLA,NTPDES,NTPLAC 
102WRITE(UNIT=NULOUT,FMT=112) SZZ0N,SZZ0M,SZZ0B,SZZ0U,SZZ0D
103WRITE(UNIT=NULOUT,FMT=113) SALBN,SALBX,SALBM,SALBG,SALBB,SALBD,&
104 & SEMIN,SEMIX,SEMIM,SEMIG,SEMIB,SEMID 
105WRITE(UNIT=NULOUT,FMT=114) SDEPN,SDEPX,SDEPD,SARGN,SARGX,SARGD,&
106 & SSABN,SSABX,SSABD,SRSMN,SRSMX,SRSMD 
107111 FORMAT(' COMMON YOMCLI',/,&
108 & ' SMASK=',F4.2,' SMANQ=',F6.0,' STHER=',F4.2,/&
109 & ' NTPMER=',I2,' NTPGLA=',I2,' NTPDES=',I2,' NTPLAC=',I2) 
110112 FORMAT(' LONGUEUR DE RUGOSITE :',/,&
111 & ' minimum    mer    banquise  villes   desert ',&
112 & /,5F9.3) 
113113 FORMAT(' ALBEDO ET EMISSIVITE :',/,&
114 & ' minimum  maximum    mer    glacier  banquise  desert ',&
115 & 2(/,6F9.3)) 
116114 FORMAT(' PROFONDEUR, % ARGILE, % SABLE, RESIS. MIN. :',/,&
117 & ' minimum  maximum   desert ',&
118 & 4(/,3F9.3)) 
119
120IF (LHOOK) CALL DR_HOOK('VAL923',1,ZHOOK_HANDLE)
121END SUBROUTINE VAL923
Note: See TracBrowser for help on using the repository browser.