Changeset 6619


Ignore:
Timestamp:
2016-05-25T14:52:53+02:00 (5 years ago)
Author:
cetlod
Message:

Minor bug correction, see ticket #1667

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/TOP_SRC/trcdta.F90

    r6607 r6619  
    199199                  DO ji = 1, jpi 
    200200                     DO jk = 1, jpk                        ! determines the intepolated T-S profiles at each (i,j) points 
    201                         zl = fsdept_n(ji,jj,jk) 
     201                        zl = gdept_n(ji,jj,jk) 
    202202                        IF(     zl < gdept_1d(1  ) ) THEN         ! above the first level of data 
    203203                           ztp(jk) =  ptrc(ji,jj,1,jl) 
     
    230230                        ik = mbkt(ji,jj)  
    231231                        IF( ik > 1 ) THEN 
    232                            zl = ( gdept_1d(ik) - fsdept_n(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 
     232                           zl = ( gdept_1d(ik) - gdept_n(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 
    233233                           ptrc(ji,jj,ik,jl) = (1.-zl) * ptrc(ji,jj,ik,jl) + zl * ptrc(ji,jj,ik-1,jl) 
    234234                        ENDIF 
Note: See TracChangeset for help on using the changeset viewer.