- Timestamp:
- 2018-11-07T18:25:49+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r9866_HPC_03_globcom/src/OCE/SBC/sbcfwb.F90
r9598 r10288 44 44 !! NEMO/OCE 4.0 , NEMO Consortium (2018) 45 45 !! $Id$ 46 !! Software governed by the CeCILL licen ce (./LICENSE)46 !! Software governed by the CeCILL license (see ./LICENSE) 47 47 !!---------------------------------------------------------------------- 48 48 CONTAINS … … 155 155 ztmsk_neg(:,:) = tmask_i(:,:) - ztmsk_pos(:,:) 156 156 ! 157 zsurf_neg = glob_sum( e1e2t(:,:)*ztmsk_neg(:,:) ) ! Area filled by <0 and >0 erp158 zsurf_pos = glob_sum( e1e2t(:,:)*ztmsk_pos(:,:) )159 157 ! ! fwf global mean (excluding ocean to ice/snow exchanges) 160 158 z_fwf = glob_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) - snwice_fmass(:,:) ) ) / area 161 159 ! 162 160 IF( z_fwf < 0._wp ) THEN ! spread out over >0 erp area to increase evaporation 161 zsurf_pos = glob_sum( e1e2t(:,:)*ztmsk_pos(:,:) ) 163 162 zsurf_tospread = zsurf_pos 164 163 ztmsk_tospread(:,:) = ztmsk_pos(:,:) 165 164 ELSE ! spread out over <0 erp area to increase precipitation 165 zsurf_neg = glob_sum( e1e2t(:,:)*ztmsk_neg(:,:) ) ! Area filled by <0 and >0 erp 166 166 zsurf_tospread = zsurf_neg 167 167 ztmsk_tospread(:,:) = ztmsk_neg(:,:)
Note: See TracChangeset
for help on using the changeset viewer.