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 |
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) ) |
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 |