- Timestamp:
- 2019-03-19T12:32:19+01:00 (5 years ago)
- Location:
- branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA
- Files:
-
- 26 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90
r10759 r10774 1258 1258 WRITE(numout,*) ' model does not use Conservative Temperature' 1259 1259 ENDIF 1260 IF(lflush) CALL flush(numout) 1260 1261 ENDIF 1261 1262 ! … … 1272 1273 IF(lwp) WRITE(numout,*) 1273 1274 IF(lwp) WRITE(numout,*) ' use of TEOS-10 equation of state (cons. temp. and abs. salinity)' 1275 IF(lwp .AND. lflush) CALL flush(numout) 1274 1276 ! 1275 1277 rdeltaS = 32._wp … … 1462 1464 IF(lwp) WRITE(numout,*) 1463 1465 IF(lwp) WRITE(numout,*) ' use of EOS-80 equation of state (pot. temp. and pract. salinity)' 1466 IF(lwp .AND. lflush) CALL flush(numout) 1464 1467 ! 1465 1468 rdeltaS = 20._wp … … 1662 1665 WRITE(numout,*) ' 2nd cabbel. coef. rn_nu = ', rn_nu 1663 1666 WRITE(numout,*) ' Caution: rn_beta0=0 incompatible with ddm parameterization ' 1667 IF(lflush) CALL flush(numout) 1664 1668 ENDIF 1665 1669 ! … … 1681 1685 IF(lwp) WRITE(numout,*) ' rau0 * rcp rau0_rcp = ', rau0_rcp 1682 1686 IF(lwp) WRITE(numout,*) ' 1. / ( rau0 * rcp ) r1_rau0_rcp = ', r1_rau0_rcp 1687 IF(lwp .AND. lflush) CALL flush(numout) 1683 1688 ! 1684 1689 END SUBROUTINE eos_init -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90
r10759 r10774 223 223 WRITE(numout,*) ' upstream scheme within muscl ln_traadv_msc_ups = ', ln_traadv_msc_ups 224 224 WRITE(numout,*) ' TVD advection scheme with zts ln_traadv_tvd_zts = ', ln_traadv_tvd_zts 225 IF(lflush) CALL flush(numout) 225 226 ENDIF 226 227 … … 260 261 IF( nadv == 7 ) WRITE(numout,*) ' TVD ZTS scheme is used' 261 262 IF( nadv == -1 ) WRITE(numout,*) ' esopa test: use all advection scheme' 263 IF(lflush) CALL flush(numout) 262 264 ENDIF 263 265 ! -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_cen2.F90
r9321 r10774 139 139 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~ ' 140 140 IF(lwp) WRITE(numout,*) 141 IF(lwp .AND. lflush) CALL flush(numout) 141 142 ! 142 143 IF( .NOT. ALLOCATED( upsmsk ) ) THEN -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_eiv.F90
r7179 r10774 97 97 IF(lwp) WRITE(numout,*) 'tra_adv_eiv : eddy induced advection on ', cdtype,' :' 98 98 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ add to velocity fields the eiv component' 99 IF(lwp .AND. lflush) CALL flush(numout) 99 100 # if defined key_diaeiv 100 101 IF( cdtype == 'TRA') THEN -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_mle.F90
r10759 r10774 306 306 WRITE(numout,*) ' =1 no MLE in case of convection ; =0 always MLE nn_conv = ', nn_conv 307 307 WRITE(numout,*) ' Density difference used to define ML for FK rn_rho_c_mle = ', rn_rho_c_mle 308 IF(lflush) CALL flush(numout) 308 309 ENDIF 309 310 ! … … 317 318 WRITE(numout,*) ' Mixed Layer Eddy parametrisation NOT used' 318 319 ENDIF 320 321 IF(lflush) CALL flush(numout) 322 319 323 ENDIF 320 324 ! … … 325 329 IF(lwp) WRITE(numout,*) ' ML buoyancy criteria = ', rb_c, ' m/s2 ' 326 330 IF(lwp) WRITE(numout,*) ' associated ML density criteria defined in zdfmxl = ', rho_c, 'kg/m3' 331 IF(lwp .AND. lflush) CALL flush(numout) 327 332 ! 328 333 IF( nn_mle == 0 ) THEN ! MLE array allocation & initialisation -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl.F90
r7771 r10774 99 99 IF(lwp) WRITE(numout,*) '~~~~~~~' 100 100 IF(lwp) WRITE(numout,*) 101 IF(lwp .AND. lflush) CALL flush(numout) 101 102 ! 102 103 ! -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl2.F90
r7179 r10774 84 84 IF(lwp) WRITE(numout,*) 'tra_adv_muscl2 : MUSCL2 advection scheme on ', cdtype 85 85 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~' 86 IF(lwp .AND. lflush) CALL flush(numout) 86 87 ENDIF 87 88 ! -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_qck.F90
r7179 r10774 101 101 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 102 102 IF(lwp) WRITE(numout,*) 103 IF(lwp .AND. lflush) CALL flush(numout) 103 104 ENDIF 104 105 l_trd = .FALSE. -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_tvd.F90
r10020 r10774 102 102 IF(lwp) WRITE(numout,*) 'tra_adv_tvd : TVD advection scheme on ', cdtype 103 103 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 104 IF(lwp .AND. lflush) CALL flush(numout) 104 105 ! 105 106 ENDIF … … 383 384 IF(lwp) WRITE(numout,*) 'tra_adv_tvd_zts : TVD ZTS advection scheme on ', cdtype 384 385 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 386 IF(lwp .AND. lflush) CALL flush(numout) 385 387 ENDIF 386 388 ! -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_ubs.F90
r7179 r10774 101 101 IF(lwp) WRITE(numout,*) 'tra_adv_ubs : horizontal UBS advection scheme on ', cdtype 102 102 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 103 IF(lwp .AND. lflush) CALL flush(numout) 103 104 ENDIF 104 105 ! -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trabbc.F90
r10759 r10774 158 158 WRITE(numout,*) ' Constant geothermal flux value rn_geoflx_cst = ', rn_geoflx_cst 159 159 WRITE(numout,*) 160 IF(lflush) CALL flush(numout) 160 161 ENDIF 161 162 … … 196 197 IF(lwp) WRITE(numout,*) ' *** no geothermal heat flux' 197 198 ENDIF 199 IF(lwp .AND. lflush) CALL flush(numout) 198 200 ! 199 201 END SUBROUTINE tra_bbc_init -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trabbl.F90
r10759 r10774 358 358 IF(lwp) WRITE(numout,*) 'trabbl:bbl : Compute bbl velocities and diffusive coefficients in ', cdtype 359 359 IF(lwp) WRITE(numout,*) '~~~~~~~~~~' 360 IF(lflush) CALL flush(numout) 360 361 ENDIF 361 362 ! !* bottom variables (T, S, alpha, beta, depth, velocity) … … 530 531 WRITE(numout,*) ' diffusive bbl coefficient rn_ahtbbl = ', rn_ahtbbl, ' m2/s' 531 532 WRITE(numout,*) ' advective bbl coefficient rn_gambbl = ', rn_gambbl, ' s' 533 IF(lflush) CALL flush(numout) 532 534 ENDIF 533 535 … … 538 540 IF( nn_bbl_adv == 1 ) WRITE(numout,*) ' * Advective BBL using upper velocity' 539 541 IF( nn_bbl_adv == 2 ) WRITE(numout,*) ' * Advective BBL using velocity = F( delta rho)' 542 IF(lflush) CALL flush(numout) 540 543 ENDIF 541 544 -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/tradmp.F90
r10759 r10774 214 214 WRITE(numout,*) ' Damping file name cn_resto = ', cn_resto 215 215 WRITE(numout,*) 216 IF(lflush) CALL flush(numout) 216 217 ENDIF 217 218 … … 227 228 CASE ( 2 ) ; IF(lwp) WRITE(numout,*) ' no tracer damping in the mixed layer' 228 229 END SELECT 230 IF(lwp .AND. lflush) CALL flush(numout) 229 231 230 232 !TG: Initialisation of dtatsd - Would it be better to have dmpdta routine -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traldf.F90
r6498 r10774 159 159 WRITE(numout,*) ' see ldf_tra_init report for lateral mixing parameters' 160 160 WRITE(numout,*) 161 IF(lflush) CALL flush(numout) 161 162 ENDIF 162 163 … … 227 228 IF( lk_esopa ) THEN 228 229 IF(lwp) WRITE(numout,*) ' esopa control: use all lateral physics options' 230 IF(lwp .AND. lflush) CALL flush(numout) 229 231 nldf = -1 230 232 ENDIF … … 238 240 IF( nldf == 2 ) WRITE(numout,*) ' bilaplacian operator' 239 241 IF( nldf == 3 ) WRITE(numout,*) ' Rotated bilaplacian' 242 IF(lflush) CALL flush(numout) 240 243 ENDIF 241 244 … … 278 281 WRITE(numout,*) 'tra:ldf_ano : lateral diffusion acting on anomalies' 279 282 WRITE(numout,*) '~~~~~~~~~~~' 283 IF(lflush) CALL flush(numout) 280 284 ENDIF 281 285 … … 292 296 zs_ref(:,:,:) = 35.0 * tmask(:,:,:) 293 297 IF(lwp) WRITE(numout,*) ' homogeneous ocean T = ', zt0, ' S = ',zs0 298 IF(lwp .AND. lflush) CALL flush(numout) 294 299 295 300 ! Initialisation of gtui/gtvi in case of no cavity … … 359 364 WRITE(numout,*) 'tra:ldf_ano : lateral diffusion acting on the full fields' 360 365 WRITE(numout,*) '~~~~~~~~~~~' 366 IF(lflush) CALL flush(numout) 361 367 ENDIF 362 368 END SUBROUTINE ldf_ano -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_bilap.F90
r7179 r10774 102 102 IF(lwp) WRITE(numout,*) 'tra_ldf_bilap : iso-level biharmonic operator on ', cdtype 103 103 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~' 104 IF(lwp .AND. lflush) CALL flush(numout) 104 105 ENDIF 105 106 ! ! =========== -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_bilapg.F90
r8280 r10774 88 88 IF(lwp) WRITE(numout,*) 'tra_ldf_bilapg : horizontal biharmonic operator in s-coordinate on ', cdtype 89 89 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~' 90 IF(lwp .AND. lflush) CALL flush(numout) 90 91 ENDIF 91 92 -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_iso.F90
r7771 r10774 137 137 IF(lwp) WRITE(numout,*) 'tra_ldf_iso : rotated laplacian diffusion operator on ', cdtype 138 138 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 139 IF(lwp .AND. lflush) CALL flush(numout) 139 140 ENDIF 140 141 ! -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_iso_grif.F90
r7179 r10774 128 128 IF(lwp) WRITE(numout,*) 'tra_ldf_iso_grif : rotated laplacian diffusion operator on ', cdtype 129 129 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 130 IF(lwp .AND. lflush) CALL flush(numout) 130 131 ALLOCATE( ah_wslp2(jpi,jpj,jpk) , zdkt3d(jpi,jpj,0:1), STAT=ierr ) 131 132 IF( lk_mpp ) CALL mpp_sum ( ierr ) -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_lap.F90
r7179 r10774 85 85 IF(lwp) WRITE(numout,*) 'tra_ldf_lap : iso-level laplacian diffusion on ', cdtype 86 86 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ ' 87 IF(lwp .AND. lflush) CALL flush(numout) 87 88 ENDIF 88 89 -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/tranpc.F90
r6486 r10774 164 164 END DO 165 165 WRITE(numout,*) 166 IF(lflush) CALL flush(numout) 166 167 ENDIF 167 168 … … 169 170 IF( jiter == 1 ) inpcc = inpcc + 1 170 171 171 IF( lp_monitor_point ) WRITE(numout, *) 'Negative N2 at ikp =',ikp,' for layer #', ilayer 172 IF( lp_monitor_point ) THEN 173 WRITE(numout, *) 'Negative N2 at ikp =',ikp,' for layer #', ilayer 174 IF(lflush) CALL flush(numout) 175 ENDIF 172 176 173 177 !! ikup is the uppermost point where mixing will start: … … 223 227 WRITE(numout,*) ' => Mean Alfa in that portion =', zalfa 224 228 WRITE(numout,*) ' => Mean Beta in that portion =', zbeta 229 IF(lflush) CALL flush(numout) 225 230 ENDIF 226 231 … … 285 290 END DO 286 291 WRITE(numout,*) 292 IF(lflush) CALL flush(numout) 287 293 ENDIF 288 294 ! … … 323 329 WRITE(numout,*) 'Exiting tra_npc , kt = ',kt,', => numb. of statically instable water-columns: ', inpcc 324 330 WRITE(numout,*) 331 IF(lflush) CALL flush(numout) 325 332 ENDIF 326 333 ! -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90
r9163 r10774 105 105 IF(lwp) WRITE(numout,*) 'tra_nxt : achieve the time stepping by Asselin filter and array swap' 106 106 IF(lwp) WRITE(numout,*) '~~~~~~~' 107 IF(lwp .AND. lflush) CALL flush(numout) 107 108 ! 108 109 rbcp = 0.25_wp * (1._wp + atfp) * (1._wp + atfp) * ( 1._wp - atfp) ! Brown & Campana parameter for semi-implicit hpg … … 243 244 IF(lwp) WRITE(numout,*) 'tra_nxt_fix : time stepping', cdtype 244 245 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 246 IF(lwp .AND. lflush) CALL flush(numout) 245 247 ENDIF 246 248 ! … … 317 319 IF(lwp) WRITE(numout,*) 'tra_nxt_vvl : time stepping', cdtype 318 320 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 321 IF(lwp .AND. lflush) CALL flush(numout) 319 322 ENDIF 320 323 ! -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traqsr.F90
r10770 r10774 119 119 IF(lwp) WRITE(numout,*) 'tra_qsr : penetration of the surface solar radiation' 120 120 IF(lwp) WRITE(numout,*) '~~~~~~~' 121 IF(lwp .AND. lflush) CALL flush(numout) 121 122 IF( .NOT.ln_traqsr ) RETURN 122 123 ENDIF … … 135 136 IF( ln_rstart .AND. & ! Restart: read in restart file 136 137 & iom_varid( numror, 'qsr_hc_b', ldstop = .FALSE. ) > 0 ) THEN 137 IF(lwp .AND. nprint >0) WRITE(numout,*) ' nit000-1 qsr tracer content forcing field red in the restart file' 138 IF(lwp .AND. nprint >0) THEN 139 WRITE(numout,*) ' nit000-1 qsr tracer content forcing field red in the restart file' 140 IF(lflush) CALL flush(numout) 141 ENDIF 138 142 zfact = 0.5e0 139 143 CALL iom_get( numror, jpdom_autoglo, 'qsr_hc_b', qsr_hc_b ) ! before heat content trend due to Qsr flux … … 369 373 & 'at it= ', kt,' date= ', ndastp 370 374 WRITE(numout,*) '~~~~' 375 IF(lflush) CALL flush(numout) 371 376 ENDIF 372 377 IF(nn_timing == 2) CALL timing_start('iom_rstput') … … 456 461 WRITE(numout,*) ' RGB & 2 bands: shortess depth of extinction rn_si0 = ', rn_si0 457 462 WRITE(numout,*) ' 2 bands: longest depth of extinction rn_si1 = ', rn_si1 463 IF(lflush) CALL flush(numout) 458 464 ENDIF 459 465 … … 487 493 IF( nqsr == 4 ) WRITE(numout,*) ' 2 bands light penetration' 488 494 IF( nqsr == 5 ) WRITE(numout,*) ' bio-model light penetration' 495 IF(lflush) CALL flush(numout) 489 496 ENDIF 490 497 ! … … 569 576 ENDIF 570 577 ! 578 IF(lwp .AND. lflush) CALL flush(numout) 579 ! 571 580 ENDIF 572 581 ! ! ---------------------------------- ! … … 578 587 IF(lwp) THEN 579 588 WRITE(numout,*) 580 IF(lwp) WRITE(numout,*) ' level of light extinction = ', nksr, ' ref depth = ', gdepw_1d(nksr+1), ' m' 589 WRITE(numout,*) ' level of light extinction = ', nksr, ' ref depth = ', gdepw_1d(nksr+1), ' m' 590 IF(lflush) CALL flush(numout) 581 591 ENDIF 582 592 ! 583 593 IF( lk_vvl ) THEN ! variable volume 584 594 IF(lwp) WRITE(numout,*) ' key_vvl: light distribution will be computed at each time step' 595 IF(lwp .AND. lflush) CALL flush(numout) 585 596 ELSE ! constant volume: computes one for all 586 597 zz0 = rn_abs * r1_rau0_rcp … … 606 617 WRITE(numout,*) 'tra_qsr_init : NO solar flux penetration' 607 618 WRITE(numout,*) '~~~~~~~~~~~~' 619 IF(lflush) CALL flush(numout) 608 620 ENDIF 609 621 ENDIF -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90
r10770 r10774 133 133 IF(lwp) WRITE(numout,*) 'tra_sbc : TRAcer Surface Boundary Condition' 134 134 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 135 IF(lwp .AND. lflush) CALL flush(numout) 135 136 ENDIF 136 137 … … 157 158 & iom_varid( numror, 'sbc_hc_b', ldstop = .FALSE. ) > 0 ) THEN 158 159 IF(lwp) WRITE(numout,*) ' nit000-1 surface tracer content forcing fields red in the restart file' 160 IF(lwp .AND. lflush) CALL flush(numout) 159 161 zfact = 0.5_wp 160 162 CALL iom_get( numror, jpdom_autoglo, 'sbc_hc_b', sbc_tsc_b(:,:,jp_tem) ) ! before heat content sbc trend … … 214 216 & 'at it= ', kt,' date= ', ndastp 215 217 WRITE(numout,*) '~~~~' 218 IF(lflush) CALL flush(numout) 216 219 ENDIF 217 220 IF(nn_timing == 2) CALL timing_start('iom_rstput') … … 262 265 & 'at it= ', kt,' date= ', ndastp 263 266 WRITE(numout,*) '~~~~' 267 IF(lflush) CALL flush(numout) 264 268 ENDIF 265 269 IF(nn_timing == 2) CALL timing_start('iom_rstput') -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf.F90
r8104 r10774 163 163 IF( nzdf == 0 ) WRITE(numout,*) ' Explicit time-splitting scheme' 164 164 IF( nzdf == 1 ) WRITE(numout,*) ' Implicit (euler backward) scheme' 165 IF(lflush) CALL flush(numout) 165 166 ENDIF 166 167 ! -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_exp.F90
r6486 r10774 100 100 IF(lwp) WRITE(numout,*) 'tra_zdf_exp : explicit vertical mixing on ', cdtype 101 101 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 102 IF(lwp .AND. lflush) CALL flush(numout) 102 103 ENDIF 103 104 -
branches/UKMO/dev_r5518_GO6_package_text_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90
r6486 r10774 99 99 IF(lwp)WRITE(numout,*) 'tra_zdf_imp : implicit vertical mixing on ', cdtype 100 100 IF(lwp)WRITE(numout,*) '~~~~~~~~~~~ ' 101 IF(lwp .AND. lflush) CALL flush(numout) 101 102 ! 102 103 IF( lk_vvl ) THEN ; r_vvl = 1._wp ! Variable volume indicator
Note: See TracChangeset
for help on using the changeset viewer.