Changeset 9166 for branches/2017/dev_merge_2017/NEMOGCM/TOOLS/NESTING/src/agrif_create_bathy.f90

Timestamp:

2017-12-22T16:45:51+01:00 (6 years ago)

Author:

jchanut

Message:

NESTING TOOLS:
Fixes to account for user defined number of ghostcells - still set to 1 to be consistent with NEMO
Set child grid bathymetry near boundaries to nearest neighbor interpolation from parent
Update Coarse grid bathymetry so that each cell volume matches child grid average

File:

• branches/2017/dev_merge_2017/NEMOGCM/TOOLS/NESTING/src/agrif_create_bathy.f90 (modified) (4 diffs)

branches/2017/dev_merge_2017/NEMOGCM/TOOLS/NESTING/src/agrif_create_bathy.f90

@@ -173 +173 @@
      CALL init_constant_bathy(G0%bathy_meter,bathy_fin_constant)
      !
-     boundary = connectionsize*irafx + nbghostcellsfine 
+     boundary = connectionsize*irafx + nbghostcellsfine + 1 
      !
      ! connection carried out by copying parent grid values for the fine points
@@ -459 +459 @@
         CALL Check_interp(G0,gdepw_ps_interp)      ! interpolation in connection zone (3 coarse grid cells)
         !
-        boundary = connectionsize*irafx + nbghostcellsfine                     
+        boundary = connectionsize*irafx + nbghostcellsfine + 1                     
         G1%gdepw_ps(1:boundary,:) = gdepw_ps_interp(1:boundary,:)
         G1%gdepw_ps(:,1:boundary) = gdepw_ps_interp(:,1:boundary)
@@ -531 +531 @@
            ! Smoothing to connect the connection zone (3 + nb_connection_pts coarse grid cells) and the interior domain
            !
-           boundary = (connectionsize+nb_connection_pts)*irafx + nbghostcellsfine 
+           boundary = (connectionsize+nb_connection_pts)*irafx + nbghostcellsfine + 1 
            CALL smooth_topo(G1%gdepw_ps(boundary:nxfin-boundary+1,boundary:nyfin-boundary+1),nbiter)
            G1%bathy_meter = G1%gdepw_ps                         
@@ -611 +611 @@
      CALL init_constant_bathy(G0%bathy_meter,bathy_fin_constant)
      !
-     boundary = connectionsize*irafx + nbghostcellsfine   
+     boundary = connectionsize*irafx + nbghostcellsfine + 1   
      !             
      G1%bathy_meter(1:boundary,:) = bathy_fin_constant(1:boundary,:)