36 |
! (latitudes, in degrees or radians, in increasing or decreasing |
! (latitudes, in degrees or radians, in increasing or decreasing |
37 |
! order, from pole to pole) |
! order, from pole to pole) |
38 |
|
|
|
LOGICAL, intent(in), optional:: interbar |
|
39 |
REAL, intent(out):: xf(:), yf(:) ! longitudes and latitudes, in rad |
REAL, intent(out):: xf(:), yf(:) ! longitudes and latitudes, in rad |
40 |
REAL, intent(inout):: champd(:, :) |
REAL, intent(inout):: champd(:, :) |
41 |
|
LOGICAL, intent(in), optional:: interbar |
42 |
|
|
43 |
! Variables locales: |
! Local: |
|
|
|
44 |
INTEGER lons, lats |
INTEGER lons, lats |
45 |
LOGICAL radianlon ! "xd" is in degrees |
LOGICAL radianlon ! "xd" is in degrees |
46 |
logical invlon ! "xd" contains longitudes between 0 and 2 pi |
logical invlon ! "xd" contains longitudes between 0 and 2 pi |
103 |
end do |
end do |
104 |
|
|
105 |
xf(i:) = xf(i:) - 2 * pi |
xf(i:) = xf(i:) - 2 * pi |
106 |
xf(:) = cshift(xf, shift=i - 1) |
xf(:) = cshift(xf, shift = i - 1) |
107 |
champd(:, :) = cshift(champd, shift=i - 1) |
champd(:, :) = cshift(champd, shift = i - 1) |
108 |
ENDIF |
ENDIF |
109 |
|
|
110 |
IF (yd(1) < yd(lats)) THEN |
IF (yd(1) < yd(lats)) THEN |