Skip to content

Commit b7f31dc

Browse files
coastwxcoastwx
authored
2024 Updates and Bug Fixes for the P-X LSM (#2023)
Pleim-Xiu LSM MODIS LCZ Compatibility & Surface Evaporation Update TYPE: bug fix & physics refinement KEYWORDS: MODIS, LCZ, P-X LSM, Latent Heat Flux SOURCE: Robert Gilliam & Jon Pleim, US EPA DESCRIPTION OF CHANGES: Problem: User indicated that the P-X LSM errored for MODIS LCZ 61 NUM_LAND_CAT configuration. Currently, we account for evaporation from transpiration, soil in both vegetated and non-veg parts, and wet leaves. But we only account for the latent heat effects on Tg from transpiration and evaporation from non-veg soil. This fix adds latent heat effects on Tg from soil in vegetated parts and from wet leaves. Fix for rare case where GRDFLX goes NaN because of a divide by zero based on a soil parameter when a water cell turns to sea ice. Solution: Logic checks in module_physics_init.F and module_sf_pxlsm.F were adjusted for 61 category inputs. P-X LSM data table, module_sf_pxlsm_data.F was updated for MODIS 61 categories. Default for LCZ 51-61 was set to MODIS urban class. We also added updates for the evaporation from vegetation and wet canopy. ISSUE: Fixes: #1965 LIST OF MODIFIED FILES: M phys/module_physics_init.F M phys/module_sf_pxlsm.F M phys/module_sf_pxlsm_data.F - Tested 61 class LCZ with PX LSM for a 1 day simulation with updated codes. Ran base MODIS 21 class scheme with same code before and after LCZ update. The results were identical after a 24 hour simulation. This confirms updates do not impact other MODIS settings in the P-X LSM. The MODIS 21 was by nature not identical to the MODIS 61, but similar enough and differences follow underlying differences in MODIS datasets. - The Jenkins tests are all passing. RELEASE NOTE: Pleim-Xiu LSM is now compatible with 61 category MODIS LCZ landuse dataset. A mode of latent heat effects on Tg from vegetated parts and from wet leaves is added to Pleim-Xiu LSM.
1 parent 9e265af commit b7f31dc

3 files changed

Lines changed: 133 additions & 44 deletions

File tree

‎phys/module_physics_init.F‎

Lines changed: 2 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -3520,7 +3520,8 @@ SUBROUTINE bl_init(STEPBL,BLDT,DT,RUBLTEN,RVBLTEN,RTHBLTEN, &
35203520
CASE (PXLSMSCHEME)
35213521
if(config_flags%num_land_cat .ne. 20 .and. config_flags%num_land_cat .ne. 21 .and. & ! MODIS
35223522
config_flags%num_land_cat .ne. 24 .and. config_flags%num_land_cat .ne. 28 .and. & ! USGS
3523-
config_flags%num_land_cat .ne. 40 .and. config_flags%num_land_cat .ne. 50 ) & ! NLCD
3523+
config_flags%num_land_cat .ne. 40 .and. config_flags%num_land_cat .ne. 50 .and. & ! NLCD
3524+
config_flags%num_land_cat .ne. 61 ) & ! MODIS W/LCZ
35243525
CALL wrf_error_fatal ( 'module_physics_init: PX LSM option requires USGS, MODIS, or NLCD' )
35253526
CALL LSMINIT(VEGFRA,SNOW,SNOWC,SNOWH,CANWAT,SMSTAV, &
35263527
SMSTOT, SFCRUNOFF,UDRUNOFF,ACSNOW, &

‎phys/module_sf_pxlsm.F‎

Lines changed: 43 additions & 30 deletions
Original file line numberDiff line numberDiff line change
@@ -370,6 +370,8 @@ SUBROUTINE pxlsm(U3D, V3D, DZ8W, QV3D,T3D,TH3D, RHO, &
370370
LAND_USE_TYPE = 'MODIS'
371371
ELSE IF (NLCAT == 21) THEN
372372
LAND_USE_TYPE = 'MODIS'
373+
ELSE IF (NLCAT == 61) THEN
374+
LAND_USE_TYPE = 'MODIS'
373375
ELSE IF (NLCAT == 24) THEN
374376
LAND_USE_TYPE = 'USGS'
375377
ELSE IF (NLCAT == 28) THEN
@@ -526,17 +528,19 @@ SUBROUTINE pxlsm(U3D, V3D, DZ8W, QV3D,T3D,TH3D, RHO, &
526528
! are converted to ice/snow for more reasonable treatment.
527529
IF( (XICE(I,J).GE.0.5) .OR. &
528530
(SST(I,J).LE.270.0.AND.XLAND(I,J).GE.1.50) ) THEN
529-
XLAND(I,J) = 1.0
530-
IFLAND = 1.0
531-
ZNT(I,J) = 0.001 ! Ice
532-
SMOIS(I,1,J) = 1.0 ! FWSAT
533-
SMOIS(I,2,J) = 1.0 ! FWSAT
534-
XICE(I,J) = 1.0
535-
ALBEDO(I,J) = 0.7
536-
SNOWC(I,J) = 1.0
537-
SNOW_FRA = 1.0
538-
VEGF_PX(I,J) = 0.0
539-
LAI(I,J) = 0.0
531+
XLAND(I,J) = 1.0
532+
IFLAND = 1.0
533+
ZNT(I,J) = 0.001
534+
SMOIS(I,1,J) = 1.0
535+
SMOIS(I,2,J) = 1.0
536+
XICE(I,J) = 1.0
537+
ALBEDO(I,J) = 0.7
538+
SNOWC(I,J) = 1.0
539+
SNOW_FRA = 1.0
540+
VEGF_PX(I,J) = 0.0
541+
LAI_PX(I,J) = 0.0
542+
LAI(I,J) = 0.0
543+
FCGSAT = 3.670
540544
ENDIF
541545
!-------------------------------------------------------------
542546

@@ -652,10 +656,15 @@ SUBROUTINE pxlsm(U3D, V3D, DZ8W, QV3D,T3D,TH3D, RHO, &
652656

653657
END DO ! Time internal PX time loop
654658

655-
IF (IFLAND .LT. 1.5) TSK(I,J) = TSLB(I,1,J) ! ATV 02/20: changed for compatibility with sst_skin = 1.
656-
CANWAT(I,J)= WR * 1000. ! convert WR back to mm for CANWAT
657-
RAW = RA(I,J) + 4.503 / USTAR
658-
QSFC(I,J) = QFX(I,J) * RAW / DENS1 + QV1
659+
IF (IFLAND .GE. 1.5) THEN
660+
TSK(I,J) = SST(I,J) ! Skin temp set to sea surface temperature for open water
661+
GRDFLX(I,J) = 0.0
662+
ELSE
663+
TSK(I,J) = TSLB(I,1,J) ! Skin temp set to 1 cm soil temperature in PX for now
664+
ENDIF
665+
CANWAT(I,J) = WR * 1000. ! convert WR back to mm for CANWAT
666+
RAW = RA(I,J) + 4.503 / USTAR
667+
QSFC(I,J) = QFX(I,J) * RAW / DENS1 + QV1
659668

660669
ENDDO ! END MIAN I LOOP
661670
ENDDO ! END MAIN J LOOP
@@ -1175,6 +1184,7 @@ SUBROUTINE SURFPX(DTPBL, IFLAND, ISNOW, NUDGEX, XICE1, SOLDN, GSW, & !in
11751184
REAL :: RSOIL, LDRY, DP ! Soil model updates - JEP 12/14
11761185
REAL :: C1MAX,ZZA,ZZB,ZDEL,ZLY,ZA,ZB,ZY2
11771186
REAL :: Rinc, Hcan ! JEP 2020
1187+
REAL :: CQ3BG,CQ3VW,CQ3VG,SIGG
11781188
11791189
!... Parameters
11801190
REAL :: ZOBS, GAMAH, BETAH, SIGF, BH, CT_SNOW, CT_IMPERV
@@ -1287,7 +1297,7 @@ SUBROUTINE SURFPX(DTPBL, IFLAND, ISNOW, NUDGEX, XICE1, SOLDN, GSW, & !in
12871297
VEGFRC, ISNOW, ISTI, IFLAND, LAI, BETAP, &
12881298
WG, W2, WR, &
12891299
RSTMIN, WWLT, WFC, RSOIL, RINC, &
1290-
EG, ER, ETR, CQ4, RS, FASS)
1300+
EG, ER, ETR, CQ4, RS, FASS, SIGG)
12911301
!--------------------------------------------------------------------
12921302
12931303
!--------------------------------------------------------------------
@@ -1330,18 +1340,17 @@ SUBROUTINE SURFPX(DTPBL, IFLAND, ISNOW, NUDGEX, XICE1, SOLDN, GSW, & !in
13301340
! IMPERVIOUS weighting scheme -- Subtract highly accurate impervious fraction from cell
13311341
! remainder is split between ground and vegetation. CT is a weighted fractional average.
13321342
! Snow CT is then applied for final heat capacity
1333-
IMF = AMAX1(0.0,IMPERV/100.0)
1334-
VEGF = (1.0 - IMF) * VEGFRC
1335-
SOILF= (1.0 - IMF) * (1.0 - VEGFRC)
1336-
CT = 1./( IMF/CT_IMPERV + VEGF/CV + SOILF/CG)
1337-
CT = 1./(SNOW_FRA/CT_SNOW + (1-SNOW_FRA)/CT)
1338-
CAPG = 1.0/CT
1339-
1343+
IMF = AMAX1(0.0,IMPERV/100.0)
1344+
VEGF = (1.0 - IMF) * VEGFRC
1345+
SOILF = (1.0 - IMF) * (1.0 - VEGFRC)
1346+
CT = 1./( IMF/CT_IMPERV + VEGF/CV + SOILF/CG)
1347+
CT = 1./(SNOW_FRA/CT_SNOW + (1-SNOW_FRA)/CT)
1348+
CAPG = 1.0/CT
13401349
SOILFLX = 2.0 * PI * TAUINV * (TG - T2)
13411350
GRDFLX = SOILFLX / CT
13421351
ENDIF
13431352
!-----------------------------------------------------------------------------------------
1344-
1353+
13451354
!--------------------------------------------------------------------
13461355
!-- ASSIMILATION --- COMPUTE SOIL MOISTURE NUDGING FROM TA2 and RH2
13471356
!-------COMPUTE ASSIMILATION COEFFICIENTS FOR ALL I
@@ -1374,7 +1383,10 @@ SUBROUTINE SURFPX(DTPBL, IFLAND, ISNOW, NUDGEX, XICE1, SOLDN, GSW, & !in
13741383
!-- Calculate the coefficients for implicit calculation of TG
13751384
CQ1 = (1.0 - 0.622 * LV * CRANKP / (r_d * TG)) * QSS
13761385
CQ2 = 0.622 * LV * QSS * CRANKP / (r_d * TG * TG)
1377-
CQ3 = DENS1 * (1.0 - VEGFRC) / (RAW + RSOIL)
1386+
CQ3BG = DENS1 * (1.0 - VEGFRC) / (RAW + RSOIL)
1387+
CQ3VW = DENS1 * VEGFRC*SIGG / RAW
1388+
CQ3VG = DENS1 * VEGFRC / (RAW + RSOIL + RINC)
1389+
CQ3 = CQ3BG + CQ3VW + CQ3VG
13781390
COEFFNP1 = 1.0 + DTPBL * CRANKP * (4.0 * EMISSI * STBOLT * TG ** 3 &
13791391
* CT + DENS1 * CPAIR / RAH * CPOT * CT + 2.0 * PI &
13801392
* TAUINV ) + DTPBL * (CT * LV * CQ2 * (CQ3 + CQ4))
@@ -1497,8 +1509,8 @@ END SUBROUTINE surfpx
14971509
SUBROUTINE QFLUX (DENS1, QV1, TA1, RG, RAW, QSS, & ! in
14981510
VEGFRC, ISNOW, ISTI, IFLAND, LAI, BETAP, & ! in
14991511
WG, W2, WR, & ! in
1500-
RSTMIN, WWLT, WFC, RSOIL, RINC, & ! in !Soil model updates - JEP 12/14
1501-
EG, ER, ETR, CQ4, RS, FASS) ! out
1512+
RSTMIN, WWLT, WFC, RSOIL, RINC, & ! in
1513+
EG, ER, ETR, CQ4, RS, FASS, SIGG ) ! out
15021514
15031515
!-------------------------------------------------------------------------
15041516
!
@@ -1535,9 +1547,10 @@ SUBROUTINE QFLUX (DENS1, QV1, TA1, RG, RAW, QSS, & ! in
15351547
!-- EG evaporation from ground (bare soil)
15361548
!-- ER evaporation from canopy
15371549
!-- ETR transpiration from vegetation
1538-
!-- CQ4
1550+
!-- CQ4 CQ4 is used for the implicit calculation of TG in SURFACE
15391551
!-- RS surface resistence
15401552
!-- FASS parameter for soil moisture nudging
1553+
!-- SIGG near ground-first layer moisture difference for evap or dew (1)
15411554
!-------------------------------------------------------------------------
15421555
!-------------------------------------------------------------------------
15431556
@@ -1554,12 +1567,12 @@ SUBROUTINE QFLUX (DENS1, QV1, TA1, RG, RAW, QSS, & ! in
15541567
REAL , INTENT(INOUT) :: BETAP, RSOIL
15551568
REAL, INTENT(IN) :: WWLT, WFC, RINC
15561569
1557-
REAL , INTENT(OUT) :: EG, ER, ETR, CQ4, RS, FASS
1570+
REAL , INTENT(OUT) :: EG, ER, ETR, CQ4, RS, FASS, SIGG
15581571
15591572
!... Local Variables
15601573
15611574
!... Real
1562-
REAL :: WRMAX, DELTA, SIGG, RADL, RADF, W2AVAIL, W2MXAV
1575+
REAL :: WRMAX, DELTA, RADL, RADF, W2AVAIL, W2MXAV
15631576
REAL :: FTOT, F1, F2, F3, F4
15641577
REAL :: FSHELT, GS, GA, FX
15651578
REAL :: PAR, F1MAX

‎phys/module_sf_pxlsm_data.F‎

Lines changed: 88 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -31,10 +31,13 @@ MODULE module_sf_pxlsm_data
3131
! 18 175. 30. 70. 50. 3.4 2.0 0.80 15. 45. wooded tundra
3232
! 19 120. 15. 40. 20. 2.4 1.0 0.40 15. 50. mixed tundra
3333
! 20 100. 10. 20. 5. 1.4 0.1 .015 25. 75. barren tundra
34+
! 21 9999. 0.1 00. 00. 0.0 0.0 0.01 8.0 08. inland lakes
35+
! 22-50 Unassigned
36+
! 51-61 150. 80. 5. 5. 2.0 0.5 0.04 11. 46. Urban LCZ 1-10 * Static urban settings *
3437
!------------------------------------------------------------------------------------
3538
!****************************************************************************************
3639
!****************************************************************************************
37-
REAL, DIMENSION(21), TARGET :: RSMIN_MODIS, Z00_MODIS, &
40+
REAL, DIMENSION(61), TARGET :: RSMIN_MODIS, Z00_MODIS, &
3841
VEG0_MODIS, VEGMN0_MODIS, &
3942
LAI0_MODIS, LAIMN0_MODIS, &
4043
SNUP0_MODIS, ALBF_MODIS, &
@@ -44,55 +47,127 @@ MODULE module_sf_pxlsm_data
4447
/ 175.0, 120.0, 175.0, 200.0, 200.0, &
4548
200.0, 200.0, 150.0, 120.0, 100.0, &
4649
200.0, 70.0, 150.0, 100.0, 9999.0, &
47-
100.0, 9999.0, 175.0, 120.0, 100.0, 9999.0 /
50+
100.0, 9999.0, 175.0, 120.0, 100.0, &
51+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
52+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
53+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
54+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
55+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
56+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
57+
150.0, 150.0, 150.0, 150.0, 150.0, &
58+
150.0, 150.0, 150.0, 150.0, 150.0, 150.0 /
4859

4960
DATA Z00_MODIS &
5061
/ 100.0, 90.0, 100.0, 100.0, 100.0, &
5162
15.0, 15.0, 25.0, 15.0, 7.0, &
5263
20.0, 10.0, 80.0, 30.0, 1.2, &
53-
5.0, 0.1, 30.0, 15.0, 10.0, 0.1 /
64+
5.0, 0.1, 30.0, 15.0, 10.0, &
65+
0.1, 9999.0, 9999.0, 9999.0, 9999.0, &
66+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
67+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
68+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
69+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
70+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
71+
80.0, 80.0, 80.0, 80.0, 80.0, &
72+
80.0, 80.0, 80.0, 80.0, 80.0, 80.0 /
5473

5574
DATA VEG0_MODIS &
5675
/ 93.0, 92.0, 60.0, 91.0, 92.0, &
5776
40.0, 20.0, 70.0, 70.0, 50.0, &
5877
65.0, 90.0, 5.0, 80.0, 0.1, &
59-
0.5, 0.0, 70.0, 40.0, 20.0, 0.0 /
78+
0.5, 0.0, 70.0, 40.0, 20.0, &
79+
0.0, 9999.0, 9999.0, 9999.0, 9999.0, &
80+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
81+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
82+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
83+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
84+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
85+
5.0, 5.0, 5.0, 5.0, 5.0, &
86+
5.0, 5.0, 5.0, 5.0, 5.0, 5.0 /
6087

6188
DATA VEGMN0_MODIS &
6289
/ 93.0, 92.0, 60.0, 91.0, 92.0, &
6390
20.0, 10.0, 60.0, 40.0, 20.0, &
6491
35.0, 20.0, 5.0, 40.0, 0.1, &
65-
0.5, 0.0, 50.0, 20.0, 5.0, 0.0 /
92+
0.5, 0.0, 50.0, 20.0, 5.0, &
93+
0.0, 9999.0, 9999.0, 9999.0, 9999.0, &
94+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
95+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
96+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
97+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
98+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
99+
5.0, 5.0, 5.0, 5.0, 5.0, &
100+
5.0, 5.0, 5.0, 5.0, 5.0, 5.0 /
66101

67102
DATA LAI0_MODIS &
68103
/ 5.5, 6.0, 3.0, 6.0, 5.5, &
69104
1.5, 1.5, 2.3, 1.5, 1.5, &
70105
2.5, 3.5, 2.0, 3.5, 0.1, &
71-
0.2, 0.0, 3.4, 2.4, 1.4, 0.0 /
106+
0.2, 0.0, 3.4, 2.4, 1.4, &
107+
0.0, 9999.0, 9999.0, 9999.0, 9999.0, &
108+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
109+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
110+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
111+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
112+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
113+
2.0, 2.0, 2.0, 2.0, 2.0, &
114+
2.0, 2.0, 2.0, 2.0, 2.0, 2.0 /
72115

73116
DATA LAIMN0_MODIS &
74117
/ 3.5, 3.5, 1.5, 2.0, 2.5, &
75118
1.0, 1.3, 2.0, 1.5, 1.5, &
76119
2.0, 1.5, 1.5, 1.5, 0.1, &
77-
0.1, 0.0, 2.0, 1.0, 0.1, 0.0 /
120+
0.1, 0.0, 2.0, 1.0, 0.1, &
121+
0.0, 9999.0, 9999.0, 9999.0, 9999.0, &
122+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
123+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
124+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
125+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
126+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
127+
1.5, 1.5, 1.5, 1.5, 1.5, &
128+
1.5, 1.5, 1.5, 1.5, 1.5, 1.5 /
78129

79130
DATA SNUP0_MODIS &
80-
/ 0.08, 0.08, 0.08, 0.08, 0.08, &
81-
0.03, 0.035, 0.03, 0.04, 0.04, &
82-
0.08, 0.04, 0.04, 0.04, 0.02, &
83-
0.02, 0.01, 0.80, 0.40, 0.015, 0.01 /
131+
/ 0.08, 0.08, 0.08, 0.08, 0.08, &
132+
0.03, 0.035, 0.03, 0.04, 0.04, &
133+
0.08, 0.04, 0.04, 0.04, 0.02, &
134+
0.02, 0.01, 0.80, 0.40, 0.015, &
135+
0.01, 9999.0, 9999.0, 9999.0, 9999.0, &
136+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
137+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
138+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
139+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
140+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
141+
0.04, 0.04, 0.04, 0.04, 0.04, &
142+
0.04, 0.04, 0.04, 0.04, 0.04, 0.04 /
84143

85144
DATA ALBF_MODIS &
86145
/ 12.0, 12.0, 14.0, 16.0, 13.0, &
87146
22.0, 20.0, 22.0, 20.0, 19.0, &
88147
17.0, 18.0, 11.0, 18.0, 60.0, &
89-
25.0, 8.0, 15.0, 15.0, 25.0, 8.0 /
148+
25.0, 8.0, 15.0, 15.0, 25.0, &
149+
8.0, 9999.0, 9999.0, 9999.0, 9999.0, &
150+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
151+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
152+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
153+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
154+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
155+
11.0, 11.0, 11.0, 11.0, 11.0, &
156+
11.0, 11.0, 11.0, 11.0, 11.0, 11.0 /
90157

91158
DATA SNOALB_MODIS &
92159
/ 30.0, 30.0, 30.0, 40.0, 35.0, &
93160
50.0, 60.0, 50.0, 50.0, 70.0, &
94161
50.0, 66.0, 46.0, 68.0, 82.0, &
95-
75.0, 8.0, 45.0, 55.0, 75.0, 8.0 /
162+
75.0, 8.0, 45.0, 55.0, 75.0, &
163+
8.0, 9999.0, 9999.0, 9999.0, 9999.0, &
164+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
165+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
166+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
167+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
168+
9999.0, 9999.0, 9999.0, 9999.0, 9999.0, &
169+
46.0, 46.0, 46.0, 46.0, 46.0, &
170+
46.0, 46.0, 46.0, 46.0, 46.0, 46.0 /
96171

97172
!****************************************************************************************
98173
!****************************************************************************************

0 commit comments

Comments
 (0)