Ignore:
Timestamp:
2019-11-08T12:26:42+01:00 (11 months ago)
Author:
mathiot
Message:

ENHANCE-02_ISF_nemo: add timing in the main isf routine + various bug fixes + cosmetic changes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfpar.F90

    r11852 r11876  
    2323   USE isfdiags , ONLY: isf_diags_flx  ! ice shelf diags subroutine 
    2424   ! 
    25    USE dom_oce  , ONLY: jpi,jpj        ! ocean space and time domain 
     25   USE dom_oce  , ONLY: jpi,jpj, bathy ! ocean space and time domain 
    2626   USE phycst   , ONLY: r1_rau0_rcp    ! physical constants 
    2727   ! 
     
    124124      END WHERE 
    125125      ! 
    126       ! compute ktop 
    127       CALL isf_tbl_ktop(ztblmin, misfkt_par) 
     126      ! ensure ztblmax <= bathy 
     127      WHERE ( ztblmax(:,:) > bathy(:,:) ) 
     128         ztblmax(:,:) = bathy(:,:) 
     129      END WHERE 
     130      ! 
     131      ! compute ktop and update ztblmin to gdepw_0(misfkt_par)  
     132      CALL isf_tbl_ktop(ztblmin, misfkt_par) !   out: misfkt_par 
     133      !                                      ! inout: ztblmin 
    128134      ! 
    129135      ! initial tbl thickness 
     
    132138      ! define iceshelf parametrisation mask 
    133139      mskisf_par = 0 
    134       WHERE ( ztblmax > 0._wp ) 
     140      WHERE ( rhisf0_tbl_par(:,:) > 0._wp ) 
    135141         mskisf_par(:,:) = 1._wp 
    136142      END WHERE 
Note: See TracChangeset for help on using the changeset viewer.