- Timestamp:
- 2020-03-13T17:33:02+01:00 (4 years ago)
- Location:
- branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN
- Files:
-
- 23 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/divcur.F90
r6487 r12555 104 104 IF(lwp) WRITE(numout,*) 'div_cur : horizontal velocity divergence and relative vorticity' 105 105 IF(lwp) WRITE(numout,*) '~~~~~~~ NOT optimal for auto-tasking case' 106 IF(lwp .AND. lflush) CALL flush(numout) 106 107 ENDIF 107 108 … … 285 286 IF(lwp) WRITE(numout,*) 'div_cur : horizontal velocity divergence and' 286 287 IF(lwp) WRITE(numout,*) '~~~~~~~ relative vorticity' 288 IF(lwp .AND. lflush) CALL flush(numout) 287 289 ENDIF 288 290 -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv.F90
r6486 r12555 108 108 READ ( numnam_cfg, namdyn_adv, IOSTAT = ios, ERR = 902 ) 109 109 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_adv in configuration namelist', lwp ) 110 IF(lwm ) WRITE ( numond, namdyn_adv )110 IF(lwm .AND. nprint > 2) WRITE ( numond, namdyn_adv ) 111 111 112 112 IF(lwp) THEN ! Namelist print … … 152 152 IF( nadv == 3 ) WRITE(numout,*) ' flux form : UBS scheme is used' 153 153 IF( nadv == -1 ) WRITE(numout,*) ' esopa test: use all advection formulation' 154 IF(lflush) CALL flush(numout) 154 155 ENDIF 155 156 ! -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_cen2.F90
r6486 r12555 66 66 WRITE(numout,*) 'dyn_adv_cen2 : 2nd order flux form momentum advection' 67 67 WRITE(numout,*) '~~~~~~~~~~~~' 68 IF(lflush) CALL flush(numout) 68 69 ENDIF 69 70 ! -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv_ubs.F90
r6486 r12555 90 90 IF(lwp) WRITE(numout,*) 'dyn_adv_ubs : UBS flux form momentum advection' 91 91 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 92 IF(lwp .AND. lflush) CALL flush(numout) 92 93 ENDIF 93 94 ! -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90
r8447 r12555 136 136 IF(lwp) THEN 137 137 WRITE(numout,*) " ! restore original density" 138 IF(lflush) CALL flush(numout) 138 139 ENDIF 139 140 rhd(:,:,:) = z_rhd_st(:,:,:) ! restore original density … … 176 177 READ ( numnam_cfg, namdyn_hpg, IOSTAT = ios, ERR = 902 ) 177 178 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_hpg in configuration namelist', lwp ) 178 IF(lwm ) WRITE ( numond, namdyn_hpg )179 IF(lwm .AND. nprint > 2) WRITE ( numond, namdyn_hpg ) 179 180 ! 180 181 IF(lwp) THEN ! Control print … … 190 191 WRITE(numout,*) ' s-coord. (Pressure Jacobian: Cubic polynomial) ln_hpg_prj = ', ln_hpg_prj 191 192 WRITE(numout,*) ' time stepping: centered (F) or semi-implicit (T) ln_dynhpg_imp = ', ln_dynhpg_imp 193 IF(lflush) CALL flush(numout) 192 194 ENDIF 193 195 ! … … 260 262 IF(lwp) WRITE(numout,*) 'dyn:hpg_zco : hydrostatic pressure gradient trend' 261 263 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ z-coordinate case ' 264 IF(lwp .AND. lflush) CALL flush(numout) 262 265 ENDIF 263 266 … … 325 328 IF(lwp) WRITE(numout,*) 'dyn:hpg_zps : hydrostatic pressure gradient trend' 326 329 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ z-coordinate with partial steps - vector optimization' 330 IF(lwp .AND. lflush) CALL flush(numout) 327 331 ENDIF 328 332 … … 423 427 IF(lwp) WRITE(numout,*) 'dyn:hpg_sco : hydrostatic pressure gradient trend' 424 428 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ s-coordinate case, OPA original scheme used' 429 IF(lwp .AND. lflush) CALL flush(numout) 425 430 ENDIF 426 431 … … 514 519 IF(lwp) WRITE(numout,*) 'dyn:hpg_isf : hydrostatic pressure gradient trend for ice shelf' 515 520 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ s-coordinate case, OPA original scheme used' 521 IF(lwp .AND. lflush) CALL flush(numout) 516 522 ENDIF 517 523 … … 768 774 IF(lwp) WRITE(numout,*) 'dyn:hpg_djc : hydrostatic pressure gradient trend' 769 775 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ s-coordinate case, density Jacobian with cubic polynomial scheme' 776 IF(lwp .AND. lflush) CALL flush(numout) 770 777 ENDIF 771 778 … … 1001 1008 IF(lwp) WRITE(numout,*) 'dyn:hpg_prj : hydrostatic pressure gradient trend' 1002 1009 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ s-coordinate case, cubic spline pressure Jacobian' 1010 IF(lwp .AND. lflush) CALL flush(numout) 1003 1011 ENDIF 1004 1012 -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynkeg.F90
r6486 r12555 88 88 IF(lwp) WRITE(numout,*) 'dyn_keg : kinetic energy gradient trend, scheme number=', kscheme 89 89 IF(lwp) WRITE(numout,*) '~~~~~~~' 90 IF(lflush) CALL flush(numout) 90 91 ENDIF 91 92 -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf.F90
r6486 r12555 100 100 IF(lwp) WRITE(numout,*) 'dyn_ldf : no lateral diffusion on momentum setup' 101 101 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 102 IF(lflush) CALL flush(numout) 102 103 ENDIF 103 104 END SELECT … … 139 140 WRITE(numout,*) ' horizontal (geopotential) ln_dynldf_hor = ', ln_dynldf_hor 140 141 WRITE(numout,*) ' iso-neutral ln_dynldf_iso = ', ln_dynldf_iso 142 IF(lflush) CALL flush(numout) 141 143 ENDIF 142 144 … … 229 231 IF( nldf == 4 ) WRITE(numout,*) ' laplacian and bilaplacian operators' 230 232 IF( nldf == 5 ) WRITE(numout,*) ' rotated laplacian and bilaplacian operators' 233 IF(lflush) CALL flush(numout) 231 234 ENDIF 232 235 ! -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_bilap.F90
r6486 r12555 90 90 WRITE(numout,*) 'dyn_ldf_bilap : iso-level bilaplacian operator' 91 91 WRITE(numout,*) '~~~~~~~~~~~~~' 92 IF(lwp .AND. lflush) CALL flush(numout) 92 93 ENDIF 93 94 -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_bilapg.F90
r8280 r12555 98 98 IF(lwp) WRITE(numout,*) 'dyn_ldf_bilapg : horizontal biharmonic operator in s-coordinate' 99 99 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~' 100 IF(lwp .AND. lflush) CALL flush(numout) 100 101 ! ! allocate dyn_ldf_bilapg arrays 101 102 IF( dyn_ldf_bilapg_alloc() /= 0 ) CALL ctl_stop('STOP', 'dyn_ldf_bilapg: failed to allocate arrays') … … 123 124 IF(lwp) WRITE(numout,*) ' max slop: u', SQRT( MAXVAL(uslp*uslp)), ' v ', SQRT(MAXVAL(vslp)), & 124 125 & ' wi', sqrt(MAXVAL(wslpi)) , ' wj', sqrt(MAXVAL(wslpj)) 126 IF(lwp .AND. lflush) CALL flush(numout) 125 127 endif 126 128 !!end -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_iso.F90
r6486 r12555 127 127 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ s-coordinate horizontal diffusive operator' 128 128 ! ! allocate dyn_ldf_bilap arrays 129 IF(lflush) CALL flush(numout) 129 130 IF( dyn_ldf_iso_alloc() /= 0 ) CALL ctl_stop('STOP', 'dyn_ldf_iso: failed to allocate arrays') 130 131 ENDIF … … 151 152 IF(lwp) WRITE(numout,*) ' max slop: u', SQRT( MAXVAL(uslp*uslp)), ' v ', SQRT(MAXVAL(vslp)), & 152 153 & ' wi', sqrt(MAXVAL(wslpi)) , ' wj', sqrt(MAXVAL(wslpj)) 154 IF(lwp .AND. lflush) CALL flush(numout) 153 155 endif 154 156 !!end -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap.F90
r6486 r12555 70 70 IF(lwp) WRITE(numout,*) 'dyn_ldf : iso-level harmonic (laplacian) operator' 71 71 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 72 IF(lflush) CALL flush(numout) 72 73 ENDIF 73 74 ! ! =============== -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynnept.F90
r6486 r12555 126 126 READ ( numnam_cfg, namdyn_nept, IOSTAT = ios, ERR = 902 ) 127 127 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_nept in configuration namelist', lwp ) 128 IF(lwm ) WRITE ( numond, namdyn_nept )128 IF(lwm .AND. nprint > 2) WRITE ( numond, namdyn_nept ) 129 129 130 130 IF(lwp) THEN ! Control print … … 145 145 WRITE(numout,*) 146 146 ENDIF 147 IF(lflush) CALL flush(numout) 147 148 ENDIF 148 149 ! … … 304 305 CALL lbc_lnk( zunep, 'U', -1.0_wp ) 305 306 CALL lbc_lnk( zvnep, 'V', -1.0_wp ) 306 WRITE(numout,*) ' zunep: min, max = ', unemin,unemax 307 WRITE(numout,*) ' zvnep: min, max = ', vnemin,vnemax 308 WRITE(numout,*) 307 IF(lwp .AND. nprint > 0) THEN 308 WRITE(numout,*) ' zunep: min, max = ', unemin,unemax 309 WRITE(numout,*) ' zvnep: min, max = ', vnemin,vnemax 310 WRITE(numout,*) 311 ENDIF 309 312 310 313 !! Compute, once and for all, the horizontal divergence (zhdivnep) … … 328 331 END DO 329 332 END DO 330 WRITE(numout,*) ' zhdivnep: min, max = ', zhdivmin,zhdivmax 331 WRITE(numout,*) ' zmrotnep: min, max = ', zmrotmin,zmrotmax 332 WRITE(numout,*) 333 334 IF(lwp .AND. nprint > 0) THEN 335 WRITE(numout,*) ' zhdivnep: min, max = ', zhdivmin,zhdivmax 336 WRITE(numout,*) ' zmrotnep: min, max = ', zmrotmin,zmrotmax 337 WRITE(numout,*) 338 IF(lflush) CALL flush(numout) 339 ENDIF 333 340 334 341 !! Deallocate temporary workspace arrays, which are all local to … … 384 391 IF(lwp) WRITE(numout,*) 'horizontal velocity divergence and' 385 392 IF(lwp) WRITE(numout,*) 'relative vorticity of Neptune flow' 393 IF(lwp .AND. lflush) CALL flush(numout) 386 394 #if defined key_noslip_accurate 387 395 !!---------------------------------------------------------------------- -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90
r6487 r12555 119 119 IF(lwp) WRITE(numout,*) 'dyn_nxt : time stepping' 120 120 IF(lwp) WRITE(numout,*) '~~~~~~~' 121 IF(lflush) CALL flush(numout) 121 122 ENDIF 122 123 -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg.F90
r6486 r12555 231 231 WRITE(numout,*) ' Free surface with time splitting lk_dynspg_ts = ', lk_dynspg_ts 232 232 WRITE(numout,*) ' Filtered free surface cst volume lk_dynspg_flt = ', lk_dynspg_flt 233 IF(lflush) CALL flush(numout) 233 234 ENDIF 234 235 … … 266 267 IF( nspg == 1 ) WRITE(numout,*) ' free surface with time splitting scheme' 267 268 IF( nspg == 2 ) WRITE(numout,*) ' filtered free surface' 269 IF(lflush) CALL flush(numout) 268 270 ENDIF 269 271 -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_exp.F90
r9583 r12555 75 75 ! 76 76 IF( lk_vvl .AND. lwp ) WRITE(numout,*) ' lk_vvl=T : spg is included in dynhpg' 77 IF(lflush) CALL flush(numout) 77 78 ENDIF 78 79 -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_flt.F90
r9321 r12555 120 120 IF(lwp) WRITE(numout,*) 'dyn_spg_flt : surface pressure gradient trend' 121 121 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ (free surface constant volume case)' 122 IF(lwp .AND. lflush) CALL flush(numout) 122 123 123 124 ! set to zero free surface specific arrays -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90
r9321 r12555 195 195 IF(lwp) WRITE(numout,*) '~~~~~~~~~~ free surface with time splitting' 196 196 IF(lwp) WRITE(numout,*) 197 IF(lflush) CALL flush(numout) 197 198 ! 198 199 IF (neuler==0) ll_init=.TRUE. … … 1088 1089 READ ( numnam_cfg, namsplit, IOSTAT = ios, ERR = 902 ) 1089 1090 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsplit in configuration namelist', lwp ) 1090 IF(lwm ) WRITE ( numond, namsplit )1091 IF(lwm .AND. nprint > 2) WRITE ( numond, namsplit ) 1091 1092 ! 1092 1093 ! ! Max courant number for ext. grav. waves … … 1161 1162 IF(lwp) WRITE(numout,*) ' Barotropic time step [s] is :', rdtbt 1162 1163 IF(lwp) WRITE(numout,*) ' Maximum Courant number is :', zcmax 1164 IF(lwp .AND. lflush) CALL flush(numout) 1163 1165 ! 1164 1166 IF ((.NOT.ln_bt_av).AND.(.NOT.ln_bt_fw)) THEN -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynvor.F90
r8280 r12555 234 234 IF(lwp) WRITE(numout,*) 'dyn:vor_ene : vorticity term: energy conserving scheme' 235 235 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 236 IF(lflush) CALL flush(numout) 236 237 ENDIF 237 238 … … 346 347 IF(lwp) WRITE(numout,*) 'dyn:vor_mix : vorticity term: mixed energy/enstrophy conserving scheme' 347 348 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 349 IF(lflush) CALL flush(numout) 348 350 ENDIF 349 351 … … 463 465 IF(lwp) WRITE(numout,*) 'dyn:vor_ens : vorticity term: enstrophy conserving scheme' 464 466 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 467 IF(lflush) CALL flush(numout) 465 468 ENDIF 466 469 … … 587 590 IF(lwp) WRITE(numout,*) 'dyn:vor_een : vorticity term: energy and enstrophy conserving scheme' 588 591 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 592 IF(lflush) CALL flush(numout) 589 593 #if ! defined key_vvl 590 594 IF( .NOT.ALLOCATED(ze3f) ) THEN … … 731 735 READ ( numnam_cfg, namdyn_vor, IOSTAT = ios, ERR = 902 ) 732 736 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_vor in configuration namelist', lwp ) 733 IF(lwm ) WRITE ( numond, namdyn_vor )737 IF(lwm .AND. nprint >2) WRITE ( numond, namdyn_vor ) 734 738 735 739 IF(lwp) THEN ! Namelist print … … 743 747 WRITE(numout,*) ' enstrophy and energy conserving scheme ln_dynvor_een = ', ln_dynvor_een 744 748 WRITE(numout,*) ' enstrophy and energy conserving scheme (old) ln_dynvor_een_old= ', ln_dynvor_een_old 749 IF(lflush) CALL flush(numout) 745 750 ENDIF 746 751 … … 798 803 IF( nvor == 3 ) WRITE(numout,*) ' vorticity scheme : energy and enstrophy conserving scheme' 799 804 IF( nvor == -1 ) WRITE(numout,*) ' esopa test: use all lateral physics options' 805 IF(lflush) CALL flush(numout) 800 806 ENDIF 801 807 ! -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynzad.F90
r6486 r12555 74 74 IF(lwp)WRITE(numout,*) 75 75 IF(lwp)WRITE(numout,*) 'dyn_zad : arakawa advection scheme' 76 IF(lwp .AND. lflush) CALL flush(numout) 76 77 ENDIF 77 78 … … 188 189 IF(lwp)WRITE(numout,*) 189 190 IF(lwp)WRITE(numout,*) 'dyn_zad_zts : arakawa advection scheme with sub-timesteps' 191 IF(lwp .AND. lflush) CALL flush(numout) 190 192 ENDIF 191 193 -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf.F90
r6486 r12555 136 136 IF( nzdf == 0 ) WRITE(numout,*) ' Explicit time-splitting scheme' 137 137 IF( nzdf == 1 ) WRITE(numout,*) ' Implicit (euler backward) scheme' 138 IF(lflush) CALL flush(numout) 138 139 ENDIF 139 140 ! -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf_exp.F90
r6486 r12555 73 73 WRITE(numout,*) 'dyn_zdf_exp : vertical momentum diffusion - explicit operator' 74 74 WRITE(numout,*) '~~~~~~~~~~~ ' 75 IF(lflush) CALL flush(numout) 75 76 ENDIF 76 77 -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf_imp.F90
r6795 r12555 81 81 IF(lwp) WRITE(numout,*) 'dyn_zdf_imp : vertical momentum diffusion implicit operator' 82 82 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ ' 83 IF(lwp .AND. lflush) CALL flush(numout) 83 84 ! 84 85 IF( lk_vvl ) THEN ; r_vvl = 1._wp ! Variable volume indicator -
branches/UKMO/dev_r5518_GO6_starthour_obsoper/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90
r8400 r12555 87 87 IF(lwp) WRITE(numout,*) 'ssh_nxt : after sea surface height' 88 88 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 89 IF(lwp .AND. lflush) CALL flush(numout) 89 90 ! 90 91 ENDIF … … 184 185 IF(lwp) WRITE(numout,*) 'wzv : now vertical velocity ' 185 186 IF(lwp) WRITE(numout,*) '~~~~~ ' 187 IF(lwp .AND. lflush) CALL flush(numout) 186 188 ! 187 189 wn(:,:,jpk) = 0._wp ! bottom boundary condition: w=0 (set once for all) … … 265 267 IF(lwp) WRITE(numout,*) 'ssh_swp : Asselin time filter and swap of sea surface height' 266 268 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 269 IF(lwp .AND. lflush) CALL flush(numout) 267 270 ENDIF 268 271
Note: See TracChangeset
for help on using the changeset viewer.