- Timestamp:
- 2010-07-29T17:05:35+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/TRP/trcadv.F90
r2030 r2034 39 39 # include "vectopt_loop_substitute.h90" 40 40 !!---------------------------------------------------------------------- 41 !! NEMO/OPA 3. 2 , LOCEAN-IPSL (2009)42 !! $Id: trcadv.F90 1601 2009-08-11 10:09:19Z ctlod $41 !! NEMO/OPA 3.3 , LOCEAN-IPSL (2010) 42 !! $Id: trcadv.F90 2024 2010-07-29 10:57:35Z cetlod $ 43 43 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 44 44 !!---------------------------------------------------------------------- … … 76 76 ! ! add the eiv transport (if necessary) 77 77 IF( lk_traldf_eiv ) CALL tra_adv_eiv( kt, zun, zvn, zwn, 'TRC' ) 78 79 SELECT CASE ( nadv ) ! compute advection trend and add it to general trend80 CASE ( 1 ) ; CALL tra_adv_cen2 ( kt , 'TRC', zun, zvn, zwn, &81 & trb, trn , tra, jptra ) ! 2nd order centered scheme82 CASE ( 2 ) ; CALL tra_adv_tvd ( kt , 'TRC', zun, zvn, zwn, &83 & trb, trn , tra, jptra ) ! TVD scheme84 CASE ( 3 ) ; CALL tra_adv_muscl ( kt , 'TRC', zun, zvn, zwn, &85 & trb, tra , jptra ) ! MUSCL scheme86 CASE ( 4 ) ; CALL tra_adv_muscl2( kt , 'TRC', zun, zvn, zwn, &87 & trb, trn , tra, jptra ) ! MUSCL2 scheme88 CASE ( 5 ) ; CALL tra_adv_ubs ( kt , 'TRC', zun, zvn, zwn, &89 & trb, trn , tra, jptra ) ! UBS scheme90 CASE ( 6 ) ; CALL tra_adv_qck ( kt , 'TRC', zun, zvn, zwn, &91 & trb, trn , tra, jptra ) ! QUICKEST scheme92 78 ! 93 CASE (-1 ) ! esopa: test all possibility with control print 94 CALL tra_adv_cen2 ( kt , 'TRC', zun, zvn, zwn, & 95 & trb, trn , tra, jptra ) 96 WRITE(charout, FMT="('adv1')") ; CALL prt_ctl_trc_info(charout) 97 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 98 ! 99 CALL tra_adv_tvd ( kt , 'TRC', zun, zvn, zwn, & 100 & trb, trn , tra, jptra ) 101 WRITE(charout, FMT="('adv2')") ; CALL prt_ctl_trc_info(charout) 102 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 103 ! 104 CALL tra_adv_muscl ( kt , 'TRC', zun, zvn, zwn, & 105 & trb, tra, jptra ) 106 WRITE(charout, FMT="('adv3')") ; CALL prt_ctl_trc_info(charout) 107 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 108 ! 109 CALL tra_adv_muscl2( kt , 'TRC', zun, zvn, zwn, & 110 & trb, trn , tra, jptra ) 111 WRITE(charout, FMT="('adv4')") ; CALL prt_ctl_trc_info(charout) 112 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 113 ! 114 CALL tra_adv_ubs ( kt , 'TRC', zun, zvn, zwn, & 115 & trb, trn , tra, jptra ) 116 WRITE(charout, FMT="('adv5')") ; CALL prt_ctl_trc_info(charout) 117 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 118 ! 119 CALL tra_adv_qck ( kt , 'TRC', zun, zvn, zwn, & 120 & trb, trn , tra, jptra ) 121 WRITE(charout, FMT="('adv6')") ; CALL prt_ctl_trc_info(charout) 122 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 123 ! 79 SELECT CASE ( nadv ) !== compute advection trend and add it to general trend ==! 80 CASE ( 1 ) ; CALL tra_adv_cen2 ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra ) ! 2nd order centered 81 CASE ( 2 ) ; CALL tra_adv_tvd ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra ) ! TVD 82 CASE ( 3 ) ; CALL tra_adv_muscl ( kt, 'TRC', zun, zvn, zwn, trb, tra, jptra ) ! MUSCL 83 CASE ( 4 ) ; CALL tra_adv_muscl2( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra ) ! MUSCL2 84 CASE ( 5 ) ; CALL tra_adv_ubs ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra ) ! UBS 85 CASE ( 6 ) ; CALL tra_adv_qck ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra ) ! QUICKEST 86 ! 87 CASE (-1 ) !== esopa: test all possibility with control print ==! 88 CALL tra_adv_cen2 ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra ) 89 WRITE(charout, FMT="('adv1')") ; CALL prt_ctl_trc_info(charout) 90 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 91 CALL tra_adv_tvd ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra ) 92 WRITE(charout, FMT="('adv2')") ; CALL prt_ctl_trc_info(charout) 93 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 94 CALL tra_adv_muscl ( kt, 'TRC', zun, zvn, zwn, trb, tra, jptra ) 95 WRITE(charout, FMT="('adv3')") ; CALL prt_ctl_trc_info(charout) 96 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 97 CALL tra_adv_muscl2( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra ) 98 WRITE(charout, FMT="('adv4')") ; CALL prt_ctl_trc_info(charout) 99 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 100 CALL tra_adv_ubs ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra ) 101 WRITE(charout, FMT="('adv5')") ; CALL prt_ctl_trc_info(charout) 102 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 103 CALL tra_adv_qck ( kt, 'TRC', zun, zvn, zwn, trb, trn, tra, jptra ) 104 WRITE(charout, FMT="('adv6')") ; CALL prt_ctl_trc_info(charout) 105 CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 106 ! 124 107 END SELECT 125 108
Note: See TracChangeset
for help on using the changeset viewer.