Ignore:
Timestamp:
2020-05-05T17:57:08+02:00 (9 months ago)
Author:
mathiot
Message:

ticket #2457: add fix suggested in ticket.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • utils/tools_ticket2457/DOMAINcfg/src/domzgr.F90

    r12414 r12871  
    895895      IF(lwp) WRITE(numout,*) '              mbathy is recomputed : bathy_level file is NOT used' 
    896896 
    897       ! compute position of the ice shelf grounding line 
    898       ! set bathy and isfdraft to 0 where grounded 
    899       IF ( ln_isfcav ) CALL zgr_isf_zspace 
    900  
    901897      ! bathymetry in level (from bathy_meter) 
    902898      ! =================== 
     
    916912      END DO 
    917913 
    918       ! Check compatibility between bathy and iceshelf draft 
    919       ! insure at least 2 wet level on the vertical under an ice shelf 
    920       ! compute misfdep and adjust isf draft if needed 
    921       IF ( ln_isfcav ) CALL zgr_isf_kspace 
     914      ! isf draft in z and level 
     915      ! ======================== 
     916      ! if under ice shelf cavities, compute misfdep for ocean point (i.e. the top level) 
     917      IF ( ln_isfcav ) THEN 
     918         ! compute position of the ice shelf grounding line 
     919         CALL zgr_isf_zspace 
     920         ! compute misfdep and adjust isf draft if needed 
     921         CALL zgr_isf_kspace 
     922      END IF 
    922923 
    923924      ! Scale factors and depth at T- and W-points 
Note: See TracChangeset for help on using the changeset viewer.