- 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.F90
r8563 r8564 295 295 DO jk = 1, nlay_i 296 296 DO ji = 1, nidx 297 ztmelts = -tmut * s _i_1d(ji,jk)297 ztmelts = -tmut * sz_i_1d(ji,jk) 298 298 ! Conversion q(S,T) -> T (second order equation) 299 299 zbbb = ( rcp - cpic ) * ztmelts + e_i_1d(ji,jk) * r1_rhoic - lfus … … 365 365 CALL tab_2d_1d( nidx, idxice(1:nidx), h_s_1d(1:nidx), h_s(:,:,kl) ) 366 366 CALL tab_2d_1d( nidx, idxice(1:nidx), t_su_1d(1:nidx), t_su(:,:,kl) ) 367 CALL tab_2d_1d( nidx, idxice(1:nidx), s m_i_1d(1:nidx), sm_i(:,:,kl) )367 CALL tab_2d_1d( nidx, idxice(1:nidx), s_i_1d(1:nidx), s_i(:,:,kl) ) 368 368 DO jk = 1, nlay_s 369 369 CALL tab_2d_1d( nidx, idxice(1:nidx), t_s_1d(1:nidx,jk), t_s(:,:,jk,kl) ) … … 373 373 CALL tab_2d_1d( nidx, idxice(1:nidx), t_i_1d(1:nidx,jk), t_i(:,:,jk,kl) ) 374 374 CALL tab_2d_1d( nidx, idxice(1:nidx), e_i_1d(1:nidx,jk), e_i(:,:,jk,kl) ) 375 CALL tab_2d_1d( nidx, idxice(1:nidx), s _i_1d(1:nidx,jk), s_i(:,:,jk,kl) )375 CALL tab_2d_1d( nidx, idxice(1:nidx), sz_i_1d(1:nidx,jk), sz_i(:,:,jk,kl) ) 376 376 END DO 377 377 ! … … 454 454 ! 455 455 ! Change thickness to volume (replaces routine ice_var_eqv2glo) 456 v_i_1d(1:nidx) 457 v_s_1d(1:nidx) 458 s mv_i_1d(1:nidx) = sm_i_1d(1:nidx) * v_i_1d(1:nidx)456 v_i_1d(1:nidx) = h_i_1d(1:nidx) * a_i_1d(1:nidx) 457 v_s_1d(1:nidx) = h_s_1d(1:nidx) * a_i_1d(1:nidx) 458 sv_i_1d(1:nidx) = s_i_1d(1:nidx) * v_i_1d(1:nidx) 459 459 460 460 CALL tab_1d_2d( nidx, idxice(1:nidx), at_i_1d(1:nidx), at_i ) … … 463 463 CALL tab_1d_2d( nidx, idxice(1:nidx), h_s_1d(1:nidx), h_s(:,:,kl) ) 464 464 CALL tab_1d_2d( nidx, idxice(1:nidx), t_su_1d(1:nidx), t_su(:,:,kl) ) 465 CALL tab_1d_2d( nidx, idxice(1:nidx), s m_i_1d(1:nidx), sm_i(:,:,kl) )465 CALL tab_1d_2d( nidx, idxice(1:nidx), s_i_1d(1:nidx), s_i(:,:,kl) ) 466 466 DO jk = 1, nlay_s 467 467 CALL tab_1d_2d( nidx, idxice(1:nidx), t_s_1d(1:nidx,jk), t_s(:,:,jk,kl) ) … … 471 471 CALL tab_1d_2d( nidx, idxice(1:nidx), t_i_1d(1:nidx,jk), t_i(:,:,jk,kl) ) 472 472 CALL tab_1d_2d( nidx, idxice(1:nidx), e_i_1d(1:nidx,jk), e_i(:,:,jk,kl) ) 473 CALL tab_1d_2d( nidx, idxice(1:nidx), s _i_1d(1:nidx,jk), s_i(:,:,jk,kl) )473 CALL tab_1d_2d( nidx, idxice(1:nidx), sz_i_1d(1:nidx,jk), sz_i(:,:,jk,kl) ) 474 474 END DO 475 475 ! … … 520 520 CALL tab_1d_2d( nidx, idxice(1:nidx), diag_fc_su_1d(1:nidx), diag_fc_su ) 521 521 ! extensive variables 522 CALL tab_1d_2d( nidx, idxice(1:nidx), v_i_1d (1:nidx), v_i(:,:,kl) )523 CALL tab_1d_2d( nidx, idxice(1:nidx), v_s_1d (1:nidx), v_s(:,:,kl) )524 CALL tab_1d_2d( nidx, idxice(1:nidx), s mv_i_1d(1:nidx), smv_i(:,:,kl) )522 CALL tab_1d_2d( nidx, idxice(1:nidx), v_i_1d (1:nidx), v_i (:,:,kl) ) 523 CALL tab_1d_2d( nidx, idxice(1:nidx), v_s_1d (1:nidx), v_s (:,:,kl) ) 524 CALL tab_1d_2d( nidx, idxice(1:nidx), sv_i_1d(1:nidx), sv_i(:,:,kl) ) 525 525 ! 526 526 END SELECT
Note: See TracChangeset
for help on using the changeset viewer.