1 |
|
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 |
✓✗ |
1 |
IF (LHOOK) CALL DR_HOOK('VAL923',0,ZHOOK_HANDLE) |
48 |
|
1 |
SMASK= 0.5_JPRB |
49 |
|
|
! Value for missing data + 1 |
50 |
|
1 |
SMANQ=-9998._JPRB |
51 |
|
|
! Land-use types for sea, ice-cap, desert, lakes |
52 |
|
1 |
NTPMER= 1 |
53 |
|
1 |
NTPGLA= 2 |
54 |
|
1 |
NTPDES= 3 |
55 |
|
1 |
NTPLAC= 5 |
56 |
|
|
! Roughness length : minimum, sea, sea-ice, urban areas, desert |
57 |
|
1 |
SZZ0N= 0.001_JPRB |
58 |
|
1 |
SZZ0M= 0.001_JPRB |
59 |
|
1 |
SZZ0B= 0.001_JPRB |
60 |
|
1 |
SZZ0U= 2.500_JPRB |
61 |
|
1 |
SZZ0D= 0.001_JPRB |
62 |
|
|
! Ration of thermal to kinetic roughness length |
63 |
|
1 |
STHER= 0.10_JPRB |
64 |
|
|
! Albedo : minimum, maximum, sea, ice-cap, sea-ice, desert |
65 |
✗✓ |
1 |
IF (LDNEW) THEN |
66 |
|
|
SALBN= 0.05_JPRB |
67 |
|
|
SALBX= 0.80_JPRB |
68 |
|
|
ELSE |
69 |
|
1 |
SALBN= 0.07_JPRB |
70 |
|
1 |
SALBX= 0.70_JPRB |
71 |
|
|
ENDIF |
72 |
|
1 |
SALBM= 0.07_JPRB |
73 |
|
1 |
SALBG= 0.75_JPRB |
74 |
|
1 |
SALBB= 0.65_JPRB |
75 |
|
1 |
SALBD= 0.10_JPRB |
76 |
|
|
! Emissivity : minimum, maximum, sea, ice-cap, sea-ice, desert |
77 |
|
1 |
SEMIN= 0.90_JPRB |
78 |
|
1 |
SEMIX= 1.00_JPRB |
79 |
|
1 |
SEMIM= 0.96_JPRB |
80 |
|
1 |
SEMIG= 0.98_JPRB |
81 |
|
1 |
SEMIB= 0.97_JPRB |
82 |
|
1 |
SEMID= 0.943_JPRB |
83 |
|
|
! Soil depth : minimum, maximum, desert |
84 |
|
1 |
SDEPN= 0.10_JPRB |
85 |
|
1 |
SDEPX= 8.00_JPRB |
86 |
|
1 |
SDEPD= 0.10_JPRB |
87 |
|
|
! Percentage of clay : minimum, maximum, desert |
88 |
|
1 |
SARGN= 3._JPRB |
89 |
|
1 |
SARGX= 58._JPRB |
90 |
|
1 |
SARGD= 3._JPRB |
91 |
|
|
! Percentage of sand : minimum, maximum, desert |
92 |
|
1 |
SSABN= 6._JPRB |
93 |
|
1 |
SSABX= 92._JPRB |
94 |
|
1 |
SSABD= 92._JPRB |
95 |
|
|
! Minimum surface resistance : minimum, maximum, desert |
96 |
|
1 |
SRSMX=5000._JPRB |
97 |
|
1 |
SRSMN= 1.0_JPRB |
98 |
|
1 |
SRSMD=5000._JPRB |
99 |
|
|
|
100 |
|
1 |
WRITE(UNIT=NULOUT,FMT=111) SMASK,SMANQ,STHER,& |
101 |
|
2 |
& NTPMER,NTPGLA,NTPDES,NTPLAC |
102 |
|
1 |
WRITE(UNIT=NULOUT,FMT=112) SZZ0N,SZZ0M,SZZ0B,SZZ0U,SZZ0D |
103 |
|
1 |
WRITE(UNIT=NULOUT,FMT=113) SALBN,SALBX,SALBM,SALBG,SALBB,SALBD,& |
104 |
|
2 |
& SEMIN,SEMIX,SEMIM,SEMIG,SEMIB,SEMID |
105 |
|
1 |
WRITE(UNIT=NULOUT,FMT=114) SDEPN,SDEPX,SDEPD,SARGN,SARGX,SARGD,& |
106 |
|
2 |
& 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 |
✓✗ |
1 |
IF (LHOOK) CALL DR_HOOK('VAL923',1,ZHOOK_HANDLE) |
121 |
|
1 |
END SUBROUTINE VAL923 |