[Met_help] update on point_stat soil verification

John Halley Gotway johnhg at rap.ucar.edu
Thu Aug 27 10:04:54 MDT 2009


Jon,

Great.  I'm glad you were able to get it working.

I'm in a workshop this week and have limited access to email, but feel
free to write with any more questions that come up, and I'll do my best to
answer.

Thanks,
John

> John,
>
> For the time being, I commented out the following sections of the file
> src/point_stat/point_stat_conf_info.cc:
>
> ~line 178:
>       // Only support vertical levels at 0m, 2m, or 10m
> // J.Case (8/26/2009) -- comment out to test soil verification for now.
> //      if(gci.lvl_type == VertLevel &&
> //         gci.lvl_1 != 0 && gci.lvl_1 != 2 && gci.lvl_1 != 10) {
> //
> //         cerr << "\n\nERROR: process_config() -> "
> //              << "Only vertical level values of 2 and 10 meters "
> //              << "are currently supported.\n\n"
> //              << flush;
> //         exit(1);
> //      }
>       gc_pd[i].set_fcst_gci(gci);
>    }
>
> And ~line 230:
>       // Only support vertical levels at 0m, 2m, or 10m
> // J.Case (8/26/2009) -- comment out to test soil verification for now.
> //      if(gci.lvl_type == VertLevel &&
> //         gci.lvl_1 != 0 && gci.lvl_1 != 2 && gci.lvl_1 != 10) {
> //
> //         cerr << "\n\nERROR: PointStatConfInfo::process_config() -> "
> //              << "Only vertical level values of 2 and 10 meters "
> //              << "are currently supported.\n\n"
> //              << flush;
> //         exit(1);
> //      }
>       gc_pd[i].set_obs_gci(gci);
>
>
> The point_stat program now appears to work on the soil moisture, soil
> temperature, 2-m temperature, and 2-m dewpoint data from the various SCAN
> sites.  I looked at the _mpr.txt file, and it depicts pairs of data as you
> advertised, matching each obs soil level with the fcst soil layer is lies
> in, including duplicates (i.e. observed 10 cm soil readings matches with
> both 0-10 cm and 10-40 cm fcst soil layers).
>
> At some point, we can clean up these "if" blocks to make it more generic.
> In the meantime, just commenting it out is sufficient for me.
>
> Cheers,
> Jonathan
>
> --------------------------------------------------------------------------------------------------
> Jonathan Case, ENSCO Inc.
> NASA Short-term Prediction Research and Transition Center (aka SPoRT
> Center)
> 320 Sparkman Drive, Room 3062
> Huntsville, AL 35805
> Voice: 256.961.7504
> Fax: 256.961.7788
> Emails: Jonathan.Case-1 at nasa.gov / case.jonathan at ensco.com
> --------------------------------------------------------------------------------------------------
>
> "Whether the weather is cold, or whether the weather is hot, we'll weather
>   the weather whether we like it or not!"
>
>
>




More information about the Met_help mailing list