/[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.12 by heimbach, Fri Jan 13 16:48:03 2006 UTC revision 1.13 by heimbach, Sat Jan 14 18:09:45 2006 UTC
# Line 36  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC Line 36  CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
36       I   infile,       I   infile,
37       I   filePrec,       I   filePrec,
38       O   arrayout,       O   arrayout,
39       I   irecord, xG, yG,       I   irecord, xG_in, yG,
40       I   lon_0, lon_inc,       I   lon_0, lon_inc,
41       I   lat_0, lat_inc,       I   lat_0, lat_inc,
42       I   nx_in, ny_in, method, mythid)       I   nx_in, ny_in, method, mythid)
# 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)        _RS      xG(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
84          _RS      threeSixtyRS
85          _RL      threeSixtyRL
86          PARAMETER ( threeSixtyRS = 360. )
87          PARAMETER ( threeSixtyRL = 360. )
88    
89  !     nomalize xg  !     nomalize xG
90        do bj=myByLo(myThid),myByHi(myThid)        do bj=myByLo(myThid),myByHi(myThid)
91         do bi=myBxLo(myThid),myBxHi(myThid)         do bi=myBxLo(myThid),myBxHi(myThid)
92          do j=1-OLy,sNy+OLy          do j=1-OLy,sNy+OLy
93           do i=1-OLx,sNx+OLx           do i=1-OLx,sNx+OLx
94            xg(i,j,bi,bj) = xg_in(i,j,bi,bj)            xG(i,j,bi,bj) = xG_in(i,j,bi,bj)
95            xg(i,j,bi,bj) = mod(xg(i,j,bi,bj),360.)            xG(i,j,bi,bj) = mod(xG(i,j,bi,bj),threeSixtyRS)
96            if ( xg(i,j,bi,bj) .lt. 0. )            if ( xG(i,j,bi,bj) .lt. 0. )
97       &     xg(i,j,bi,bj) = 360.+xg(i,j,bi,bj)       &     xG(i,j,bi,bj) = threeSixtyRS+xG(i,j,bi,bj)
98           enddo           enddo
99          enddo          enddo
100         enddo         enddo
# Line 111  C     _BEGIN_MASTER( myThid ) Line 115  C     _BEGIN_MASTER( myThid )
115  C setup input grid  C setup input grid
116         do i=-1,nx_in+2         do i=-1,nx_in+2
117          x_in(i) = lon_0 + (i-1.)*lon_inc          x_in(i) = lon_0 + (i-1.)*lon_inc
118          x_in(i) = mod(x_in(i),360.)          x_in(i) = mod(x_in(i),threeSixtyRL)
119          if ( x_in(i) .lt. 0. )          if ( x_in(i) .lt. 0. )
120       &   x_in(i) = 360.+x_in(i)       &   x_in(i) = threeSixtyRL+x_in(i)
121         enddo         enddo
122    
123         y_in(0) = lat_0 - lat_inc(1)         y_in(0) = lat_0 - lat_inc(1)
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22