Changeset 9124 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/BDY
- Timestamp:
- 2017-12-19T09:26:25+01:00 (7 years ago)
- Location:
- branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/BDY
- Files:
-
- 10 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90
r9019 r9124 87 87 !!--------------------------------------------------------------------------- 88 88 ! 89 IF( nn_timing == 1) CALL timing_start('bdy_dta')89 IF( ln_timing ) CALL timing_start('bdy_dta') 90 90 ! 91 91 ! Initialise data arrays once for all from initial conditions where required … … 395 395 ENDIF 396 396 ! 397 IF( nn_timing == 1 )CALL timing_stop('bdy_dta')397 IF( ln_timing ) CALL timing_stop('bdy_dta') 398 398 ! 399 399 END SUBROUTINE bdy_dta … … 441 441 !!--------------------------------------------------------------------------- 442 442 ! 443 IF( nn_timing == 1 ) CALL timing_start('bdy_dta_init')444 !445 443 IF(lwp) WRITE(numout,*) 446 444 IF(lwp) WRITE(numout,*) 'bdy_dta_ini : initialization of data at the open boundaries' … … 827 825 END DO ! jbdy 828 826 ! 829 IF( nn_timing == 1 ) CALL timing_stop('bdy_dta_init')830 !831 827 END SUBROUTINE bdy_dta_init 832 828 -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/BDY/bdydyn.F90
r7646 r9124 17 17 !!---------------------------------------------------------------------- 18 18 USE wrk_nemo ! Memory Allocation 19 USE timing ! Timing20 19 USE oce ! ocean dynamics and tracers 21 20 USE dom_oce ! ocean space and time domain … … 53 52 REAL(wp), POINTER, DIMENSION(:,:) :: pua2d, pva2d ! after barotropic velocities 54 53 !!---------------------------------------------------------------------- 55 !56 IF( nn_timing == 1 ) CALL timing_start('bdy_dyn')57 54 ! 58 55 ll_dyn2d = .true. … … 129 126 CALL wrk_dealloc( jpi,jpj, pua2d, pva2d ) 130 127 ! 131 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn')132 !133 128 END SUBROUTINE bdy_dyn 134 129 -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/BDY/bdydyn2d.F90
r9023 r9124 14 14 !! bdy_ssh : Duplicate sea level across open boundaries 15 15 !!---------------------------------------------------------------------- 16 USE timing ! Timing17 16 USE oce ! ocean dynamics and tracers 18 17 USE dom_oce ! ocean space and time domain … … 95 94 !!---------------------------------------------------------------------- 96 95 ! 97 IF( nn_timing == 1 ) CALL timing_start('bdy_dyn2d_frs')98 !99 96 igrd = 2 ! Relaxation of zonal velocity 100 97 DO jb = 1, idx%nblen(igrd) … … 115 112 CALL lbc_bdy_lnk( pva2d, 'V', -1., ib_bdy) ! Boundary points should be updated 116 113 ! 117 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn2d_frs')118 !119 120 114 END SUBROUTINE bdy_dyn2d_frs 121 115 … … 153 147 REAL(wp) :: zflag, z1_2 ! " " 154 148 !!---------------------------------------------------------------------- 155 156 IF( nn_timing == 1 ) CALL timing_start('bdy_dyn2d_fla')157 149 158 150 z1_2 = 0.5_wp … … 218 210 CALL lbc_bdy_lnk( pva2d, 'V', -1., ib_bdy ) ! 219 211 ! 220 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn2d_fla')221 !222 212 END SUBROUTINE bdy_dyn2d_fla 223 213 … … 244 234 INTEGER :: ii, ij, iibm1, ijbm1 ! indices 245 235 !!---------------------------------------------------------------------- 246 247 IF( nn_timing == 1 ) CALL timing_start('bdy_dyn2d_orlanski')248 236 ! 249 237 igrd = 2 ! Orlanski bc on u-velocity; … … 258 246 CALL lbc_bdy_lnk( pva2d, 'V', -1., ib_bdy ) ! 259 247 ! 260 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn2d_orlanski')261 !262 248 END SUBROUTINE bdy_dyn2d_orlanski 249 263 250 264 251 SUBROUTINE bdy_ssh( zssh ) -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/BDY/bdydyn3d.F90
r9090 r9124 80 80 !!---------------------------------------------------------------------- 81 81 ! 82 IF( nn_timing == 1 ) CALL timing_start('bdy_dyn3d_spe')83 !84 82 igrd = 2 ! Relaxation of zonal velocity 85 83 DO jb = 1, idx%nblenrim(igrd) … … 104 102 IF( kt == nit000 ) CLOSE( unit = 102 ) 105 103 ! 106 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn3d_spe')107 !108 104 END SUBROUTINE bdy_dyn3d_spe 105 109 106 110 107 SUBROUTINE bdy_dyn3d_zgrad( idx, dta, kt , ib_bdy ) … … 126 123 !!---------------------------------------------------------------------- 127 124 ! 128 IF( nn_timing == 1 ) CALL timing_start('bdy_dyn3d_zgrad')129 !130 125 igrd = 2 ! Copying tangential velocity into bdy points 131 126 DO jb = 1, idx%nblenrim(igrd) … … 152 147 CALL lbc_bdy_lnk( va, 'V', -1., ib_bdy ) 153 148 ! 154 IF( kt .eq. nit000 ) CLOSE( unit = 102 ) 155 156 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn3d_zgrad') 157 149 IF( kt == nit000 ) CLOSE( unit = 102 ) 150 ! 158 151 END SUBROUTINE bdy_dyn3d_zgrad 152 159 153 160 154 SUBROUTINE bdy_dyn3d_zro( idx, dta, kt, ib_bdy ) … … 174 168 REAL(wp) :: zwgt ! boundary weight 175 169 !!---------------------------------------------------------------------- 176 !177 IF( nn_timing == 1 ) CALL timing_start('bdy_dyn3d_zro')178 170 ! 179 171 igrd = 2 ! Everything is at T-points here … … 199 191 IF( kt == nit000 ) CLOSE( unit = 102 ) 200 192 ! 201 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn3d_zro')202 !203 193 END SUBROUTINE bdy_dyn3d_zro 204 194 … … 225 215 !!---------------------------------------------------------------------- 226 216 ! 227 IF( nn_timing == 1 ) CALL timing_start('bdy_dyn3d_frs')228 !229 217 igrd = 2 ! Relaxation of zonal velocity 230 218 DO jb = 1, idx%nblen(igrd) … … 250 238 ! 251 239 IF( kt == nit000 ) CLOSE( unit = 102 ) 252 !253 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn3d_frs')254 240 ! 255 241 END SUBROUTINE bdy_dyn3d_frs … … 273 259 INTEGER :: jb, igrd ! dummy loop indices 274 260 !!---------------------------------------------------------------------- 275 276 IF( nn_timing == 1 ) CALL timing_start('bdy_dyn3d_orlanski')277 261 ! 278 262 !! Note that at this stage the ub and ua arrays contain the baroclinic velocities. … … 289 273 CALL lbc_bdy_lnk( va, 'V', -1., ib_bdy ) 290 274 ! 291 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn3d_orlanski')292 !293 275 END SUBROUTINE bdy_dyn3d_orlanski 294 276 … … 309 291 !!---------------------------------------------------------------------- 310 292 ! 311 IF( nn_timing == 1) CALL timing_start('bdy_dyn3d_dmp')293 IF( ln_timing ) CALL timing_start('bdy_dyn3d_dmp') 312 294 ! 313 295 DO ib_bdy=1, nb_bdy … … 339 321 CALL lbc_lnk_multi( ua, 'U', -1., va, 'V', -1. ) ! Boundary points should be updated 340 322 ! 341 IF( nn_timing == 1) CALL timing_stop('bdy_dyn3d_dmp')323 IF( ln_timing ) CALL timing_stop('bdy_dyn3d_dmp') 342 324 ! 343 325 END SUBROUTINE bdy_dyn3d_dmp 326 344 327 345 328 SUBROUTINE bdy_dyn3d_nmn( idx, ib_bdy ) … … 357 340 INTEGER :: jb, igrd ! dummy loop indices 358 341 !!---------------------------------------------------------------------- 359 360 IF( nn_timing == 1 ) CALL timing_start('bdy_dyn3d_nmn')361 342 ! 362 343 !! Note that at this stage the ub and ua arrays contain the baroclinic velocities. … … 373 354 CALL lbc_bdy_lnk( va, 'V', -1., ib_bdy ) 374 355 ! 375 IF( nn_timing == 1 ) CALL timing_stop('bdy_dyn3d_nmn')376 !377 356 END SUBROUTINE bdy_dyn3d_nmn 378 357 -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/BDY/bdyice.F90
r9019 r9124 58 58 !!---------------------------------------------------------------------- 59 59 ! 60 IF( ln_timing ) CALL timing_start('bdy_ice') 61 ! 60 62 CALL ice_var_glo2eqv 61 63 ! … … 63 65 ! 64 66 SELECT CASE( cn_ice_lim(ib_bdy) ) 65 CASE('none') 66 CYCLE 67 CASE('frs') 68 CALL bdy_ice_frs( idx_bdy(ib_bdy), dta_bdy(ib_bdy), kt, ib_bdy ) 67 CASE('none') ; CYCLE 68 CASE('frs' ) ; CALL bdy_ice_frs( idx_bdy(ib_bdy), dta_bdy(ib_bdy), kt, ib_bdy ) 69 69 CASE DEFAULT 70 70 CALL ctl_stop( 'bdy_ice : unrecognised option for open boundaries for ice fields' ) … … 73 73 END DO 74 74 ! 75 CALL ice_var_zapsmall 76 CALL ice_var_agg(1) 75 CALL ice_var_zapsmall 76 CALL ice_var_agg(1) 77 ! 77 78 IF( ln_icectl ) CALL ice_prt( kt, iiceprt, jiceprt, 1, ' - ice thermo bdy - ' ) 79 IF( ln_timing ) CALL timing_stop('bdy_ice') 78 80 ! 79 81 END SUBROUTINE bdy_ice … … 102 104 REAL(wp) :: ztmelts, zdh 103 105 !!------------------------------------------------------------------------------ 104 !105 IF( ln_timing ) CALL timing_start('bdy_ice_frs')106 106 ! 107 107 jgrd = 1 ! Everything is at T-points here … … 255 255 IF( jpl > 1 ) CALL ice_itd_reb( kt ) 256 256 ! 257 IF( ln_timing ) CALL timing_stop('bdy_ice_frs')258 !259 257 END SUBROUTINE bdy_ice_frs 260 258 … … 277 275 REAL(wp) :: zmsk1, zmsk2, zflag 278 276 !!------------------------------------------------------------------------------ 279 !280 IF( ln_timing ) CALL timing_start('bdy_ice_dyn')281 277 ! 282 278 DO ib_bdy=1, nb_bdy … … 355 351 END DO 356 352 ! 357 IF( ln_timing ) CALL timing_stop('bdy_ice_dyn')358 !359 353 END SUBROUTINE bdy_ice_dyn 360 354 -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90
r9105 r9124 29 29 USE lib_mpp ! for mpp_sum 30 30 USE iom ! I/O 31 USE timing ! Timing32 31 33 32 IMPLICIT NONE … … 74 73 INTEGER :: ios ! Local integer output status for namelist read 75 74 !!---------------------------------------------------------------------- 76 !77 IF( nn_timing == 1 ) CALL timing_start('bdy_init')78 75 79 76 ! ------------------------ … … 112 109 ! 113 110 ENDIF 114 !115 IF( nn_timing == 1 ) CALL timing_stop('bdy_init')116 111 ! 117 112 END SUBROUTINE bdy_init … … 156 151 INTEGER :: ios ! Local integer output status for namelist read 157 152 !!---------------------------------------------------------------------- 158 !159 IF( nn_timing == 1 ) CALL timing_start('bdy_segs')160 153 ! 161 154 cgrid = (/'t','u','v'/) … … 1314 1307 IF( nb_bdy>0 ) DEALLOCATE( nbidta, nbjdta, nbrdta ) 1315 1308 ! 1316 IF( nn_timing == 1 ) CALL timing_stop('bdy_segs')1317 !1318 1309 END SUBROUTINE bdy_segs 1319 1310 -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/BDY/bdylib.F90
r7646 r9124 18 18 USE in_out_manager ! 19 19 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 20 USE timing ! Timing21 20 22 21 IMPLICIT NONE … … 51 50 !!---------------------------------------------------------------------- 52 51 ! 53 IF( nn_timing == 1 ) CALL timing_start('bdy_frs')54 !55 52 igrd = 1 ! Everything is at T-points here 56 53 DO ib = 1, idx%nblen(igrd) … … 63 60 END DO 64 61 ! 65 IF( nn_timing == 1 ) CALL timing_stop('bdy_frs')66 !67 62 END SUBROUTINE bdy_frs 63 68 64 69 65 SUBROUTINE bdy_spe( idx, pta, dta ) … … 82 78 INTEGER :: ii, ij ! 2D addresses 83 79 !!---------------------------------------------------------------------- 84 !85 IF( nn_timing == 1 ) CALL timing_start('bdy_spe')86 80 ! 87 81 igrd = 1 ! Everything is at T-points here … … 94 88 END DO 95 89 ! 96 IF( nn_timing == 1 ) CALL timing_stop('bdy_spe')97 !98 90 END SUBROUTINE bdy_spe 91 99 92 100 93 SUBROUTINE bdy_orl( idx, ptb, pta, dta, ll_npo ) … … 115 108 !!---------------------------------------------------------------------- 116 109 ! 117 IF( nn_timing == 1 ) CALL timing_start('bdy_orl')118 !119 110 igrd = 1 ! Everything is at T-points here 120 111 ! 121 112 CALL bdy_orlanski_3d( idx, igrd, ptb(:,:,:), pta(:,:,:), dta, ll_npo ) 122 113 ! 123 IF( nn_timing == 1 ) CALL timing_stop('bdy_orl')124 !125 114 END SUBROUTINE bdy_orl 115 126 116 127 117 SUBROUTINE bdy_orlanski_2d( idx, igrd, phib, phia, phi_ext, ll_npo ) … … 161 151 REAL(wp), POINTER, DIMENSION(:,:) :: pe_ydif ! scale factors for y-derivatives 162 152 !!---------------------------------------------------------------------- 163 !164 IF( nn_timing == 1 ) CALL timing_start('bdy_orlanski_2d')165 153 ! 166 154 ! ----------------------------------! … … 279 267 END DO 280 268 ! 281 IF( nn_timing == 1 ) CALL timing_stop('bdy_orlanski_2d')282 !283 269 END SUBROUTINE bdy_orlanski_2d 284 270 … … 320 306 REAL(wp), POINTER, DIMENSION(:,:) :: pe_ydif ! scale factors for y-derivatives 321 307 !!---------------------------------------------------------------------- 322 !323 IF( nn_timing == 1 ) CALL timing_start('bdy_orlanski_3d')324 308 ! 325 309 ! ----------------------------------! … … 441 425 END DO 442 426 ! 443 IF( nn_timing == 1 ) CALL timing_stop('bdy_orlanski_3d')444 !445 427 END SUBROUTINE bdy_orlanski_3d 446 428 … … 462 444 INTEGER :: ii, ij, ip, jp ! 2D addresses 463 445 !!---------------------------------------------------------------------- 464 !!465 IF( nn_timing == 1 ) CALL timing_start('bdy_nmn')466 446 ! 467 447 SELECT CASE(igrd) … … 514 494 END DO 515 495 ! 516 IF( nn_timing == 1 ) CALL timing_stop('bdy_nmn')517 !518 496 END SUBROUTINE bdy_nmn 519 497 -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/BDY/bdytides.F90
r7646 r9124 26 26 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 27 27 USE wrk_nemo ! Memory allocation 28 USE timing ! timing29 28 30 29 IMPLICIT NONE … … 83 82 NAMELIST/nambdy_tide/filtide, ln_bdytide_2ddta, ln_bdytide_conj 84 83 !!---------------------------------------------------------------------- 85 !86 IF( nn_timing == 1 ) CALL timing_start('bdytide_init')87 84 ! 88 85 IF (nb_bdy>0) THEN … … 270 267 END DO ! loop on ib_bdy 271 268 ! 272 IF( nn_timing == 1 ) CALL timing_stop('bdytide_init')273 !274 269 END SUBROUTINE bdytide_init 275 270 … … 300 295 REAL(wp), DIMENSION(jpmax_harmo) :: z_sist, z_cost 301 296 !!---------------------------------------------------------------------- 302 !303 IF( nn_timing == 1 ) CALL timing_start('bdytide_update')304 297 ! 305 298 ilen0(1) = SIZE(td%ssh(:,1,1)) … … 363 356 END DO 364 357 ! 365 IF( nn_timing == 1 ) CALL timing_stop('bdytide_update')366 !367 358 END SUBROUTINE bdytide_update 368 359 … … 391 382 REAL(wp) :: z_arg, z_sarg, zramp, zoff, z_cost, z_sist 392 383 !!---------------------------------------------------------------------- 393 !394 IF( nn_timing == 1 ) CALL timing_start('bdy_dta_tides')395 384 ! 396 385 lk_first_btstp=.TRUE. … … 485 474 END IF 486 475 END DO 487 !488 IF( nn_timing == 1 ) CALL timing_stop('bdy_dta_tides')489 476 ! 490 477 END SUBROUTINE bdy_dta_tides -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/BDY/bdytra.F90
r7646 r9124 78 78 END SUBROUTINE bdy_tra 79 79 80 80 81 SUBROUTINE bdy_rnf( idx, pta, jpa ) 81 82 !!---------------------------------------------------------------------- … … 96 97 !!---------------------------------------------------------------------- 97 98 ! 98 IF( nn_timing == 1 ) CALL timing_start('bdy_rnf')99 !100 99 igrd = 1 ! Everything is at T-points here 101 100 DO ib = 1, idx%nblenrim(igrd) … … 110 109 END DO 111 110 ! 112 IF( nn_timing == 1 ) CALL timing_stop('bdy_rnf')113 !114 111 END SUBROUTINE bdy_rnf 112 115 113 116 114 SUBROUTINE bdy_tra_dmp( kt ) … … 130 128 !!---------------------------------------------------------------------- 131 129 ! 132 IF( nn_timing == 1) CALL timing_start('bdy_tra_dmp')130 IF( ln_timing ) CALL timing_start('bdy_tra_dmp') 133 131 ! 134 132 DO ib_bdy = 1, nb_bdy … … 149 147 END DO 150 148 ! 151 IF( nn_timing == 1) CALL timing_stop('bdy_tra_dmp')149 IF( ln_timing ) CALL timing_stop('bdy_tra_dmp') 152 150 ! 153 151 END SUBROUTINE bdy_tra_dmp -
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/BDY/bdyvol.F90
r7646 r9124 19 19 USE in_out_manager ! I/O manager 20 20 USE lib_mpp ! for mppsum 21 USE timing ! Timing22 21 USE lib_fortran ! Fortran routines library 23 22 … … 73 72 TYPE(OBC_INDEX), POINTER :: idx 74 73 !!----------------------------------------------------------------------------- 75 !76 IF( nn_timing == 1 ) CALL timing_start('bdy_vol')77 74 ! 78 75 IF( ln_vol ) THEN … … 165 162 END IF 166 163 ! 167 IF( nn_timing == 1 ) CALL timing_stop('bdy_vol')168 !169 164 END IF ! ln_vol 170 165 !
Note: See TracChangeset
for help on using the changeset viewer.