1 | SUBROUTINE VAL923(LDNEW) |
---|
2 | |
---|
3 | USE PARKIND1 ,ONLY : JPIM ,JPRB |
---|
4 | USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
---|
5 | |
---|
6 | USE 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 |
---|
14 | USE 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 | |
---|
41 | IMPLICIT NONE |
---|
42 | |
---|
43 | LOGICAL ,INTENT(IN) :: LDNEW |
---|
44 | REAL(KIND=JPRB) :: ZHOOK_HANDLE |
---|
45 | |
---|
46 | ! Threshold defining the mask |
---|
47 | IF (LHOOK) CALL DR_HOOK('VAL923',0,ZHOOK_HANDLE) |
---|
48 | SMASK= 0.5_JPRB |
---|
49 | ! Value for missing data + 1 |
---|
50 | SMANQ=-9998._JPRB |
---|
51 | ! Land-use types for sea, ice-cap, desert, lakes |
---|
52 | NTPMER= 1 |
---|
53 | NTPGLA= 2 |
---|
54 | NTPDES= 3 |
---|
55 | NTPLAC= 5 |
---|
56 | ! Roughness length : minimum, sea, sea-ice, urban areas, desert |
---|
57 | SZZ0N= 0.001_JPRB |
---|
58 | SZZ0M= 0.001_JPRB |
---|
59 | SZZ0B= 0.001_JPRB |
---|
60 | SZZ0U= 2.500_JPRB |
---|
61 | SZZ0D= 0.001_JPRB |
---|
62 | ! Ration of thermal to kinetic roughness length |
---|
63 | STHER= 0.10_JPRB |
---|
64 | ! Albedo : minimum, maximum, sea, ice-cap, sea-ice, desert |
---|
65 | IF (LDNEW) THEN |
---|
66 | SALBN= 0.05_JPRB |
---|
67 | SALBX= 0.80_JPRB |
---|
68 | ELSE |
---|
69 | SALBN= 0.07_JPRB |
---|
70 | SALBX= 0.70_JPRB |
---|
71 | ENDIF |
---|
72 | SALBM= 0.07_JPRB |
---|
73 | SALBG= 0.75_JPRB |
---|
74 | SALBB= 0.65_JPRB |
---|
75 | SALBD= 0.10_JPRB |
---|
76 | ! Emissivity : minimum, maximum, sea, ice-cap, sea-ice, desert |
---|
77 | SEMIN= 0.90_JPRB |
---|
78 | SEMIX= 1.00_JPRB |
---|
79 | SEMIM= 0.96_JPRB |
---|
80 | SEMIG= 0.98_JPRB |
---|
81 | SEMIB= 0.97_JPRB |
---|
82 | SEMID= 0.943_JPRB |
---|
83 | ! Soil depth : minimum, maximum, desert |
---|
84 | SDEPN= 0.10_JPRB |
---|
85 | SDEPX= 8.00_JPRB |
---|
86 | SDEPD= 0.10_JPRB |
---|
87 | ! Percentage of clay : minimum, maximum, desert |
---|
88 | SARGN= 3._JPRB |
---|
89 | SARGX= 58._JPRB |
---|
90 | SARGD= 3._JPRB |
---|
91 | ! Percentage of sand : minimum, maximum, desert |
---|
92 | SSABN= 6._JPRB |
---|
93 | SSABX= 92._JPRB |
---|
94 | SSABD= 92._JPRB |
---|
95 | ! Minimum surface resistance : minimum, maximum, desert |
---|
96 | SRSMX=5000._JPRB |
---|
97 | SRSMN= 1.0_JPRB |
---|
98 | SRSMD=5000._JPRB |
---|
99 | |
---|
100 | WRITE(UNIT=NULOUT,FMT=111) SMASK,SMANQ,STHER,& |
---|
101 | & NTPMER,NTPGLA,NTPDES,NTPLAC |
---|
102 | WRITE(UNIT=NULOUT,FMT=112) SZZ0N,SZZ0M,SZZ0B,SZZ0U,SZZ0D |
---|
103 | WRITE(UNIT=NULOUT,FMT=113) SALBN,SALBX,SALBM,SALBG,SALBB,SALBD,& |
---|
104 | & SEMIN,SEMIX,SEMIM,SEMIG,SEMIB,SEMID |
---|
105 | WRITE(UNIT=NULOUT,FMT=114) SDEPN,SDEPX,SDEPD,SARGN,SARGX,SARGD,& |
---|
106 | & SSABN,SSABX,SSABD,SRSMN,SRSMX,SRSMD |
---|
107 | 111 FORMAT(' COMMON YOMCLI',/,& |
---|
108 | & ' SMASK=',F4.2,' SMANQ=',F6.0,' STHER=',F4.2,/& |
---|
109 | & ' NTPMER=',I2,' NTPGLA=',I2,' NTPDES=',I2,' NTPLAC=',I2) |
---|
110 | 112 FORMAT(' LONGUEUR DE RUGOSITE :',/,& |
---|
111 | & ' minimum mer banquise villes desert ',& |
---|
112 | & /,5F9.3) |
---|
113 | 113 FORMAT(' ALBEDO ET EMISSIVITE :',/,& |
---|
114 | & ' minimum maximum mer glacier banquise desert ',& |
---|
115 | & 2(/,6F9.3)) |
---|
116 | 114 FORMAT(' PROFONDEUR, % ARGILE, % SABLE, RESIS. MIN. :',/,& |
---|
117 | & ' minimum maximum desert ',& |
---|
118 | & 4(/,3F9.3)) |
---|
119 | |
---|
120 | IF (LHOOK) CALL DR_HOOK('VAL923',1,ZHOOK_HANDLE) |
---|
121 | END SUBROUTINE VAL923 |
---|