Changeset 10499 for NEMO/trunk/src/OCE/DIA

Timestamp:

2019-01-10T16:12:24+01:00 (6 years ago)

Author:

deazer

Message:

Fix ticket #2154

Location:

NEMO/trunk/src/OCE/DIA

Files:

• dia25h.F90 (modified) (2 diffs)
• diatmb.F90 (modified) (2 diffs)

NEMO/trunk/src/OCE/DIA/dia25h.F90

@@ -13 +13 @@
    USE in_out_manager  ! I/O units
    USE iom             ! I/0 library
+   USE wet_dry
 
    IMPLICIT NONE
@@ -211 +212 @@
          CALL iom_put( "salin25h", zw3d  )   ! salinity
          zw2d(:,:) = sshn_25h(:,:)*tmask(:,:,1) + zmdi*(1.0-tmask(:,:,1))
-         CALL iom_put( "ssh25h", zw2d )   ! sea surface 
+         IF( ll_wd ) THEN
+            CALL iom_put( "ssh25h", zw2d+ssh_ref )   ! sea surface 
+         ELSE
+            CALL iom_put( "ssh25h", zw2d )   ! sea surface
+         ENDIF
          ! Write velocities (instantaneous)
          zw3d(:,:,:) = un_25h(:,:,:)*umask(:,:,:) + zmdi*(1.0-umask(:,:,:))

NEMO/trunk/src/OCE/DIA/diatmb.F90

@@ -12 +12 @@
    USE in_out_manager  ! I/O units
    USE iom             ! I/0 library
+   USE wet_dry
 
    IMPLICIT NONE
@@ -108 +109 @@
       CALL dia_calctmb( tsn(:,:,:,jp_tem), zwtmb )
       !ssh already output but here we output it masked
-      CALL iom_put( "sshnmasked", sshn(:,:)*tmask(:,:,1) + zmdi*(1.0 - tmask(:,:,1)) )
+      IF( ll_wd ) THEN
+         CALL iom_put( "sshnmasked", (sshn(:,:)+ssh_ref)*tmask(:,:,1) + zmdi*(1.0 - tmask(:,:,1)) )
+      ELSE
+         CALL iom_put( "sshnmasked", sshn(:,:)*tmask(:,:,1) + zmdi*(1.0 - tmask(:,:,1)) )
+      ENDIF
+
       CALL iom_put( "top_temp"  , zwtmb(:,:,1) )    ! tmb Temperature
       CALL iom_put( "mid_temp"  , zwtmb(:,:,2) )    ! tmb Temperature