diff --git a/lis/metforcing/gefs/read_gefs_operational.F90 b/lis/metforcing/gefs/read_gefs_operational.F90 index fa1fe7d94..42c9e8fd2 100644 --- a/lis/metforcing/gefs/read_gefs_operational.F90 +++ b/lis/metforcing/gefs/read_gefs_operational.F90 @@ -329,9 +329,10 @@ subroutine read_gefs_operational(n, m, findex, order, filename, ferror) call gefs_shift_longitude( gefs_struc(n)%nc, gefs_struc(n)%nr, & numpts, gefs_grib_data ) - pcp_flag = .true. ! Spatially interp GEFS forcing field to LIS domain: + pcp_flag = .true. call interp_gefs(n, findex, gefs_grib_data, pcp_flag, varfield ) + pcp_flag = .false. do r=1,LIS_rc%lnr(n) do c=1,LIS_rc%lnc(n) diff --git a/lis/metforcing/gefs/timeinterp_gefs.F90 b/lis/metforcing/gefs/timeinterp_gefs.F90 index 6cb3a4fe3..610c77702 100644 --- a/lis/metforcing/gefs/timeinterp_gefs.F90 +++ b/lis/metforcing/gefs/timeinterp_gefs.F90 @@ -76,7 +76,7 @@ subroutine timeinterp_gefs(n,findex) call LIS_time2date(btime,bdoy,gmt2,byr,bmo,bda,bhr,bmn) - !== Interpolate data in time + !== Interpolate data in time ! Check if bookend times differ, else stop ... if( (gefs_struc(n)%fcsttime2-gefs_struc(n)%fcsttime1)==0 ) then @@ -130,7 +130,6 @@ subroutine timeinterp_gefs(n,findex) call ESMF_FieldGet(swdField,localDE=0,farrayPtr=swd,rc=status) call LIS_verify(status) - zdoy=LIS_rc%doy do t=1,LIS_rc%ntiles(n)/LIS_rc%nensem(n) do m=1,gefs_struc(n)%max_ens_members do k=1,mfactor @@ -138,6 +137,7 @@ subroutine timeinterp_gefs(n,findex) index1 = LIS_domain(n)%tile(tid)%index ! Compute and apply zenith angle weights + zdoy=LIS_rc%doy call zterp( 0, LIS_domain(n)%grid(index1)%lat, & LIS_domain(n)%grid(index1)%lon, gmt1, gmt2, & LIS_rc%gmt,zdoy,zw1,zw2,czb,cze,czm,LIS_rc ) diff --git a/lis/testcases/metforcing/gefs/MODEL_OUTPUT_LIST.TBL b/lis/testcases/metforcing/gefs/MODEL_OUTPUT_LIST.TBL old mode 100755 new mode 100644 diff --git a/lis/testcases/metforcing/gefs/README b/lis/testcases/metforcing/gefs/README old mode 100755 new mode 100644 diff --git a/lis/testcases/metforcing/gefs/input.ctl b/lis/testcases/metforcing/gefs/input.ctl old mode 100755 new mode 100644 diff --git a/lis/testcases/metforcing/gefs/input_testcase.ctl b/lis/testcases/metforcing/gefs/input_testcase.ctl old mode 100755 new mode 100644 diff --git a/lis/testcases/metforcing/gefs/ldt.config b/lis/testcases/metforcing/gefs/ldt.config old mode 100755 new mode 100644 diff --git a/lis/testcases/metforcing/gefs/lis.config b/lis/testcases/metforcing/gefs/lis.config old mode 100755 new mode 100644 diff --git a/lis/testcases/metforcing/gefs/output.ctl b/lis/testcases/metforcing/gefs/output.ctl old mode 100755 new mode 100644 index b465a4196..cc48988ec --- a/lis/testcases/metforcing/gefs/output.ctl +++ b/lis/testcases/metforcing/gefs/output.ctl @@ -1,4 +1,4 @@ -DSET ./OUTPUT/SURFACEMODEL/%y4%m2/LIS_HIST_%y4%m2%d2%h200.d01.nc +DSET ^OUTPUT/SURFACEMODEL/%y4%m2/LIS_HIST_%y4%m2%d2%h200.d01.nc DTYPE netcdf OPTIONS template UNDEF -9999 diff --git a/lis/testcases/metforcing/gefs/testcase.ctl b/lis/testcases/metforcing/gefs/testcase.ctl old mode 100755 new mode 100644 index a9339b1ec..ce879cbd2 --- a/lis/testcases/metforcing/gefs/testcase.ctl +++ b/lis/testcases/metforcing/gefs/testcase.ctl @@ -1,4 +1,4 @@ -DSET ./TARGET_OUTPUT/SURFACEMODEL/%y4%m2/LIS_HIST_%y4%m2%d2%h200.d01.nc +DSET ^TARGET_OUTPUT/SURFACEMODEL/%y4%m2/LIS_HIST_%y4%m2%d2%h200.d01.nc DTYPE netcdf OPTIONS template UNDEF -9999