关于CLM type的四级数据结构

刘振

大家好，请问有谁知道clm的四级数据结构中的坐标转换问题吗？在lnd/src/clm4_0/biogeophys目录下的HydrologyMod.F90中，我找到了土壤单位体积含水率的公式：
    do j = 1, nlevgrnd
       do fc = 1, num_nolakec
          c = filter_nolakec(fc)
          if ((ctype(c) == icol_sunwall .or. ctype(c) == icol_shadewall &
               .or. ctype(c) == icol_roof) .and. j > nlevurb) then
          else
            endwb(c) = endwb(c) + h2osoi_ice(c,j) + h2osoi_liq(c,j)
            h2osoi_vol(c,j) = h2osoi_liq(c,j)/(dz(c,j)*denh2o) + h2osoi_ice(c,j)/(dz(c,j)*denice)
          end if
       end do
    end do
integer, intent(in) :: num_nolakec                 ! number of column non-lake points in column filter
integer, intent(in) :: filter_nolakec(ubc-lbc+1)   ! column filter for non-lake points

其中这两个变量的含义我不是很理解。
还有个问题就是在最后模式的结果中h2osoi都是经纬度的设置的。从h2osoi_vol(c,j)到h2osoi_vol（lon,lat,nlevgrnd）是利用什么函数转化的吗？谢谢