- Timestamp:
- 2018-07-13T09:28:50+02:00 (6 years ago)
- Location:
- NEMO/branches/2018/dev_r9838_ENHANCE04_RK3
- Files:
-
- 1 edited
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r9838_ENHANCE04_RK3/src/ICE/icedyn_rhg_evp.F90
r9660 r9939 114 114 INTEGER :: jter ! local integers 115 115 ! 116 REAL(wp) :: zrhoco ! r au0 * rn_cio116 REAL(wp) :: zrhoco ! rho0 * rn_cio 117 117 REAL(wp) :: zdtevp, z1_dtevp ! time step for subcycling 118 118 REAL(wp) :: ecc2, z1_ecc2 ! square of yield ellipse eccenticity … … 221 221 ! 1) define some variables and initialize arrays 222 222 !------------------------------------------------------------------------------! 223 zrhoco = r au0 * rn_cio223 zrhoco = rho0 * rn_cio 224 224 225 225 ! ecc2: square of yield ellipse eccenticrity … … 271 271 zintb = REAL( nn_fsbc + 1 ) / REAL( nn_fsbc ) * 0.5_wp 272 272 ! 273 zpice(:,:) = ssh_m(:,:) + ( zintn * snwice_mass(:,:) + zintb * snwice_mass_b(:,:) ) * r1_r au0273 zpice(:,:) = ssh_m(:,:) + ( zintn * snwice_mass(:,:) + zintb * snwice_mass_b(:,:) ) * r1_rho0 274 274 ! 275 275 ELSE !== non-embedded sea ice: use ocean surface for slope calculation ==! … … 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) … … 799 799 zfac_y = 0.5 * v_ice(ji,jj) * e1v(ji,jj) * rswitch 800 800 801 zdiag_xmtrp_ice(ji,jj) = rhoi c * zfac_x * ( vt_i(ji+1,jj) + vt_i(ji,jj) )! ice mass transport, X-component802 zdiag_ymtrp_ice(ji,jj) = rhoi c * zfac_y * ( vt_i(ji,jj+1) + vt_i(ji,jj) )! '' Y- ''803 804 zdiag_xmtrp_snw(ji,jj) = rhos n * zfac_x * ( vt_s(ji+1,jj) + vt_s(ji,jj) )! snow mass transport, X-component805 zdiag_ymtrp_snw(ji,jj) = rhos n * zfac_y * ( vt_s(ji,jj+1) + vt_s(ji,jj) )! '' Y- ''806 807 zdiag_xatrp(ji,jj) = zfac_x * ( at_i(ji+1,jj) + at_i(ji,jj) ) ! area transport, X-component808 zdiag_yatrp(ji,jj) = zfac_y * ( at_i(ji,jj+1) + at_i(ji,jj) ) ! '' Y- ''801 zdiag_xmtrp_ice(ji,jj) = rhoi * zfac_x * ( vt_i(ji+1,jj) + vt_i(ji,jj) ) ! ice mass transport, X-component 802 zdiag_ymtrp_ice(ji,jj) = rhoi * zfac_y * ( vt_i(ji,jj+1) + vt_i(ji,jj) ) ! '' Y- '' 803 804 zdiag_xmtrp_snw(ji,jj) = rhos * zfac_x * ( vt_s(ji+1,jj) + vt_s(ji,jj) ) ! snow mass transport, X-component 805 zdiag_ymtrp_snw(ji,jj) = rhos * zfac_y * ( vt_s(ji,jj+1) + vt_s(ji,jj) ) ! '' Y- '' 806 807 zdiag_xatrp(ji,jj) = zfac_x * ( at_i(ji+1,jj) + at_i(ji,jj) ) ! area transport, X-component 808 zdiag_yatrp(ji,jj) = zfac_y * ( at_i(ji,jj+1) + at_i(ji,jj) ) ! '' Y- '' 809 809 810 810 END DO
Note: See TracChangeset
for help on using the changeset viewer.