有用过read_wrf_nc.f的请帮我看一看，谢谢

sunshinezhpb

想用read_wrf_nc.f改下垫面，修改之后改LU_INDEX成功了，但是对于LANDUSEF这个四维变量怎么都没改成，还请大家帮忙看一下错在何处：
定义real,    dimension(dim1,dim2,20,dim3)        ::  data_real1 (我用的modis数据只分了20类）
elseif ( var == 'LU_INDEX')then
    do i = 1, 88
     do j = 1, 88
      if(data_real(i,j,1) == 12)then
           data_real(i,j,1)=13
        endif
     enddo
    enddo

  elseif ( var == 'LANDUSEF'.and.var1 == 'LU_INDEX') then
    do k =1, 20
     do i = 1, 88
      do j = 1, 88
          if (data_real(i,j,1) == 12)then
           data_real1(i,j,12,1) = 0
           data_real1(i,j,13,1) = 1
          endif
        enddo
     enddo
    enddo
改LU_INDEX成功，改LANDUSEF则提示 Variable given was not one of above - so no change will be made to any variables.然后我又改了下，
elseif ( var == 'LANDUSEF') then
    do k =1, 20
     do i = 1, 88
      do j = 1, 88
           data_real1(i,j,12,1) = 0
           data_real1(i,j,13,1) = 1
        enddo
     enddo
    enddo
这次fortran提示我已经改动,但我用ncdump提出LANDUSEF仍然没变，求解LANDUSEF这种四维的怎么改？

月是故乡明

楼主好，提示是否修改时，我选择yes，为什么就屏幕就一直不停刷出“y”，然后就没结果了呢

damien1987

楼主是怎么改的？用ncdump提取出LU_INDEX.txt然后呢？手动修改还是写代码实现的？

紫罗兰2013

楼主该下垫面时对哪个文件进行操作？