Changeset 11789


Ignore:
Timestamp:
2019-10-24T17:20:17+02:00 (11 months ago)
Author:
dancopsey
Message:

Add missing divide by icea fraction to sea ice fluxes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0_fix_cpl_oce_only/src/OCE/SBC/sbccpl.F90

    r11657 r11789  
    20512051      CASE ('coupled') 
    20522052         IF (ln_scale_ice_fluxes) THEN 
    2053             qml_ice(:,:,:) = frcv(jpr_topm)%z3(:,:,:) * a_i_last_couple(:,:,:) 
    2054             qcn_ice(:,:,:) = frcv(jpr_botm)%z3(:,:,:) * a_i_last_couple(:,:,:) 
     2053            WHERE( a_i(:,:,:) > 0.0_wp ) qml_ice(:,:,:) = frcv(jpr_topm)%z3(:,:,:) * a_i_last_couple(:,:,:) / a_i(:,:,:) 
     2054            WHERE( a_i(:,:,:) <= 0.0_wp ) qml_ice(:,:,:) = 0.0_wp 
     2055            WHERE( a_i(:,:,:) > 0.0_wp ) qcn_ice(:,:,:) = frcv(jpr_botm)%z3(:,:,:) * a_i_last_couple(:,:,:) / a_i(:,:,:) 
     2056            WHERE( a_i(:,:,:) <= 0.0_wp ) qcn_ice(:,:,:) = 0.0_wp 
    20552057         ELSE 
    20562058            qml_ice(:,:,:) = frcv(jpr_topm)%z3(:,:,:) 
Note: See TracChangeset for help on using the changeset viewer.