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 10843 for NEMO/branches/UKMO/dev_r10037_GPU/src/TOP/trcdta.F90 – NEMO

Ignore:
Timestamp:
2019-04-05T16:01:32+02:00 (5 years ago)
Author:
andmirek
Message:

ticket #2197 merge with dev_r9950_GO8_package at 10320

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/dev_r10037_GPU/src/TOP/trcdta.F90

    r9950 r10843  
    4040 
    4141   !!---------------------------------------------------------------------- 
    42    !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
     42   !! NEMO/TOP 4.0 , NEMO Consortium (2018) 
    4343   !! $Id$ 
    44    !! Software governed by the CeCILL licence     (./LICENSE) 
     44   !! Software governed by the CeCILL license (see ./LICENSE) 
    4545   !!---------------------------------------------------------------------- 
    4646CONTAINS 
     
    217217         ELSE                                !==   z- or zps- coordinate   ==! 
    218218            ! zps-coordinate (partial steps) interpolation at the last ocean level 
    219             IF( ln_zps ) THEN 
    220                DO jj = 1, jpj 
    221                   DO ji = 1, jpi 
    222                      ik = mbkt(ji,jj)  
    223                      IF( ik > 1 ) THEN 
    224                         zl = ( gdept_1d(ik) - gdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 
    225                         ptrcdta(ji,jj,ik) = (1.-zl) * ptrcdta(ji,jj,ik) + zl * ptrcdta(ji,jj,ik-1) 
    226                      ENDIF 
    227                      ik = mikt(ji,jj) 
    228                      IF( ik > 1 ) THEN 
    229                         zl = ( gdept_0(ji,jj,ik) - gdept_1d(ik) ) / ( gdept_1d(ik+1) - gdept_1d(ik) ) 
    230                         ptrcdta(ji,jj,ik) = (1.-zl) * ptrcdta(ji,jj,ik) + zl * ptrcdta(ji,jj,ik+1) 
    231                      ENDIF 
    232                   END DO 
    233               END DO 
    234             ENDIF 
     219!            IF( ln_zps ) THEN 
     220!               DO jj = 1, jpj 
     221!                  DO ji = 1, jpi 
     222!                     ik = mbkt(ji,jj)  
     223!                     IF( ik > 1 ) THEN 
     224!                        zl = ( gdept_1d(ik) - gdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 
     225!                        ptrcdta(ji,jj,ik) = (1.-zl) * ptrcdta(ji,jj,ik) + zl * ptrcdta(ji,jj,ik-1) 
     226!                     ENDIF 
     227!                     ik = mikt(ji,jj) 
     228!                     IF( ik > 1 ) THEN 
     229!                        zl = ( gdept_0(ji,jj,ik) - gdept_1d(ik) ) / ( gdept_1d(ik+1) - gdept_1d(ik) ) 
     230!                        ptrcdta(ji,jj,ik) = (1.-zl) * ptrcdta(ji,jj,ik) + zl * ptrcdta(ji,jj,ik+1) 
     231!                     ENDIF 
     232!                  END DO 
     233!              END DO 
     234!            ENDIF 
    235235            ! 
    236236         ENDIF 
Note: See TracChangeset for help on using the changeset viewer.