Changeset 2877 for branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/compute_sections.f90
- Timestamp:
- 2011-09-29T18:54:19+02:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/compute_sections.f90
r2860 r2877 301 301 DO ji=1,jpi 302 302 coord_t=COORD_SECTION(glamf(ji,jj),gphif(ji,jj)) 303 zdistFirst = distance(coord_t,coordFirst) 304 zdistLast = distance(coord_t,coordLast) 303 !cbr zdistFirst = distance(coord_t,coordFirst) 304 zdistFirst = distance2(coord_t,coordFirst) 305 !cbr zdistLast = distance(coord_t,coordLast) 306 zdistLast = distance2(coord_t,coordLast) 305 307 IF( zdistFirst .LT. zdistmesh .AND. zdistFirst .LT. zdistante )THEN 306 308 sec%listPoint(1) = POINT_SECTION(ji,jj) … … 391 393 392 394 !initialize distnew value (with distance between section's extremities) 393 zdistnew = distance(coordFirst,coordLast,sec%ll_date_line) 395 !cbr zdistnew = distance(coordFirst,coordLast,sec%ll_date_line) 396 zdistnew = distance2(coordFirst,coordLast,sec%ll_date_line) 394 397 prevPoint = POINT_SECTION(0,0) 395 398 jseg = 1 … … 707 710 708 711 !now compute new slopeSection with ij-coordinates of first and last point 712 sec%slopeSection = 0 ! default value 709 713 IF( sec%nb_point .ne. 0 )THEN 710 714 IF ( sec%listPoint(sec%nb_point)%I .NE. sec%listPoint(1)%I ) THEN
Note: See TracChangeset
for help on using the changeset viewer.