Changeset 8868 for branches/2017/dev_METO_2017/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

Timestamp:

2017-12-01T09:43:23+01:00 (6 years ago)

Author:

timgraham

Message:

Merged dev_r8789_sbc into branch

File:

• branches/2017/dev_METO_2017/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90 (modified) (4 diffs)

branches/2017/dev_METO_2017/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

@@ -20 +20 @@
    !!   dyn_asm_inc    : Apply the dynamic (u and v) increments
    !!   ssh_asm_inc    : Apply the SSH increment
+   !!   ssh_asm_div    : Apply divergence associated with SSH increment
    !!   seaice_asm_inc : Apply the seaice increment
    !!----------------------------------------------------------------------
@@ -48 +49 @@
    PUBLIC   dyn_asm_inc    !: Apply the dynamic (u and v) increments
    PUBLIC   ssh_asm_inc    !: Apply the SSH increment
+   PUBLIC   ssh_asm_div    !: Apply the SSH divergence
    PUBLIC   seaice_asm_inc !: Apply the seaice increment
 
@@ -768 +770 @@
             ENDIF
             !
+#if defined key_asminc
+         ELSE IF( kt == nitiaufin_r+1 ) THEN
+            !
+            ssh_iau(:,:) = 0._wp
+            !
+#endif
          ENDIF
          !                          !-----------------------------------------
@@ -792 +800 @@
    END SUBROUTINE ssh_asm_inc
 
+   SUBROUTINE ssh_asm_div( kt, phdivn )
+      !!----------------------------------------------------------------------
+      !!                  ***  ROUTINE ssh_asm_div  ***
+      !!
+      !! ** Purpose :   ssh increment with z* is incorporated via a correction of the local divergence          
+      !!                across all the water column
+      !!
+      !! ** Method  :
+      !!                CAUTION : sshiau is positive (inflow) decreasing the
+      !!                          divergence and expressed in m/s
+      !!
+      !! ** Action  :   phdivn   decreased by the ssh increment
+      !!----------------------------------------------------------------------
+      INTEGER, INTENT(IN) :: kt                               ! ocean time-step index
+      REAL(wp), DIMENSION(:,:,:), INTENT(inout) ::   phdivn   ! horizontal divergence
+      !!
+      INTEGER  ::   jk                                        ! dummy loop index
+      REAL(wp), DIMENSION(:,:)  , POINTER       ::   ztim     ! local array
+      !!----------------------------------------------------------------------
+      ! 
+#if defined key_asminc
+      CALL ssh_asm_inc( kt ) !==   (calculate increments)
+      !
+      IF( ln_linssh ) THEN 
+         phdivn(:,:,1) = phdivn(:,:,1) - ssh_iau(:,:) / e3t_n(:,:,1) * tmask(:,:,1)
+      ELSE 
+         CALL wrk_alloc( jpi,jpj, ztim)
+         ztim(:,:) = ssh_iau(:,:) / ( ht_n(:,:) + 1.0 - ssmask(:,:) )
+         DO jk = 1, jpkm1                                 
+            phdivn(:,:,jk) = phdivn(:,:,jk) - ztim(:,:) * tmask(:,:,jk) 
+         END DO
+         !
+         CALL wrk_dealloc( jpi,jpj, ztim)
+      ENDIF
+#endif
+      !
+   END SUBROUTINE ssh_asm_div
 
    SUBROUTINE seaice_asm_inc( kt, kindic )