Changeset 2948
- Timestamp:
- 2011-10-17T18:21:39+02:00 (13 years ago)
- Location:
- branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/compute_sections.f90
r2947 r2948 51 51 nmesh ! number of cells in processor domain 52 52 INTEGER :: itest , jtest ! dummy integer 53 REAL(wp),SAVE :: zdistmesh ! Taller cell of the mesh in ocean 53 REAL(wp),SAVE :: zdistmesh ! Taller cell of the mesh in ocean in degrees 54 54 REAL(wp) :: & 55 55 zdistEst , zdistNorth , zdistWest , zdistSouth , &! temporary scalars … … 106 106 107 107 IF( jsec == 1 )THEN 108 !Find the taller cell of the mesh in ocean (used in compsec) 108 !Find the taller cell of the mesh in ocean (used in compsec), 109 !unit= degrees 109 110 zdistmesh=0. 110 111 DO jj=2,jpj … … 161 162 !intersection section/[AB]=? 162 163 coordTemp = intersec(sec,coord_a,coord_b) !compute intersection 163 coordTemp%lon = coordTemp%lon !* zmask(ji,jj)-9999.*(1-zmask(ji,jj))164 coordTemp%lat = coordTemp%lat !* zmask(ji,jj)-9999.*(1-zmask(ji,jj))165 164 166 165 IF( coordTemp%lon .NE. -9999 )THEN … … 474 473 !------------------------------------------------------------------ 475 474 ! BE CAREFUL! When the domain crosses the date line (ll_date_domain): 475 ! 476 476 ! When we will compute distances between W/E/S/N points and endingPoint, 477 477 ! we have to check if theses 4 lines crosses the date line 478 478 ! (test: ABS(coordTemp%lon - coordLast%lon).GT. 180) 479 ! If it's true,we have to add 360° to coordLast%long to compute the 480 ! correct distance through the date line and not through the complementary 481 ! in the mesh. 479 ! 480 ! If on of this 4 lines crosses the date line, we have to add 360° its 481 ! extremities longitudes to compute the correct distance through the date line 482 ! and not around the Earth 482 483 483 484 ! c.1 compute distWest: distance between west point and endingPoint -
branches/2011/dev_r2802_MERCATOR10_diadct/NEMOGCM/TOOLS/SECTIONS_DIADCT/src/sections_tools.f90
r2947 r2948 303 303 !3. Compute intersection (zX,zY) ! 304 304 !================================! 305 IF( za1 == 1.e20 ) THEN ! Case X=c ste305 IF( za1 == 1.e20 ) THEN ! Case X=constant 306 306 zX = sec%coordSec(1)%lon 307 307 zY = za2 * zX + zb2 308 ELSE IF ( za2 == 1.e20) THEN ! Case X=c ste308 ELSE IF ( za2 == 1.e20) THEN ! Case X=constante 309 309 zX = coord_b%lon 310 310 zY = za1 * zX + zb1 … … 387 387 !! 388 388 !! ** Purpose : organize a list point by latitude , and after by longitude 389 !! 390 !! History 389 !! using "numerical recipies " routine. 391 390 !!--------------------------------------------------------------------- 392 391 !! * arguments
Note: See TracChangeset
for help on using the changeset viewer.