- Timestamp:
- 2019-09-18T16:11:52+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10973_AGRIF-01_jchanut_small_jpi_jpj/src/OCE/step.F90
r10364 r11564 112 112 IF( ln_tide ) CALL sbc_tide( kstp ) ! update tide potential 113 113 IF( ln_apr_dyn ) CALL sbc_apr ( kstp ) ! atmospheric pressure (NB: call before bdy_dta which needs ssh_ib) 114 IF( ln_bdy ) CALL bdy_dta ( kstp, time_offset=+1 ) ! update dynamic & tracer data at open boundaries114 IF( ln_bdy ) CALL bdy_dta ( kstp, kt_offset = +1 ) ! update dynamic & tracer data at open boundaries 115 115 CALL sbc ( kstp ) ! Sea Boundary Condition (including sea-ice) 116 116 … … 165 165 CALL eos ( tsn, rhd, rhop, gdept_n(:,:,:) ) ! now in situ density for hpg computation 166 166 167 !!jc: fs simplification168 !!jc: lines below are useless if ln_linssh=F. Keep them here (which maintains a bug if ln_linssh=T and ln_zps=T, cf ticket #1636)169 !! but ensures reproductible results170 !! with previous versions using split-explicit free surface171 IF( ln_zps .AND. .NOT. ln_isfcav ) &172 & CALL zps_hde ( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps: before horizontal gradient173 & rhd, gru , grv ) ! of t, s, rd at the last ocean level174 IF( ln_zps .AND. ln_isfcav ) &175 & CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv, gtui, gtvi, & ! Partial steps for top cell (ISF)176 & rhd, gru , grv , grui, grvi ) ! of t, s, rd at the first ocean level177 !!jc: fs simplification178 167 179 168 ua(:,:,:) = 0._wp ! set dynamics trends to zero … … 198 187 CALL div_hor ( kstp ) ! Horizontal divergence (2nd call in time-split case) 199 188 IF(.NOT.ln_linssh) CALL dom_vvl_sf_nxt( kstp, kcall=2 ) ! after vertical scale factors (update depth average component) 189 ENDIF 190 CALL dyn_zdf ( kstp ) ! vertical diffusion 191 192 IF( ln_dynspg_ts ) THEN 200 193 CALL wzv ( kstp ) ! now cross-level velocity 201 194 IF( ln_zad_Aimp ) CALL wAimp ( kstp ) ! Adaptive-implicit vertical advection partitioning 202 195 ENDIF 203 204 CALL dyn_zdf ( kstp ) ! vertical diffusion205 196 206 197 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 212 203 ! diagnostics and outputs 213 204 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 214 IF( l k_floats ) CALL flo_stp ( kstp ) ! drifting Floats205 IF( ln_floats ) CALL flo_stp ( kstp ) ! drifting Floats 215 206 IF( ln_diacfl ) CALL dia_cfl ( kstp ) ! Courant number diagnostics 216 207 IF( lk_diahth ) CALL dia_hth ( kstp ) ! Thermocline depth (20 degres isotherm depth) 217 IF( l k_diadct ) CALL dia_dct ( kstp ) ! Transports208 IF( ln_diadct ) CALL dia_dct ( kstp ) ! Transports 218 209 CALL dia_ar5 ( kstp ) ! ar5 diag 219 IF( l k_diaharm ) CALL dia_harm( kstp ) ! Tidal harmonic analysis210 IF( ln_diaharm ) CALL dia_harm( kstp ) ! Tidal harmonic analysis 220 211 CALL dia_wri ( kstp ) ! ocean model: outputs 221 212 !
Note: See TracChangeset
for help on using the changeset viewer.