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.
isfnxt.F90 in NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF – NEMO

source: NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfnxt.F90 @ 11541

Last change on this file since 11541 was 11541, checked in by mathiot, 3 years ago

ENHANCE-02_ISF: simplify use of ln_isf, add extra comments + minor changes (ticket #2142)

File size: 4.0 KB
Line 
1MODULE isfnxt
2   !!=========================================================================
3   !!                       ***  MODULE  isfnxt  ***
4   !! Ice shelf update: compute the dynnxt ice shelf contribution
5   !!=========================================================================
6   !! History :  OPA  !  2019-09  (P. Mathiot)  Original code
7   !!-------------------------------------------------------------------------
8 
9   !!-------------------------------------------------------------------------
10   !!   isfnxt       : aplly correction need for the ice shelf
11   !!-------------------------------------------------------------------------
12
13   USE isf
14   USE phycst
15   USE dom_oce
16   USE in_out_manager
17
18   IMPLICIT NONE
19
20   PRIVATE
21
22   PUBLIC isf_dynnxt ! isf_tranxt
23
24CONTAINS
25
26   SUBROUTINE isf_dynnxt ( kt, pcoef )
27      !!--------------------------------------------------------------------
28      !!                  ***  ROUTINE isf_dynnxt  ***
29      !!
30      !! ** Purpose : compute the ice shelf volume filter correction for cavity, param, ice sheet coupling case
31      !!
32      !!--------------------------------------------------------------------
33      !!-------------------------- OUT -------------------------------------
34      INTEGER ,                     INTENT(in   ) :: kt
35      !
36      REAL(wp),                     INTENT(in   ) :: pcoef           ! atfp * rdt
37      !!-------------------------- IN  -------------------------------------
38      !!--------------------------------------------------------------------
39      INTEGER :: jk  ! loop index
40      !!--------------------------------------------------------------------
41      !
42      ! ice shelf cavity
43      IF ( ln_isfcav_mlt ) CALL isf_dynnxt_mlt(misfkt_cav, misfkb_cav, rhisf_tbl_cav, rfrac_tbl_cav, fwfisf_cav, fwfisf_cav_b, pcoef)
44      !
45      ! ice shelf parametrised
46      IF ( ln_isfpar_mlt ) CALL isf_dynnxt_mlt(misfkt_par, misfkb_par, rhisf_tbl_par, rfrac_tbl_par, fwfisf_par, fwfisf_par_b, pcoef)
47      !
48      IF ( ln_isfcpl .AND. ln_rstart .AND. kt == nit000+1 ) THEN
49         DO jk = 1, jpkm1
50            e3t_b(:,:,jk) =   e3t_b(:,:,jk) - pcoef * risfcpl_vol(:,:,jk) * r1_e1e2t(:,:)
51         END DO
52      END IF
53      !
54   END SUBROUTINE isf_dynnxt
55
56   SUBROUTINE isf_dynnxt_mlt ( ktop, kbot, phtbl, pfrac, pfwf, pfwf_b, pcoef )
57      !!--------------------------------------------------------------------
58      !!                  ***  ROUTINE isf_dynnxt_mlt  ***
59      !!
60      !! ** Purpose : compute the ice shelf volume filter correction for cavity or param
61      !!
62      !!--------------------------------------------------------------------
63      !!-------------------------- OUT -------------------------------------
64      !!-------------------------- IN  -------------------------------------
65      INTEGER , DIMENSION(jpi,jpj), INTENT(in   ) :: ktop , kbot     ! top and bottom level of tbl
66      REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: pfrac, phtbl    ! fraction of bottom cell included in tbl, tbl thickness
67      REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: pfwf , pfwf_b   ! now/before fwf
68      REAL(wp),                     INTENT(in   ) :: pcoef           ! atfp * rdt * r1_rau0
69      !!----------------------------------------------------------------------
70      INTEGER :: ji,jj,jk
71      REAL(wp), DIMENSION(jpi,jpj) :: zfwfinc
72      !!----------------------------------------------------------------------
73      !
74      ! compute fwf conservation correction
75      zfwfinc(:,:) = pcoef * ( pfwf_b(:,:) - pfwf(:,:) ) / phtbl(:,:) * r1_rau0
76      !
77      ! add the increment in the tbl
78      DO jk = 1, jpkm1
79         DO jj = 1, jpj
80            DO ji = 1, jpi
81               IF( ktop(ji,jj) <= jk .AND. jk < kbot(ji,jj)  ) THEN
82                  e3t_b(ji,jj,jk) = e3t_b(ji,jj,jk) - zfwfinc(ji,jj) * e3t_n(ji,jj,jk)
83               ELSEIF ( jk == kbot(ji,jj) ) THEN
84                  e3t_b(ji,jj,jk) = e3t_b(ji,jj,jk) - zfwfinc(ji,jj) * e3t_n(ji,jj,jk) * pfrac(ji,jj)
85               ENDIF
86            END DO
87         END DO
88      END DO
89      !
90   END SUBROUTINE isf_dynnxt_mlt
91
92END MODULE isfnxt
Note: See TracBrowser for help on using the repository browser.