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 10726 for NEMO/releases/release-4.0/src/OCE/ICB/icbclv.F90 – NEMO

Ignore:
Timestamp:
2019-02-27T16:06:35+01:00 (5 years ago)
Author:
mathiot
Message:

changes related to bug fixes described in tickets #2228, #2229, #2238 and #1595 in NEMO 4.0

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/release-4.0/src/OCE/ICB/icbclv.F90

    r10425 r10726  
    2424   USE icbdia         ! iceberg diagnostics 
    2525   USE icbutl         ! iceberg utility routines 
     26   USE icb_oce        ! iceberg parameters  
    2627 
    2728   IMPLICIT NONE 
     
    5758      ! this assumes that input is given as equivalent water flux so that pure water density is appropriate 
    5859 
    59       zfact = ( (1000._wp)**3 / ( NINT(rday) * nyear_len(1) ) ) * 850._wp 
    60       berg_grid%calving(:,:) = src_calving(:,:) * tmask_i(:,:) * zfact 
     60      zfact = ( (1000._wp)**3 / ( NINT(rday) * nyear_len(1) ) ) * rn_rho_bergs 
     61      berg_grid%calving(:,:) = src_calving(:,:) * zfact * tmask_i(:,:) * tmask(:,:,1) 
    6162 
    6263      ! Heat in units of W/m2, and mask (just in case) 
    63       berg_grid%calving_hflx(:,:) = src_calving_hflx(:,:) * tmask_i(:,:) 
     64      berg_grid%calving_hflx(:,:) = src_calving_hflx(:,:) * tmask_i(:,:) * tmask(:,:,1) 
    6465 
    6566      IF( ll_first_call .AND. .NOT. l_restarted_bergs ) THEN      ! This is a hack to simplify initialization 
Note: See TracChangeset for help on using the changeset viewer.