Changeset 2948 for branches/2011/dev_r2802_MERCATOR10_diadct

Timestamp:

2011-10-17T18:21:39+02:00 (13 years ago)

Author:

cbricaud

Message:

corrections after review

Location:

branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src

Files:

• compute_sections.f90 (modified) (4 diffs)
• sections_tools.f90 (modified) (2 diffs)

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

@@ -51 +51 @@
         nmesh                    ! number of cells in processor domain
      INTEGER :: itest , jtest    ! dummy integer
-     REAL(wp),SAVE :: zdistmesh  ! Taller cell of the mesh in ocean
+     REAL(wp),SAVE :: zdistmesh  ! Taller cell of the mesh in ocean in degrees
      REAL(wp)      :: &
         zdistEst   , zdistNorth , zdistWest , zdistSouth ,  &! temporary scalars
@@ -106 +106 @@
 
      IF( jsec == 1 )THEN
-        !Find the taller cell of the mesh in ocean (used in compsec)
+        !Find the taller cell of the mesh in ocean (used in compsec), 
+        !unit= degrees 
         zdistmesh=0.
         DO jj=2,jpj
@@ -161 +162 @@
               !intersection section/[AB]=?
               coordTemp     = intersec(sec,coord_a,coord_b)     !compute intersection 
-              coordTemp%lon = coordTemp%lon !* zmask(ji,jj)-9999.*(1-zmask(ji,jj))
-              coordTemp%lat = coordTemp%lat !* zmask(ji,jj)-9999.*(1-zmask(ji,jj))
  
               IF( coordTemp%lon .NE. -9999 )THEN
@@ -474 +473 @@
               !------------------------------------------------------------------
               ! BE CAREFUL! When the domain crosses the date line (ll_date_domain):
+              !
               ! When we will compute distances between W/E/S/N points and endingPoint,
               ! we have to check if theses 4 lines crosses the date line
               ! (test: ABS(coordTemp%lon - coordLast%lon).GT. 180)
-              ! If it's true,we have to add 360° to coordLast%long to compute the 
-              ! correct distance through the date line and not through the complementary
-              ! in the mesh.
+              !
+              ! If on of this 4 lines  crosses the date line, we have to add 360° its 
+              ! extremities longitudes to compute the correct distance through the date line
+              ! and not around the Earth
        
               ! c.1 compute distWest: distance between west point and endingPoint

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

@@ -303 +303 @@
           !3. Compute intersection (zX,zY) !
           !================================!
-          IF( za1 == 1.e20 ) THEN ! Case X=cste 
+          IF( za1 == 1.e20 ) THEN ! Case X=constant
               zX = sec%coordSec(1)%lon
               zY = za2 * zX + zb2
-          ELSE IF ( za2 == 1.e20) THEN ! Case X=cste
+          ELSE IF ( za2 == 1.e20) THEN ! Case X=constante
               zX = coord_b%lon
               zY = za1 * zX + zb1
@@ -387 +387 @@
      !!
      !! ** Purpose : organize a list point by latitude , and after by longitude
-     !!
-     !!     History
+     !!              using "numerical recipies " routine.
      !!---------------------------------------------------------------------
      !! * arguments