- Timestamp:
- 2019-11-13T12:40:52+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/ENHANCE-02_ISF_nemo/tests/ISOMIP+/MY_SRC/sbcfwb.F90
r11889 r11896 105 105 ! 106 106 IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN 107 y_fwfnow(1) = local_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) - snwice_fmass(:,:) ) ) 108 CALL mpp_delay_sum( 'sbcfwb', 'fwb', y_fwfnow(:), z_fwfprv(:), kt == nitend - nn_fsbc + 1 ) 109 z_fwfprv(1) = z_fwfprv(1) / area 110 zcoef = z_fwfprv(1) * rcp 111 emp(:,:) = emp(:,:) - z_fwfprv(1) * tmask(:,:,1) 112 qns(:,:) = qns(:,:) + zcoef * sst_m(:,:) * tmask(:,:,1) ! account for change to the heat budget due to fw correction 113 ENDIF 114 ! 115 CASE ( 4 ) !== global mean fwf set to zero (ISOMIP case) ==! 116 ! 117 IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN 107 118 z_fwf = glob_sum( 'sbcfwb', e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf_cav(:,:) + fwfisf_par(:,:) - snwice_fmass(:,:) ) ) 108 119 z_fwf = z_fwf / area 109 120 zcoef = z_fwf * rcp 110 emp(:,:) = emp(:,:) - z_fwf * tmask(:,:,1)111 qns(:,:) = qns(:,:) + zcoef * sst_m(:,:) * tmask(:,:,1) ! account for change to the heat budget due to fw correction112 sfx(:,:) = sfx(:,:) + z_fwf * sss_m(:,:) * tmask(:,:,1)121 emp(:,:) = emp(:,:) - z_fwf * tmask(:,:,1) ! (Eq. 34 AD2015) 122 qns(:,:) = qns(:,:) + zcoef * ( -1.9 ) * tmask(:,:,1) ! (Eq. 35 AD2015) ! could be sst_m if we don't want any bouyancy fluxes 123 sfx(:,:) = sfx(:,:) + z_fwf * ( 33.8 ) * tmask(:,:,1) ! (Eq. 36 AD2015) ! could be sss_m if we don't want any bouyancy fluxes 113 124 ENDIF 114 125 !
Note: See TracChangeset
for help on using the changeset viewer.