- Timestamp:
- 2011-03-15T16:27:46+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90
r2636 r2690 50 50 !! * Substitutions 51 51 # include "domzgr_substitute.h90" 52 !!---------------------------------------------------------------------- ---52 !!---------------------------------------------------------------------- 53 53 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 54 54 !! $Id$ 55 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 56 !!------------------------------------------------------------------------- 57 55 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 56 !!---------------------------------------------------------------------- 58 57 CONTAINS 59 58 … … 92 91 !! un,vn now horizontal velocity of next time-step 93 92 !!---------------------------------------------------------------------- 94 USE oce, ONLY : ze3u_f => ta ! use ta as 3D workspace 95 USE oce, ONLY : ze3v_f => sa ! use sa as 3D workspace 96 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 97 USE wrk_nemo, ONLY: zs_t => wrk_2d_1, zs_u_1 => wrk_2d_2, & 98 zs_v_1 => wrk_2d_3 93 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 94 USE oce , ONLY: ze3u_f => ta , ze3v_f => sa ! (ta,sa) used as 3D workspace 95 USE wrk_nemo, ONLY: zs_t => wrk_2d_1 , zs_u_1 => wrk_2d_2 , zs_v_1 => wrk_2d_3 96 ! 99 97 INTEGER, INTENT( in ) :: kt ! ocean time-step index 100 ! !98 ! 101 99 INTEGER :: ji, jj, jk ! dummy loop indices 102 100 #if ! defined key_dynspg_flt 103 101 REAL(wp) :: z2dt ! temporary scalar 104 102 #endif 105 REAL(wp) :: zue3a , zue3n , zue3b ! temporary scalar 106 REAL(wp) :: zve3a , zve3n , zve3b ! - - 107 REAL(wp) :: zuf , zvf ! - - 108 REAL(wp) :: zec ! - - 109 REAL(wp) :: zv_t_ij , zv_t_ip1j ! - - 110 REAL(wp) :: zv_t_ijp1 ! - - 103 REAL(wp) :: zue3a, zue3n, zue3b, zuf ! local scalars 104 REAL(wp) :: zve3a, zve3n, zve3b, zvf ! - - 105 REAL(wp) :: zec, zv_t_ij, zv_t_ip1j, zv_t_ijp1 111 106 !!---------------------------------------------------------------------- 112 107 … … 163 158 CALL obc_dyn( kt ) 164 159 ! 165 IF ( lk_dynspg_exp .OR. lk_dynspg_ts) THEN160 IF( .NOT. lk_dynspg_flt ) THEN 166 161 ! Flather boundary condition : - Update sea surface height on each open boundary 167 ! sshn (= after ssh ) for explicit case 168 ! sshn_b (= after ssha_b) for time-splitting case 162 ! sshn (= after ssh ) for explicit case (lk_dynspg_exp=T) 163 ! sshn_b (= after ssha_b) for time-splitting case (lk_dynspg_ts=T) 169 164 ! - Correct the barotropic velocities 170 165 CALL obc_dyn_bt( kt ) … … 180 175 # elif defined key_bdy 181 176 ! !* BDY open boundaries 182 IF( lk_dynspg_exp .OR. lk_dynspg_ts ) THEN ! except for filtered option 183 CALL bdy_dyn_frs( kt ) 184 ENDIF 177 IF( .NOT. lk_dynspg_flt ) CALL bdy_dyn_frs( kt ) 185 178 # endif 186 179 ! … … 325 318 & tab3d_2=vn, clinfo2=' Vn: ' , mask2=vmask ) 326 319 ! 327 IF( wrk_not_released(2, 1,2,3) ) CALL ctl_stop('dyn_nxt: failed to release workspace arrays .')320 IF( wrk_not_released(2, 1,2,3) ) CALL ctl_stop('dyn_nxt: failed to release workspace arrays') 328 321 ! 329 322 END SUBROUTINE dyn_nxt
Note: See TracChangeset
for help on using the changeset viewer.