- Timestamp:
- 2019-12-05T13:18:21+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/UKMO_MERGE_2019/src/OCE/DYN/sshwzv.F90
r11822 r12068 19 19 !!---------------------------------------------------------------------- 20 20 USE oce ! ocean dynamics and tracers variables 21 USE isf ! ice shelf 22 USE isfutils 21 23 USE dom_oce ! ocean space and time domain variables 22 24 USE sbc_oce ! surface boundary condition: ocean … … 256 258 zcoef = atfp * rdt * r1_rau0 257 259 pssh(:,:,Kmm) = pssh(:,:,Kmm) - zcoef * ( emp_b(:,:) - emp (:,:) & 258 & - rnf_b(:,:) + rnf (:,:) & 259 & + fwfisf_b(:,:) - fwfisf(:,:) ) * ssmask(:,:) 260 & - rnf_b(:,:) + rnf (:,:) & 261 & + fwfisf_cav_b(:,:) - fwfisf_cav(:,:) & 262 & + fwfisf_par_b(:,:) - fwfisf_par(:,:) ) * ssmask(:,:) 263 264 ! ice sheet coupling 265 IF ( ln_isf .AND. ln_isfcpl .AND. kt == nit000+1) pssh(:,:,Kbb) = pssh(:,:,Kbb) - atfp * rdt * ( risfcpl_ssh(:,:) - 0.0 ) * ssmask(:,:) 266 260 267 ENDIF 261 268 ENDIF … … 350 357 zCu = MAX( Cu_adv(ji,jj,jk) , Cu_adv(ji,jj,jk-1) ) 351 358 ! alt: 352 ! IF ( w n(ji,jj,jk) > 0._wp ) THEN359 ! IF ( ww(ji,jj,jk) > 0._wp ) THEN 353 360 ! zCu = Cu_adv(ji,jj,jk) 354 361 ! ELSE
Note: See TracChangeset
for help on using the changeset viewer.