- Timestamp:
- 2018-12-01T12:33:31+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/OCE/SBC/sbcfwb.F90
r10314 r10365 154 154 WHERE( erp < 0._wp ) ztmsk_pos = 0._wp 155 155 ztmsk_neg(:,:) = tmask_i(:,:) - ztmsk_pos(:,:) 156 !157 zsurf_neg = glob_sum( 'sbcfwb', e1e2t(:,:)*ztmsk_neg(:,:) ) ! Area filled by <0 and >0 erp158 zsurf_pos = glob_sum( 'sbcfwb', e1e2t(:,:)*ztmsk_pos(:,:) )159 156 ! ! fwf global mean (excluding ocean to ice/snow exchanges) 160 157 z_fwf = glob_sum( 'sbcfwb', e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) - snwice_fmass(:,:) ) ) / area 161 158 ! 162 159 IF( z_fwf < 0._wp ) THEN ! spread out over >0 erp area to increase evaporation 163 zsurf_tospread = zsurf_pos 164 ztmsk_tospread(:,:) = ztmsk_pos(:,:) 160 zsurf_pos = glob_sum( 'sbcfwb', e1e2t(:,:)*ztmsk_pos(:,:) ) 161 zsurf_tospread = zsurf_pos 162 ztmsk_tospread(:,:) = ztmsk_pos(:,:) 165 163 ELSE ! spread out over <0 erp area to increase precipitation 166 zsurf_tospread = zsurf_neg 167 ztmsk_tospread(:,:) = ztmsk_neg(:,:) 164 zsurf_neg = glob_sum( 'sbcfwb', e1e2t(:,:)*ztmsk_neg(:,:) ) ! Area filled by <0 and >0 erp 165 zsurf_tospread = zsurf_neg 166 ztmsk_tospread(:,:) = ztmsk_neg(:,:) 168 167 ENDIF 169 168 !
Note: See TracChangeset
for help on using the changeset viewer.