Changeset 15745 for NEMO/branches/UKMO/tools_r4.0-HEAD_dev_DMP_TOOLS/DMP_TOOLS/src/coast_dist.F90

Timestamp:

2022-03-08T18:55:08+01:00 (2 years ago)

Author:

dbruciaferri

Message:

modifications for restoring upper 300 m to EN4

File:

• NEMO/branches/UKMO/tools_r4.0-HEAD_dev_DMP_TOOLS/DMP_TOOLS/src/coast_dist.F90 (modified) (2 diffs)

NEMO/branches/UKMO/tools_r4.0-HEAD_dev_DMP_TOOLS/DMP_TOOLS/src/coast_dist.F90

@@ -9 +9 @@
    CONTAINS
 
-   SUBROUTINE coast_dist_weight( presto ) 
+   SUBROUTINE coast_dist_weight( presto, zdct ) 
       !!----------------------------------------------------------------------
       !!                 *** ROUTINE coast_dist_weight ***
@@ -20 +20 @@
       IMPLICIT NONE
       REAL(wp), DIMENSION(jpi,jpj), INTENT( inout ) :: presto
-      REAL(wp), DIMENSION(jpi,jpj) :: zdct
-      REAL(wp) :: zinfl = 1000.e3_wp  ! Distance of influence of coast line (could be
-                                  ! a namelist setting)
-      INTEGER :: jj, ji           ! dummy loop indices
+      REAL(wp), DIMENSION(jpi,jpj), INTENT( in    ) :: zdct
+      REAL(wp) :: zinfl = 200.e3_wp  ! Distance of influence of coast line (could be
+                                     ! a namelist setting)
+      REAL(wp) :: rdct               ! coastal distance factor
+      INTEGER  :: jj, ji             ! dummy loop indices
       
 
-      CALL cofdis( zdct )
+      zinfl = pn_cst * 1000._wp ! in m
+
+      !CALL cofdis( zdct )
       DO jj = 1, jpj
          DO ji = 1, jpi
-            zdct(ji,jj) = MIN( zinfl, zdct(ji,jj) )
-            presto(ji,jj) = presto(ji, jj) * 0.5_wp * (  1._wp - COS( rpi*zdct(ji,jj)/zinfl) )
+            !zdct(ji,jj) = MIN( zinfl, zdct(ji,jj) )
+            rdct = MIN( zinfl, zdct(ji,jj) ) / zinfl
+            !presto(ji,jj) = presto(ji, jj) * 0.5_wp * (  1._wp - COS( rpi*zdct(ji,jj)/zinfl) )
+            presto(ji,jj) = presto(ji, jj) * 0.5_wp * (  1._wp - COS( rpi*rdct) )
          END DO
       END DO