Changeset 15745 for NEMO/branches/UKMO/tools_r4.0-HEAD_dev_DMP_TOOLS/DMP_TOOLS/src/make_dmp_file.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/make_dmp_file.F90 (modified) (5 diffs)

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

@@ -49 +49 @@
 
   ALLOCATE( resto(jpi, jpj) )
+  ALLOCATE( dcst(jpi, jpj) )
 
   !Create output file
@@ -57 +58 @@
   !Calculate surface and bottom damping coefficients
   zsdmp = 1._wp / ( pn_surf * rday )
-  zbdmp = 1._wp / ( pn_bot  * rday )
+  IF (pn_bot > 0.) THEN
+     zbdmp = 1._wp / ( pn_bot  * rday )
+  ELSE
+     zbdmp = 0.
+  ENDIF
+
+  ! Calculate distance from the coast
+  IF (ln_coast) CALL cofdis( dcst ) 
 
   !Loop through levels and read in tmask for each level as starting point for
@@ -74 +82 @@
            DO jj = 1, jpj
               DO ji = 1, jpi
-                 resto(ji,jj) = tmask(ji, jj) * (zbdmp + (zsdmp-zbdmp) * EXP(-gdept(ji,jj)/pn_dep))
+                 IF ( ln_exp ) THEN
+                    resto(ji,jj) = tmask(ji, jj) * (zbdmp + (zsdmp-zbdmp) * EXP(-gdept(ji,jj)/pn_dep))
+                 ELSE
+                    resto(ji,jj) = tmask(ji, jj) * (zbdmp + (zsdmp-zbdmp) / &
+                       &           (1._wp + EXP(rn_k * (gdept(ji,jj) - pn_dep) / 100._wp)))
+                 ENDIF 
               END DO
            END DO
@@ -94 +107 @@
            IF (ln_coast) THEN
               ! Reduce damping in vicinity of coastlines
-              CALL coast_dist_weight(resto)
+              CALL coast_dist_weight(resto, dcst)
            ENDIF
         ENDIF
@@ -109 +122 @@
         !Any user modifications can be added in the custom module
         IF ( ln_custom ) THEN
-              CALL custom_resto( resto )
+              WHERE (resto(:,:) <= 1.e-07) resto(:,:) = 0._wp
+              !CALL custom_resto( resto )
         ENDIF
      ENDIF