Changeset 5965 for branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv.F90
- Timestamp:
- 2015-12-01T16:35:30+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO14.5_SST_BIAS_CORRECTION/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv.F90
r4624 r5965 5 5 !!============================================================================== 6 6 !! History : 1.0 ! 2006-11 (G. Madec) Original code 7 !! 3.3 ! 2010-10 (C. Ethe, G. Madec) reorganisation of initialisation phase 7 !! 3.3 ! 2010-10 (C. Ethe, G. Madec) reorganisation of initialisation phase 8 !! 3.6 ! 2015-05 (N. Ducousso, G. Madec) add Hollingsworth scheme as an option 8 9 !!---------------------------------------------------------------------- 9 10 … … 17 18 USE dynkeg ! kinetic energy gradient (dyn_keg routine) 18 19 USE dynzad ! vertical advection (dyn_zad routine) 20 ! 19 21 USE in_out_manager ! I/O manager 20 22 USE lib_mpp ! MPP library … … 25 27 26 28 PUBLIC dyn_adv ! routine called by step module 27 PUBLIC dyn_adv_init ! routine called by opa module29 PUBLIC dyn_adv_init ! routine called by opa module 28 30 31 ! !* namdyn_adv namelist * 29 32 LOGICAL, PUBLIC :: ln_dynadv_vec !: vector form flag 33 INTEGER, PUBLIC :: nn_dynkeg !: scheme of kinetic energy gradient: =0 C2 ; =1 Hollingsworth 30 34 LOGICAL, PUBLIC :: ln_dynadv_cen2 !: flux form - 2nd order centered scheme flag 31 35 LOGICAL, PUBLIC :: ln_dynadv_ubs !: flux form - 3rd order UBS scheme flag 36 LOGICAL, PUBLIC :: ln_dynzad_zts !: vertical advection with sub-timestepping (requires vector form) 32 37 33 38 INTEGER :: nadv ! choice of the formulation and scheme for the advection … … 37 42 # include "vectopt_loop_substitute.h90" 38 43 !!---------------------------------------------------------------------- 39 !! NEMO/OPA 3. 3 , NEMO Consortium (2010)44 !! NEMO/OPA 3.6 , NEMO Consortium (2015) 40 45 !! $Id$ 41 46 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) … … 62 67 SELECT CASE ( nadv ) ! compute advection trend and add it to general trend 63 68 CASE ( 0 ) 64 CALL dyn_keg ( kt ) ! vector form : horizontal gradient of kinetic energy 65 CALL dyn_zad ( kt ) ! vector form : vertical advection 66 CASE ( 1 ) 67 CALL dyn_adv_cen2( kt ) ! 2nd order centered scheme 68 CASE ( 2 ) 69 CALL dyn_adv_ubs ( kt ) ! 3rd order UBS scheme 69 CALL dyn_keg ( kt, nn_dynkeg ) ! vector form : horizontal gradient of kinetic energy 70 CALL dyn_zad ( kt ) ! vector form : vertical advection 71 CASE ( 1 ) 72 CALL dyn_keg ( kt, nn_dynkeg ) ! vector form : horizontal gradient of kinetic energy 73 CALL dyn_zad_zts ( kt ) ! vector form : vertical advection with sub-timestepping 74 CASE ( 2 ) 75 CALL dyn_adv_cen2( kt ) ! 2nd order centered scheme 76 CASE ( 3 ) 77 CALL dyn_adv_ubs ( kt ) ! 3rd order UBS scheme 70 78 ! 71 CASE (-1 ) ! esopa: test all possibility with control print72 CALL dyn_keg ( kt )79 CASE (-1 ) ! esopa: test all possibility with control print 80 CALL dyn_keg ( kt, nn_dynkeg ) 73 81 CALL dyn_zad ( kt ) 74 82 CALL dyn_adv_cen2( kt ) … … 88 96 !! momentum advection formulation & scheme and set nadv 89 97 !!---------------------------------------------------------------------- 90 INTEGER :: ioptio 91 INTEGER :: ios ! Local integer output status for namelist read 92 !! 93 NAMELIST/namdyn_adv/ ln_dynadv_vec, ln_dynadv_cen2 , ln_dynadv_ubs 98 INTEGER :: ioptio, ios ! Local integer 99 ! 100 NAMELIST/namdyn_adv/ ln_dynadv_vec, nn_dynkeg, ln_dynadv_cen2 , ln_dynadv_ubs, ln_dynzad_zts 94 101 !!---------------------------------------------------------------------- 95 102 ! 96 103 REWIND( numnam_ref ) ! Namelist namdyn_adv in reference namelist : Momentum advection scheme 97 104 READ ( numnam_ref, namdyn_adv, IOSTAT = ios, ERR = 901) … … 108 115 WRITE(numout,*) '~~~~~~~~~~~' 109 116 WRITE(numout,*) ' Namelist namdyn_adv : chose a advection formulation & scheme for momentum' 110 WRITE(numout,*) ' Vector/flux form (T/F) ln_dynadv_vec = ', ln_dynadv_vec 111 WRITE(numout,*) ' 2nd order centred advection scheme ln_dynadv_cen2 = ', ln_dynadv_cen2 112 WRITE(numout,*) ' 3rd order UBS advection scheme ln_dynadv_ubs = ', ln_dynadv_ubs 117 WRITE(numout,*) ' Vector/flux form (T/F) ln_dynadv_vec = ', ln_dynadv_vec 118 WRITE(numout,*) ' = 0 standard scheme ; =1 Hollingsworth scheme nn_dynkeg = ', nn_dynkeg 119 WRITE(numout,*) ' 2nd order centred advection scheme ln_dynadv_cen2 = ', ln_dynadv_cen2 120 WRITE(numout,*) ' 3rd order UBS advection scheme ln_dynadv_ubs = ', ln_dynadv_ubs 121 WRITE(numout,*) ' Sub timestepping of vertical advection ln_dynzad_zts = ', ln_dynzad_zts 113 122 ENDIF 114 123 … … 120 129 121 130 IF( ioptio /= 1 ) CALL ctl_stop( 'Choose ONE advection scheme in namelist namdyn_adv' ) 131 IF( ln_dynzad_zts .AND. .NOT. ln_dynadv_vec ) & 132 CALL ctl_stop( 'Sub timestepping of vertical advection requires vector form; set ln_dynadv_vec = .TRUE.' ) 133 IF( nn_dynkeg /= nkeg_C2 .AND. nn_dynkeg /= nkeg_HW ) & 134 CALL ctl_stop( 'KEG scheme wrong value of nn_dynkeg' ) 122 135 123 136 ! ! Set nadv 124 137 IF( ln_dynadv_vec ) nadv = 0 125 IF( ln_dynadv_cen2 ) nadv = 1 126 IF( ln_dynadv_ubs ) nadv = 2 138 IF( ln_dynzad_zts ) nadv = 1 139 IF( ln_dynadv_cen2 ) nadv = 2 140 IF( ln_dynadv_ubs ) nadv = 3 127 141 IF( lk_esopa ) nadv = -1 128 142 129 143 IF(lwp) THEN ! Print the choice 130 144 WRITE(numout,*) 131 IF( nadv == 0 ) WRITE(numout,*) ' vector form : keg + zad + vor is used' 132 IF( nadv == 1 ) WRITE(numout,*) ' flux form : 2nd order scheme is used' 133 IF( nadv == 2 ) WRITE(numout,*) ' flux form : UBS scheme is used' 145 IF( nadv == 0 ) WRITE(numout,*) ' vector form : keg + zad + vor is used' 146 IF( nadv == 1 ) WRITE(numout,*) ' vector form : keg + zad_zts + vor is used' 147 IF( nadv == 0 .OR. nadv == 1 ) THEN 148 IF( nn_dynkeg == nkeg_C2 ) WRITE(numout,*) 'with Centered standard keg scheme' 149 IF( nn_dynkeg == nkeg_HW ) WRITE(numout,*) 'with Hollingsworth keg scheme' 150 ENDIF 151 IF( nadv == 2 ) WRITE(numout,*) ' flux form : 2nd order scheme is used' 152 IF( nadv == 3 ) WRITE(numout,*) ' flux form : UBS scheme is used' 134 153 IF( nadv == -1 ) WRITE(numout,*) ' esopa test: use all advection formulation' 135 154 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.