Changeset 3294 for trunk/NEMOGCM/NEMO/OPA_SRC/step.F90
- Timestamp:
- 2012-01-28T17:44:18+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/step.F90
r2715 r3294 23 23 !! 3.3 ! 2010-05 (K. Mogensen, A. Weaver, M. Martin, D. Lea) Assimilation interface 24 24 !! - ! 2010-10 (C. Ethe, G. Madec) reorganisation of initialisation phase + merge TRC-TRA 25 !! 3.4 ! 2011-04 (G. Madec, C. Ethe) Merge of dtatem and dtasal 25 26 !!---------------------------------------------------------------------- 26 27 … … 35 36 USE agrif_opa_sponge ! Momemtum and tracers sponges 36 37 #endif 37 USE asminc ! assimilation increments (tra_asm_inc, dyn_asm_inc routines)38 38 39 39 IMPLICIT NONE … … 94 94 ! Update data, open boundaries, surface boundary condition (including sea-ice) 95 95 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 96 IF( lk_dtatem ) CALL dta_tem( kstp ) ! update 3D temperature data97 IF( lk_dtasal ) CALL dta_sal( kstp ) ! update 3D salinity data98 96 CALL sbc ( kstp ) ! Sea Boundary Condition (including sea-ice) 97 IF( lk_tide ) CALL sbc_tide( kstp ) 99 98 IF( lk_obc ) CALL obc_dta( kstp ) ! update dynamic and tracer data at open boundaries 100 99 IF( lk_obc ) CALL obc_rad( kstp ) ! compute phase velocities at open boundaries 101 IF( lk_bdy ) CALL bdy_dta _frs( kstp ) ! update dynamic and tracer data for FRS conditions (BDY)100 IF( lk_bdy ) CALL bdy_dta( kstp, time_offset=+1 ) ! update dynamic and tracer data at open boundaries 102 101 103 102 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 107 106 108 107 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 109 ! Ocean physics update (ua, va, t a,sa used as workspace)108 ! Ocean physics update (ua, va, tsa used as workspace) 110 109 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 111 110 CALL bn2( tsb, rn2b ) ! before Brunt-Vaisala frequency … … 158 157 159 158 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 160 ! diagnostics and outputs (ua, va, t a,sa used as workspace)159 ! diagnostics and outputs (ua, va, tsa used as workspace) 161 160 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 162 161 IF( lk_floats ) CALL flo_stp( kstp ) ! drifting Floats … … 164 163 IF( lk_diafwb ) CALL dia_fwb( kstp ) ! Fresh water budget diagnostics 165 164 IF( ln_diaptr ) CALL dia_ptr( kstp ) ! Poleward TRansports diagnostics 165 IF( lk_diadct ) CALL dia_dct( kstp ) ! Transports 166 166 IF( lk_diaar5 ) CALL dia_ar5( kstp ) ! ar5 diag 167 IF( lk_diaharm ) CALL dia_harm( kstp ) ! Tidal harmonic analysis 167 168 CALL dia_wri( kstp ) ! ocean model: outputs 168 169 … … 185 186 IF( ln_trabbc ) CALL tra_bbc ( kstp ) ! bottom heat flux 186 187 IF( lk_trabbl ) CALL tra_bbl ( kstp ) ! advective (and/or diffusive) bottom boundary layer scheme 187 IF( l k_tradmp ) CALL tra_dmp ( kstp ) ! internal damping trends188 IF( ln_tradmp ) CALL tra_dmp ( kstp ) ! internal damping trends 188 189 CALL tra_adv ( kstp ) ! horizontal & vertical advection 189 190 IF( lk_zdfkpp ) CALL tra_kpp ( kstp ) ! KPP non-local tracer fluxes 190 191 CALL tra_ldf ( kstp ) ! lateral mixing 191 192 #if defined key_agrif 192 CALL tra_unswap193 193 IF(.NOT. Agrif_Root()) CALL Agrif_Sponge_tra ! tracers sponge 194 CALL tra_swap195 194 #endif 196 195 CALL tra_zdf ( kstp ) ! vertical mixing and after tracer fields … … 210 209 CALL tra_nxt( kstp ) ! tracer fields at next time step 211 210 ENDIF 212 CALL tra_unswap ! udate T & S 3D arrays (to be suppressed) 213 214 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 215 ! Dynamics (ta, sa used as workspace) 211 212 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 213 ! Dynamics (tsa used as workspace) 216 214 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 217 215 ua(:,:,:) = 0.e0 ! set dynamics trends to zero … … 220 218 IF( ln_asmiau .AND. & 221 219 & ln_dyninc ) CALL dyn_asm_inc( kstp ) ! apply dynamics assimilation increment 220 IF( ln_neptsimp ) CALL dyn_nept_cor( kstp ) ! subtract Neptune velocities (simplified) 222 221 CALL dyn_adv( kstp ) ! advection (vector or flux form) 223 222 CALL dyn_vor( kstp ) ! vorticity term including Coriolis 224 223 CALL dyn_ldf( kstp ) ! lateral mixing 224 IF( ln_neptsimp ) CALL dyn_nept_cor( kstp ) ! add Neptune velocities (simplified) 225 225 #if defined key_agrif 226 226 IF(.NOT. Agrif_Root()) CALL Agrif_Sponge_dyn ! momemtum sponge … … 250 250 251 251 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 252 ! Trends (ua, va, t a,sa used as workspace)252 ! Trends (ua, va, tsa used as workspace) 253 253 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 254 254 IF( nstop == 0 ) THEN … … 264 264 IF( lk_cpl ) CALL sbc_cpl_snd( kstp ) ! coupled mode : field exchanges 265 265 ! 266 IF( nn_timing == 1 .AND. kstp == nit000 ) CALL timing_reset 266 267 ! 267 268 END SUBROUTINE stp
Note: See TracChangeset
for help on using the changeset viewer.