Changeset 2878 for branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/compute_sections.f90

Timestamp:

2011-09-30T08:26:28+02:00 (13 years ago)

Author:

cbricaud

Message:

 

File:

• branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/compute_sections.f90 (modified) (5 diffs)

branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/compute_sections.f90

@@ -2 +2 @@
    !!=====================================================================
    !!                       ***  MODULE  diadct  ***
-   !! Ocean diagnostics: Compute the transport trough a sec.
+   !! Ocean diagnostics: Compute the transport trough a section
    !!
+   !! History: 2011: Clement Bricaud, Mercator-Ocean
    !!
-   !! History: 2011: cbricaud Mercator-Ocean
    !!===============================================================
    !! * Modules used
@@ -301 +301 @@
            DO ji=1,jpi  
               coord_t=COORD_SECTION(glamf(ji,jj),gphif(ji,jj))
-!cbr              zdistFirst = distance(coord_t,coordFirst)
               zdistFirst = distance2(coord_t,coordFirst)
-!cbr              zdistLast = distance(coord_t,coordLast)
               zdistLast = distance2(coord_t,coordLast)
               IF( zdistFirst .LT. zdistmesh .AND. zdistFirst .LT. zdistante )THEN
@@ -370 +368 @@
 
      IF( sec%nb_point .ne. 0 )THEN
-     !IF( nb_inmesh .NE. 0 )THEN
 
         !The serie of mesh's points that form the section will 'link' 
@@ -393 +390 @@
 
         !initialize distnew value (with distance between section's extremities) 
-!cbr        zdistnew  = distance(coordFirst,coordLast,sec%ll_date_line) 
         zdistnew  = distance2(coordFirst,coordLast,sec%ll_date_line) 
         prevPoint  = POINT_SECTION(0,0)
@@ -596 +592 @@
                  CALL write_debug(jsec,cltmp)
                    
-                 ! on passe à l'ancienne methode: le point parmies les 4 pts (NSWE)  qui se rapproche
-                 ! le + du dernier pt
-                 !----------------------------- 
                  !be carreful! we can't go backward.
-
                  zdistNorth = zdistNorth2    ; zdistSouth = zdistSouth2
                  zdistEst   = zdistEst2      ; zdistWest  = zdistWest2