61 USE modd_data_cover_par
, ONLY : jpcover
80 CHARACTER(LEN=6),
DIMENSION(8) :: YDATA_TOWN
81 CHARACTER(LEN=6),
DIMENSION(9) :: YDATA_MATE
82 CHARACTER(LEN=6) :: YSTRING6
84 CHARACTER(LEN=200):: YFMT
87 REAL,
DIMENSION(JPCOVER) :: ZDATA_SVF_ROAD,ZDATA_SVF_WALL
88 REAL,
DIMENSION(SIZE(XDATA_HC_ROOF,1),SIZE(XDATA_HC_ROOF,2)) :: ZDATA_HC_ROOF
89 REAL(KIND=JPRB) :: ZHOOK_HANDLE
92 IF (
lhook)
CALL dr_hook(
'WRITE_COVER_TEX_TEB',0,zhook_handle)
112 WHERE (zdata_hc_roof.EQ.
xundef)
113 zdata_hc_roof=
xundef*1.e+6
115 WHERE (zdata_hc_road.EQ.
xundef)
116 zdata_hc_road=
xundef*1.e+6
118 WHERE (zdata_hc_wall.EQ.
xundef)
119 zdata_hc_wall=
xundef*1.e+6
128 IF (
clang==
'EN')
THEN 129 WRITE(
ntex,*)
'{\bf town parameters (1)}\' 131 WRITE(NTEX,*) '{\bf param\`etres de ville(1)}
133 WRITE(
ntex,*)
'\medskip\' 134 WRITE(NTEX,*) '\begin{tabular}{||r|l||c|c|c|c|c|c|c||}
' 135 WRITE(NTEX,*) '\
hline' 136 WRITE(NTEX,*) '\
hline' 137 WRITE(NTEX,*) '&&$z_{0_{town}}$&bld frac.&$h$&h/l&{h/w} $^\star$&
' 138 WRITE(NTEX,*) '$\alpha_{_{roof}}$&
' 139 WRITE(NTEX,*) '$\epsilon_{_{roof}}$
' 140 WRITE(NTEX,*) '\
hline' 141 WRITE(NTEX,*) '\
hline' 146 .AND..NE.
IF (XDATA_TOWN(I)>0. XDATA_Z0_TOWN(I)XUNDEF) THEN 148 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(XDATA_Z0_TOWN(I)),'.
',DEC(XDATA_Z0_TOWN(I)),'' 149 WRITE(YSTRING6, FMT=YFMT) XDATA_Z0_TOWN(I) 150 YDATA_TOWN(1) = YSTRING6 151 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(XDATA_BLD(I)),'.
',DEC(XDATA_BLD(I)),'' 152 WRITE(YSTRING6, FMT=YFMT) XDATA_BLD(I) 153 YDATA_TOWN(2) = YSTRING6 154 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(XDATA_BLD_HEIGHT(I)),'.
',DEC(XDATA_BLD_HEIGHT(I)),'' 155 WRITE(YSTRING6, FMT=YFMT) XDATA_BLD_HEIGHT(I) 156 YDATA_TOWN(3) = YSTRING6 157 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(XDATA_WALL_O_HOR(I)),'.
',DEC(XDATA_WALL_O_HOR(I)),'' 158 WRITE(YSTRING6, FMT=YFMT) XDATA_WALL_O_HOR(I) 159 YDATA_TOWN(4) = YSTRING6 160 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(XDATA_CAN_HW_RATIO(I)),'.
',DEC(XDATA_CAN_HW_RATIO(I)),'' 161 WRITE(YSTRING6, FMT=YFMT) XDATA_CAN_HW_RATIO(I) 162 YDATA_TOWN(5) = YSTRING6 163 WRITE(YSTRING6, FMT='(f3.2)
') XDATA_ALB_ROOF(I) 164 YDATA_TOWN(6) = YSTRING6 165 WRITE(YSTRING6, FMT='(f3.2)
') XDATA_EMIS_ROOF(I) 166 YDATA_TOWN(7) = YSTRING6 169 I,' &
',CNAME(I),' &
', YDATA_TOWN(1),' &
',YDATA_TOWN(2),' &
', & 170 YDATA_TOWN(3),' &
', YDATA_TOWN(4),' &
', & 171 YDATA_TOWN(5),' &
', YDATA_TOWN(6),' &
', & 174 WRITE(NTEX,*) '\
hline' 176 CALL HLINE(NTEX,GLINE,I) 179 WRITE(NTEX,*) '\end{tabular}
' 181 WRITE(NTEX,*) '\bigskip
' 190 IF (CLANG=='en
') THEN 191 WRITE(NTEX,*) '{\bf town parameters(2)}
193 WRITE(
ntex,*)
'{\bf param\`etres de ville (2)}\' 195 WRITE(NTEX,*) '\begin{tabular}{||r|l||c|c|c|c|c|c||}
' 196 WRITE(NTEX,*) '\
hline' 197 WRITE(NTEX,*) '\
hline' 198 WRITE(NTEX,*) '&&$\alpha_{_{road}}$&
' 199 WRITE(NTEX,*) '$\epsilon_{_{road}}$&$\psi_{s_{road}}$ $^\dagger$&
' 200 WRITE(NTEX,*) '$\alpha_{_{wall}}$&
' 201 WRITE(NTEX,*) '$\epsilon_{_{wall}}$&$\psi_{s_{wall}}$ $^\ddagger$
' 202 WRITE(NTEX,*) '\
hline' 203 WRITE(NTEX,*) '\
hline' 208 .AND..NE.
IF (XDATA_TOWN(I)>0. XDATA_Z0_TOWN(I)XUNDEF) THEN 210 WRITE(YSTRING6, FMT='(f3.2)
') XDATA_ALB_ROAD(I) 211 YDATA_TOWN(1) = YSTRING6 212 WRITE(YSTRING6, FMT='(f3.2)
') XDATA_EMIS_ROAD(I) 213 YDATA_TOWN(2) = YSTRING6 214 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(ZDATA_SVF_ROAD(I)),'.
',DEC(ZDATA_SVF_ROAD(I)),'' 215 WRITE(YSTRING6, FMT=YFMT) ZDATA_SVF_ROAD(I) 216 YDATA_TOWN(3) = YSTRING6 217 WRITE(YSTRING6, FMT='(f3.2)
') XDATA_ALB_WALL(I) 218 YDATA_TOWN(4) = YSTRING6 219 WRITE(YSTRING6, FMT='(f3.2)
') XDATA_EMIS_WALL(I) 220 YDATA_TOWN(5) = YSTRING6 221 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(ZDATA_SVF_WALL(I)),'.
',DEC(ZDATA_SVF_WALL(I)),'' 222 WRITE(YSTRING6, FMT='(f3.2)
') ZDATA_SVF_WALL(I) 223 YDATA_TOWN(6) = YSTRING6 226 I,' &
',CNAME(I),' &
', YDATA_TOWN(1),' &
',YDATA_TOWN(2),' &
', & 227 YDATA_TOWN(3),' &
', YDATA_TOWN(4),' &
', & 228 YDATA_TOWN(5),' &
', YDATA_TOWN(6),' ' 230 WRITE(NTEX,*) '\
hline' 232 CALL HLINE(NTEX,GLINE,I) 235 WRITE(NTEX,*) '\end{tabular}
' 237 WRITE(NTEX,*) '\bigskip
' 239 WRITE(NTEX,*) '\bigskip
' 244 !------------------------------------------------------------------------------- 250 IF (CLANG=='en
') THEN 251 WRITE(NTEX,*) '$^\star$ {h/w} is computed from the vegetation fraction
' 252 WRITE(NTEX,*) 'fraction of artificial area, and h/l
' 254 WRITE(NTEX,*) "$^\star$ {h/W} est calcul\'e \`a partir de la fraction de v\'eg\'etation," 255 WRITE(NTEX,*) 'de la fraction de batiments dans la surface construite
' 258 WRITE(NTEX,*) '$^\dagger$ $\psi_{s_{road}} = \sqrt{\frac{h^2}{w^2}+1}-
' 260 WRITE(NTEX,*) '$^\ddagger$ $\psi_{s_{wall}} = \frac{1}{2}\frac{w}{h}\left
' 264 WRITE(NTEX,*) '\bigskip
' 265 IF (CLANG=='en
') THEN 266 WRITE(NTEX,*) '{\bf roof thermal parameters }
' 268 WRITE(NTEX,*) '{\bf param\`etres thermiques des toits}
' 270 WRITE(NTEX,*) '\medskip
271 WRITE(
ntex,*)
'\begin{tabular}{||r|l||c|c|c|c|c|c|c|c|c||}' 272 WRITE(
ntex,*)
'\hline' 273 WRITE(
ntex,*)
'\hline' 274 WRITE(
ntex,*)
'&&$C_{_{roof}}^{(1)}$&$C_{_{roof}}^{(2)}$&$C_{_{roof}}^{(3)}$&' 275 WRITE(
ntex,*)
'$\lambda_{_{roof}}^{(1)}$&$\lambda_{_{roof}}^{(2)}$&$\lambda_{_{roof}}^{(3)}$&' 276 WRITE(
ntex,*)
'$D_{_{roof}}^{(1)}$&$D_{_{roof}}^{(2)}$&$D_{_{roof}}^{(3)}$\\' 277 WRITE(
ntex,*)
'\hline' 278 WRITE(
ntex,*)
'\hline' 287 WRITE(yfmt,
'(A2,I1,A1,I1,A1)')
'(F',
nb(zdata_hc_roof(i,1)*1.e-6),
'.'')' 288 WRITE(ystring6, fmt=yfmt) zdata_hc_roof(i,1)*1.e-6
289 ydata_mate(1) = ystring6
290 WRITE(yfmt,
'(A2,I1,A1,I1,A1)')
'(F',
nb(zdata_hc_roof(i,2)*1.e-6),
'.'')' 291 WRITE(ystring6, fmt=yfmt) zdata_hc_roof(i,2)*1.e-6
292 ydata_mate(2) = ystring6
293 WRITE(yfmt,
'(A2,I1,A1,I1,A1)')
'(F',
nb(zdata_hc_roof(i,3)*1.e-6),
'.'')' 294 WRITE(ystring6, fmt=yfmt) zdata_hc_roof(i,3)*1.e-6
295 ydata_mate(3) = ystring6
298 ydata_mate(4) = ystring6
301 ydata_mate(5) = ystring6
304 ydata_mate(6) = ystring6
307 ydata_mate(7) = ystring6
310 ydata_mate(8) = ystring6
313 ydata_mate(9) = ystring6
317 i,
' &',
cname(i),
' &', ydata_mate(1),
' &',ydata_mate(2),
' &', &
318 ydata_mate(3),
' &', ydata_mate(4),
' &', &
319 ydata_mate(5),
' &', ydata_mate(6),
' &', &
320 ydata_mate(7),
' &', ydata_mate(8),
' &',ydata_mate(9),
' \\' 321 WRITE(
ntex,*)
'\hline' 327 WRITE(
ntex,*)
'\end{tabular}' 331 WRITE(
ntex,*)
'\bigskip' 337 IF (
clang==
'EN')
THEN 338 WRITE(
ntex,*)
'{\bf road thermal parameters}\\' 340 WRITE(
ntex,*)
'{\bf param\`etres thermiques des rues}\\' 342 WRITE(
ntex,*)
'\medskip\' 343 WRITE(NTEX,*) '\begin{tabular}{||r|l||c|c|c|c|c|c|c|c|c||}
' 344 WRITE(NTEX,*) '\
hline' 345 WRITE(NTEX,*) '\
hline' 346 WRITE(NTEX,*) '&&$c_{_{road}}^{(1)}$&$c_{_{road}}^{(2)}$&$c_{_{road}}^
' 347 WRITE(NTEX,*) '$\lambda_{_{road}}^{(1)}$&$\lambda_{_{road}}^{(2)}$&$\lambda_
' 348 WRITE(NTEX,*) '$d_{_{road}}^{(1)}$&$d_{_{road}}^{(2)}$&$d_{_{road}}^{(
' 349 WRITE(NTEX,*) '\
hline' 350 WRITE(NTEX,*) '\
hline' 355 .AND..NE.
IF (XDATA_TOWN(I)>0. XDATA_Z0_TOWN(I)XUNDEF) THEN 357 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(ZDATA_HC_ROAD(I,1)*1.E-6),'',DEC(ZDATA_HC_ROAD(I,1)*1.E-6),'' 358 WRITE(YSTRING6, FMT=YFMT) ZDATA_HC_ROAD(I,1)*1.E-6 359 YDATA_MATE(1) = YSTRING6 360 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(ZDATA_HC_ROAD(I,2)*1.E-6),'',DEC(ZDATA_HC_ROAD(I,2)*1.E-6),'' 361 WRITE(YSTRING6, FMT=YFMT) ZDATA_HC_ROAD(I,2)*1.E-6 362 YDATA_MATE(2) = YSTRING6 363 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(ZDATA_HC_ROAD(I,3)*1.E-6),'',DEC(ZDATA_HC_ROAD(I,3)*1.E-6),'' 364 WRITE(YSTRING6, FMT=YFMT) ZDATA_HC_ROAD(I,3)*1.E-6 365 YDATA_MATE(3) = YSTRING6 366 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(XDATA_TC_ROAD(I,1)),'.
',DEC(XDATA_TC_ROAD(I,1)),'' 367 WRITE(YSTRING6, FMT=YFMT) XDATA_TC_ROAD(I,1) 368 YDATA_MATE(4) = YSTRING6 369 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(XDATA_TC_ROAD(I,2)),'.
',DEC(XDATA_TC_ROAD(I,2)),'' 370 WRITE(YSTRING6, FMT=YFMT) XDATA_TC_ROAD(I,2) 371 YDATA_MATE(5) = YSTRING6 372 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(XDATA_TC_ROAD(I,3)),'.
',DEC(XDATA_TC_ROAD(I,3)),'' 373 WRITE(YSTRING6, FMT=YFMT) XDATA_TC_ROAD(I,3) 374 YDATA_MATE(6) = YSTRING6 375 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(XDATA_D_ROAD(I,1)),'.
',DEC(XDATA_D_ROAD(I,1)),'' 376 WRITE(YSTRING6, FMT=YFMT) XDATA_D_ROAD(I,1) 377 YDATA_MATE(7) = YSTRING6 378 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(XDATA_D_ROAD(I,2)),'.
',DEC(XDATA_D_ROAD(I,2)),'' 379 WRITE(YSTRING6, FMT=YFMT) XDATA_D_ROAD(I,2) 380 YDATA_MATE(8) = YSTRING6 381 WRITE(YFMT,'(a2,i1,a1,i1,a1)
') '(f
',NB(XDATA_D_ROAD(I,3)),'.
',DEC(XDATA_D_ROAD(I,3)),'' 382 WRITE(YSTRING6, FMT=YFMT) XDATA_D_ROAD(I,3) 383 YDATA_MATE(9) = YSTRING6 387 I,' &
',CNAME(I),' &
', YDATA_MATE(1),' &
',YDATA_MATE(2),' &
', & 388 YDATA_MATE(3),' &
', YDATA_MATE(4),' &
', & 389 YDATA_MATE(5),' &
', YDATA_MATE(6),' &
', & 390 YDATA_MATE(7),' &
', YDATA_MATE(8),' &
',YDATA_MATE(9),' ' 391 WRITE(NTEX,*) '\
hline' 394 CALL HLINE(NTEX,GLINE,I) 397 WRITE(NTEX,*) '\end{tabular}
' 401 WRITE(NTEX,*) '\bigskip
' 407 IF (CLANG=='en
') THEN 408 WRITE(NTEX,*) '{\bf wall thermal parameters}
' 410 WRITE(NTEX,*) '{\bf param\`etres thermiques des murs}
' 412 WRITE(NTEX,*) '\medskip
413 WRITE(
ntex,*)
'\begin{tabular}{||r|l||c|c|c|c|c|c|c|c|c||}' 414 WRITE(
ntex,*)
'\hline' 415 WRITE(
ntex,*)
'\hline' 416 WRITE(
ntex,*)
'&&$C_{_{wall}}^{(1)}$&$C_{_{wall}}^{(2)}$&$C_{_{wall}}^{(3)}$&' 417 WRITE(
ntex,*)
'$\lambda_{_{wall}}^{(1)}$&$\lambda_{_{wall}}^{(2)}$&$\lambda_{_{wall}}^{(3)}$&' 418 WRITE(
ntex,*)
'$D_{_{wall}}^{(1)}$&$D_{_{wall}}^{(2)}$&$D_{_{wall}}^{(3)}$\\' 419 WRITE(
ntex,*)
'\hline' 420 WRITE(
ntex,*)
'\hline' 427 WRITE(yfmt,
'(A2,I1,A1,I1,A1)')
'(F',
nb(zdata_hc_wall(i,1)*1.e-6),
'.'')' 428 WRITE(ystring6, fmt=yfmt) zdata_hc_wall(i,1)*1.e-6
429 ydata_mate(1) = ystring6
430 WRITE(yfmt,
'(A2,I1,A1,I1,A1)')
'(F',
nb(zdata_hc_wall(i,2)*1.e-6),
'.'')' 431 WRITE(ystring6, fmt=yfmt) zdata_hc_wall(i,2)*1.e-6
432 ydata_mate(2) = ystring6
433 WRITE(yfmt,
'(A2,I1,A1,I1,A1)')
'(F',
nb(zdata_hc_wall(i,3)*1.e-6),
'.'')' 434 WRITE(ystring6, fmt=yfmt) zdata_hc_wall(i,3)*1.e-6
435 ydata_mate(3) = ystring6
438 ydata_mate(4) = ystring6
441 ydata_mate(5) = ystring6
444 ydata_mate(6) = ystring6
447 ydata_mate(7) = ystring6
450 ydata_mate(8) = ystring6
453 ydata_mate(9) = ystring6
457 i,
' &',
cname(i),
' &', ydata_mate(1),
' &',ydata_mate(2),
' &', &
458 ydata_mate(3),
' &', ydata_mate(4),
' &', &
459 ydata_mate(5),
' &', ydata_mate(6),
' &', &
460 ydata_mate(7),
' &', ydata_mate(8),
' &',ydata_mate(9),
' \\' 461 WRITE(
ntex,*)
'\hline' 467 WRITE(
ntex,*)
'\end{tabular}' 471 WRITE(
ntex,*)
'\clearpage' 474 IF (
lhook)
CALL dr_hook(
'WRITE_COVER_TEX_TEB',1,zhook_handle)
real, dimension(:,:), allocatable xdata_hc_roof
real, dimension(:), allocatable xdata_alb_road
real, dimension(:), allocatable xdata_emis_wall
real, dimension(:,:), allocatable xdata_d_roof
real, dimension(:), allocatable xdata_bld
subroutine hline(KTEX, GLINE, I)
subroutine write_cover_tex_teb
integer function nb(PX, KMAX)
real, dimension(:,:), allocatable xdata_tc_road
real, dimension(:), allocatable xdata_wall_o_hor
real, dimension(:,:), allocatable xdata_hc_wall
real, dimension(:), allocatable xdata_alb_wall
real, dimension(:,:), allocatable xdata_tc_roof
real, dimension(:,:), allocatable xdata_d_road
real, dimension(:), allocatable xdata_emis_road
integer function dec(PX, KMAX)
real, dimension(:), allocatable xdata_can_hw_ratio
real, dimension(:), allocatable xdata_bld_height
real, dimension(:,:), allocatable xdata_tc_wall
character(len=60), dimension(:), allocatable cname
real, dimension(:), allocatable xdata_z0_town
real, dimension(:), allocatable xdata_emis_roof
real, dimension(:), allocatable xdata_alb_roof
real, dimension(:,:), allocatable xdata_hc_road
real, dimension(:), allocatable xdata_town
real, dimension(:,:), allocatable xdata_d_wall
integer, parameter nlines