[已解决] NCL如何画全球网格数据的陆地或者海洋平均？

LeslieTse

齐萃妮 发表于 2023-6-12 13:14
楼主您好，我也遇到了相似的问题，想要计算陆地的平均，请问您是怎么解决的呢

  贴一段代码吧，供你参考。需要同样分辨率的landfraction变量来把海洋区域的格点设为缺测。
  yyyymm1      = cd_calendar(tim1, -1)
  ymStrt1      = 185101
  ymLast1      = 201412
  iStrt1       = ind(yyyymm1.eq.ymStrt1)
  iLast1       = ind(yyyymm1.eq.ymLast1)
  T_1          = b1->TREFHT(iStrt1:iLast1,:,:)
  T_1          = T_1 - 273.15
  lf_1         = c1->LANDFRAC(iStrt1:iLast1,:,:)
  T_1_l        = mask(T_1, lf_1, 1)
  T_1_o        = mask(T_1, lf_1, 0)
  copy_VarMeta(T_1, T_1_l)
  copy_VarMeta(T_1, T_1_o)
  T_1_y        = month_to_annual(T_1, 1)
  T_1_ly       = month_to_annual(T_1_l, 1)
  T_1_oy       = month_to_annual(T_1_o, 1)
  T_1_y_ave    = wgt_areaave_Wrap(T_1_y, clat1, 1.0, 1)
  T_1_ly_ave   = wgt_areaave_Wrap(T_1_ly, cos(T_1_ly&lat*rad1), 1.0, 1)
  T_1_oy_ave   = wgt_areaave_Wrap(T_1_oy, cos(T_1_oy&lat*rad1), 1.0, 1)