- 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/OPA_SRC/BDY/bdyice.F90
r8563 r8564 174 174 ! 175 175 ! Ice salinity, age, temperature 176 s m_i(ji,jj,jl) = rswitch * rn_ice_sal(ib_bdy) + ( 1.0 - rswitch ) * rn_simin176 s_i (ji,jj,jl) = rswitch * rn_ice_sal(ib_bdy) + ( 1.0 - rswitch ) * rn_simin 177 177 oa_i(ji,jj,jl) = rswitch * rn_ice_age(ib_bdy) * a_i(ji,jj,jl) 178 178 t_su(ji,jj,jl) = rswitch * rn_ice_tem(ib_bdy) + ( 1.0 - rswitch ) * rn_ice_tem(ib_bdy) … … 181 181 END DO 182 182 DO jk = 1, nlay_i 183 t_i (ji,jj,jk,jl) = rswitch * rn_ice_tem(ib_bdy) + ( 1.0 - rswitch ) * rt0184 s _i(ji,jj,jk,jl) = rswitch * rn_ice_sal(ib_bdy) + ( 1.0 - rswitch ) * rn_simin183 t_i (ji,jj,jk,jl) = rswitch * rn_ice_tem(ib_bdy) + ( 1.0 - rswitch ) * rt0 184 sz_i(ji,jj,jk,jl) = rswitch * rn_ice_sal(ib_bdy) + ( 1.0 - rswitch ) * rn_simin 185 185 END DO 186 186 ! … … 188 188 ! 189 189 ! Ice salinity, age, temperature 190 s m_i(ji,jj,jl) = rswitch * sm_i(ii,ij,jl) + ( 1.0 - rswitch ) * rn_simin190 s_i (ji,jj,jl) = rswitch * s_i (ii,ij,jl) + ( 1.0 - rswitch ) * rn_simin 191 191 oa_i(ji,jj,jl) = rswitch * oa_i(ii,ij,jl) 192 192 t_su(ji,jj,jl) = rswitch * t_su(ii,ij,jl) + ( 1.0 - rswitch ) * rt0 … … 195 195 END DO 196 196 DO jk = 1, nlay_i 197 t_i (ji,jj,jk,jl) = rswitch * t_i(ii,ij,jk,jl) + ( 1.0 - rswitch ) * rt0198 s _i(ji,jj,jk,jl) = rswitch * s_i(ii,ij,jk,jl) + ( 1.0 - rswitch ) * rn_simin197 t_i (ji,jj,jk,jl) = rswitch * t_i (ii,ij,jk,jl) + ( 1.0 - rswitch ) * rt0 198 sz_i(ji,jj,jk,jl) = rswitch * sz_i(ii,ij,jk,jl) + ( 1.0 - rswitch ) * rn_simin 199 199 END DO 200 200 ! … … 202 202 ! 203 203 IF( nn_icesal == 1 ) THEN ! constant salinity : overwrite rn_icesal 204 s m_i(ji,jj ,jl) = rn_icesal205 s _i(ji,jj,:,jl) = rn_icesal204 s_i (ji,jj ,jl) = rn_icesal 205 sz_i(ji,jj,:,jl) = rn_icesal 206 206 ENDIF 207 207 ! 208 208 ! contents 209 s mv_i(ji,jj,jl) = MIN( sm_i(ji,jj,jl) , sss_m(ji,jj) ) * v_i(ji,jj,jl)209 sv_i(ji,jj,jl) = MIN( s_i(ji,jj,jl) , sss_m(ji,jj) ) * v_i(ji,jj,jl) 210 210 DO jk = 1, nlay_s 211 211 ! Snow energy of melting … … 215 215 END DO 216 216 DO jk = 1, nlay_i 217 ztmelts = - tmut * s _i(ji,jj,jk,jl) + rt0 !Melting temperature in K217 ztmelts = - tmut * sz_i(ji,jj,jk,jl) + rt0 !Melting temperature in K 218 218 ! heat content per unit volume 219 219 e_i(ji,jj,jk,jl) = rswitch * rhoic * & … … 227 227 END DO 228 228 ! 229 CALL lbc_bdy_lnk( 230 CALL lbc_bdy_lnk( 231 CALL lbc_bdy_lnk( 232 CALL lbc_bdy_lnk( 233 CALL lbc_bdy_lnk( 234 ! 235 CALL lbc_bdy_lnk( s mv_i(:,:,jl), 'T', 1., ib_bdy )236 CALL lbc_bdy_lnk( s m_i(:,:,jl), 'T', 1., ib_bdy )237 CALL lbc_bdy_lnk( 238 CALL lbc_bdy_lnk( 229 CALL lbc_bdy_lnk( a_i(:,:,jl), 'T', 1., ib_bdy ) 230 CALL lbc_bdy_lnk( h_i(:,:,jl), 'T', 1., ib_bdy ) 231 CALL lbc_bdy_lnk( h_s(:,:,jl), 'T', 1., ib_bdy ) 232 CALL lbc_bdy_lnk( v_i(:,:,jl), 'T', 1., ib_bdy ) 233 CALL lbc_bdy_lnk( v_s(:,:,jl), 'T', 1., ib_bdy ) 234 ! 235 CALL lbc_bdy_lnk( sv_i(:,:,jl), 'T', 1., ib_bdy ) 236 CALL lbc_bdy_lnk( s_i(:,:,jl), 'T', 1., ib_bdy ) 237 CALL lbc_bdy_lnk( oa_i(:,:,jl), 'T', 1., ib_bdy ) 238 CALL lbc_bdy_lnk( t_su(:,:,jl), 'T', 1., ib_bdy ) 239 239 DO jk = 1, nlay_s 240 240 CALL lbc_bdy_lnk(t_s(:,:,jk,jl), 'T', 1., ib_bdy )
Note: See TracChangeset
for help on using the changeset viewer.