- Timestamp:
- 2015-11-20T09:39:06+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/OPA_SRC/step.F90
r5620 r5901 24 24 !! - ! 2010-10 (C. Ethe, G. Madec) reorganisation of initialisation phase + merge TRC-TRA 25 25 !! 3.4 ! 2011-04 (G. Madec, C. Ethe) Merge of dtatem and dtasal 26 !! ! 2012-07 (J. Simeon, G. Madec, C. Ethe) Online coarsening of outputs 27 !! 3.7 ! 2014-04 (F. Roquet, G. Madec) New equations of state 26 !! 3.6 ! 2012-07 (J. Simeon, G. Madec. C. Ethe) Online coarsening of outputs 27 !! 3.6 ! 2014-04 (F. Roquet, G. Madec) New equations of state 28 !! 3.7 ! 2014-10 (G. Madec) LDF simplication 29 !! - ! 2014-12 (G. Madec) remove KPP scheme 28 30 !!---------------------------------------------------------------------- 29 31 … … 37 39 PRIVATE 38 40 39 PUBLIC stp ! called by opa.F9041 PUBLIC stp ! called by nemogcm.F90 40 42 41 43 !! * Substitutions … … 43 45 !!gm # include "zdfddm_substitute.h90" 44 46 !!---------------------------------------------------------------------- 45 !! NEMO/OPA 3.7 , NEMO Consortium (201 4)47 !! NEMO/OPA 3.7 , NEMO Consortium (2015) 46 48 !! $Id$ 47 49 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) … … 50 52 51 53 #if defined key_agrif 52 SUBROUTINE stp( )54 RECURSIVE SUBROUTINE stp( ) 53 55 INTEGER :: kstp ! ocean time-step index 54 56 #else … … 69 71 !! -5- Compute the momentum trends 70 72 !! -6- Update the horizontal velocity 71 !! -7- Compute the diagnostics variables (rd,N2, div,cur,w)73 !! -7- Compute the diagnostics variables (rd,N2, hdiv,w) 72 74 !! -8- Outputs and diagnostics 73 75 !!---------------------------------------------------------------------- … … 76 78 INTEGER :: kcall ! optional integer argument (dom_vvl_sf_nxt) 77 79 !! --------------------------------------------------------------------- 78 79 80 #if defined key_agrif 80 81 kstp = nit000 + Agrif_Nb_Step() 81 ! IF ( Agrif_Root() .and. lwp) Write(*,*) '---' 82 ! IF (lwp) Write(*,*) 'Grid Number',Agrif_Fixed(),' time step ',kstp 83 IF ( kstp == (nit000 + 1) ) lk_agrif_fstep = .FALSE. 82 IF( lk_agrif_debug ) THEN 83 IF( Agrif_Root() .and. lwp) WRITE(*,*) '---' 84 IF(lwp) WRITE(*,*) 'Grid Number', Agrif_Fixed(),' time step ', kstp, 'int tstep', Agrif_NbStepint() 85 ENDIF 86 IF( kstp == nit000 + 1 ) lk_agrif_fstep = .FALSE. 84 87 # if defined key_iomput 85 88 IF( Agrif_Nbstepint() == 0 ) CALL iom_swap( cxios_context ) 86 89 # endif 87 90 #endif 88 indic = 0 ! reset to no error condition 89 IF( kstp == nit000 ) THEN 90 ! must be done after nemo_init for AGRIF+XIOS+OASIS 91 CALL iom_init( cxios_context ) ! iom_put initialization 92 IF( ln_crs ) CALL iom_init( TRIM(cxios_context)//"_crs" ) ! initialize context for coarse grid 93 ENDIF 94 91 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 92 ! update I/O and calendar 93 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 94 indic = 0 ! reset to no error condition 95 96 IF( kstp == nit000 ) THEN ! initialize IOM context (must be done after nemo_init for AGRIF+XIOS+OASIS) 97 CALL iom_init( cxios_context ) ! for model grid (including passible AGRIF zoom) 98 IF( ln_crs ) CALL iom_init( TRIM(cxios_context)//"_crs" ) ! for coarse grid 99 ENDIF 95 100 IF( kstp /= nit000 ) CALL day( kstp ) ! Calendar (day was already called at nit000 in day_init) 96 CALL iom_setkt( kstp - nit000 + 1, cxios_context ) ! tell iom we are at time step kstp 97 IF( ln_crs ) CALL iom_setkt( kstp - nit000 + 1, TRIM(cxios_context)//"_crs" ) ! tell iom we are at time step kstp 98 99 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 100 ! Update data, open boundaries, surface boundary condition (including sea-ice) 101 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 102 IF( lk_tide ) CALL sbc_tide( kstp ) 103 IF( lk_bdy ) THEN 104 IF( ln_apr_dyn) CALL sbc_apr( kstp ) ! bdy_dta needs ssh_ib 105 CALL bdy_dta ( kstp, time_offset=+1 ) ! update dynamic & tracer data at open boundaries 106 ENDIF 107 CALL sbc ( kstp ) ! Sea Boundary Condition (including sea-ice) 108 ! clem: moved here for bdy ice purpose 101 CALL iom_setkt( kstp - nit000 + 1, cxios_context ) ! tell IOM we are at time step kstp 102 IF( ln_crs ) CALL iom_setkt( kstp - nit000 + 1, TRIM(cxios_context)//"_crs" ) ! tell IOM we are at time step kstp 103 104 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 105 ! Update external forcing (tides, open boundaries, and surface boundary condition (including sea-ice) 106 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 107 IF( lk_tide ) CALL sbc_tide( kstp ) ! update tide potential 108 IF( ln_apr_dyn ) CALL sbc_apr ( kstp ) ! atmospheric pressure (NB: call before bdy_dta which needs ssh_ib) 109 IF( lk_bdy ) CALL bdy_dta ( kstp, time_offset=+1 ) ! update dynamic & tracer data at open boundaries 110 CALL sbc ( kstp ) ! Sea Boundary Condition (including sea-ice) 111 109 112 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 110 113 ! Update stochastic parameters and random T/S fluctuations 111 114 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 112 CALL sto_par( kstp ) ! Stochastic parameters115 CALL sto_par( kstp ) ! Stochastic parameters 113 116 114 117 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 120 123 CALL bn2 ( tsb, rab_b, rn2b ) ! before Brunt-Vaisala frequency 121 124 CALL bn2 ( tsn, rab_n, rn2 ) ! now Brunt-Vaisala frequency 125 122 126 ! 123 127 ! VERTICAL PHYSICS … … 127 131 IF( lk_zdftke ) CALL zdf_tke( kstp ) ! TKE closure scheme for Kz 128 132 IF( lk_zdfgls ) CALL zdf_gls( kstp ) ! GLS closure scheme for Kz 129 IF( lk_zdfkpp ) CALL zdf_kpp( kstp ) ! KPP closure scheme for Kz130 133 IF( lk_zdfcst ) THEN ! Constant Kz (reset avt, avm[uv] to the background value) 131 134 avt (:,:,:) = rn_avt0 * wmask (:,:,:) … … 133 136 avmv(:,:,:) = rn_avm0 * wvmask(:,:,:) 134 137 ENDIF 138 135 139 IF( ln_rnf_mouth ) THEN ! increase diffusivity at rivers mouths 136 DO jk = 2, nkrnf ; avt(:,:,jk) = avt(:,:,jk) + 2. e0* rn_avt_rnf * rnfmsk(:,:) * tmask(:,:,jk) ; END DO140 DO jk = 2, nkrnf ; avt(:,:,jk) = avt(:,:,jk) + 2._wp * rn_avt_rnf * rnfmsk(:,:) * tmask(:,:,jk) ; END DO 137 141 ENDIF 138 142 IF( ln_zdfevd ) CALL zdf_evd( kstp ) ! enhanced vertical eddy diffusivity … … 140 144 IF( lk_zdftmx ) CALL zdf_tmx( kstp ) ! tidal vertical mixing 141 145 142 IF( lk_zdfddm .AND. .NOT. lk_zdfkpp ) & 143 & CALL zdf_ddm( kstp ) ! double diffusive mixing 146 IF( lk_zdfddm ) CALL zdf_ddm( kstp ) ! double diffusive mixing 144 147 145 148 CALL zdf_mxl( kstp ) ! mixed layer depth … … 151 154 ! LATERAL PHYSICS 152 155 ! 153 IF( lk_ldfslp ) THEN ! slope of lateral mixing 156 IF( l_ldfslp ) THEN ! slope of lateral mixing 157 !!gm : why this here ???? 154 158 IF(ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations 159 !!gm 155 160 CALL eos( tsb, rhd, gdept_0(:,:,:) ) ! before in situ density 161 156 162 IF( ln_zps .AND. .NOT. ln_isfcav) & 157 163 & CALL zps_hde ( kstp, jpts, tsb, gtsu, gtsv, & ! Partial steps: before horizontal gradient 158 164 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 165 159 166 IF( ln_zps .AND. ln_isfcav) & 160 & CALL zps_hde_isf( kstp, jpts, tsb, gtsu, gtsv, & ! Partial steps for top cell (ISF)167 & CALL zps_hde_isf( kstp, jpts, tsb, gtsu, gtsv, gtui, gtvi, & ! Partial steps for top cell (ISF) 161 168 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & 162 & gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the first ocean level 163 IF( ln_traldf_grif ) THEN ! before slope for Griffies operator 164 CALL ldf_slp_grif( kstp ) 165 ELSE 166 CALL ldf_slp( kstp, rhd, rn2b ) ! before slope for Madec operator 169 & grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the first ocean level 170 171 IF( ln_traldf_triad ) THEN 172 CALL ldf_slp_triad( kstp ) ! before slope for triad operator 173 ELSE 174 CALL ldf_slp ( kstp, rhd, rn2b ) ! before slope for standard operator 167 175 ENDIF 168 176 ENDIF 169 #if defined key_traldf_c2d 170 IF( lk_traldf_eiv ) CALL ldf_eiv( kstp ) ! eddy induced velocity coefficient 171 #endif 172 #if defined key_traldf_c3d && defined key_traldf_smag 173 CALL ldf_tra_smag( kstp ) ! eddy induced velocity coefficient 174 # endif 175 #if defined key_dynldf_c3d && defined key_dynldf_smag 176 CALL ldf_dyn_smag( kstp ) ! eddy induced velocity coefficient 177 # endif 178 179 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 180 ! Ocean dynamics : hdiv, rot, ssh, e3, wn 181 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 182 CALL ssh_nxt ( kstp ) ! after ssh (includes call to div_cur) 177 ! ! eddy diffusivity coeff. and/or eiv coeff. 178 IF( l_ldftra_time .OR. l_ldfeiv_time ) CALL ldf_tra( kstp ) 179 180 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 181 ! Ocean dynamics : hdiv, ssh, e3, wn 182 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 183 CALL ssh_nxt ( kstp ) ! after ssh (includes call to div_hor) 183 184 IF( lk_vvl ) CALL dom_vvl_sf_nxt( kstp ) ! after vertical scale factors 184 185 CALL wzv ( kstp ) ! now cross-level velocity … … 188 189 ! Note that the computation of vertical velocity above, hence "after" sea level 189 190 ! is necessary to compute momentum advection for the rhs of barotropic loop: 191 !!gm : why also here ???? 190 192 IF(ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations 193 !!gm 191 194 CALL eos ( tsn, rhd, rhop, fsdept_n(:,:,:) ) ! now in situ density for hpg computation 192 IF( ln_zps .AND. .NOT. ln_isfcav) & 193 & CALL zps_hde ( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps: before horizontal gradient 194 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 195 IF( ln_zps .AND. ln_isfcav) & 196 & CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps for top cell (ISF) 195 196 IF( ln_zps .AND. .NOT. ln_isfcav) & ! Partial steps: bottom before horizontal gradient 197 & CALL zps_hde ( kstp, jpts, tsn, gtsu, gtsv, & ! of t, s, rd at the last ocean level 198 & rhd, gru , grv ) 199 IF( ln_zps .AND. ln_isfcav) & ! Partial steps: top & bottom before horizontal gradient 200 & CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv, gtui, gtvi, & 197 201 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & 198 & gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the last ocean level199 200 ua(:,:,:) = 0. e0! set dynamics trends to zero201 va(:,:,:) = 0. e0202 IF( l n_asmiau .AND. &202 & grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) 203 204 ua(:,:,:) = 0._wp ! set dynamics trends to zero 205 va(:,:,:) = 0._wp 206 IF( lk_asminc .AND. ln_asmiau .AND. & 203 207 & ln_dyninc ) CALL dyn_asm_inc ( kstp ) ! apply dynamics assimilation increment 204 IF( ln_neptsimp ) CALL dyn_nept_cor ( kstp ) ! subtract Neptune velocities (simplified)205 208 IF( lk_bdy ) CALL bdy_dyn3d_dmp( kstp ) ! bdy damping trends 206 209 CALL dyn_adv ( kstp ) ! advection (vector or flux form) 207 210 CALL dyn_vor ( kstp ) ! vorticity term including Coriolis 208 211 CALL dyn_ldf ( kstp ) ! lateral mixing 209 IF( ln_neptsimp ) CALL dyn_nept_cor ( kstp ) ! add Neptune velocities (simplified)210 212 #if defined key_agrif 211 213 IF(.NOT. Agrif_Root()) CALL Agrif_Sponge_dyn ! momentum sponge … … 217 219 va_sv(:,:,:) = va(:,:,:) 218 220 219 CALL div_ cur( kstp ) ! Horizontal divergence & Relative vorticity(2nd call in time-split case)221 CALL div_hor( kstp ) ! Horizontal divergence (2nd call in time-split case) 220 222 IF( lk_vvl ) CALL dom_vvl_sf_nxt( kstp, kcall=2 ) ! after vertical scale factors (update depth average component) 221 223 CALL wzv ( kstp ) ! now cross-level velocity … … 225 227 ! diagnostics and outputs (ua, va, tsa used as workspace) 226 228 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 227 IF( lk_floats ) 228 IF( lk_diahth ) 229 IF( .NOT. ln_cpl )CALL dia_fwb( kstp ) ! Fresh water budget diagnostics230 IF( lk_diadct ) 231 IF( lk_diaar5 ) 232 IF( lk_diaharm ) 233 234 ! 235 IF( ln_crs ) 229 IF( lk_floats ) CALL flo_stp( kstp ) ! drifting Floats 230 IF( lk_diahth ) CALL dia_hth( kstp ) ! Thermocline depth (20 degres isotherm depth) 231 IF(.NOT.ln_cpl ) CALL dia_fwb( kstp ) ! Fresh water budget diagnostics 232 IF( lk_diadct ) CALL dia_dct( kstp ) ! Transports 233 IF( lk_diaar5 ) CALL dia_ar5( kstp ) ! ar5 diag 234 IF( lk_diaharm ) CALL dia_harm( kstp ) ! Tidal harmonic analysis 235 CALL dia_wri( kstp ) ! ocean model: outputs 236 ! 237 IF( ln_crs ) CALL crs_fld( kstp ) ! ocean model: online field coarsening & output 236 238 237 239 #if defined key_top … … 242 244 #endif 243 245 244 245 246 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 246 247 ! Active tracers (ua, va used as workspace) 247 248 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 248 tsa(:,:,:,:) = 0. e0! set tracer trends to zero249 250 IF( l n_asmiau .AND. &249 tsa(:,:,:,:) = 0._wp ! set tracer trends to zero 250 251 IF( lk_asminc .AND. ln_asmiau .AND. & 251 252 & ln_trainc ) CALL tra_asm_inc( kstp ) ! apply tracer assimilation increment 252 253 CALL tra_sbc ( kstp ) ! surface boundary condition … … 257 258 IF( lk_bdy ) CALL bdy_tra_dmp( kstp ) ! bdy damping trends 258 259 CALL tra_adv ( kstp ) ! horizontal & vertical advection 259 IF( lk_zdfkpp ) CALL tra_kpp ( kstp ) ! KPP non-local tracer fluxes260 260 CALL tra_ldf ( kstp ) ! lateral mixing 261 261 262 !!gm : why CALL to dia_ptr has been moved here??? (use trends info?) 262 263 IF( ln_diaptr ) CALL dia_ptr ! Poleward adv/ldf TRansports diagnostics 264 !!gm 263 265 264 266 #if defined key_agrif … … 270 272 IF( ln_zdfnpc ) CALL tra_npc( kstp ) ! update after fields by non-penetrative convection 271 273 CALL tra_nxt( kstp ) ! tracer fields at next time step 274 !!gm : why again a call to sto_pts ??? 272 275 IF( ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations 276 !!gm 273 277 CALL eos ( tsa, rhd, rhop, fsdept_n(:,:,:) ) ! Time-filtered in situ density for hpg computation 274 278 IF( ln_zps .AND. .NOT. ln_isfcav) & … … 276 280 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 277 281 IF( ln_zps .AND. ln_isfcav) & 278 & CALL zps_hde_isf( kstp, jpts, tsa, gtsu, gtsv, & ! Partial steps for top cell (ISF)282 & CALL zps_hde_isf( kstp, jpts, tsa, gtsu, gtsv, gtui, gtvi, & ! Partial steps for top/bottom cells 279 283 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & 280 & gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the last ocean level284 & grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) 281 285 ELSE ! centered hpg (eos then time stepping) 282 286 IF ( .NOT. lk_dynspg_ts ) THEN ! eos already called in time-split case 287 !!gm : why again a call to sto_pts ??? 283 288 IF( ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations 289 !!gm 284 290 CALL eos ( tsn, rhd, rhop, fsdept_n(:,:,:) ) ! now in situ density for hpg computation 285 291 IF( ln_zps .AND. .NOT. ln_isfcav) & 286 & CALL zps_hde ( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps: b efore horizontal gradient287 & rhd, gru , grv ) ! of t, s, rd at the last ocean level292 & CALL zps_hde ( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps: bottom before horizontal gradient 293 & rhd, gru , grv ) ! of t, s, rd at the last ocean level 288 294 IF( ln_zps .AND. ln_isfcav) & 289 & CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps for top cell (ISF)295 & CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv, gtui, gtvi, & ! Partial steps for top/bottom cells 290 296 & rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv , & 291 & g tui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the last ocean level297 & grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi ) ! of t, s, rd at the last ocean level 292 298 ENDIF 293 299 IF( ln_zdfnpc ) CALL tra_npc( kstp ) ! update after fields by non-penetrative convection … … 303 309 ua(:,:,:) = ua_sv(:,:,:) 304 310 va(:,:,:) = va_sv(:,:,:) 305 ! Revert now divergence and rotational to previously computed ones306 !(needed because of the time swap in div_cur, at the beginning of each time step)307 hdivn(:,:,:) = hdivb(:,:,:)308 rotn(:,:,:) = rotb(:,:,:)309 311 310 312 CALL dyn_bfr( kstp ) ! bottom friction 311 313 CALL dyn_zdf( kstp ) ! vertical diffusion 312 314 ELSE 313 ua(:,:,:) = 0. e0! set dynamics trends to zero314 va(:,:,:) = 0. e0315 316 IF( l n_asmiau .AND. &315 ua(:,:,:) = 0._wp ! set dynamics trends to zero 316 va(:,:,:) = 0._wp 317 318 IF( lk_asminc .AND. ln_asmiau .AND. & 317 319 & ln_dyninc ) CALL dyn_asm_inc( kstp ) ! apply dynamics assimilation increment 318 320 IF( ln_bkgwri ) CALL asm_bkg_wri( kstp ) ! output background fields 319 IF( ln_neptsimp ) CALL dyn_nept_cor( kstp ) ! subtract Neptune velocities (simplified)320 321 IF( lk_bdy ) CALL bdy_dyn3d_dmp(kstp ) ! bdy damping trends 321 322 CALL dyn_adv( kstp ) ! advection (vector or flux form) 322 323 CALL dyn_vor( kstp ) ! vorticity term including Coriolis 323 324 CALL dyn_ldf( kstp ) ! lateral mixing 324 IF( ln_neptsimp ) CALL dyn_nept_cor( kstp ) ! add Neptune velocities (simplified)325 325 #if defined key_agrif 326 326 IF(.NOT. Agrif_Root()) CALL Agrif_Sponge_dyn ! momemtum sponge … … 335 335 CALL ssh_swp( kstp ) ! swap of sea surface height 336 336 IF( lk_vvl ) CALL dom_vvl_sf_swp( kstp ) ! swap of vertical scale factors 337 337 ! 338 339 !!gm : This does not only concern the dynamics ==>>> add a new title 340 !!gm2: why ouput restart before AGRIF update? 341 IF( lrst_oce ) CALL rst_write( kstp ) ! write output ocean restart file 342 343 #if defined key_agrif 344 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 345 ! AGRIF 346 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 347 CALL Agrif_Integrate_ChildGrids( stp ) 348 349 IF ( Agrif_NbStepint().EQ.0 ) THEN 350 CALL Agrif_Update_Tra() ! Update active tracers 351 CALL Agrif_Update_Dyn() ! Update momentum 352 ENDIF 353 #endif 338 354 IF( ln_diahsb ) CALL dia_hsb( kstp ) ! - ML - global conservation diagnostics 339 355 IF( lk_diaobs ) CALL dia_obs( kstp ) ! obs-minus-model (assimilation) diagnostics (call after dynamics update) 340 356 341 357 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 342 ! Control and restarts358 ! Control 343 359 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 344 360 CALL stp_ctl( kstp, indic ) … … 350 366 CALL iom_close( numror ) ! close input ocean restart file 351 367 IF(lwm) CALL FLUSH ( numond ) ! flush output namelist oce 352 IF( lwm.AND.numoni /= -1 ) CALL FLUSH ( numoni ) ! flush output namelist ice353 ENDIF354 IF( lrst_oce ) CALL rst_write ( kstp ) ! write output ocean restart file368 IF(lwm.AND.numoni /= -1 ) & 369 & CALL FLUSH ( numoni ) ! flush output namelist ice (if exist) 370 ENDIF 355 371 356 372 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 357 373 ! Coupled mode 358 374 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 375 !!gm why lk_oasis and not lk_cpl ???? 359 376 IF( lk_oasis ) CALL sbc_cpl_snd( kstp ) ! coupled mode : field exchanges 360 377 !
Note: See TracChangeset
for help on using the changeset viewer.
