Changeset 15004 for NEMO/trunk/tests
- Timestamp:
- 2021-06-16T12:33:18+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/tests/ISOMIP+/MY_SRC/sbcfwb.F90
r14995 r15004 118 118 ! 119 119 IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN 120 y_fwfnow(1) = local_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) +fwfisf_par(:,:) - snwice_fmass(:,:) ) )120 y_fwfnow(1) = local_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) - fwfisf_cav(:,:) - fwfisf_par(:,:) - snwice_fmass(:,:) ) ) 121 121 CALL mpp_delay_sum( 'sbcfwb', 'fwb', y_fwfnow(:), z_fwfprv(:), kt == nitend - nn_fsbc + 1 ) 122 122 z_fwfprv(1) = z_fwfprv(1) / area … … 132 132 ! 133 133 IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN 134 z_fwf = glob_sum( 'sbcfwb', e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) +fwfisf_par(:,:) - snwice_fmass(:,:) ) )134 z_fwf = glob_sum( 'sbcfwb', e1e2t(:,:) * ( emp(:,:) - rnf(:,:) - fwfisf_cav(:,:) - fwfisf_par(:,:) - snwice_fmass(:,:) ) ) 135 135 ! 136 136 ! correction for ice sheet coupling testing (ie remove the excess through the surface) … … 201 201 ztmsk_neg(:,:) = tmask_i(:,:) - ztmsk_pos(:,:) 202 202 ! ! fwf global mean (excluding ocean to ice/snow exchanges) 203 z_fwf = glob_sum( 'sbcfwb', e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) +fwfisf_par(:,:) - snwice_fmass(:,:) ) ) / area203 z_fwf = glob_sum( 'sbcfwb', e1e2t(:,:) * ( emp(:,:) - rnf(:,:) - fwfisf_cav(:,:) - fwfisf_par(:,:) - snwice_fmass(:,:) ) ) / area 204 204 ! 205 205 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.