source: NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/ISF/isfdynatf.F90 @ 12616

Last change on this file since 12616 was 12616, checked in by techene, 11 months ago

all: add e3 substitute (sometimes it requires to add ze3t/u/v/w) and limit precompiled files lines to about 130 character, OCE/ASM/asminc.F90, OCE/DOM/domzgr_substitute.h90, OCE/ISF/isfcpl.F90, OCE/SBC/sbcice_cice, OCE/CRS/crsini.F90 : add key_LF

File size: 4.2 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_rau0         ! physical constant
16   USE dom_oce                        ! time and space domain
17
18   USE in_out_manager
19
20   IMPLICIT NONE
21
22   PRIVATE
23
24   PUBLIC isf_dynatf
25   !! * Substitutions
26#  include "do_loop_substitute.h90"
27#  include "domzgr_substitute.h90"
28
29CONTAINS
30
31   SUBROUTINE isf_dynatf ( kt, Kmm, pe3t_f, pcoef )
32      !!--------------------------------------------------------------------
33      !!                  ***  ROUTINE isf_dynatf  ***
34      !!
35      !! ** Purpose : compute the ice shelf volume filter correction for cavity, param, ice sheet coupling case
36      !!
37      !!-------------------------- OUT -------------------------------------
38      INTEGER                         , INTENT(in   ) :: kt       ! ocean time step
39      INTEGER                         , INTENT(in   ) :: Kmm      ! ocean time level index
40      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) :: pe3t_f   ! time filtered scale factor to be corrected
41      !
42      REAL(wp)                        , INTENT(in   ) :: pcoef    ! atfp * rdt * r1_rau0
43      !!--------------------------------------------------------------------
44      INTEGER :: jk  ! loop index
45      !!--------------------------------------------------------------------
46      !
47      ! ice shelf cavity
48      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)
49      !
50      ! ice shelf parametrised
51      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)
52      !
53      IF ( ln_isfcpl .AND. ln_rstart .AND. kt == nit000+1 ) THEN
54         DO jk = 1, jpkm1
55            pe3t_f(:,:,jk) =   pe3t_f(:,:,jk) - pcoef * risfcpl_vol(:,:,jk) * r1_e1e2t(:,:)
56         END DO
57      END IF
58      !
59   END SUBROUTINE isf_dynatf
60
61   SUBROUTINE isf_dynatf_mlt ( Kmm, pe3t_f, ktop, kbot, phtbl, pfrac, pfwf, pfwf_b, pcoef )
62      !!--------------------------------------------------------------------
63      !!                  ***  ROUTINE isf_dynatf_mlt  ***
64      !!
65      !! ** Purpose : compute the ice shelf volume filter correction for cavity or param
66      !!
67      !!-------------------------- IN  -------------------------------------
68      INTEGER                         , INTENT(in   ) :: Kmm             ! ocean time level index
69      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) :: pe3t_f          ! time-filtered scale factor to be corrected
70      INTEGER , DIMENSION(jpi,jpj)    , INTENT(in   ) :: ktop , kbot     ! top and bottom level of tbl
71      REAL(wp), DIMENSION(jpi,jpj)    , INTENT(in   ) :: pfrac, phtbl    ! fraction of bottom cell included in tbl, tbl thickness
72      REAL(wp), DIMENSION(jpi,jpj)    , INTENT(in   ) :: pfwf , pfwf_b   ! now/before fwf
73      REAL(wp),                         INTENT(in   ) :: pcoef           ! atfp * rdt * r1_rau0
74      !!----------------------------------------------------------------------
75      INTEGER :: ji,jj,jk
76      REAL(wp), DIMENSION(jpi,jpj) :: zfwfinc
77      !!----------------------------------------------------------------------
78      !
79      ! compute fwf conservation correction
80      zfwfinc(:,:) = pcoef * ( pfwf_b(:,:) - pfwf(:,:) ) / ( ht(:,:) + 1._wp - ssmask(:,:) ) * r1_rau0
81      !
82      ! add the increment
83      DO jk = 1, jpkm1
84         pe3t_f(:,:,jk) = pe3t_f(:,:,jk) - tmask(:,:,jk) * zfwfinc(:,:)   &
85            &                              * e3t(:,:,jk,Kmm)
86      END DO
87      !
88   END SUBROUTINE isf_dynatf_mlt
89
90END MODULE isfdynatf
Note: See TracBrowser for help on using the repository browser.