Changeset 15334 for NEMO/trunk/src/ICE/icecor.F90
- Timestamp:
- 2021-10-05T23:18:34+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ICE/icecor.F90
r14997 r15334 53 53 INTEGER, INTENT(in) :: kn ! 1 = after dyn ; 2 = after thermo 54 54 ! 55 INTEGER :: ji, jj, j k, jl! dummy loop indices55 INTEGER :: ji, jj, jl ! dummy loop indices 56 56 REAL(wp) :: zsal, zzc 57 57 !!---------------------------------------------------------------------- … … 99 99 END DO 100 100 ENDIF 101 101 ! 102 102 IF( kn /= 0 ) THEN ! no zapsmall if kn=0 (for bdy for instance) because we do not want ice-ocean exchanges (wfx,sfx,hfx) 103 103 ! otherwise conservation diags will fail … … 105 105 CALL ice_var_zapsmall ! Zap small values ! 106 106 ! !----------------------------------------------------- 107 ENDIF108 ! !-----------------------------------------------------109 IF( kn == 2 ) THEN ! Ice drift case: Corrections to avoid wrong values !110 DO_2D( 0, 0, 0, 0 ) !-----------------------------------------------------111 IF ( at_i(ji,jj) == 0._wp ) THEN ! what to do if there is no ice112 IF ( at_i(ji+1,jj) == 0._wp ) u_ice(ji ,jj) = 0._wp ! right side113 IF ( at_i(ji-1,jj) == 0._wp ) u_ice(ji-1,jj) = 0._wp ! left side114 IF ( at_i(ji,jj+1) == 0._wp ) v_ice(ji,jj ) = 0._wp ! upper side115 IF ( at_i(ji,jj-1) == 0._wp ) v_ice(ji,jj-1) = 0._wp ! bottom side116 ENDIF117 END_2D118 CALL lbc_lnk( 'icecor', u_ice, 'U', -1.0_wp, v_ice, 'V', -1.0_wp )119 107 ENDIF 120 108 !
Note: See TracChangeset
for help on using the changeset viewer.