- Timestamp:
- 2017-09-27T11:11:01+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd_sal.F90
r8563 r8564 58 58 INTEGER :: ji, jk ! dummy loop indices 59 59 REAL(wp) :: iflush, igravdr ! local scalars 60 REAL(wp) :: zs_sni, zs m_i_gd, zsm_i_fl, zsm_i_si, zsm_i_bg ! local scalars60 REAL(wp) :: zs_sni, zs_i_gd, zs_i_fl, zs_i_si, zs_i_bg ! local scalars 61 61 REAL(wp) :: z1_time_gd, z1_time_fl 62 62 !!--------------------------------------------------------------------- … … 77 77 IF( h_i_1d(ji) > 0._wp ) THEN 78 78 zs_sni = sss_1d(ji) * ( rhoic - rhosn ) * r1_rhoic ! Salinity of snow ice 79 zs m_i_si = ( zs_sni - sm_i_1d(ji) ) * dh_snowice(ji) / h_i_1d(ji) ! snow-ice80 zs m_i_bg = ( s_i_new(ji) - sm_i_1d(ji) ) * MAX( 0._wp, dh_i_bott(ji) ) / h_i_1d(ji) ! bottom growth79 zs_i_si = ( zs_sni - s_i_1d(ji) ) * dh_snowice(ji) / h_i_1d(ji) ! snow-ice 80 zs_i_bg = ( s_i_new(ji) - s_i_1d(ji) ) * MAX( 0._wp, dh_i_bott(ji) ) / h_i_1d(ji) ! bottom growth 81 81 ! Update salinity (nb: salt flux already included in icethd_dh) 82 s m_i_1d(ji) = sm_i_1d(ji) + zsm_i_bg + zsm_i_si82 s_i_1d(ji) = s_i_1d(ji) + zs_i_bg + zs_i_si 83 83 ENDIF 84 84 … … 90 90 igravdr = MAX( 0._wp , SIGN( 1._wp , t_bo_1d(ji) - t_su_1d(ji) ) ) ! =1 if t_su < t_bo 91 91 92 zs m_i_gd = - igravdr * MAX( sm_i_1d(ji) - rn_sal_gd , 0._wp ) * z1_time_gd ! gravity drainage93 zs m_i_fl = - iflush * MAX( sm_i_1d(ji) - rn_sal_fl , 0._wp ) * z1_time_fl ! flushing92 zs_i_gd = - igravdr * MAX( s_i_1d(ji) - rn_sal_gd , 0._wp ) * z1_time_gd ! gravity drainage 93 zs_i_fl = - iflush * MAX( s_i_1d(ji) - rn_sal_fl , 0._wp ) * z1_time_fl ! flushing 94 94 95 95 ! Update salinity 96 s m_i_1d(ji) = sm_i_1d(ji) + zsm_i_fl + zsm_i_gd96 s_i_1d(ji) = s_i_1d(ji) + zs_i_fl + zs_i_gd 97 97 98 98 ! Salt flux 99 sfx_bri_1d(ji) = sfx_bri_1d(ji) - rhoic * a_i_1d(ji) * h_i_1d(ji) * ( zs m_i_fl + zsm_i_gd ) * r1_rdtice99 sfx_bri_1d(ji) = sfx_bri_1d(ji) - rhoic * a_i_1d(ji) * h_i_1d(ji) * ( zs_i_fl + zs_i_gd ) * r1_rdtice 100 100 ENDIF 101 101 END DO
Note: See TracChangeset
for help on using the changeset viewer.