- Timestamp:
- 2019-08-22T11:49:08+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/dev_r10037_GPU/src/ICE/icedyn_rhg_evp.F90
r10843 r11467 107 107 !! Kimmritz et al., Ocean Modelling 2016 & 2017 108 108 !!------------------------------------------------------------------- 109 USE scoce, ONLY : z1_e1t0 => scr2D1, z1_e2t0 => scr2D2, & ! scale factors 110 zp_delt => scr2D3, & ! P/delta at T points 111 zbeta => scr2D4, & ! beta coef from Kimmritz 2017 112 ! 113 zdt_m => scr2D5, & ! (dt / ice-snow_mass) on T points 114 zaU => scr2D6 , zaV => scr2D7, & ! ice fraction on U/V points 115 zmU_t => scr2D8, zmV_t => scr2D9, & ! (ice-snow_mass / dt) on U/V points 116 zmf => scr2D10, & ! coriolis parameter at T points 117 zTauU_ia => scr2D11, ztauV_ia => scr2D12, & ! ice-atm. stress at U-V points 118 zspgU => scr2D13, zspgV => scr2D14, & ! surface pressure gradient at U/V points 119 v_oceU => scr2D15, u_oceV => scr2D16, v_iceU => scr2D17, u_iceV => scr2D18, & ! ocean/ice u/v component on V/U points 120 zfU => scr2D19 , zfV => scr2D20, & ! internal stresses 121 ! 122 zds => scr2D21, & ! shear 123 zs1 => scr2D22, zs2 => scr2D23, zs12 => scr2D38, & ! stress tensor components 124 zu_ice => scr2D24, zv_ice => scr2D25, zresr => scr2D26, & ! check convergence 125 zpice => scr2D27, & ! array used for the calculation of ice surface slope: 126 ! ! ocean surface (ssh_m) if ice is not embedded 127 ! ! ice top surface if ice is embedded 128 zCorx => scr2D28, zCory => scr2D29, & ! Coriolis stress array 129 ztaux_oi => scr2D30, ztauy_oi => scr2D31, & ! Ocean-to-ice stress array 130 ! 131 zswitchU => scr2D32, zswitchV => scr2D33, & ! dummy arrays 132 zmaskU => scr2D34, zmaskV => scr2D35, & ! mask for ice presence 133 zfmask => scr2D36, zwf => scr2D37 ! mask at F points for the ice 134 135 136 109 137 INTEGER , INTENT(in ) :: kt ! time step 110 138 REAL(wp), DIMENSION(:,:), INTENT(inout) :: pstress1_i, pstress2_i, pstress12_i ! … … 127 155 REAL(wp) :: zshear, zdum1, zdum2 128 156 ! 129 REAL(wp), DIMENSION(jpi,jpj) :: z1_e1t0, z1_e2t0 ! scale factors130 REAL(wp), DIMENSION(jpi,jpj) :: zp_delt ! P/delta at T points131 REAL(wp), DIMENSION(jpi,jpj) :: zbeta ! beta coef from Kimmritz 2017132 !133 REAL(wp), DIMENSION(jpi,jpj) :: zdt_m ! (dt / ice-snow_mass) on T points134 REAL(wp), DIMENSION(jpi,jpj) :: zaU , zaV ! ice fraction on U/V points135 REAL(wp), DIMENSION(jpi,jpj) :: zmU_t, zmV_t ! (ice-snow_mass / dt) on U/V points136 REAL(wp), DIMENSION(jpi,jpj) :: zmf ! coriolis parameter at T points137 REAL(wp), DIMENSION(jpi,jpj) :: zTauU_ia , ztauV_ia ! ice-atm. stress at U-V points138 REAL(wp), DIMENSION(jpi,jpj) :: zspgU , zspgV ! surface pressure gradient at U/V points139 REAL(wp), DIMENSION(jpi,jpj) :: v_oceU, u_oceV, v_iceU, u_iceV ! ocean/ice u/v component on V/U points140 REAL(wp), DIMENSION(jpi,jpj) :: zfU , zfV ! internal stresses141 !142 REAL(wp), DIMENSION(jpi,jpj) :: zds ! shear143 REAL(wp), DIMENSION(jpi,jpj) :: zs1, zs2, zs12 ! stress tensor components144 REAL(wp), DIMENSION(jpi,jpj) :: zu_ice, zv_ice, zresr ! check convergence145 REAL(wp), DIMENSION(jpi,jpj) :: zpice ! array used for the calculation of ice surface slope:146 ! ! ocean surface (ssh_m) if ice is not embedded147 ! ! ice top surface if ice is embedded148 REAL(wp), DIMENSION(jpi,jpj) :: zCorx, zCory ! Coriolis stress array149 REAL(wp), DIMENSION(jpi,jpj) :: ztaux_oi, ztauy_oi ! Ocean-to-ice stress array150 !151 REAL(wp), DIMENSION(jpi,jpj) :: zswitchU, zswitchV ! dummy arrays152 REAL(wp), DIMENSION(jpi,jpj) :: zmaskU, zmaskV ! mask for ice presence153 REAL(wp), DIMENSION(jpi,jpj) :: zfmask, zwf ! mask at F points for the ice154 155 157 REAL(wp), PARAMETER :: zepsi = 1.0e-20_wp ! tolerance parameter 156 158 REAL(wp), PARAMETER :: zmmin = 1._wp ! ice mass (kg/m2) below which ice velocity equals ocean velocity
Note: See TracChangeset
for help on using the changeset viewer.