Changeset 11475 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/OCE/DOM
- Timestamp:
- 2019-08-28T09:22:35+02:00 (5 years ago)
- Location:
- NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/OCE/DOM
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/OCE/DOM/dom_oce.F90
r11053 r11475 12 12 !! 3.7 ! 2015-11 (G. Madec) introduce surface and scale factor ratio 13 13 !! - ! 2015-11 (G. Madec, A. Coward) time varying zgr by default 14 !! 4.1 ! 2019-08 (A. Coward, D. Storkey) rename prognostic variables in preparation for new time scheme. 14 15 !!---------------------------------------------------------------------- 15 16 -
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps_rewrite_time_filterswap/src/OCE/DOM/domvvl.F90
r11053 r11475 13 13 !! dom_vvl_init : define initial vertical scale factors, depths and column thickness 14 14 !! dom_vvl_sf_nxt : Compute next vertical scale factors 15 !! dom_vvl_sf_ swp: Swap vertical scale factors and update the vertical grid15 !! dom_vvl_sf_update : Swap vertical scale factors and update the vertical grid 16 16 !! dom_vvl_interpol : Interpolate vertical scale factors from one grid point to another 17 17 !! dom_vvl_rst : read/write restart file … … 37 37 PUBLIC dom_vvl_init ! called by domain.F90 38 38 PUBLIC dom_vvl_sf_nxt ! called by step.F90 39 PUBLIC dom_vvl_sf_ swp! called by step.F9039 PUBLIC dom_vvl_sf_update ! called by step.F90 40 40 PUBLIC dom_vvl_interpol ! called by dynnxt.F90 41 41 … … 566 566 567 567 568 SUBROUTINE dom_vvl_sf_ swp( kt, Kbb, Kmm, Kaa )569 !!---------------------------------------------------------------------- 570 !! *** ROUTINE dom_vvl_sf_ swp***568 SUBROUTINE dom_vvl_sf_update( kt, Kbb, Kmm, Kaa ) 569 !!---------------------------------------------------------------------- 570 !! *** ROUTINE dom_vvl_sf_update *** 571 571 !! 572 !! ** Purpose : compute time filter and swap of scale factors572 !! ** Purpose : for z tilde case: compute time filter and swap of scale factors 573 573 !! compute all depths and related variables for next time step 574 574 !! write outputs and restart file 575 575 !! 576 !! ** Method : - swap of e3t with trick for volume/tracer conservation 576 !! ** Method : - swap of e3t with trick for volume/tracer conservation (ONLY FOR Z TILDE CASE) 577 577 !! - reconstruct scale factor at other grid points (interpolate) 578 578 !! - recompute depths and water height fields 579 579 !! 580 !! ** Action : - e3t_(b/n), tilde_e3t_(b/n) and e3(u/v)_nready for next time step580 !! ** Action : - tilde_e3t_(b/n) ready for next time step 581 581 !! - Recompute: 582 582 !! e3(u/v)_b … … 599 599 IF( ln_linssh ) RETURN ! No calculation in linear free surface 600 600 ! 601 IF( ln_timing ) CALL timing_start('dom_vvl_sf_ swp')601 IF( ln_timing ) CALL timing_start('dom_vvl_sf_update') 602 602 ! 603 603 IF( kt == nit000 ) THEN 604 604 IF(lwp) WRITE(numout,*) 605 IF(lwp) WRITE(numout,*) 'dom_vvl_sf_ swp : - time filter and swap of scale factors'606 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~ - interpolate scale factors and compute depths for next time step'605 IF(lwp) WRITE(numout,*) 'dom_vvl_sf_update : - interpolate scale factors and compute depths for next time step' 606 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~~~' 607 607 ENDIF 608 608 ! … … 624 624 ! Horizontal scale factor interpolations 625 625 ! -------------------------------------- 626 ! - ML - e3u(:,:,:,Kbb) and e3v(:,:,:,Kbb) are al lready computed in dynnxt626 ! - ML - e3u(:,:,:,Kbb) and e3v(:,:,:,Kbb) are already computed in dynnxt 627 627 ! - JC - hu(:,:,:,Kbb), hv(:,:,:,:,Kbb), hur_b, hvr_b also 628 628 … … 667 667 IF( lrst_oce ) CALL dom_vvl_rst( kt, Kbb, Kmm, 'WRITE' ) 668 668 ! 669 IF( ln_timing ) CALL timing_stop('dom_vvl_sf_ swp')670 ! 671 END SUBROUTINE dom_vvl_sf_ swp669 IF( ln_timing ) CALL timing_stop('dom_vvl_sf_update') 670 ! 671 END SUBROUTINE dom_vvl_sf_update 672 672 673 673
Note: See TracChangeset
for help on using the changeset viewer.