Changeset 13353


Ignore:
Timestamp:
2020-07-29T13:06:44+02:00 (13 days ago)
Author:
dancopsey
Message:

Correctly scale ice fluxes by ice fractions from last coupling point.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.1_penetrating_solar/src/OCE/SBC/sbccpl.F90

    r13347 r13353  
    20442044         !      ! ===> here we receive the qtr_ice_top array from the coupler 
    20452045         CASE ('coupled') 
    2046             qtr_ice_top(:,:,:) = frcv(jpr_qtr)%z3(:,:,:) 
     2046            IF (ln_scale_ice_fluxes) THEN 
     2047               WHERE( a_i(:,:,:) > 0.0_wp ) qtr_ice_top(:,:,:) = frcv(jpr_qtr)%z3(:,:,:) * a_i_last_couple(:,:,:) / a_i(:,:,:) 
     2048               WHERE( a_i(:,:,:) <= 0.0_wp ) qtr_ice_top(:,:,:) = 0.0_wp 
     2049            ELSE 
     2050               qtr_ice_top(:,:,:) = frcv(jpr_qtr)%z3(:,:,:) 
     2051            ENDIF 
    20472052 
    20482053         !      if we are not getting this data from the coupler then assume zero (fully opaque ice) 
Note: See TracChangeset for help on using the changeset viewer.