Ignore:
Timestamp:
2013-11-20T18:04:09+01:00 (8 years ago)
Author:
cetlod
Message:

dev_MERGE_2013 : changes to improve compilation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/C1D/dtauvd.F90

    r4245 r4294  
    174174            DO ji = 1, jpi                ! determines the interpolated U & V current profiles at each (i,j) point 
    175175               DO jk = 1, jpk 
    176                   zl = fsdept_0(ji,jj,jk) 
    177                   IF    ( zl < gdept_0(1  ) ) THEN          ! extrapolate above the first level of data 
     176                  zl = fsdept(ji,jj,jk) 
     177                  IF    ( zl < gdept_1d(1  ) ) THEN          ! extrapolate above the first level of data 
    178178                     zup(jk) =  puvd(ji,jj,1    ,1) 
    179179                     zvp(jk) =  puvd(ji,jj,1    ,2) 
    180                   ELSEIF( zl > gdept_0(jpk) ) THEN          ! extrapolate below the last level of data 
     180                  ELSEIF( zl > gdept_1d(jpk) ) THEN          ! extrapolate below the last level of data 
    181181                     zup(jk) =  puvd(ji,jj,jpkm1,1) 
    182182                     zvp(jk) =  puvd(ji,jj,jpkm1,2) 
    183183                  ELSE                                      ! inbetween : vertical interpolation between jkk & jkk+1 
    184184                     DO jkk = 1, jpkm1                      ! when  gdept(jkk) < zl < gdept(jkk+1) 
    185                         IF( (zl-gdept_0(jkk)) * (zl-gdept_0(jkk+1)) <= 0._wp ) THEN 
    186                            zi = ( zl - gdept_0(jkk) ) / (gdept_0(jkk+1)-gdept_0(jkk)) 
     185                        IF( (zl-gdept_1d(jkk)) * (zl-gdept_1d(jkk+1)) <= 0._wp ) THEN 
     186                           zi = ( zl - gdept_1d(jkk) ) / (gdept_1d(jkk+1)-gdept_1d(jkk)) 
    187187                           zup(jk) = puvd(ji,jj,jkk,1) + ( puvd(ji,jj,jkk+1,1 ) - puvd(ji,jj,jkk,1) ) * zi  
    188188                           zvp(jk) = puvd(ji,jj,jkk,2) + ( puvd(ji,jj,jkk+1,2 ) - puvd(ji,jj,jkk,2) ) * zi 
     
    212212                  ik = mbkt(ji,jj)  
    213213                  IF( ik > 1 ) THEN 
    214                      zl = ( gdept_0(ik) - fsdept_0(ji,jj,ik) ) / ( gdept_0(ik) - gdept_0(ik-1) ) 
     214                     zl = ( gdept_1d(ik) - gdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 
    215215                     puvd(ji,jj,ik,1) = (1.-zl) * puvd(ji,jj,ik,1) + zl * puvd(ji,jj,ik-1,1) 
    216216                     puvd(ji,jj,ik,2) = (1.-zl) * puvd(ji,jj,ik,2) + zl * puvd(ji,jj,ik-1,2) 
Note: See TracChangeset for help on using the changeset viewer.