Updated for large domain case and neg interpolated landuse

noaa-oar-arl · Jul 12, 2022 · 7bc87a8 · 7bc87a8
1 parent aabc189
commit 7bc87a8
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 0 deletions.
diff --git a/parallel/src/rdfv3.f90 b/parallel/src/rdfv3.f90
@@ -1362,6 +1362,7 @@ END SUBROUTINE windrotation
     ENDIF
     call myinterp(dum2d,met_nx,met_ny,atmp,xindex,yindex,ncols_x,nrows_x,1)
     landuse(1:ncols_x,1:nrows_x) = atmp(1:ncols_x,1:nrows_x)
+    where(landuse.lt.0) landuse=0.0 !case of larger domains with fv3 water type/edge
     WRITE (*,f6100) 'vtype ', landuse(lprt_metx, lprt_mety), 'category'
   ELSE
     WRITE (*,f9400) TRIM(pname), 'vtype', TRIM(nf90_strerror(rcode))

diff --git a/serial/src/rdfv3.f90 b/serial/src/rdfv3.f90
@@ -1347,6 +1347,7 @@ END SUBROUTINE windrotation
     ENDIF
     call myinterp(dum2d,met_nx,met_ny,atmp,xindex,yindex,ncols_x,nrows_x,1)
     landuse(1:ncols_x,1:nrows_x) = atmp(1:ncols_x,1:nrows_x)
+    where(landuse.lt.0) landuse=0.0 !case of larger domains with fv3 water type/edge
     WRITE (*,f6100) 'vtype ', landuse(lprt_metx, lprt_mety), 'category'
   ELSE
     WRITE (*,f9400) TRIM(pname), 'vtype', TRIM(nf90_strerror(rcode))