Changeset 11466
- Timestamp:
- 2019-08-21T09:18:21+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/OCE/DYN/sshwzv.F90
r11426 r11466 226 226 !! ** Purpose : achieve the sea surface height time stepping by 227 227 !! applying Asselin time filter and swapping the arrays 228 !! ssh(:,:,Kaa) already computed in ssh_nxt228 !! pssh(:,:,Kaa) already computed in ssh_nxt 229 229 !! 230 230 !! ** Method : - apply Asselin time fiter to now ssh (excluding the forcing 231 231 !! from the filter, see Leclair and Madec 2010) and swap : 232 !! ssh(:,:,Kmm) = ssh(:,:,Kaa) + atfp * ( ssh(:,:,Kbb) -2 ssh(:,:,Kmm) +ssh(:,:,Kaa) )232 !! pssh(:,:,Kmm) = pssh(:,:,Kaa) + atfp * ( pssh(:,:,Kbb) -2 pssh(:,:,Kmm) + pssh(:,:,Kaa) ) 233 233 !! - atfp * rdt * ( emp_b - emp ) / rau0 234 !! ssh(:,:,Kmm) =ssh(:,:,Kaa)235 !! 236 !! ** action : - ssh(:,:,Kbb),ssh(:,:,Kmm) : before & now sea surface height234 !! pssh(:,:,Kmm) = pssh(:,:,Kaa) 235 !! 236 !! ** action : - pssh(:,:,Kbb), pssh(:,:,Kmm) : before & now sea surface height 237 237 !! ready for the next time step 238 238 !! … … 256 256 IF ( .NOT.( neuler == 0 .AND. kt == nit000 ) ) THEN ! Only do time filtering for leapfrog timesteps 257 257 ! ! filtered "now" field 258 ssh(:,:,Kmm) = ssh(:,:,Kmm) + atfp * ( ssh(:,:,Kbb) - 2 * ssh(:,:,Kmm) +ssh(:,:,Kaa) )258 pssh(:,:,Kmm) = pssh(:,:,Kmm) + atfp * ( pssh(:,:,Kbb) - 2 * pssh(:,:,Kmm) + pssh(:,:,Kaa) ) 259 259 IF( .NOT.ln_linssh ) THEN ! "now" <-- with forcing removed 260 260 zcoef = atfp * rdt * r1_rau0 261 ssh(:,:,Kmm) =ssh(:,:,Kmm) - zcoef * ( emp_b(:,:) - emp (:,:) &261 pssh(:,:,Kmm) = pssh(:,:,Kmm) - zcoef * ( emp_b(:,:) - emp (:,:) & 262 262 & - rnf_b(:,:) + rnf (:,:) & 263 263 & + fwfisf_b(:,:) - fwfisf(:,:) ) * ssmask(:,:) … … 265 265 ENDIF 266 266 ! 267 IF(ln_ctl) CALL prt_ctl( tab2d_1= ssh(:,:,Kmm), clinfo1='ssh(:,:,Kmm) - : ', mask1=tmask )267 IF(ln_ctl) CALL prt_ctl( tab2d_1=pssh(:,:,Kmm), clinfo1=' pssh(:,:,Kmm) - : ', mask1=tmask ) 268 268 ! 269 269 IF( ln_timing ) CALL timing_stop('ssh_atf')
Note: See TracChangeset
for help on using the changeset viewer.