Changeset 8543


Ignore:
Timestamp:
2017-09-19T12:27:07+02:00 (3 years ago)
Author:
cetlod
Message:

bugfix on trcdta to be consistent with dtatsd

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/TOP_SRC/trcdta.F90

    r6688 r8543  
    200200               DO ji = 1, jpi 
    201201                  DO jk = 1, jpk                        ! determines the intepolated T-S profiles at each (i,j) points 
    202                      zl = fsdept_n(ji,jj,jk) 
     202                     zl = gdept_0(ji,jj,jk) 
    203203                     IF(     zl < gdept_1d(1  ) ) THEN         ! above the first level of data 
    204204                        ztp(jk) = ztrcdta(ji,jj,1) 
     
    229229                     ik = mbkt(ji,jj)  
    230230                     IF( ik > 1 ) THEN 
    231                         zl = ( gdept_1d(ik) - fsdept_n(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 
     231                        zl = ( gdept_1d(ik) - gdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 
    232232                        ztrcdta(ji,jj,ik) = (1.-zl) * ztrcdta(ji,jj,ik) + zl * ztrcdta(ji,jj,ik-1) 
    233233                     ENDIF 
    234234                     ik = mikt(ji,jj) 
    235235                     IF( ik > 1 ) THEN 
    236                         zl = ( fsdept_n(ji,jj,ik) - gdept_1d(ik) ) / ( gdept_1d(ik+1) - gdept_1d(ik) ) 
     236                        zl = ( gdept_0(ji,jj,ik) - gdept_1d(ik) ) / ( gdept_1d(ik+1) - gdept_1d(ik) ) 
    237237                        ztrcdta(ji,jj,ik) = (1.-zl) * ztrcdta(ji,jj,ik) + zl * ztrcdta(ji,jj,ik+1) 
    238238                     ENDIF 
Note: See TracChangeset for help on using the changeset viewer.