- Timestamp:
- 2015-06-04T16:12:19+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r5134_UKMO4_CF_compliance/NEMOGCM/NEMO/OPA_SRC/step.F90
r5120 r5350 106 106 107 107 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 108 ! Update stochastic parameters and random T/S fluctuations 109 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 110 CALL sto_par( kstp ) ! Stochastic parameters 111 112 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 108 113 ! Ocean physics update (ua, va, tsa used as workspace) 109 114 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< … … 145 150 ! 146 151 IF( lk_ldfslp ) THEN ! slope of lateral mixing 152 IF(ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations 147 153 CALL eos( tsb, rhd, gdept_0(:,:,:) ) ! before in situ density 148 154 IF( ln_zps .AND. .NOT. ln_isfcav) & … … 180 186 ! Note that the computation of vertical velocity above, hence "after" sea level 181 187 ! is necessary to compute momentum advection for the rhs of barotropic loop: 188 IF(ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations 182 189 CALL eos ( tsn, rhd, rhop, fsdept_n(:,:,:) ) ! now in situ density for hpg computation 183 190 IF( ln_zps .AND. .NOT. ln_isfcav) & … … 216 223 ! diagnostics and outputs (ua, va, tsa used as workspace) 217 224 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 218 IF( lk_floats ) CALL flo_stp( kstp ) ! drifting Floats 219 IF( lk_diahth ) CALL dia_hth( kstp ) ! Thermocline depth (20 degres isotherm depth) 220 IF( .NOT. lk_cpl ) CALL dia_fwb( kstp ) ! Fresh water budget diagnostics 221 IF( ln_diaptr ) CALL dia_ptr( kstp ) ! Poleward TRansports diagnostics 222 IF( lk_diadct ) CALL dia_dct( kstp ) ! Transports 223 IF( lk_diaar5 ) CALL dia_ar5( kstp ) ! ar5 diag 224 IF( lk_diaharm ) CALL dia_harm( kstp ) ! Tidal harmonic analysis 225 CALL dia_wri( kstp ) ! ocean model: outputs 226 ! 227 IF( ln_crs ) CALL crs_fld( kstp ) ! ocean model: online field coarsening & output 228 225 IF( lk_floats ) CALL flo_stp( kstp ) ! drifting Floats 226 IF( lk_diahth ) CALL dia_hth( kstp ) ! Thermocline depth (20 degres isotherm depth) 227 IF( .NOT. lk_cpl ) CALL dia_fwb( kstp ) ! Fresh water budget diagnostics 228 IF( lk_diadct ) CALL dia_dct( kstp ) ! Transports 229 IF( lk_diaar5 ) CALL dia_ar5( kstp ) ! ar5 diag 230 IF( lk_diaharm ) CALL dia_harm( kstp ) ! Tidal harmonic analysis 231 CALL dia_wri( kstp ) ! ocean model: outputs 232 ! 233 IF( ln_crs ) CALL crs_fld( kstp ) ! ocean model: online field coarsening & output 229 234 230 235 #if defined key_top … … 252 257 IF( lk_zdfkpp ) CALL tra_kpp ( kstp ) ! KPP non-local tracer fluxes 253 258 CALL tra_ldf ( kstp ) ! lateral mixing 259 260 IF( ln_diaptr ) CALL dia_ptr ! Poleward adv/ldf TRansports diagnostics 261 254 262 #if defined key_agrif 255 263 IF(.NOT. Agrif_Root()) CALL Agrif_Sponge_tra ! tracers sponge … … 260 268 IF( ln_zdfnpc ) CALL tra_npc( kstp ) ! update after fields by non-penetrative convection 261 269 CALL tra_nxt( kstp ) ! tracer fields at next time step 270 IF( ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations 262 271 CALL eos ( tsa, rhd, rhop, fsdept_n(:,:,:) ) ! Time-filtered in situ density for hpg computation 263 272 IF( ln_zps .AND. .NOT. ln_isfcav) & … … 270 279 ELSE ! centered hpg (eos then time stepping) 271 280 IF ( .NOT. lk_dynspg_ts ) THEN ! eos already called in time-split case 281 IF( ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations 272 282 CALL eos ( tsn, rhd, rhop, fsdept_n(:,:,:) ) ! now in situ density for hpg computation 273 283 IF( ln_zps .AND. .NOT. ln_isfcav) & … … 338 348 CALL iom_close( numror ) ! close input ocean restart file 339 349 IF(lwm) CALL FLUSH ( numond ) ! flush output namelist oce 340 IF( lwm.AND.numoni /= -1 ) CALL FLUSH ( numoni ) ! flush output namelist ice 350 IF( lwm.AND.numoni /= -1 ) CALL FLUSH ( numoni ) ! flush output namelist ice 341 351 ENDIF 342 352 IF( lrst_oce ) CALL rst_write ( kstp ) ! write output ocean restart file
Note: See TracChangeset
for help on using the changeset viewer.