- Timestamp:
- 2016-07-19T10:38:35+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5549_BDY_ZEROGRAD/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv.F90
r5322 r6808 39 39 40 40 !! * Substitutions 41 # include "domzgr_substitute.h90"42 41 # include "vectopt_loop_substitute.h90" 43 42 !!---------------------------------------------------------------------- … … 76 75 CASE ( 3 ) 77 76 CALL dyn_adv_ubs ( kt ) ! 3rd order UBS scheme 78 !79 CASE (-1 ) ! esopa: test all possibility with control print80 CALL dyn_keg ( kt, nn_dynkeg )81 CALL dyn_zad ( kt )82 CALL dyn_adv_cen2( kt )83 CALL dyn_adv_ubs ( kt )84 77 END SELECT 85 78 ! … … 104 97 READ ( numnam_ref, namdyn_adv, IOSTAT = ios, ERR = 901) 105 98 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_adv in reference namelist', lwp ) 106 99 ! 107 100 REWIND( numnam_cfg ) ! Namelist namdyn_adv in configuration namelist : Momentum advection scheme 108 101 READ ( numnam_cfg, namdyn_adv, IOSTAT = ios, ERR = 902 ) … … 126 119 IF( ln_dynadv_cen2 ) ioptio = ioptio + 1 127 120 IF( ln_dynadv_ubs ) ioptio = ioptio + 1 128 IF( lk_esopa ) ioptio = 1129 121 130 122 IF( ioptio /= 1 ) CALL ctl_stop( 'Choose ONE advection scheme in namelist namdyn_adv' ) … … 139 131 IF( ln_dynadv_cen2 ) nadv = 2 140 132 IF( ln_dynadv_ubs ) nadv = 3 141 IF( lk_esopa ) nadv = -1142 133 143 134 IF(lwp) THEN ! Print the choice … … 151 142 IF( nadv == 2 ) WRITE(numout,*) ' flux form : 2nd order scheme is used' 152 143 IF( nadv == 3 ) WRITE(numout,*) ' flux form : UBS scheme is used' 153 IF( nadv == -1 ) WRITE(numout,*) ' esopa test: use all advection formulation'154 144 ENDIF 155 145 !
Note: See TracChangeset
for help on using the changeset viewer.