Changeset 9935 for NEMO/trunk/src/ICE/icedyn_rhg_evp.F90
- Timestamp:
- 2018-07-12T16:12:48+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ICE/icedyn_rhg_evp.F90
r9866 r9935 285 285 286 286 ! Ice/snow mass at U-V points 287 zm1 = ( rhos n * vt_s(ji ,jj ) + rhoic* vt_i(ji ,jj ) )288 zm2 = ( rhos n * vt_s(ji+1,jj ) + rhoic* vt_i(ji+1,jj ) )289 zm3 = ( rhos n * vt_s(ji ,jj+1) + rhoic* vt_i(ji ,jj+1) )287 zm1 = ( rhos * vt_s(ji ,jj ) + rhoi * vt_i(ji ,jj ) ) 288 zm2 = ( rhos * vt_s(ji+1,jj ) + rhoi * vt_i(ji+1,jj ) ) 289 zm3 = ( rhos * vt_s(ji ,jj+1) + rhoi * vt_i(ji ,jj+1) ) 290 290 zmassU = 0.5_wp * ( zm1 * e1e2t(ji,jj) + zm2 * e1e2t(ji+1,jj) ) * r1_e1e2u(ji,jj) * umask(ji,jj,1) 291 291 zmassV = 0.5_wp * ( zm1 * e1e2t(ji,jj) + zm3 * e1e2t(ji,jj+1) ) * r1_e1e2v(ji,jj) * vmask(ji,jj,1) … … 795 795 zfac_y = 0.5 * v_ice(ji,jj) * e1v(ji,jj) * rswitch 796 796 797 zdiag_xmtrp_ice(ji,jj) = rhoi c* zfac_x * ( vt_i(ji+1,jj) + vt_i(ji,jj) ) ! ice mass transport, X-component798 zdiag_ymtrp_ice(ji,jj) = rhoi c* zfac_y * ( vt_i(ji,jj+1) + vt_i(ji,jj) ) ! '' Y- ''799 800 zdiag_xmtrp_snw(ji,jj) = rhos n* zfac_x * ( vt_s(ji+1,jj) + vt_s(ji,jj) ) ! snow mass transport, X-component801 zdiag_ymtrp_snw(ji,jj) = rhos n* zfac_y * ( vt_s(ji,jj+1) + vt_s(ji,jj) ) ! '' Y- ''802 803 zdiag_xatrp(ji,jj) = zfac_x * ( at_i(ji+1,jj) + at_i(ji,jj) ) 804 zdiag_yatrp(ji,jj) = zfac_y * ( at_i(ji,jj+1) + at_i(ji,jj) ) 797 zdiag_xmtrp_ice(ji,jj) = rhoi * zfac_x * ( vt_i(ji+1,jj) + vt_i(ji,jj) ) ! ice mass transport, X-component 798 zdiag_ymtrp_ice(ji,jj) = rhoi * zfac_y * ( vt_i(ji,jj+1) + vt_i(ji,jj) ) ! '' Y- '' 799 800 zdiag_xmtrp_snw(ji,jj) = rhos * zfac_x * ( vt_s(ji+1,jj) + vt_s(ji,jj) ) ! snow mass transport, X-component 801 zdiag_ymtrp_snw(ji,jj) = rhos * zfac_y * ( vt_s(ji,jj+1) + vt_s(ji,jj) ) ! '' Y- '' 802 803 zdiag_xatrp(ji,jj) = zfac_x * ( at_i(ji+1,jj) + at_i(ji,jj) ) ! area transport, X-component 804 zdiag_yatrp(ji,jj) = zfac_y * ( at_i(ji,jj+1) + at_i(ji,jj) ) ! '' Y- '' 805 805 806 806 END DO
Note: See TracChangeset
for help on using the changeset viewer.