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

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

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

revision 1.16 by dimitri, Fri Jan 25 02:43:19 2008 UTC revision 1.17 by mlosch, Fri Jan 25 16:02:56 2008 UTC
# Line 120  c     get record numbers and interpolati Line 120  c     get record numbers and interpolati
120       &        )       &        )
121    
122            if ( first ) then            if ( first ) then
123              if (useExfYearlyFields) then  
124  C     Complete filenames with YR or _YEAR extension              call exf_GetYearlyFieldName(
125                il = ilnblnk( uvecfile )       I         useExfYearlyFields, twoDigitYear, uvecperiod, year0,
126                if (twoDigitYear) then       I         uvecfile,
127                  if (year0.ge.2000) then       O         uvecfile0,
128                    write(uvecfile0(1:128),'(a,i2.2)')       I         mytime, myiter, mythid )
129       &                uvecfile(1:il),year0-2000              call exf_GetYearlyFieldName(
130                  else       I         useExfYearlyFields, twoDigitYear, vvecperiod, year0,
131                    write(uvecfile0(1:128),'(a,i2.2)')       I         vvecfile,
132       &                uvecfile(1:il),year0-1900       O         vvecfile0,
133                  endif       I         mytime, myiter, mythid )
134                else  
                 write(uvecfile0(1:128),'(2a,i4.4)')  
      &              uvecfile(1:il),'_',year0  
               endif  
               il = ilnblnk( vvecfile )  
               if (twoDigitYear) then  
                 if (year0.ge.2000) then  
                   write(vvecfile0(1:128),'(a,i2.2)')  
      &                vvecfile(1:il),year0-2000  
                 else  
                   write(vvecfile0(1:128),'(a,i2.2)')  
      &                vvecfile(1:il),year0-1900  
                 endif  
               else  
                 write(vvecfile0(1:128),'(2a,i4.4)')  
      &              vvecfile(1:il),'_',year0  
               endif  
             else  
               uvecfile0 = uvecfile  
               vvecfile0 = vvecfile  
             endif  
135  c     scalar interpolation to (xC,yC) locations  c     scalar interpolation to (xC,yC) locations
136              interp_method=12              interp_method=12
137              call exf_interp( uvecfile0, exf_iprec              call exf_interp( uvecfile0, exf_iprec
# Line 219  c     apply mask Line 199  c     apply mask
199            if (( first ) .or. ( changed )) then            if (( first ) .or. ( changed )) then
200              call exf_SwapFFields( uvec0, uvec1, mythid )              call exf_SwapFFields( uvec0, uvec1, mythid )
201              call exf_SwapFFields( vvec0, vvec1, mythid )              call exf_SwapFFields( vvec0, vvec1, mythid )
202              if (useExfYearlyFields) then  
203  C     Complete filenames with YR or _YEAR extension              call exf_GetYearlyFieldName(
204                il = ilnblnk( uvecfile )       I         useExfYearlyFields, twoDigitYear, uvecperiod, year1,
205                if (twoDigitYear) then       I         uvecfile,
206                  if (year1.ge.2000) then       O         uvecfile1,
207                    write(uvecfile1(1:128),'(a,i2.2)')       I         mytime, myiter, mythid )
208       &                uvecfile(1:il),year1-2000              call exf_GetYearlyFieldName(
209                  else       I         useExfYearlyFields, twoDigitYear, vvecperiod, year1,
210                    write(uvecfile1(1:128),'(a,i2.2)')       I         vvecfile,
211       &                uvecfile(1:il),year1-1900       O         vvecfile1,
212                  endif       I         mytime, myiter, mythid )
213                else  
                 write(uvecfile1(1:128),'(2a,i4.4)')  
      &              uvecfile(1:il),'_',year1  
               endif  
               il = ilnblnk( vvecfile )  
               if (twoDigitYear) then  
                 if (year1.ge.2000) then  
                   write(vvecfile1(1:128),'(a,i2.2)')  
      &                vvecfile(1:il),year1-2000  
                 else  
                   write(vvecfile1(1:128),'(a,i2.2)')  
      &                vvecfile(1:il),year1-1900  
                 endif  
               else  
                 write(vvecfile1(1:128),'(2a,i4.4)')  
      &              vvecfile(1:il),'_',year1  
               endif  
             else  
               uvecfile1 = uvecfile  
               vvecfile1 = vvecfile  
             endif  
214  c     scalar interpolation to (xC,yC) locations  c     scalar interpolation to (xC,yC) locations
215              interp_method=12              interp_method=12
216              call exf_interp( uvecfile1, exf_iprec              call exf_interp( uvecfile1, exf_iprec
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22