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 – 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
Files:
14 edited

Legend:

Unmodified
Added
Removed

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DOM/dom_oce.F90

    r2239 r2240  
    120120   LOGICAL, PUBLIC ::   ln_zps     =  .FALSE.   !: z-coordinate - partial step 
    121121   LOGICAL, PUBLIC ::   ln_sco     =  .FALSE.   !: s-coordinate or hybrid z-s coordinate 
    122 #if defined key_zco 
    123    LOGICAL, PUBLIC, PARAMETER ::   lk_zco = .TRUE.    !: z-coordinate flag (1D arrays) 
    124 #else 
    125    LOGICAL, PUBLIC, PARAMETER ::   lk_zco = .FALSE.   !: z-coordinate flag (3D arrays) 
    126122 
    127123   !! All coordinates 
    … …  
    133129   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   e3vw            !: analytical vertical scale factors at  VW-- 
    134130   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   e3w   , e3uw    !:                                        W--UW  points (m) 
    135 #endif 
    136131#if defined key_vvl 
    137132   LOGICAL, PUBLIC, PARAMETER ::   lk_vvl = .TRUE.    !: variable grid flag 

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DOM/domvvl.F90

    r2236 r2240  
    6262      ENDIF 
    6363 
    64       IF( lk_zco )   CALL ctl_stop( 'dom_vvl : key_zco is incompatible with variable volume option key_vvl') 
    6564 
    6665      fsdept(:,:,:) = gdept (:,:,:) 

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DOM/domwri.F90

    r2004 r2240  
    168168      CALL iom_rstput( 0, 0, inum4, 'mbathy', zprt, ktype = jp_i2 ) 
    169169             
    170 #if ! defined key_zco 
    171170      IF( ln_sco ) THEN                                         ! s-coordinate 
    172171         CALL iom_rstput( 0, 0, inum4, 'hbatt', hbatt )         !    ! depth 
    … …  
    235234      ENDIF 
    236235       
    237 #endif 
    238236       
    239237      IF( ln_zco ) THEN 

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DOM/domzgr.F90

    r2236 r2240  
    7474      !!              - vertical coordinate (gdep., e3.) depending on the  
    7575      !!                coordinate chosen : 
    76       !!                   ln_zco=T   z-coordinate   (forced if lk_zco) 
     76      !!                   ln_zco=T   z-coordinate    
    7777      !!                   ln_zps=T   z-coordinate with partial steps 
    7878      !!                   ln_zco=T   s-coordinate  
    … …  
    103103      IF( ln_sco ) ioptio = ioptio + 1 
    104104      IF ( ioptio /= 1 )   CALL ctl_stop( ' none or several vertical coordinate options used' ) 
    105       IF( lk_zco ) THEN 
    106           IF(lwp) WRITE(numout,*) '          z-coordinate with reduced incore memory requirement' 
    107           IF( ln_zps .OR. ln_sco )   CALL ctl_stop( ' reduced memory with zps or sco option is impossible' ) 
    108       ENDIF 
    109105 
    110106      ! Build the vertical coordinate system 
    … …  
    679675      !! ** Purpose :   define the z-coordinate system 
    680676      !! 
    681       !! ** Method  :   set 3D coord. arrays to reference 1D array if lk_zco=T 
     677      !! ** Method  :   set 3D coord. arrays to reference 1D array  
    682678      !!---------------------------------------------------------------------- 
    683679      INTEGER  ::   jk 
    684680      !!---------------------------------------------------------------------- 
    685681      ! 
    686       IF( .NOT.lk_zco ) THEN 
    687          DO jk = 1, jpk 
    688             fsdept(:,:,jk) = gdept_0(jk) 
    689             fsdepw(:,:,jk) = gdepw_0(jk) 
    690             fsde3w(:,:,jk) = gdepw_0(jk) 
    691             fse3t (:,:,jk) = e3t_0(jk) 
    692             fse3u (:,:,jk) = e3t_0(jk) 
    693             fse3v (:,:,jk) = e3t_0(jk) 
    694             fse3f (:,:,jk) = e3t_0(jk) 
    695             fse3w (:,:,jk) = e3w_0(jk) 
    696             fse3uw(:,:,jk) = e3w_0(jk) 
    697             fse3vw(:,:,jk) = e3w_0(jk) 
    698          END DO 
    699       ENDIF 
     682      DO jk = 1, jpk 
     683         fsdept(:,:,jk) = gdept_0(jk) 
     684         fsdepw(:,:,jk) = gdepw_0(jk) 
     685         fsde3w(:,:,jk) = gdepw_0(jk) 
     686         fse3t (:,:,jk) = e3t_0(jk) 
     687         fse3u (:,:,jk) = e3t_0(jk) 
     688         fse3v (:,:,jk) = e3t_0(jk) 
     689         fse3f (:,:,jk) = e3t_0(jk) 
     690         fse3w (:,:,jk) = e3w_0(jk) 
     691         fse3uw(:,:,jk) = e3w_0(jk) 
     692         fse3vw(:,:,jk) = e3w_0(jk) 
     693      END DO 
    700694      ! 
    701695   END SUBROUTINE zgr_zco 
    702696 
    703 #if defined key_zco 
    704    !!---------------------------------------------------------------------- 
    705    !!   'key_zco' :                                              "pure" zco (gdep & e3 are 1D arrays) 
    706    !!---------------------------------------------------------------------- 
    707    SUBROUTINE zgr_zps      ! Empty routine 
    708    END SUBROUTINE zgr_zps 
    709    SUBROUTINE zgr_sco      ! Empty routine 
    710    END SUBROUTINE zgr_sco 
    711  
    712 #else 
    713697   !!---------------------------------------------------------------------- 
    714698   !!   Default option :                      zco, zps and/or sco available (gedp & e3 are 3D arrays) 
    … …  
    15581542   END SUBROUTINE zgr_sco 
    15591543 
    1560 #endif 
    15611544 
    15621545   !!====================================================================== 

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DOM/domzgr_substitute.h90

    r2236 r2240  
    88   !!            3.1  !  2009-02  (G. Madec, M. Leclair)  pure z* coordinate 
    99   !!---------------------------------------------------------------------- 
    10 #if defined key_zco 
    11 ! reference for pure z-coordinate (1D - no i,j and time dependency) 
    12 #   define  fsdept_0(i,j,k)  gdept_0(k) 
    13 #   define  fsdepw_0(i,j,k)  gdepw_0(k) 
    14 #   define  fsde3w_0(i,j,k)  gdepw_0(k) 
    15 #   define  fse3t_0(i,j,k)   e3t_0(k) 
    16 #   define  fse3u_0(i,j,k)   e3t_0(k) 
    17 #   define  fse3v_0(i,j,k)   e3t_0(k) 
    18 #   define  fse3f_0(i,j,k)   e3t_0(k) 
    19 #   define  fse3w_0(i,j,k)   e3w_0(k) 
    20 #   define  fse3uw_0(i,j,k)  e3w_0(k) 
    21 #   define  fse3vw_0(i,j,k)  e3w_0(k) 
    22 #else 
    2310! reference for s- or zps-coordinate (3D no time dependency) 
    2411#   define  fsdept_0(i,j,k)  gdept(i,j,k) 
    … …  
    3219#   define  fse3uw_0(i,j,k)  e3uw(i,j,k) 
    3320#   define  fse3vw_0(i,j,k)  e3vw(i,j,k) 
    34 #endif 
    3521#if defined key_vvl 
    3622! s* or z*-coordinate (3D + time dependency) + use of additional now arrays (..._1) 

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DYN/divcur.F90

    r2236 r2240  
    5454      !!         - compute the now divergence given by : 
    5555      !!         hdivn = 1/(e1t*e2t*e3t) ( di[e2u*e3u un] + dj[e1v*e3v vn] ) 
    56       !!      Note: if lk_zco=T, e3u=e3v=e3t, they are simplified in the 
    5756      !!      above expression 
    5857      !!         - apply lateral boundary conditions on hdivn  
    … …  
    112111         DO jj = 2, jpjm1 
    113112            DO ji = fs_2, fs_jpim1   ! vector opt. 
    114 #if defined key_zco 
    115                hdivn(ji,jj,jk) = (  e2u(ji,jj) * un(ji,jj,jk) - e2u(ji-1,jj  ) * un(ji-1,jj  ,jk)      & 
    116                   &               + e1v(ji,jj) * vn(ji,jj,jk) - e1v(ji  ,jj-1) * vn(ji  ,jj-1,jk)  )   & 
    117                   &            / ( e1t(ji,jj) * e2t(ji,jj) ) 
    118 #else 
    119113               hdivn(ji,jj,jk) =   & 
    120114                  (  e2u(ji,jj)*fse3u(ji,jj,jk) * un(ji,jj,jk) - e2u(ji-1,jj  )*fse3u(ji-1,jj  ,jk) * un(ji-1,jj  ,jk)       & 
    121115                   + e1v(ji,jj)*fse3v(ji,jj,jk) * vn(ji,jj,jk) - e1v(ji  ,jj-1)*fse3v(ji  ,jj-1,jk) * vn(ji  ,jj-1,jk)  )    & 
    122116                  / ( e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) ) 
    123 #endif 
    124117            END DO 
    125118         END DO 
    … …  
    266259      !!      - compute the now divergence given by : 
    267260      !!         hdivn = 1/(e1t*e2t*e3t) ( di[e2u*e3u un] + dj[e1v*e3v vn] ) 
    268       !!      Note: if lk_zco=T, e3u=e3v=e3t, they are simplified in the  
    269261      !!      above expression 
    270262      !!      - apply lateral boundary conditions on hdivn  
    … …  
    319311         DO jj = 2, jpjm1 
    320312            DO ji = fs_2, fs_jpim1   ! vector opt. 
    321 #if defined key_zco 
    322                hdivn(ji,jj,jk) = (  e2u(ji,jj) * un(ji,jj,jk) - e2u(ji-1,jj  ) * un(ji-1,jj  ,jk)      & 
    323                   &               + e1v(ji,jj) * vn(ji,jj,jk) - e1v(ji  ,jj-1) * vn(ji  ,jj-1,jk)  )   &  
    324                   / ( e1t(ji,jj) * e2t(ji,jj) ) 
    325 #else 
    326313               hdivn(ji,jj,jk) =   & 
    327314                  (  e2u(ji,jj)*fse3u(ji,jj,jk) * un(ji,jj,jk) - e2u(ji-1,jj  )*fse3u(ji-1,jj  ,jk) * un(ji-1,jj  ,jk)       & 
    328315                   + e1v(ji,jj)*fse3v(ji,jj,jk) * vn(ji,jj,jk) - e1v(ji  ,jj-1)*fse3v(ji  ,jj-1,jk) * vn(ji  ,jj-1,jk)  )    & 
    329316                  / ( e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) ) 
    330 #endif 
    331317            END DO   
    332318         END DO   

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DYN/dynldf_bilap.F90

    r1152 r2240  
    162162               zuf(ji,jj,jk) = fmask(ji,jj,jk) * (  zcv(ji+1,jj  ) - zcv(ji,jj)      & 
    163163                  &                            - zcu(ji  ,jj+1) + zcu(ji,jj)  )   & 
    164 #if defined key_zco 
    165                   &                         / ( e1f(ji,jj)*e2f(ji,jj) ) 
    166 #else 
    167164                  &       * fse3f(ji,jj,jk) / ( e1f(ji,jj)*e2f(ji,jj) ) 
    168 #endif 
    169165            END DO   
    170166         END DO   
    … …  
    173169         DO jj = 1, jpjm1 
    174170            DO ji = 1, fs_jpim1   ! vector opt. 
    175 #if defined key_zco 
    176                zlu(ji,jj,jk) = e2u(ji,jj) * zlu(ji,jj,jk) 
    177                zlv(ji,jj,jk) = e1v(ji,jj) * zlv(ji,jj,jk) 
    178 #else 
    179171               zlu(ji,jj,jk) = e2u(ji,jj) * fse3u(ji,jj,jk) * zlu(ji,jj,jk) 
    180172               zlv(ji,jj,jk) = e1v(ji,jj) * fse3v(ji,jj,jk) * zlv(ji,jj,jk) 
    181 #endif 
    182173            END DO 
    183174         END DO 
    … …  
    186177         DO jj = 2, jpj 
    187178            DO ji = fs_2, jpi   ! vector opt. 
    188 #if defined key_zco 
    189                zbt = e1t(ji,jj) * e2t(ji,jj) 
    190 #else 
    191179               zbt = e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) 
    192 #endif 
    193180               zut(ji,jj,jk) = (  zlu(ji,jj,jk) - zlu(ji-1,jj  ,jk)   & 
    194181                  &             + zlv(ji,jj,jk) - zlv(ji  ,jj-1,jk) ) / zbt 
    … …  
    210197         DO jj = 2, jpjm1 
    211198            DO ji = fs_2, fs_jpim1   ! vector opt. 
    212 #if defined key_zco 
    213                ze2u = e2u(ji,jj) 
    214                ze2v = e1v(ji,jj) 
    215 #else 
    216199               ze2u = e2u(ji,jj) * fse3u(ji,jj,jk) 
    217200               ze2v = e1v(ji,jj) * fse3v(ji,jj,jk) 
    218 #endif 
    219201               ! horizontal biharmonic diffusive trends 
    220202               zua = - ( zuf(ji  ,jj,jk) - zuf(ji,jj-1,jk) ) / ze2u   & 

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DYN/dynldf_lap.F90

    r1954 r2240  
    5050      !!         difu = 1/e1u di[ahmt hdivb] - 1/(e2u*e3u) dj-1[e3f ahmf rotb] 
    5151      !!         difv = 1/e2v dj[ahmt hdivb] + 1/(e1v*e3v) di-1[e3f ahmf rotb] 
    52       !!      If lk_zco=T, e3f=e3u=e3v, the vertical scale factor are simplified 
    5352      !!      in the rotational part of the diffusion. 
    5453      !!      Add this before trend to the general trend (ua,va): 
    … …  
    8786         DO jj = 2, jpjm1 
    8887            DO ji = fs_2, fs_jpim1   ! vector opt. 
    89 #if defined key_zco 
    90                ! horizontal diffusive trends 
    91                ze2u = rotb (ji,jj,jk)*fsahmf(ji,jj,jk) 
    92                ze1v = hdivb(ji,jj,jk)*fsahmt(ji,jj,jk) 
    93                zua = - ( ze2u - rotb (ji,jj-1,jk)*fsahmf(ji,jj-1,jk)                   ) / e2u(ji,jj)   & 
    94                      + ( hdivb(ji+1,jj,jk)*fsahmt(ji+1,jj,jk) - ze1v                   ) / e1u(ji,jj) 
    95  
    96                zva = + ( ze2u - rotb (ji-1,jj,jk)*fsahmf(ji-1,jj,jk)                   ) / e1v(ji,jj)   & 
    97                      + ( hdivb(ji,jj+1,jk)*fsahmt(ji,jj+1,jk) - ze1v                   ) / e2v(ji,jj) 
    98 #else 
    9988               ze2u = rotb (ji,jj,jk)*fsahmf(ji,jj,jk)*fse3f(ji,jj,jk) 
    10089               ze1v = hdivb(ji,jj,jk)*fsahmt(ji,jj,jk) 
    … …  
    10594               zva = + ( ze2u - rotb (ji-1,jj,jk)*fsahmf(ji-1,jj,jk)*fse3f(ji-1,jj,jk) ) / ( e1v(ji,jj) * fse3v(ji,jj,jk) )   & 
    10695                     + ( hdivb(ji,jj+1,jk)*fsahmt(ji,jj+1,jk) - ze1v                   ) / e2v(ji,jj) 
    107 #endif 
    10896 
    10997               ! add it to the general momentum trends 

  • branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90

    r2236 r2240  
    148148            END DO 
    149149         ELSE                            ! partial steps or s-ccordinate 
    150 # if defined key_zco 
    151             zc = gdept_0(jpkm1) 
    152 # else 
    153150            zc = MAXVAL( fsdept(:,:,jpkm1) ) 
    154 # endif 
    155151            IF( lk_mpp )   CALL mpp_max( zc )   ! max over the global domain 
    156152 

  • 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.