id,summary,reporter,owner,description,type,status,priority,milestone,component,version,resolution,keywords,cc 247,allow for using WRF output with lakes,ahilboll,,"There is a hard-coded check for the value of `num_land_cat` in the WRF output, and FLEXPART-WRF would only recognize values 24 (for USGS) and 20 (for MODIS). However, when one uses lakes in WPS geogrid, the value of `num_land_cat` will be 28 or 21 for USGS and MODIS land use data, respectively. This patch modifies the check accordingly. Also, now a warning message is printed if the value of `num_land_cat` is not recognized. {{{ 1 file changed, 3 insertions(+), 2 deletions(-) readwind.f90 | 5 +++-- modified readwind.f90 @@ -1125,7 +1125,7 @@ enddo enddo enddo - if (num_land_cat.eq.24) then ! Assume USGS landuse data + if (num_land_cat.eq.24 .or. num_land_cat.eq.28) then ! Assume USGS landuse data do j = 0, nymin1 do i = 0, nxmin1 k = nint(landuse_wrf(i,j)) ! Safely convert element to integer (nearest) @@ -1174,7 +1174,7 @@ endselect enddo enddo - elseif (num_land_cat.eq.20) then ! Assume MODIS data + elseif (num_land_cat.eq.20 .or. num_land_cat.eq.21) then ! Assume MODIS data do j = 0, nymin1 do i = 0, nxmin1 k = nint(landuse_wrf(i,j)) ! Safely convert element to integer (nearest) @@ -1219,6 +1219,7 @@ enddo enddo endif + write(*,*) 'WARNING: unknown num_land_cat in WRF files:', num_land_cat endif ! lu_option.eq.1 ! snow depth }}} ",Defect,new,major,,FP other,FLEXPART-WRF,,,