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

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

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

revision 1.19 by heimbach, Wed Apr 18 13:24:28 2007 UTC revision 1.20 by heimbach, Wed Apr 18 15:34:40 2007 UTC
# Line 147  C     Complete filename with YR or _YEAR Line 147  C     Complete filename with YR or _YEAR
147                 genfile0 = genfile                 genfile0 = genfile
148              endif              endif
149    
   
150  #ifdef USE_EXF_INTERPOLATION  #ifdef USE_EXF_INTERPOLATION
151              call exf_interp( genfile0, exf_iprec              call exf_interp( genfile0, exf_iprec
152       &           , gen1, count0, gen_xout, gen_yout       &           , gen1, count0, gen_xout, gen_yout
# Line 215  c     Loop over tiles. Line 214  c     Loop over tiles.
214             do j = 1,sny             do j = 1,sny
215              do i = 1,snx              do i = 1,snx
216  c     Interpolate linearly onto the  time.  c     Interpolate linearly onto the  time.
 cph             if ( genfile .EQ. climsstfile ) then  
 cph                if (gen0(i,j,bi,bj) .lt. climtempfreeze) then  
 cph                   gen0(i,j,bi,bj) = climtempfreeze  
 cph                endif  
 cph                if (gen1(i,j,bi,bj) .lt. climtempfreeze) then  
 cph                   gen1(i,j,bi,bj) = climtempfreeze  
 cph                endif  
 cph             endif  
217               genfld(i,j,bi,bj) = exf_inscal_gen * (               genfld(i,j,bi,bj) = exf_inscal_gen * (
218       &                          fac * gen0(i,j,bi,bj) +       &                          fac * gen0(i,j,bi,bj) +
219       &              (exf_one - fac) * gen1(i,j,bi,bj) )       &              (exf_one - fac) * gen1(i,j,bi,bj) )
# Line 230  cph             endif Line 221  cph             endif
221       &            genfld(i,j,bi,bj) -       &            genfld(i,j,bi,bj) -
222       &            exf_inscal_gen * ( genremove_intercept +       &            exf_inscal_gen * ( genremove_intercept +
223       &            genremove_slope*(mytime-starttime) )       &            genremove_slope*(mytime-starttime) )
                   enddo  
                enddo  
224              enddo              enddo
225               enddo
226              enddo
227           enddo           enddo
228    
229        endif        endif
Legend:
Removed from v.1.19  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22