Changeset 9049 for branches/2017
- Timestamp:
- 2017-12-14T14:55:17+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icedyn_rhg_evp.F90
r9019 r9049 112 112 INTEGER :: jter ! local integers 113 113 ! 114 REAL(wp) :: zrhoco 115 REAL(wp) :: zdtevp, z1_dtevp 116 REAL(wp) :: ecc2, z1_ecc2 117 REAL(wp) :: zalph1, z1_alph1, zalph2, z1_alph2 118 REAL(wp) :: zm1, zm2, zm3, zmassU, zmassV 119 REAL(wp) :: zdelta, zp_delf, zds2, zdt, zdt2, zdiv, zdiv2 120 REAL(wp) :: zTauO, zTauB, zTauE, zvel 121 ! 122 REAL(wp) :: zresm 123 REAL(wp) :: zintb, zintn 114 REAL(wp) :: zrhoco ! rau0 * rn_cio 115 REAL(wp) :: zdtevp, z1_dtevp ! time step for subcycling 116 REAL(wp) :: ecc2, z1_ecc2 ! square of yield ellipse eccenticity 117 REAL(wp) :: zalph1, z1_alph1, zalph2, z1_alph2 ! alpha coef from Bouillon 2009 or Kimmritz 2017 118 REAL(wp) :: zm1, zm2, zm3, zmassU, zmassV ! ice/snow mass 119 REAL(wp) :: zdelta, zp_delf, zds2, zdt, zdt2, zdiv, zdiv2 ! temporary scalars 120 REAL(wp) :: zTauO, zTauB, zTauE, zvel ! temporary scalars 121 ! 122 REAL(wp) :: zresm ! Maximal error on ice velocity 123 REAL(wp) :: zintb, zintn ! dummy argument 124 124 REAL(wp) :: zfac_x, zfac_y 125 125 REAL(wp) :: zshear, zdum1, zdum2 … … 142 142 REAL(wp), DIMENSION(jpi,jpj) :: zu_ice, zv_ice, zresr ! check convergence 143 143 REAL(wp), DIMENSION(jpi,jpj) :: zpice ! array used for the calculation of ice surface slope: 144 ! !ocean surface (ssh_m) if ice is not embedded145 ! !ice top surface if ice is embedded144 ! ! ocean surface (ssh_m) if ice is not embedded 145 ! ! ice top surface if ice is embedded 146 146 REAL(wp), DIMENSION(jpi,jpj) :: zCorx, zCory ! Coriolis stress array 147 147 REAL(wp), DIMENSION(jpi,jpj) :: ztaux_oi, ztauy_oi ! Ocean-to-ice stress array … … 810 810 END DO 811 811 812 CALL lbc_lnk_multi( zdiag_sig1 , 'T', 1., zdiag_sig2 , 'T', 1., & 813 & zdiag_dssh_dx, 'U', -1., zdiag_dssh_dy, 'V', -1., & 814 & zdiag_corstrx, 'U', -1., zdiag_corstry, 'V', -1., & 815 & zdiag_intstrx, 'U', -1., zdiag_intstry, 'V', -1. ) 816 817 CALL lbc_lnk_multi( zdiag_utau_oi, 'U', -1., zdiag_vtau_oi, 'V', -1. ) 818 819 CALL lbc_lnk_multi( zdiag_xmtrp_ice, 'U', -1., zdiag_xmtrp_snw, 'U', -1., & 820 & zdiag_xatrp , 'U', -1., zdiag_ymtrp_ice, 'V', -1., & 821 & zdiag_ymtrp_snw, 'V', -1., zdiag_yatrp , 'V', -1. ) 812 CALL lbc_lnk_multi( zdiag_sig1 , 'T', 1., zdiag_sig2 , 'T', 1., & 813 & zdiag_dssh_dx, 'U', -1., zdiag_dssh_dy, 'V', -1., & 814 & zdiag_corstrx, 'U', -1., zdiag_corstry, 'V', -1., & 815 & zdiag_intstrx, 'U', -1., zdiag_intstry, 'V', -1. ) 816 817 CALL lbc_lnk_multi( zdiag_utau_oi , 'U', -1., zdiag_vtau_oi , 'V', -1., & 818 & zdiag_xmtrp_ice, 'U', -1., zdiag_xmtrp_snw, 'U', -1., & 819 & zdiag_xatrp , 'U', -1., zdiag_ymtrp_ice, 'V', -1., & 820 & zdiag_ymtrp_snw, 'V', -1., zdiag_yatrp , 'V', -1. ) 822 821 823 822 IF( iom_use('normstr' ) ) CALL iom_put( 'normstr' , zdiag_sig1(:,:) ) ! Normal stress
Note: See TracChangeset
for help on using the changeset viewer.