3 & LBi, UBi, UBj, Istr, Iend, Jrec, &
48 logical,
intent(out) :: update
50 integer,
intent(in) :: ng, model, ifield
51 integer,
intent(in) :: LBi, UBi, UBj, Istr, Iend, Jrec
53 real(r8),
intent(in) :: Finp(LBi:UBi,UBj,2)
55 real(r8),
intent(out) :: Fout(LBi:UBi,UBj)
61 integer :: Tindex, i, it1, it2, j
63 real(dp) :: SecScale, fac, fac1, fac2
71 lonerec=
linfo(3,ifield,ng)
72 tindex=
iinfo(8,ifield,ng)
82 fac1=anint((
tintrp(it2,ifield,ng)-
time(ng))*secscale,dp)
83 fac2=anint((
time(ng)-
tintrp(it1,ifield,ng))*secscale,dp)
90 fout(i,j)=finp(i,j,tindex)
96 ELSE IF (((fac1*fac2).ge.0.0_dp).and.(fac1+fac2).gt.0.0_dp)
THEN
97 fac=1.0_dp/(fac1+fac2)
102 fout(i,j)=fac1*finp(i,j,it1)+fac2*finp(i,j,it2)
125 10
FORMAT (/,
' SET_NGFLD - current model time', &
126 &
' exceeds ending value for variable: ',a, &
127 & /,14x,
'TDAYS = ',f15.4, &
128 & /,14x,
'Data Tmin = ',f15.4,2x,
'Data Tmax = ',f15.4, &
129 & /,14x,
'Data Tstr = ',f15.4,2x,
'Data Tend = ',f15.4, &
130 & /,14x,
'TINTRP1 = ',f15.4,2x,
'TINTRP2 = ',f15.4, &
131 & /,14x,
'FAC1 = ',f15.4,2x,
'FAC2 = ',f15.4)
logical, dimension(:,:,:), allocatable linfo
real(dp), dimension(:,:,:), allocatable tintrp
real(dp), dimension(:,:,:), allocatable finfo
character(len=maxlen), dimension(6, 0:nv) vname
integer, dimension(:,:,:), allocatable iinfo
real(dp), dimension(:), allocatable dt
logical, dimension(:), allocatable synchro_flag
real(dp), dimension(:), allocatable tdays
real(dp), parameter sec2day
real(dp), dimension(:), allocatable time
subroutine set_ngfld(ng, model, ifield, lbi, ubi, ubj, istr, iend, jrec, finp, fout, update)