- Timestamp:
- 2015-07-21T13:25:36+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/NEMO/OPA_SRC/step.F90
r5189 r5621 83 83 IF ( kstp == (nit000 + 1) ) lk_agrif_fstep = .FALSE. 84 84 # if defined key_iomput 85 IF( Agrif_Nbstepint() == 0 ) CALL iom_swap( "nemo")85 IF( Agrif_Nbstepint() == 0 ) CALL iom_swap( cxios_context ) 86 86 # endif 87 87 #endif 88 88 indic = 0 ! reset to no error condition 89 89 IF( kstp == nit000 ) THEN 90 CALL iom_init( "nemo" ) ! iom_put initialization (must be done after nemo_init for AGRIF+XIOS+OASIS) 91 IF( ln_crs ) CALL iom_init( "nemo_crs" ) ! initialize context for coarse grid 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 92 93 ENDIF 93 94 94 95 IF( kstp /= nit000 ) CALL day( kstp ) ! Calendar (day was already called at nit000 in day_init) 95 CALL iom_setkt( kstp - nit000 + 1, "nemo" ) ! say to iom thatwe are at time step kstp96 IF( ln_crs ) CALL iom_setkt( kstp - nit000 + 1, "nemo_crs" ) ! say to iom thatwe are at time step kstp96 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 97 98 98 99 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 100 101 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 101 102 IF( lk_tide ) CALL sbc_tide( kstp ) 102 IF( lk_bdy ) CALL bdy_dta ( kstp, time_offset=+1 ) ! update dynamic & tracer data at open boundaries 103 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 104 107 CALL sbc ( kstp ) ! Sea Boundary Condition (including sea-ice) 105 108 ! clem: moved here for bdy ice purpose 109 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 110 ! Update stochastic parameters and random T/S fluctuations 111 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 112 CALL sto_par( kstp ) ! Stochastic parameters 106 113 107 114 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> … … 145 152 ! 146 153 IF( lk_ldfslp ) THEN ! slope of lateral mixing 154 IF(ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations 147 155 CALL eos( tsb, rhd, gdept_0(:,:,:) ) ! before in situ density 148 156 IF( ln_zps .AND. .NOT. ln_isfcav) & … … 161 169 IF( lk_traldf_eiv ) CALL ldf_eiv( kstp ) ! eddy induced velocity coefficient 162 170 #endif 163 #if defined key_traldf_c3d && key_traldf_smag171 #if defined key_traldf_c3d && defined key_traldf_smag 164 172 CALL ldf_tra_smag( kstp ) ! eddy induced velocity coefficient 165 173 # endif 166 #if defined key_dynldf_c3d && key_dynldf_smag174 #if defined key_dynldf_c3d && defined key_dynldf_smag 167 175 CALL ldf_dyn_smag( kstp ) ! eddy induced velocity coefficient 168 176 # endif … … 179 187 ! Note that the computation of vertical velocity above, hence "after" sea level 180 188 ! is necessary to compute momentum advection for the rhs of barotropic loop: 189 IF(ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations 181 190 CALL eos ( tsn, rhd, rhop, fsdept_n(:,:,:) ) ! now in situ density for hpg computation 182 191 IF( ln_zps .AND. .NOT. ln_isfcav ) & … … 216 225 IF( lk_floats ) CALL flo_stp( kstp ) ! drifting Floats 217 226 IF( lk_diahth ) CALL dia_hth( kstp ) ! Thermocline depth (20 degres isotherm depth) 218 IF( .NOT. l k_cpl ) CALL dia_fwb( kstp ) ! Fresh water budget diagnostics227 IF( .NOT. ln_cpl ) CALL dia_fwb( kstp ) ! Fresh water budget diagnostics 219 228 IF( lk_diadct ) CALL dia_dct( kstp ) ! Transports 220 229 IF( lk_diaar5 ) CALL dia_ar5( kstp ) ! ar5 diag … … 259 268 IF( ln_zdfnpc ) CALL tra_npc( kstp ) ! update after fields by non-penetrative convection 260 269 CALL tra_nxt( kstp ) ! tracer fields at next time step 270 IF( ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations 261 271 CALL eos ( tsa, rhd, rhop, fsdept_n(:,:,:) ) ! Time-filtered in situ density for hpg computation 262 272 IF( ln_zps .AND. .NOT. ln_isfcav) & … … 268 278 ELSE ! centered hpg (eos then time stepping) 269 279 IF ( .NOT. lk_dynspg_ts ) THEN ! eos already called in time-split case 280 IF( ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations 270 281 CALL eos ( tsn, rhd, rhop, fsdept_n(:,:,:) ) ! now in situ density for hpg computation 271 282 IF( ln_zps .AND. .NOT. ln_isfcav) & … … 342 353 ! Coupled mode 343 354 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 344 IF( lk_ cpl) CALL sbc_cpl_snd( kstp ) ! coupled mode : field exchanges355 IF( lk_oasis ) CALL sbc_cpl_snd( kstp ) ! coupled mode : field exchanges 345 356 ! 346 357 #if defined key_iomput 347 358 IF( kstp == nitend .OR. indic < 0 ) THEN 348 CALL iom_context_finalize( "nemo") ! needed for XIOS+AGRIF349 IF( ln_crs ) CALL iom_context_finalize( "nemo_crs" ) !359 CALL iom_context_finalize( cxios_context ) ! needed for XIOS+AGRIF 360 IF( ln_crs ) CALL iom_context_finalize( trim(cxios_context)//"_crs" ) ! 350 361 ENDIF 351 362 #endif
Note: See TracChangeset
for help on using the changeset viewer.