/[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.13 by heimbach, Sat Jan 14 18:09:45 2006 UTC revision 1.14 by jmc, Sat Jan 14 20:47:43 2006 UTC
# Line 71  C subroutine variables Line 71  C subroutine variables
71        integer       method, mythid        integer       method, mythid
72    
73  C local variables  C local variables
74        real*8   ne_fac,nw_fac,se_fac,sw_fac  c     real*8   ne_fac,nw_fac,se_fac,sw_fac
75        integer  e_ind(snx,sny),w_ind(snx,sny)        integer  e_ind(snx,sny),w_ind(snx,sny)
76        integer  n_ind(snx,sny),s_ind(snx,sny)        integer  n_ind(snx,sny),s_ind(snx,sny)
77        real*8   px_ind(4), py_ind(4), ew_val(4)        real*8   px_ind(4), py_ind(4), ew_val(4)
# Line 82  C local variables Line 82  C local variables
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        _RS      threeSixtyRS
       _RL      threeSixtyRL  
85        PARAMETER ( threeSixtyRS = 360. )        PARAMETER ( threeSixtyRS = 360. )
       PARAMETER ( threeSixtyRL = 360. )  
86    
87  !     nomalize xG  C     put xG in interval [ lon_0 , lon_0+360 [
88        do bj=myByLo(myThid),myByHi(myThid)        do bj=myByLo(myThid),myByHi(myThid)
89         do bi=myBxLo(myThid),myBxHi(myThid)         do bi=myBxLo(myThid),myBxHi(myThid)
90          do j=1-OLy,sNy+OLy          do j=1-OLy,sNy+OLy
91           do i=1-OLx,sNx+OLx           do i=1-OLx,sNx+OLx
92            xG(i,j,bi,bj) = xG_in(i,j,bi,bj)            xG(i,j,bi,bj) = xG_in(i,j,bi,bj)-lon_0
93            xG(i,j,bi,bj) = mod(xG(i,j,bi,bj),threeSixtyRS)       &                  + threeSixtyRS*2.
94            if ( xG(i,j,bi,bj) .lt. 0. )            xG(i,j,bi,bj) = lon_0+mod(xG(i,j,bi,bj),threeSixtyRS)
      &     xG(i,j,bi,bj) = threeSixtyRS+xG(i,j,bi,bj)  
95           enddo           enddo
96          enddo          enddo
97         enddo         enddo
# Line 114  C     _BEGIN_MASTER( myThid ) Line 111  C     _BEGIN_MASTER( myThid )
111    
112  C setup input grid  C setup input grid
113         do i=-1,nx_in+2         do i=-1,nx_in+2
114          x_in(i) = lon_0 + (i-1.)*lon_inc          x_in(i) = lon_0 + (i-1)*lon_inc
         x_in(i) = mod(x_in(i),threeSixtyRL)  
         if ( x_in(i) .lt. 0. )  
      &   x_in(i) = threeSixtyRL+x_in(i)  
115         enddo         enddo
116    
117         y_in(0) = lat_0 - lat_inc(1)         y_in(0) = lat_0 - lat_inc(1)
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22