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

Ignore:
Timestamp:
2019-08-05T19:14:52+02:00 (5 years ago)
Author:
mathiot
Message:

ENHANCE-02_ISF_nemo : add comments, renaming file (AGRIF), add isfload module (ticket #2142)

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfnxt.F90

    r11395 r11403  
    11MODULE 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   !!   dyn_nxt       : obtain the next (after) horizontal velocity 
     11   !!------------------------------------------------------------------------- 
    212 
    3 USE isf 
    4 USE isfutils 
     13   USE isf 
     14   USE dom_oce 
     15   USE in_out_manager 
    516 
    6 USE dom_oce 
    7 USE in_out_manager 
     17   IMPLICIT NONE 
    818 
    9 IMPLICIT NONE 
     19   PRIVATE 
    1020 
    11 PRIVATE 
    12  
    13 PUBLIC isf_dynnxt !, isf_tranxt_mlt, isf_dynnxt_cpl, isf_tranxt_cpl 
     21   PUBLIC isf_dynnxt ! isf_tranxt 
    1422 
    1523CONTAINS 
    1624 
    1725   SUBROUTINE isf_dynnxt ( pcoef ) 
    18       !!---------------------------------------------------------------------- 
     26      !!-------------------------------------------------------------------- 
     27      !!                  ***  ROUTINE isf_dynnxt  *** 
     28      !! 
     29      !! ** Purpose : compute the ice shelf volume filter correction for cavity, param, ice sheet coupling case 
     30      !! 
     31      !!-------------------------------------------------------------------- 
     32      !!-------------------------- OUT ------------------------------------- 
    1933      REAL(wp),                     INTENT(in   ) :: pcoef           ! atfp * rdt * r1_rau0 
    20       !!---------------------------------------------------------------------- 
     34      !!-------------------------- IN  ------------------------------------- 
     35      !!-------------------------------------------------------------------- 
     36      !!-------------------------------------------------------------------- 
    2137      ! 
    2238      ! ice shelf cavity 
    … …  
    3248 
    3349   SUBROUTINE isf_dynnxt_mlt ( ktop, kbot, phtbl, pfrac, pfwf, pfwf_b, pcoef ) 
    34    !atfp * rdt * r1_rau0 
    35       !!---------------------------------------------------------------------- 
     50      !!-------------------------------------------------------------------- 
     51      !!                  ***  ROUTINE isf_dynnxt_mlt  *** 
     52      !! 
     53      !! ** Purpose : compute the ice shelf volume filter correction for cavity or param 
     54      !! 
     55      !!-------------------------------------------------------------------- 
     56      !!-------------------------- OUT ------------------------------------- 
     57      !!-------------------------- IN  ------------------------------------- 
    3658      INTEGER , DIMENSION(jpi,jpj), INTENT(in   ) :: ktop , kbot     ! top and bottom level of tbl 
    3759      REAL(wp), DIMENSION(jpi,jpj), INTENT(in   ) :: pfrac, phtbl    ! fraction of bottom cell included in tbl, tbl thickness 
    … …  
    4668      zfwfinc(:,:) = pcoef * ( pfwf_b(:,:) - pfwf(:,:) ) / phtbl(:,:) 
    4769      ! 
    48       ! add the in depth increment 
     70      ! add the increment in the tbl 
    4971      DO jk = 1, jpkm1 
    5072         DO jj = 1, jpj 
Note: See TracChangeset for help on using the changeset viewer.