New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 2240 for branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/TRA – NEMO

Ignore:
Timestamp:
2010-10-13T10:48:48+02:00 (14 years ago)
Author:
cetlod
Message:

Suppression of key_zco everywhere in the code

Location:
branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/TRA
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/TRA/traadv_muscl.F90

    r2236 r2240  
    180180         !                                 ! "Poleward" heat and salt transports (contribution of upstream fluxes) 
    181181         IF( cdtype == 'TRA' .AND. ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN   
    182             IF( lk_zco ) THEN 
    183                DO jk = 1, jpkm1 
    184                   DO jj = 2, jpjm1 
    185                      DO ji = fs_2, fs_jpim1   ! vector opt. 
    186                        zwy(ji,jj,jk) = zwy(ji,jj,jk) * fse3v(ji,jj,jk)                  
    187                      END DO 
    188                   END DO 
    189                END DO 
    190             ENDIF 
    191182            IF( jn == jp_tem )  pht_adv(:) = ptr_vj( zwy(:,:,:) ) 
    192183            IF( jn == jp_sal )  pst_adv(:) = ptr_vj( zwy(:,:,:) ) 

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/TRA/traadv_muscl2.F90

    r2236 r2240  
    202202         !                                 ! "Poleward" heat and salt transports (contribution of upstream fluxes) 
    203203         IF( cdtype == 'TRA' .AND. ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN 
    204             IF( lk_zco ) THEN 
    205                DO jk = 1, jpkm1 
    206                   DO jj = 2, jpjm1 
    207                      DO ji = fs_2, fs_jpim1   ! vector opt. 
    208                        zwy(ji,jj,jk) = zwy(ji,jj,jk) * fse3v(ji,jj,jk) 
    209                      END DO 
    210                   END DO 
    211                END DO 
    212             ENDIF 
    213204            IF( jn == jp_tem )  pht_adv(:) = ptr_vj( zwy(:,:,:) ) 
    214205            IF( jn == jp_sal )  pst_adv(:) = ptr_vj( zwy(:,:,:) ) 

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/TRA/traadv_ubs.F90

    r2236 r2240  
    184184         !                                 ! "Poleward" heat and salt transports (contribution of upstream fluxes) 
    185185         IF( cdtype == 'TRA' .AND. ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN   
    186             IF( lk_zco ) THEN 
    187                DO jk = 1, jpkm1 
    188                   DO jj = 2, jpjm1 
    189                      DO ji = fs_2, fs_jpim1   ! vector opt. 
    190                        zwy(ji,jj,jk) = zwy(ji,jj,jk) * fse3v(ji,jj,jk)                  
    191                      END DO 
    192                   END DO 
    193                END DO 
    194             ENDIF 
    195186            IF( jn == jp_tem )  pht_adv(:) = ptr_vj( zwy(:,:,:) ) 
    196187            IF( jn == jp_sal )  pst_adv(:) = ptr_vj( zwy(:,:,:) ) 

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/TRA/traldf_lap.F90

    r2104 r2240  
    5757      !!          difft = 1/(e1t*e2t*e3t) {  di-1[ aht e2u*e3u/e1u di(tb) ] 
    5858      !!                                   + dj-1[ aht e1v*e3v/e2v dj(tb) ] } 
    59       !!     Note: key_zco defined, the e3t=e3u=e3v, the trend becomes:   
    60       !!          difft = 1/(e1t*e2t) {  di-1[ aht e2u/e1u di(tb) ] 
    61       !!                               + dj-1[ aht e1v/e2v dj(tb) ] } 
    6259      !!      Add this trend to the general tracer trend pta : 
    6360      !!          pta = pta + difft 

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/TRA/tranxt.F90

    r2236 r2240  
    143143      ! Leap-Frog + Asselin filter time stepping 
    144144      IF( lk_vvl )   THEN   ;   CALL tra_nxt_vvl( kt, 'TRA', tsb, tsn, tsa, jpts )  ! variable volume level (vvl)      
    145       ELSE                  ;   CALL tra_nxt_fix( kt,       tsb, tsn, tsa, jpts )  ! fixed    volume level  
     145      ELSE                  ;   CALL tra_nxt_fix( kt, 'TRA', tsb, tsn, tsa, jpts )  ! fixed    volume level  
    146146      ENDIF 
    147147 
    … …  
    172172 
    173173 
    174    SUBROUTINE tra_nxt_fix( kt, ptb, ptn, pta, kjpt ) 
     174   SUBROUTINE tra_nxt_fix( kt, cdtype, ptb, ptn, pta, kjpt ) 
    175175      !!---------------------------------------------------------------------- 
    176176      !!                   ***  ROUTINE tra_nxt_fix  *** 
    … …  
    195195      !!              - (ta,sa) time averaged (t,s)   (ln_dynhpg_imp = T) 
    196196      !!---------------------------------------------------------------------- 
    197       INTEGER , INTENT(in   )                               ::  kt       ! ocean time-step index 
    198       INTEGER , INTENT(in   )                               ::  kjpt     ! number of tracers 
    199       REAL(wp), INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt)  ::  ptb      ! before tracer fields 
    200       REAL(wp), INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt)  ::  ptn      ! now tracer fields 
    201       REAL(wp), INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt)  ::  pta      ! tracer trend 
    202       !! 
     197      INTEGER         , INTENT(in   )                               ::  kt       ! ocean time-step index 
     198      CHARACTER(len=3), INTENT(in   )                               ::  cdtype   ! =TRA or TRC (tracer indicator) 
     199      INTEGER         , INTENT(in   )                               ::  kjpt     ! number of tracers 
     200      REAL(wp)        , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt)  ::  ptb      ! before tracer fields 
     201      REAL(wp)        , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt)  ::  ptn      ! now tracer fields 
     202      REAL(wp)        , INTENT(inout), DIMENSION(jpi,jpj,jpk,kjpt)  ::  pta      ! tracer trend 
    203203      INTEGER  :: ji, jj, jk, jn   ! dummy loop indices 
    204204      REAL(wp) :: ztd, ztm         ! temporary scalars 
    … …  
    222222      ELSE                                           ! general case (Leapfrog + Asselin filter 
    223223         ! 
    224          !                                           ! ----------------------- ! 
    225          IF( ln_dynhpg_imp ) THEN                    ! semi-implicite hpg case ! 
    226             !                                        ! ----------------------- ! 
     224         !                                                     ! ----------------------- ! 
     225         IF( ln_dynhpg_imp .AND. cdtype == 'TRA' ) THEN        ! semi-implicite hpg case ! 
     226            !                                                  ! ----------------------- ! 
    227227            DO jn = 1, kjpt 
    228228               DO jk = 1, jpkm1 
    … …  
    317317      ELSE                                           ! general case (Leapfrog + Asselin filter) 
    318318         ! 
    319          !                                           ! ----------------------- ! 
    320          IF( ln_dynhpg_imp ) THEN                    ! semi-implicite hpg case ! 
    321             !                                        ! ----------------------- ! 
     319         !                                                     ! ----------------------- ! 
     320         IF( ln_dynhpg_imp .AND. cdtype == 'TRA' ) THEN        ! semi-implicite hpg case ! 
     321            !                                                  ! ----------------------- ! 
    322322            DO jn = 1, kjpt                           
    323323               DO jk = 1, jpkm1 
Note: See TracChangeset for help on using the changeset viewer.