Changeset 52 for trunk/NEMO/OPA_SRC
- Timestamp:
- 2004-04-22T11:47:06+02:00 (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DYN/dynvor.F90
r32 r52 7 7 8 8 !!---------------------------------------------------------------------- 9 !! dyn_vor_enstrophy : enstrophy conserving scheme(ln_dynvor_ens=T)10 !! dyn_vor_energy : energy conserving scheme(ln_dynvor_ene=T)11 !! dyn_vor_mixed :mixed enstrophy/energy conserving (ln_dynvor_mix=T)12 !! dyn_vor_ctl 9 !! dyn_vor_enstrophy: enstrophy conserving scheme (ln_dynvor_ens=T) 10 !! dyn_vor_energy : energy conserving scheme (ln_dynvor_ene=T) 11 !! dyn_vor_mixed : mixed enstrophy/energy conserving (ln_dynvor_mix=T) 12 !! dyn_vor_ctl : control of the different vorticity option 13 13 !!---------------------------------------------------------------------- 14 14 !! * Modules used 15 USE oce 16 USE dom_oce 17 USE in_out_manager 15 USE oce ! ocean dynamics and tracers 16 USE dom_oce ! ocean space and time domain 17 USE in_out_manager ! I/O manager 18 18 USE trddyn_oce ! ocean momentum trends 19 19 … … 25 25 PUBLIC dyn_vor_energy ! routine called by step.F90 26 26 PUBLIC dyn_vor_mixed ! routine called by step.F90 27 PUBLIC dyn_vor_ctl ! routine called by step.F90 27 28 28 29 !! * Shared module variables … … 42 43 SUBROUTINE dyn_vor_energy( kt ) 43 44 !!---------------------------------------------------------------------- 44 !! *** ROUTINE dyn_vor ***45 !! *** ROUTINE dyn_vor_energy *** 45 46 !! 46 47 !! ** Purpose : Compute the now total vorticity trend and add it to … … 87 88 !!---------------------------------------------------------------------- 88 89 89 IF( kt == nit000 ) CALL dyn_vor_ctl ! Check options 90 IF( kt == nit000 ) THEN 91 IF(lwp) WRITE(numout,*) 92 IF(lwp) WRITE(numout,*) 'dyn_vor_energy : vorticity term: energy conserving scheme' 93 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~' 94 ENDIF 90 95 91 96 ! Local constant initialization … … 140 145 ! ! =============== 141 146 142 IF( l_ctl .AND. lwp) THEN ! print sum trends (used for debugging)147 IF(l_ctl) THEN ! print sum trends (used for debugging) 143 148 zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 144 149 zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) … … 198 203 !!---------------------------------------------------------------------- 199 204 200 IF( kt == nit000 ) CALL dyn_vor_ctl ! Check options 205 IF( kt == nit000 ) THEN 206 IF(lwp) WRITE(numout,*) 207 IF(lwp) WRITE(numout,*) 'dyn_vor_mixed : vorticity term: mixed energy/enstrophy conserving scheme' 208 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~' 209 ENDIF 201 210 202 211 ! Local constant initialization … … 251 260 ! ! =============== 252 261 253 IF( l_ctl .AND. lwp) THEN ! print sum trends (used for debugging)262 IF(l_ctl) THEN ! print sum trends (used for debugging) 254 263 zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 255 264 zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) … … 263 272 SUBROUTINE dyn_vor_enstrophy( kt ) 264 273 !!---------------------------------------------------------------------- 265 !! 274 !! *** ROUTINE dyn_vor_enstrophy *** 266 275 !! 267 276 !! ** Purpose : Compute the now total vorticity trend and add it to … … 310 319 !!---------------------------------------------------------------------- 311 320 312 IF( kt == nit000 ) CALL dyn_vor_ctl ! Check options 321 IF( kt == nit000 ) THEN 322 IF(lwp) WRITE(numout,*) 323 IF(lwp) WRITE(numout,*) 'dyn_vor_enstrophy : vorticity term: enstrophy conserving scheme' 324 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~~~' 325 ENDIF 313 326 314 327 ! Local constant initialization … … 381 394 ! ! =============== 382 395 383 IF( l_ctl .AND. lwp) THEN ! print sum trends (used for debugging)396 IF(l_ctl) THEN ! print sum trends (used for debugging) 384 397 zua = SUM( ua(2:jpim1,2:jpjm1,1:jpkm1) * umask(2:jpim1,2:jpjm1,1:jpkm1) ) 385 398 zva = SUM( va(2:jpim1,2:jpjm1,1:jpkm1) * vmask(2:jpim1,2:jpjm1,1:jpkm1) ) … … 414 427 ! Control of vorticity scheme options 415 428 ! ----------------------------------- 416 IF( ln_dynvor_ens .AND. lwp ) THEN 429 ! Control print 430 IF(lwp) THEN 417 431 WRITE(numout,*) 418 WRITE(numout,*) 'dyn_vor_enstrophy : vorticity term : enstrophy conserving scheme' 419 WRITE(numout,*) '~~~~~~~~~~~~~~~~~~' 432 WRITE(numout,*) 'dyn_vor_ctl : vorticity term : read namelist and control the consistency' 433 WRITE(numout,*) '~~~~~~~~~~~' 434 WRITE(numout,*) ' Namelist nam_dynvor : oice of the vorticity term scheme' 435 WRITE(numout,*) ' enstrophy conserving scheme ln_dynvor_ens = ', ln_dynvor_ens 436 WRITE(numout,*) ' energy conserving scheme ln_dynvor_ene = ', ln_dynvor_ene 437 WRITE(numout,*) ' mixed enstrophy/energy conserving scheme ln_dynvor_mix = ', ln_dynvor_mix 438 ENDIF 439 440 IF( ln_dynvor_ens ) THEN 441 IF(lwp) WRITE(numout,*) 442 IF(lwp) WRITE(numout,*) ' vorticity term : enstrophy conserving scheme' 420 443 ioptio = ioptio + 1 421 444 ENDIF 422 IF( ln_dynvor_ene .AND. lwp ) THEN 423 WRITE(numout,*) 424 WRITE(numout,*) 'dyn_vor_energy : vorticity term : energy conserving scheme' 425 WRITE(numout,*) '~~~~~~~~~~~~~~ ' 445 IF( ln_dynvor_ene ) THEN 446 IF(lwp) WRITE(numout,*) 447 IF(lwp) WRITE(numout,*) ' vorticity term : energy conserving scheme' 426 448 ioptio = ioptio + 1 427 449 ENDIF 428 IF( ln_dynvor_mix .AND. lwp ) THEN 429 WRITE(numout,*) 430 WRITE(numout,*) 'dyn_vor_mixed : vorticity term : mixed enstrophy/energy conserving scheme' 431 WRITE(numout,*) '~~~~~~~~~~~~~ ' 450 IF( ln_dynvor_mix ) THEN 451 IF(lwp) WRITE(numout,*) 452 IF(lwp) WRITE(numout,*) ' vorticity term : mixed enstrophy/energy conserving scheme' 432 453 ioptio = ioptio + 1 433 454 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.