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 2909 – NEMO

Changeset 2909


Ignore:
Timestamp:
2011-10-13T07:59:18+02:00 (13 years ago)
Author:
cbricaud
Message:

correction

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90

    r2908 r2909  
    897897        ENDIF 
    898898        !temperature classes transports 
    899         IF( ( sec%ztem(jcl+1) .NE. 99._wp     ) .AND. & 
     899        IF( ( sec%ztem(jcl) .NE. 99._wp     ) .AND. & 
    900900            ( sec%ztem(jcl+1) .NE. 99._wp     )       ) THEN 
    901901           classe = 'T       ' 
    … …  
    991991  !!K-1 | ptab(I,J,K-1) |                  |       interpolation between zbis and ptab(I+1,J,K)   
    992992  !!    |     .         |                  | 
    993   !!    |     .         |                  |       Z= ( 0.5*zet2*zbis + 0.5*zet1*ptab(I+1,J,K) )/(0.5*zet2+0.5*zet1)  
     993  !!    |     .         |                  |       interp = ( 0.5*zet2*zbis + 0.5*zet1*ptab(I+1,J,K) )/(0.5*zet2+0.5*zet1)  
    994994  !!    |     .         |                  | 
    995995  !!  ------------------------------------------ 
    … …  
    997997  !!    |     .         |                  | 
    998998  !!    |     .         |                  | 
    999   !!K   |    zbis...... U ..ptab(I+1,J,K)  | 
     999  !!K   |    zbis.......U...ptab(I+1,J,K)  | 
    10001000  !!    |     .         |                  | 
    10011001  !!    | ptab(I,J,K)   |                  | 
    … …  
    10091009  !!   ====>  s-coordinate 
    10101010  !!      
    1011   !!    |                |                  |   1. Compute distance between T1 and U points 
    1012   !!    |                |                  |      Compute distance between T2 and U points 
     1011  !!    |                |                  |   1. Compute distance between T1 and U points: SQRT( zdep1^2 + (0.5 * zet1 )^2 
     1012  !!    |                |                  |      Compute distance between T2 and U points: SQRT( zdep2^2 + (0.5 * zet2 )^2 
    10131013  !!    |                | ptab(I+1,J,K)    |  
    10141014  !!    |                |      T2          |   2. Interpolation between  T1 and T2 values at U point  
    … …  
    10631063     zdep2 = fsdept(ii2,ij2,kk) - zdepu 
    10641064 
     1065     !weights 
    10651066     zwgt1 = SQRT( ( 0.5 * zet1 ) * ( 0.5 * zet1 ) + ( zdep1 * zdep1 ) ) 
    10661067     zwgt2 = SQRT( ( 0.5 * zet2 ) * ( 0.5 * zet2 ) + ( zdep2 * zdep2 ) ) 
Note: See TracChangeset for help on using the changeset viewer.