New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
isfdynatf.F90 in NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/ISF – NEMO

source: NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/ISF/isfdynatf.F90 @ 13248

Last change on this file since 13248 was 13248, checked in by francesca, 4 years ago

dev_r12558_HPC-08_epico_Extra_Halo: merge with trunk@13237, see #2366

File size: 4.3 KB
Line 
1MODULE isfdynatf
2   !!=========================================================================
3   !!                       ***  MODULE  isfnxt  ***
4   !! Ice shelf update: compute the dynatf ice shelf contribution
5   !!=========================================================================
6   !! History :  OPA  !  2019-09  (P. Mathiot)  Original code
7   !!-------------------------------------------------------------------------
8 
9   !!-------------------------------------------------------------------------
10   !!   isfnxt       : apply correction needed for the ice shelf to ensure conservation
11   !!-------------------------------------------------------------------------
12
13   USE isf_oce
14
15   USE phycst , ONLY: r1_rho0         ! physical constant
16   USE dom_oce                        ! time and space domain
17   USE oce, ONLY : ssh                ! sea-surface height !!st needed for substitution
18
19   USE in_out_manager
20
21   IMPLICIT NONE
22
23   PRIVATE
24
25   PUBLIC isf_dynatf
26   !! * Substitutions
27#  include "do_loop_substitute.h90"
28#  include "domzgr_substitute.h90"
29
30CONTAINS
31
32   SUBROUTINE isf_dynatf ( kt, Kmm, pe3t_f, pcoef )
33      !!--------------------------------------------------------------------
34      !!                  ***  ROUTINE isf_dynatf  ***
35      !!
36      !! ** Purpose : compute the ice shelf volume filter correction for cavity, param, ice sheet coupling case
37      !!
38      !!-------------------------- OUT -------------------------------------
39      INTEGER                         , INTENT(in   ) :: kt       ! ocean time step
40      INTEGER                         , INTENT(in   ) :: Kmm      ! ocean time level index
41      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) :: pe3t_f   ! time filtered scale factor to be corrected
42      !
43      REAL(wp)                        , INTENT(in   ) :: pcoef    ! rn_atfp * rn_Dt * r1_rho0
44      !!--------------------------------------------------------------------
45      INTEGER :: jk  ! loop index
46      !!--------------------------------------------------------------------
47      !
48      ! ice shelf cavity
49      IF ( ln_isfcav_mlt ) CALL isf_dynatf_mlt(Kmm, pe3t_f, misfkt_cav, misfkb_cav, rhisf_tbl_cav, rfrac_tbl_cav, fwfisf_cav, fwfisf_cav_b, pcoef)
50      !
51      ! ice shelf parametrised
52      IF ( ln_isfpar_mlt ) CALL isf_dynatf_mlt(Kmm, pe3t_f, misfkt_par, misfkb_par, rhisf_tbl_par, rfrac_tbl_par, fwfisf_par, fwfisf_par_b, pcoef)
53      !
54      IF ( ln_isfcpl .AND. ln_rstart .AND. kt == nit000+1 ) THEN
55         DO jk = 1, jpkm1
56            pe3t_f(:,:,jk) =   pe3t_f(:,:,jk) - pcoef * risfcpl_vol(:,:,jk) * r1_e1e2t(:,:)
57         END DO
58      END IF
59      !
60   END SUBROUTINE isf_dynatf
61
62   SUBROUTINE isf_dynatf_mlt ( Kmm, pe3t_f, ktop, kbot, phtbl, pfrac, pfwf, pfwf_b, pcoef )
63      !!--------------------------------------------------------------------
64      !!                  ***  ROUTINE isf_dynatf_mlt  ***
65      !!
66      !! ** Purpose : compute the ice shelf volume filter correction for cavity or param
67      !!
68      !!-------------------------- IN  -------------------------------------
69      INTEGER                         , INTENT(in   ) :: Kmm             ! ocean time level index
70      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) :: pe3t_f          ! time-filtered scale factor to be corrected
71      INTEGER , DIMENSION(jpi,jpj)    , INTENT(in   ) :: ktop , kbot     ! top and bottom level of tbl
72      REAL(wp), DIMENSION(jpi,jpj)    , INTENT(in   ) :: pfrac, phtbl    ! fraction of bottom cell included in tbl, tbl thickness
73      REAL(wp), DIMENSION(jpi,jpj)    , INTENT(in   ) :: pfwf , pfwf_b   ! now/before fwf
74      REAL(wp),                         INTENT(in   ) :: pcoef           ! rn_atfp * rn_Dt * r1_rho0
75      !!----------------------------------------------------------------------
76      INTEGER :: ji,jj,jk
77      REAL(wp), DIMENSION(jpi,jpj) :: zfwfinc
78      !!----------------------------------------------------------------------
79      !
80      ! compute fwf conservation correction
81      zfwfinc(:,:) = pcoef * ( pfwf_b(:,:) - pfwf(:,:) ) / ( ht(:,:) + 1._wp - ssmask(:,:) ) * r1_rho0
82      !
83      ! add the increment
84      DO jk = 1, jpkm1
85         pe3t_f(:,:,jk) = pe3t_f(:,:,jk) - tmask(:,:,jk) * zfwfinc(:,:)   &
86            &                              * e3t(:,:,jk,Kmm)
87      END DO
88      !
89   END SUBROUTINE isf_dynatf_mlt
90
91END MODULE isfdynatf
Note: See TracBrowser for help on using the repository browser.