- Timestamp:
- 2016-07-19T10:38:35+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5549_BDY_ZEROGRAD/NEMOGCM/NEMO/OPA_SRC/SBC/sbcfwb.F90
r5120 r6808 12 12 13 13 !!---------------------------------------------------------------------- 14 !! sbc_fwb : freshwater budget for global ocean configurations15 !! in free surface and forced mode16 !!----------------------------------------------------------------------17 USE oce ! ocean dynamics and tracers18 USE dom_oce ! ocean space and time domain19 USE sbc_oce ! surface ocean boundary condition20 USE phycst ! physical constants21 USE sbc rnf ! ocean runoffs22 USE sbc isf ! ice shelf melting contribution23 USE sbcssr ! SS damping terms24 USE in_out_manager 25 USE lib_mpp 26 USE wrk_nemo 27 USE timing 28 USE lbclnk 29 USE lib_fortran 14 !! sbc_fwb : freshwater budget for global ocean configurations (free surface & forced mode) 15 !!---------------------------------------------------------------------- 16 USE oce ! ocean dynamics and tracers 17 USE dom_oce ! ocean space and time domain 18 USE sbc_oce ! surface ocean boundary condition 19 USE phycst ! physical constants 20 USE sbcrnf ! ocean runoffs 21 USE sbcisf ! ice shelf melting contribution 22 USE sbcssr ! Sea-Surface damping terms 23 ! 24 USE in_out_manager ! I/O manager 25 USE lib_mpp ! distribued memory computing library 26 USE wrk_nemo ! work arrays 27 USE timing ! Timing 28 USE lbclnk ! ocean lateral boundary conditions 29 USE lib_fortran ! 30 30 31 31 IMPLICIT NONE … … 40 40 41 41 !! * Substitutions 42 # include "domzgr_substitute.h90"43 42 # include "vectopt_loop_substitute.h90" 44 43 !!---------------------------------------------------------------------- … … 108 107 ! 109 108 IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN 110 z_fwf = glob_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + rdivisf * fwfisf(:,:) -snwice_fmass(:,:) ) ) / area ! sum over the global domain109 z_fwf = glob_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) - snwice_fmass(:,:) ) ) / area ! sum over the global domain 111 110 zcoef = z_fwf * rcp 112 111 emp(:,:) = emp(:,:) - z_fwf * tmask(:,:,1) … … 129 128 ENDIF 130 129 ! ! Update fwfold if new year start 131 ikty = 365 * 86400 / rdt tra(1)!!bug use of 365 days leap year or 360d year !!!!!!!130 ikty = 365 * 86400 / rdt !!bug use of 365 days leap year or 360d year !!!!!!! 132 131 IF( MOD( kt, ikty ) == 0 ) THEN 133 132 a_fwb_b = a_fwb ! mean sea level taking into account the ice+snow … … 162 161 zsurf_pos = glob_sum( e1e2t(:,:)*ztmsk_pos(:,:) ) 163 162 ! ! fwf global mean (excluding ocean to ice/snow exchanges) 164 z_fwf = glob_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + rdivisf *fwfisf(:,:) - snwice_fmass(:,:) ) ) / area163 z_fwf = glob_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) - snwice_fmass(:,:) ) ) / area 165 164 ! 166 165 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.