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

Timestamp:

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

Author:

cbricaud

Message:

corrections after review

File:

• branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/compute_sections.f90 (modified) (4 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