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

Ignore:
Timestamp:
2019-09-12T18:41:17+02:00 (5 years ago)
Author:
mathiot
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r11521 r11541  
    88   
    99   !!------------------------------------------------------------------------- 
    10    !!   dyn_nxt       : obtain the next (after) horizontal velocity 
     10   !!   isfnxt       : aplly correction need for the ice shelf 
    1111   !!------------------------------------------------------------------------- 
    1212 
    1313   USE isf 
     14   USE phycst 
    1415   USE dom_oce 
    1516   USE in_out_manager 
    … …  
    2324CONTAINS 
    2425 
    25    SUBROUTINE isf_dynnxt ( pcoef ) 
     26   SUBROUTINE isf_dynnxt ( kt, pcoef ) 
    2627      !!-------------------------------------------------------------------- 
    2728      !!                  ***  ROUTINE isf_dynnxt  *** 
    … …  
    3132      !!-------------------------------------------------------------------- 
    3233      !!-------------------------- OUT ------------------------------------- 
    33       REAL(wp),                     INTENT(in   ) :: pcoef           ! atfp * rdt * r1_rau0 
     34      INTEGER ,                     INTENT(in   ) :: kt 
     35      ! 
     36      REAL(wp),                     INTENT(in   ) :: pcoef           ! atfp * rdt 
    3437      !!-------------------------- IN  ------------------------------------- 
    3538      !!-------------------------------------------------------------------- 
     39      INTEGER :: jk  ! loop index 
    3640      !!-------------------------------------------------------------------- 
    3741      ! 
    … …  
    4145      ! ice shelf parametrised 
    4246      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 
    4353      ! 
    4454   END SUBROUTINE isf_dynnxt 
    … …  
    6373      ! 
    6474      ! compute fwf conservation correction 
    65       zfwfinc(:,:) = pcoef * ( pfwf_b(:,:) - pfwf(:,:) ) / phtbl(:,:) 
     75      zfwfinc(:,:) = pcoef * ( pfwf_b(:,:) - pfwf(:,:) ) / phtbl(:,:) * r1_rau0 
    6676      ! 
    6777      ! add the increment in the tbl 
Note: See TracChangeset for help on using the changeset viewer.