/[MITgcm]/MITgcm/pkg/exf/exf_interp.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_interp.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.11 by cnh, Tue Nov 8 15:53:41 2005 UTC revision 1.12 by heimbach, Fri Jan 13 16:48:03 2006 UTC
# Line 64  C subroutine variables Line 64  C subroutine variables
64        character*(*) infile        character*(*) infile
65        integer       filePrec, irecord, nx_in, ny_in        integer       filePrec, irecord, nx_in, ny_in
66        _RL           arrayout(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL           arrayout(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
67        _RS           xG      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS           xG_in      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
68        _RS           yG      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS           yG      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
69        _RL           lon_0, lon_inc        _RL           lon_0, lon_inc
70        _RL           lat_0, lat_inc(ny_in-1)        _RL           lat_0, lat_inc(ny_in-1)
# Line 80  C local variables Line 80  C local variables
80        real*4   arrayin(-1:nx_in+2 ,      -1:ny_in+2)        real*4   arrayin(-1:nx_in+2 ,      -1:ny_in+2)
81        real*8   x_in   (-1:nx_in+2), y_in(-1:ny_in+2)        real*8   x_in   (-1:nx_in+2), y_in(-1:ny_in+2)
82        integer  i, j, k, l, js, bi, bj, sp, interp_unit        integer  i, j, k, l, js, bi, bj, sp, interp_unit
83          _RS      xg(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
84    
85    !     nomalize xg
86          do bj=myByLo(myThid),myByHi(myThid)
87           do bi=myBxLo(myThid),myBxHi(myThid)
88            do j=1-OLy,sNy+OLy
89             do i=1-OLx,sNx+OLx
90              xg(i,j,bi,bj) = xg_in(i,j,bi,bj)
91              xg(i,j,bi,bj) = mod(xg(i,j,bi,bj),360.)
92              if ( xg(i,j,bi,bj) .lt. 0. )
93         &     xg(i,j,bi,bj) = 360.+xg(i,j,bi,bj)
94             enddo
95            enddo
96           enddo
97          enddo
98    
99         call exf_interp_read(         call exf_interp_read(
100       I   infile,       I   infile,
# Line 96  C     _BEGIN_MASTER( myThid ) Line 111  C     _BEGIN_MASTER( myThid )
111  C setup input grid  C setup input grid
112         do i=-1,nx_in+2         do i=-1,nx_in+2
113          x_in(i) = lon_0 + (i-1.)*lon_inc          x_in(i) = lon_0 + (i-1.)*lon_inc
114            x_in(i) = mod(x_in(i),360.)
115            if ( x_in(i) .lt. 0. )
116         &   x_in(i) = 360.+x_in(i)
117         enddo         enddo
118    
119         y_in(0) = lat_0 - lat_inc(1)         y_in(0) = lat_0 - lat_inc(1)
120         y_in(-1)= lat_0 - 2.*lat_inc(1)         y_in(-1)= lat_0 - 2.*lat_inc(1)
121         y_in(1) = lat_0         y_in(1) = lat_0
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22