- Timestamp:
- 2019-12-05T18:41:39+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/UKMO_MERGE_2019/tests/ISOMIP+/MY_SRC/sbcfwb.F90
r11931 r12077 17 17 USE dom_oce ! ocean space and time domain 18 18 USE sbc_oce ! surface ocean boundary condition 19 USE isf_oce ! ice shelf melting contribution 19 20 USE sbc_ice , ONLY : snwice_mass, snwice_mass_b, snwice_fmass 20 21 USE phycst ! physical constants 21 22 USE sbcrnf ! ocean runoffs 22 USE isf ! ice shelf melting contribution23 23 USE sbcssr ! Sea-Surface damping terms 24 USE tradmp !25 24 ! 26 25 USE in_out_manager ! I/O manager … … 49 48 CONTAINS 50 49 51 SUBROUTINE sbc_fwb( kt, kn_fwb, kn_fsbc )50 SUBROUTINE sbc_fwb( kt, kn_fwb, kn_fsbc, Kmm ) 52 51 !!--------------------------------------------------------------------- 53 52 !! *** ROUTINE sbc_fwb *** … … 66 65 INTEGER, INTENT( in ) :: kn_fsbc ! 67 66 INTEGER, INTENT( in ) :: kn_fwb ! ocean time-step index 67 INTEGER, INTENT( in ) :: Kmm ! ocean time level index 68 68 ! 69 69 INTEGER :: inum, ikty, iyear ! local integers … … 157 157 a_fwb_b = a_fwb ! mean sea level taking into account the ice+snow 158 158 ! sum over the global domain 159 a_fwb = glob_sum( 'sbcfwb', e1e2t(:,:) * ( ssh n(:,:) + snwice_mass(:,:) * r1_rau0 ) )159 a_fwb = glob_sum( 'sbcfwb', e1e2t(:,:) * ( ssh(:,:,Kmm) + snwice_mass(:,:) * r1_rau0 ) ) 160 160 a_fwb = a_fwb * 1.e+3 / ( area * rday * 365. ) ! convert in Kg/m3/s = mm/s 161 161 !!gm ! !!bug 365d year
Note: See TracChangeset
for help on using the changeset viewer.