Changeset 9753 for utils/tools/NESTING/src/agrif_create_bathy.f90

Timestamp:

2018-06-06T15:28:47+02:00 (6 years ago)

Author:

jchanut

Message:

Modify arithmetic and median bathymetry average near land so that coastline is more realistic

File:

• utils/tools/NESTING/src/agrif_create_bathy.f90 (modified) (2 diffs)

utils/tools/NESTING/src/agrif_create_bathy.f90

@@ -335 +335 @@
                  WHERE( vardep(:,:) .GT. 0. )  mask_oce = 1
 
-                 IF( SUM(mask_oce) == 0 ) THEN
+!                 IF( SUM(mask_oce) == 0 ) THEN
+                 IF( SUM(mask_oce) < 0.5*(nxhr*nyhr) ) THEN
                     G1%bathy_meter(ji,jj) = 0.
                  ELSE
@@ -562 +563 @@
             END DO
         END DO
-
-        G1%gdepw_ps(:,:) = (1.-G1%wgt(:,:)) * G1%gdepw_ps(:,:) + gdepw_ps_interp(:,:)*G1%wgt(:,:)
+        IF (.NOT.identical_grids) THEN
+           G1%gdepw_ps(:,:) = (1.-G1%wgt(:,:)) * G1%gdepw_ps(:,:) + gdepw_ps_interp(:,:)*G1%wgt(:,:)
+        ENDIF
 
         G1%bathy_meter = G1%gdepw_ps