Menu

ncap2 masking

mark
2018-02-22
2018-02-22
  • mark

    mark - 2018-02-22

    Hello,
    I'd like to create/subset a variable subject to a mask. ncap2 executed as

    ncap2 -O -v -S fcst.nco diag_viirs_ges.2015081000.nc4 tmp.nc4

    %cat fcst.nco
    where(use_flag > 0) Fcst=Observation-Obs_Minus_Forecast_unadjusted;

    produces error

    ncap2: WARNING prs_cls::ncap_var_init(): Unable to find variable Fcst in diag_viirs_ges.2015081000.nc4 or tmp.nc4
    ncap2: ERROR out(): Could not read var Fcst

    How that should be executed instead? Do I need to assign Fcst dimensions first?

    Thanks,
    Mark

    netcdf diag_viirs_ges.2015081000 {
    dimensions:
    nchans = 11 ;
    nobs = UNLIMITED ; // (63811 currently)
    ...
    int use_flag(nchans) ;
    float Observation(nobs) ;
    float Obs_Minus_Forecast_unadjusted(nobs) ;

     
    • Charlie Zender

      Charlie Zender - 2018-02-22

      Please carefully read the the ncap2 where documentation.
      Fcst must already exist. And the mask, use_flag, should share some dimensions with the LHS and RHS variables.

       

Log in to post a comment.