Changeset 4499


Ignore:
Timestamp:
2014-02-18T16:14:31+01:00 (7 years ago)
Author:
timgraham
Message:

Modified advection schemes as in #1243 to prevent trdtrc being called when key_trdtra is set but key_trdtrc is not set.

Location:
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/TRA
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_cen2.F90

    r3718 r4499  
    160160            avmb(4) =  2.5 * avmb(4)      ;      avtb(4) =  2.5 * avtb(4) 
    161161         ENDIF 
    162          ! 
    163          l_trd = .FALSE. 
    164          IF( ( cdtype == 'TRA' .AND. l_trdtra ) .OR. ( cdtype == 'TRC' .AND. l_trdtrc ) )   l_trd = .TRUE. 
    165162      ENDIF 
     163      ! 
     164      l_trd = .FALSE. 
     165      IF( ( cdtype == 'TRA' .AND. l_trdtra ) .OR. ( cdtype == 'TRC' .AND. l_trdtrc ) )   l_trd = .TRUE. 
    166166      ! 
    167167      ! Upstream / centered scheme indicator 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl.F90

    r3718 r4499  
    115115         ENDIF 
    116116         ! 
    117          ! 
    118          l_trd = .FALSE. 
    119          IF( ( cdtype == 'TRA' .AND. l_trdtra ) .OR. ( cdtype == 'TRC' .AND. l_trdtrc ) ) l_trd = .TRUE. 
    120117 
    121118         ! 
     
    137134         ! 
    138135      ENDIF  
     136      ! 
     137      l_trd = .FALSE. 
     138      IF( ( cdtype == 'TRA' .AND. l_trdtra ) .OR. ( cdtype == 'TRC' .AND. l_trdtrc ) ) l_trd = .TRUE. 
     139       
    139140      !                                                     ! =========== 
    140141      DO jn = 1, kjpt                                       ! tracer loop 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl2.F90

    r3625 r4499  
    8888         IF(lwp) WRITE(numout,*) 'tra_adv_muscl2 : MUSCL2 advection scheme on ', cdtype 
    8989         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~~' 
    90          ! 
    91          l_trd = .FALSE. 
    92          IF( ( cdtype == 'TRA' .AND. l_trdtra ) .OR. ( cdtype == 'TRC' .AND. l_trdtrc ) )   l_trd = .TRUE. 
    9390      ENDIF 
     91      ! 
     92      l_trd = .FALSE. 
     93      IF( ( cdtype == 'TRA' .AND. l_trdtra ) .OR. ( cdtype == 'TRC' .AND. l_trdtrc ) )   l_trd = .TRUE. 
    9494 
    9595      !                                                          ! =========== 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_qck.F90

    r3625 r4499  
    102102         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
    103103         IF(lwp) WRITE(numout,*) 
    104          ! 
    105          l_trd = .FALSE. 
    106          IF( ( cdtype == 'TRA' .AND. l_trdtra ) .OR. ( cdtype == 'TRC' .AND. l_trdtrc ) ) l_trd = .TRUE. 
    107104      ENDIF 
     105      ! 
     106      l_trd = .FALSE. 
     107      IF( ( cdtype == 'TRA' .AND. l_trdtra ) .OR. ( cdtype == 'TRC' .AND. l_trdtrc ) ) l_trd = .TRUE. 
    108108 
    109109      ! I. The horizontal fluxes are computed with the QUICKEST + ULTIMATE scheme 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_tvd.F90

    r3625 r4499  
    9393         IF(lwp) WRITE(numout,*) 'tra_adv_tvd : TVD advection scheme on ', cdtype 
    9494         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
    95          ! 
    96          l_trd = .FALSE. 
    97          IF( ( cdtype == 'TRA' .AND. l_trdtra ) .OR. ( cdtype == 'TRC' .AND. l_trdtrc ) ) l_trd = .TRUE. 
    9895      ENDIF 
     96      ! 
     97      l_trd = .FALSE. 
     98      IF( ( cdtype == 'TRA' .AND. l_trdtra ) .OR. ( cdtype == 'TRC' .AND. l_trdtrc ) ) l_trd = .TRUE. 
    9999      ! 
    100100      IF( l_trd )  THEN 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_ubs.F90

    r3787 r4499  
    103103         IF(lwp) WRITE(numout,*) 'tra_adv_ubs :  horizontal UBS advection scheme on ', cdtype 
    104104         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
    105          ! 
    106          l_trd = .FALSE. 
    107          IF( ( cdtype == 'TRA' .AND. l_trdtra ) .OR. ( cdtype == 'TRC' .AND. l_trdtrc ) ) l_trd = .TRUE. 
    108105      ENDIF 
     106      ! 
     107      l_trd = .FALSE. 
     108      IF( ( cdtype == 'TRA' .AND. l_trdtra ) .OR. ( cdtype == 'TRC' .AND. l_trdtrc ) ) l_trd = .TRUE. 
    109109      ! 
    110110      !                                                          ! =========== 
Note: See TracChangeset for help on using the changeset viewer.