198 |
xcoord((bi-1)*sNx+I) = xc(I,1,bi,1) |
xcoord((bi-1)*sNx+I) = xc(I,1,bi,1) |
199 |
ENDDO |
ENDDO |
200 |
ENDDO |
ENDDO |
201 |
CALL WRITE_1D_R8( xcoord, Nx, INDEX_I,'xcoord = ', |
CALL WRITE_1D_R8( xcoord, sNx*nSx, INDEX_I,'xcoord = ', |
202 |
&' /* P-point X coordinate ( m - cartesian, degrees - spherical ) */') |
&' /* P-point X coordinate ( m - cartesian, degrees - spherical ) */') |
203 |
DO bj=1,nSy |
DO bj=1,nSy |
204 |
DO J=1,sNy |
DO J=1,sNy |
205 |
ycoord((bj-1)*sNy+J) = yc(1,J,1,bj) |
ycoord((bj-1)*sNy+J) = yc(1,J,1,bj) |
206 |
ENDDO |
ENDDO |
207 |
ENDDO |
ENDDO |
208 |
CALL WRITE_1D_R8( ycoord, Ny, INDEX_J,'ycoord = ', |
CALL WRITE_1D_R8( ycoord, sNy*nSy, INDEX_J,'ycoord = ', |
209 |
&' /* P-point Y coordinate ( m - cartesian, degrees - spherical ) */') |
&' /* P-point Y coordinate ( m - cartesian, degrees - spherical ) */') |
210 |
DO K=1,Nz |
DO K=1,Nz |
211 |
zcoord(K) = zc(K) |
zcoord(K) = zc(K) |