Changeset 789 for trunk/NEMO/OPA_SRC/TRA/traadv.F90
- Timestamp:
- 2008-01-11T19:04:56+01:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/TRA/traadv.F90
r719 r789 14 14 USE dom_oce ! ocean space and time domain 15 15 USE traadv_cen2 ! 2nd order centered scheme (tra_adv_cen2 routine) 16 USE traadv_cen2_jki ! 2nd order centered scheme (tra_adv_cen2 routine)17 16 USE traadv_tvd ! TVD scheme (tra_adv_tvd routine) 18 17 USE traadv_muscl ! MUSCL scheme (tra_adv_muscl routine) … … 88 87 89 88 SELECT CASE ( nadv ) ! compute advection trend and add it to general trend 90 CASE ( 0 ) ; CALL tra_adv_cen2 ( kt, zun, zvn, zwn ) ! 2nd order centered scheme k-j-i loops 91 CASE ( 1 ) ; CALL tra_adv_cen2_jki( kt, zun, zvn, zwn ) ! 2nd order centered scheme 89 CASE ( 1 ) ; CALL tra_adv_cen2 ( kt, zun, zvn, zwn ) ! 2nd order centered scheme 92 90 CASE ( 2 ) ; CALL tra_adv_tvd ( kt, zun, zvn, zwn ) ! TVD scheme 93 91 CASE ( 3 ) ; CALL tra_adv_muscl ( kt, zun, zvn, zwn ) ! MUSCL scheme … … 99 97 CALL tra_adv_cen2 ( kt, zun, zvn, zwn ) 100 98 CALL prt_ctl( tab3d_1=ta, clinfo1=' adv0 - Ta: ', mask1=tmask, & 101 & tab3d_2=sa, clinfo2= ' Sa: ', mask2=tmask, clinfo3='tra' )102 CALL tra_adv_cen2_jki( kt, zun, zvn, zwn )103 CALL prt_ctl( tab3d_1=ta, clinfo1=' adv1 - Ta: ', mask1=tmask, &104 99 & tab3d_2=sa, clinfo2= ' Sa: ', mask2=tmask, clinfo3='tra' ) 105 100 CALL tra_adv_tvd ( kt, zun, zvn, zwn ) … … 171 166 172 167 ! ! Set nadv 173 IF( ln_traadv_cen2 ) nadv = 0174 #if defined key_mpp_omp175 168 IF( ln_traadv_cen2 ) nadv = 1 176 #endif177 169 IF( ln_traadv_tvd ) nadv = 2 178 170 IF( ln_traadv_muscl ) nadv = 3 … … 184 176 IF(lwp) THEN ! Print the choice 185 177 WRITE(numout,*) 186 IF( nadv == 0 ) WRITE(numout,*) ' 2nd order scheme is used' 187 IF( nadv == 1 ) WRITE(numout,*) ' 2nd order scheme is usedi, k-j-i case' 178 IF( nadv == 1 ) WRITE(numout,*) ' 2nd order scheme is used' 188 179 IF( nadv == 2 ) WRITE(numout,*) ' TVD scheme is used' 189 180 IF( nadv == 3 ) WRITE(numout,*) ' MUSCL scheme is used' 190 181 IF( nadv == 4 ) WRITE(numout,*) ' MUSCL2 scheme is used' 191 182 IF( nadv == 5 ) WRITE(numout,*) ' UBS scheme is used' 192 IF( nadv == 6 ) WRITE(numout,*) ' PPM scheme is used' 193 IF( nadv == 7 ) WRITE(numout,*) ' QUICKEST scheme is used' 183 IF( nadv == 6 ) WRITE(numout,*) ' QUICKEST scheme is used' 194 184 IF( nadv == -1 ) WRITE(numout,*) ' esopa test: use all advection scheme' 195 185 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.