Changeset 12761 for NEMO/branches/2020
- Timestamp:
- 2020-04-17T10:28:12+02:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DYN/sshwzv.F90
r12724 r12761 229 229 230 230 231 SUBROUTINE ssh_atf( kt, Kbb, Kmm, Kaa, pssh )231 SUBROUTINE ssh_atf( kt, Kbb, Kmm, Kaa, pssh, pssh_f ) 232 232 !!---------------------------------------------------------------------- 233 233 !! *** ROUTINE ssh_atf *** … … 244 244 !! Reference : Leclair, M., and G. Madec, 2009, Ocean Modelling. 245 245 !!---------------------------------------------------------------------- 246 INTEGER , INTENT(in ) :: kt ! ocean time-step index 247 INTEGER , INTENT(in ) :: Kbb, Kmm, Kaa ! ocean time level indices 248 REAL(wp), DIMENSION(jpi,jpj,jpt), INTENT(inout) :: pssh ! SSH field 246 INTEGER , INTENT(in ) :: kt ! ocean time-step index 247 INTEGER , INTENT(in ) :: Kbb, Kmm, Kaa ! ocean time level indices 248 REAL(wp), DIMENSION(jpi,jpj,jpt) , TARGET, INTENT(inout) :: pssh ! SSH field 249 REAL(wp), DIMENSION(jpi,jpj ), OPTIONAL, TARGET, INTENT( out) :: pssh_f ! filtered SSH field 249 250 ! 250 251 REAL(wp) :: zcoef ! local scalar 252 REAL(wp), POINTER, DIMENSION(:,:) :: zssh ! pointer for filtered SSH 251 253 !!---------------------------------------------------------------------- 252 254 ! … … 260 262 ! !== Euler time-stepping: no filter, just swap ==! 261 263 IF ( .NOT.( l_1st_euler ) ) THEN ! Only do time filtering for leapfrog timesteps 264 IF( PRESENT( pssh_f ) ) THEN ; zssh => pssh_f 265 ELSE ; zssh => pssh(:,:,Kmm) 266 ENDIF 262 267 ! ! filtered "now" field 263 268 pssh(:,:,Kmm) = pssh(:,:,Kmm) + rn_atfp * ( pssh(:,:,Kbb) - 2 * pssh(:,:,Kmm) + pssh(:,:,Kaa) ) … … 281 286 END SUBROUTINE ssh_atf 282 287 288 283 289 SUBROUTINE wAimp( kt, Kmm ) 284 290 !!---------------------------------------------------------------------- -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/oce.F90
r12724 r12761 32 32 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: Cu_adv !: vertical Courant number (adaptive-implicit) 33 33 34 !! free surface ! before ! now ! after !35 !! ------------ ! fields ! fields ! fields !34 !! free surface 35 !! ------------ 36 36 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: ssh, uu_b, vv_b !: SSH [m] and barotropic velocities [m/s] 37 37 !REAL(wp), PUBLIC, ALLOCATABLE, SAVE, TARGET, DIMENSION(:,:,:) :: ssh !: SSH [m] 38 !REAL(wp), PUBLIC, ALLOCATABLE, SAVE , DIMENSION(:,:,:) :: uu_b, vv_b !: barotropic velocities [m/s] 39 38 40 !! Arrays at barotropic time step: ! befbefore! before ! now ! after ! 39 41 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ubb_e , ub_e , un_e , ua_e !: u-external velocity -
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/stpMLF.F90
r12732 r12761 91 91 INTEGER :: kcall ! optional integer argument (dom_vvl_sf_nxt) 92 92 !!st patch 93 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgdept 93 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgdept 94 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zssh_f 94 95 !! --------------------------------------------------------------------- 95 96 #if defined key_agrif … … 248 249 ! Passive Tracer Model 249 250 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 250 CALL trc_stp ( kstp, Nbb, Nnn, Nrhs, Naa ) ! time-stepping 251 ALLOCATE( zssh_f(jpi,jpj) ) 252 CALL ssh_atf ( kstp, Nbb, Nnn, Naa , ssh, zssh_f ) ! time filtering of "now" sea surface height 253 CALL dom_qco_r3c( zssh_f, r3t_f, r3u_f, r3v_f ) ! "now" ssh/h_0 ratio from filtrered ssh 254 ! 255 CALL trc_stp ( kstp, Nbb, Nnn, Nrhs, Naa ) ! time-stepping 256 DEALLOCATE( zssh_f ) 251 257 #endif 252 258 … … 297 303 CALL finalize_sbc ( kstp, Nbb, Naa, uu, vv, ts ) ! boundary condifions 298 304 CALL ssh_atf ( kstp, Nbb, Nnn, Naa, ssh ) ! time filtering of "now" sea surface height 299 CALL dom_qco_r3c 305 CALL dom_qco_r3c ( ssh(:,:,Nnn), r3t_f, r3u_f, r3v_f ) ! "now" ssh/h_0 ratio from filtrered ssh 300 306 CALL tra_atf_qco ( kstp, Nbb, Nnn, Naa, ts ) ! time filtering of "now" tracer arrays 301 307 CALL dyn_atf_qco ( kstp, Nbb, Nnn, Naa, uu, vv ) ! time filtering of "now" velocities and scale factors
Note: See TracChangeset
for help on using the changeset viewer.