- Timestamp:
- 2019-12-05T13:18:21+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/UKMO_MERGE_2019/src/OCE/SBC/sbcfwb.F90
r11027 r12068 20 20 USE phycst ! physical constants 21 21 USE sbcrnf ! ocean runoffs 22 USE sbcisf! ice shelf melting contribution22 USE isf ! ice shelf melting contribution 23 23 USE sbcssr ! Sea-Surface damping terms 24 24 ! … … 105 105 ! 106 106 IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN 107 y_fwfnow(1) = local_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf (:,:) - snwice_fmass(:,:) ) )107 y_fwfnow(1) = local_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) - snwice_fmass(:,:) ) ) 108 108 CALL mpp_delay_sum( 'sbcfwb', 'fwb', y_fwfnow(:), z_fwfprv(:), kt == nitend - nn_fsbc + 1 ) 109 109 z_fwfprv(1) = z_fwfprv(1) / area … … 160 160 ztmsk_neg(:,:) = tmask_i(:,:) - ztmsk_pos(:,:) 161 161 ! ! fwf global mean (excluding ocean to ice/snow exchanges) 162 z_fwf = glob_sum( 'sbcfwb', e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf (:,:) - snwice_fmass(:,:) ) ) / area162 z_fwf = glob_sum( 'sbcfwb', e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) - snwice_fmass(:,:) ) ) / area 163 163 ! 164 164 IF( z_fwf < 0._wp ) THEN ! spread out over >0 erp area to increase evaporation
Note: See TracChangeset
for help on using the changeset viewer.