Changeset 12377 for NEMO/trunk/src/OCE/SBC/sbcfwb.F90
- Timestamp:
- 2020-02-12T15:39:06+01:00 (4 years ago)
- Location:
- NEMO/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEAD ext/AGRIF5 ^/vendors/AGRIF/dev_r11615_ENHANCE-04_namelists_as_internalfiles_agrif@HEAD ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL
-
- Property svn:externals
-
NEMO/trunk/src/OCE/SBC/sbcfwb.F90
r10570 r12377 17 17 USE dom_oce ! ocean space and time domain 18 18 USE sbc_oce ! surface ocean boundary condition 19 USE isf_oce , ONLY : fwfisf_cav, fwfisf_par ! 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 sbcisf ! ice shelf melting contribution23 23 USE sbcssr ! Sea-Surface damping terms 24 24 ! … … 39 39 REAL(wp) :: area ! global mean ocean surface (interior domain) 40 40 41 !! * Substitutions42 # include "vectopt_loop_substitute.h90"43 41 !!---------------------------------------------------------------------- 44 42 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 48 46 CONTAINS 49 47 50 SUBROUTINE sbc_fwb( kt, kn_fwb, kn_fsbc )48 SUBROUTINE sbc_fwb( kt, kn_fwb, kn_fsbc, Kmm ) 51 49 !!--------------------------------------------------------------------- 52 50 !! *** ROUTINE sbc_fwb *** … … 65 63 INTEGER, INTENT( in ) :: kn_fsbc ! 66 64 INTEGER, INTENT( in ) :: kn_fwb ! ocean time-step index 65 INTEGER, INTENT( in ) :: Kmm ! ocean time level index 67 66 ! 68 67 INTEGER :: inum, ikty, iyear ! local integers … … 104 103 ! 105 104 IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN 106 y_fwfnow(1) = local_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf (:,:) - snwice_fmass(:,:) ) )105 y_fwfnow(1) = local_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) - snwice_fmass(:,:) ) ) 107 106 CALL mpp_delay_sum( 'sbcfwb', 'fwb', y_fwfnow(:), z_fwfprv(:), kt == nitend - nn_fsbc + 1 ) 108 107 z_fwfprv(1) = z_fwfprv(1) / area … … 131 130 a_fwb_b = a_fwb ! mean sea level taking into account the ice+snow 132 131 ! sum over the global domain 133 a_fwb = glob_sum( 'sbcfwb', e1e2t(:,:) * ( ssh n(:,:) + snwice_mass(:,:) * r1_rau0 ) )132 a_fwb = glob_sum( 'sbcfwb', e1e2t(:,:) * ( ssh(:,:,Kmm) + snwice_mass(:,:) * r1_rau0 ) ) 134 133 a_fwb = a_fwb * 1.e+3 / ( area * rday * 365. ) ! convert in Kg/m3/s = mm/s 135 134 !!gm ! !!bug 365d year … … 159 158 ztmsk_neg(:,:) = tmask_i(:,:) - ztmsk_pos(:,:) 160 159 ! ! fwf global mean (excluding ocean to ice/snow exchanges) 161 z_fwf = glob_sum( 'sbcfwb', e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf (:,:) - snwice_fmass(:,:) ) ) / area160 z_fwf = glob_sum( 'sbcfwb', e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) - snwice_fmass(:,:) ) ) / area 162 161 ! 163 162 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.