Ignore:
Timestamp:
2011-09-30T08:26:28+02:00 (9 years ago)
Author:
cbricaud
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/compute_sections.f90

    r2877 r2878  
    22   !!===================================================================== 
    33   !!                       ***  MODULE  diadct  *** 
    4    !! Ocean diagnostics: Compute the transport trough a sec. 
     4   !! Ocean diagnostics: Compute the transport trough a section 
    55   !! 
     6   !! History: 2011: Clement Bricaud, Mercator-Ocean 
    67   !! 
    7    !! History: 2011: cbricaud Mercator-Ocean 
    88   !!=============================================================== 
    99   !! * Modules used 
     
    301301           DO ji=1,jpi   
    302302              coord_t=COORD_SECTION(glamf(ji,jj),gphif(ji,jj)) 
    303 !cbr              zdistFirst = distance(coord_t,coordFirst) 
    304303              zdistFirst = distance2(coord_t,coordFirst) 
    305 !cbr              zdistLast = distance(coord_t,coordLast) 
    306304              zdistLast = distance2(coord_t,coordLast) 
    307305              IF( zdistFirst .LT. zdistmesh .AND. zdistFirst .LT. zdistante )THEN 
     
    370368 
    371369     IF( sec%nb_point .ne. 0 )THEN 
    372      !IF( nb_inmesh .NE. 0 )THEN 
    373370 
    374371        !The serie of mesh's points that form the section will 'link'  
     
    393390 
    394391        !initialize distnew value (with distance between section's extremities)  
    395 !cbr        zdistnew  = distance(coordFirst,coordLast,sec%ll_date_line)  
    396392        zdistnew  = distance2(coordFirst,coordLast,sec%ll_date_line)  
    397393        prevPoint  = POINT_SECTION(0,0) 
     
    596592                 CALL write_debug(jsec,cltmp) 
    597593                    
    598                  ! on passe à l'ancienne methode: le point parmies les 4 pts (NSWE)  qui se rapproche 
    599                  ! le + du dernier pt 
    600                  !-----------------------------  
    601594                 !be carreful! we can't go backward. 
    602  
    603595                 zdistNorth = zdistNorth2    ; zdistSouth = zdistSouth2 
    604596                 zdistEst   = zdistEst2      ; zdistWest  = zdistWest2  
Note: See TracChangeset for help on using the changeset viewer.