fixing issues with ERA5 climate

[tslin2]

#56
Namelist Vtable, z changed to TERRAIN
call latlon_to_ij() don’t need to +360.0 for ERA5 data, if longitude is positive.

[cenlinhe]

This fix to deal with ERA-5 looks good to me and I approve it.

[CharlesZheZhang]

Thanks Cenlin and Tzu-Shun. As we discussed this afternoon, I agree to merge this pull request to develop branch.

[tslin2]

if geo longitude is negative, the longitude of precipitation needs to add 360 degree in the subroutine interp_rainfall_nearest_neighbor

hrldas/hrldas/HRLDAS_forcing/create_forcing.F

Line 2734 in 9df967d

call latlon_to_ij(datastruct%proj, geo_em%lat(ii,jj), geo_em%lon(ii,jj), x, y)

revise to

    if (forcing_type=="ERA5") then
       if (geo_em%lon(ii,jj) < 0) then
          call latlon_to_ij(datastruct%proj, geo_em%lat(ii,jj), geo_em%lon(ii,jj)+360.0, x, y)
       else
          call latlon_to_ij(datastruct%proj, geo_em%lat(ii,jj), geo_em%lon(ii,jj), x, y)
       endif
    else
       call latlon_to_ij(datastruct%proj, geo_em%lat(ii,jj), geo_em%lon(ii,jj), x, y)
    endif

[tslin2]

Another issue with the soil thickness vertical interpolation
The problem is due to the var(:,:,k) in the following line

hrldas/hrldas/IO_code/module_hrldas_netcdf_io.F

Line 1909 in 9df967d

var(:,:,k) = (src(:,ktop,:)*fraction) + (src(:,kbottom,:)*(1.0-fraction))

changing to var(:,k,:) works fine.