- Timestamp:
- 2020-02-24T14:00:21+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/src/OCE/DYN/sshwzv.F90
r12424 r12443 88 88 ENDIF 89 89 ! 90 zcoef = 0.5_wp * r1_r au090 zcoef = 0.5_wp * r1_rho0 91 91 92 92 ! !------------------------------! … … 222 222 !! ** Method : - apply Asselin time fiter to now ssh (excluding the forcing 223 223 !! from the filter, see Leclair and Madec 2010) and swap : 224 !! pssh(:,:,Kmm) = pssh(:,:,Kaa) + atfp * ( pssh(:,:,Kbb) -2 pssh(:,:,Kmm) + pssh(:,:,Kaa) )225 !! - atfp * rn_Dt * ( emp_b - emp ) / rau0224 !! pssh(:,:,Kmm) = pssh(:,:,Kaa) + rn_atfp * ( pssh(:,:,Kbb) -2 pssh(:,:,Kmm) + pssh(:,:,Kaa) ) 225 !! - rn_atfp * rn_Dt * ( emp_b - emp ) / rho0 226 226 !! 227 227 !! ** action : - pssh(:,:,Kmm) time filtered … … 246 246 IF ( .NOT.( l_1st_euler ) ) THEN ! Only do time filtering for leapfrog timesteps 247 247 ! ! filtered "now" field 248 pssh(:,:,Kmm) = pssh(:,:,Kmm) + atfp * ( pssh(:,:,Kbb) - 2 * pssh(:,:,Kmm) + pssh(:,:,Kaa) )248 pssh(:,:,Kmm) = pssh(:,:,Kmm) + rn_atfp * ( pssh(:,:,Kbb) - 2 * pssh(:,:,Kmm) + pssh(:,:,Kaa) ) 249 249 IF( .NOT.ln_linssh ) THEN ! "now" <-- with forcing removed 250 zcoef = atfp * rn_Dt * r1_rau0250 zcoef = rn_atfp * rn_Dt * r1_rho0 251 251 pssh(:,:,Kmm) = pssh(:,:,Kmm) - zcoef * ( emp_b(:,:) - emp (:,:) & 252 252 & - rnf_b(:,:) + rnf (:,:) & … … 255 255 256 256 ! ice sheet coupling 257 IF ( ln_isf .AND. ln_isfcpl .AND. kt == nit000+1) pssh(:,:,Kbb) = pssh(:,:,Kbb) - atfp * rn_Dt * ( risfcpl_ssh(:,:) - 0.0 ) * ssmask(:,:)257 IF ( ln_isf .AND. ln_isfcpl .AND. kt == nit000+1) pssh(:,:,Kbb) = pssh(:,:,Kbb) - rn_atfp * rn_Dt * ( risfcpl_ssh(:,:) - 0.0 ) * ssmask(:,:) 258 258 259 259 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.