From 7bc87a865b3497c498a08b99da1bac070b8900fa Mon Sep 17 00:00:00 2001 From: pcampbe8 Date: Tue, 12 Jul 2022 16:44:09 -0400 Subject: [PATCH] Updated for large domain case and neg interpolated landuse --- parallel/src/rdfv3.f90 | 1 + serial/src/rdfv3.f90 | 1 + 2 files changed, 2 insertions(+) diff --git a/parallel/src/rdfv3.f90 b/parallel/src/rdfv3.f90 index f1649b6..64636db 100644 --- 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 index b55b7bd..f6c5b01 100644 --- 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))