Changeset 10710


Ignore:
Timestamp:
2019-02-21T12:22:22+01:00 (18 months ago)
Author:
mathiot
Message:

branch and changes related to ticket #1595

Location:
NEMO/branches/2019/branch_ticket1595
Files:
1 edited
1 copied

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/branch_ticket1595/src/OCE/ICB/icbclv.F90

    r10425 r10710  
    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.