[Met_help] [rt.rap.ucar.edu #81581] History for displacement of objects in the MODE outpur PS figure

Julie Prestopnik via RT met_help at ucar.edu
Mon Sep 18 08:42:00 MDT 2017


----------------------------------------------------------------
  Initial Request
----------------------------------------------------------------

Hi,
Upon running MODE, i encounter a weird phenomenon:
Before running MODE, forcast data (in wrfout* format) and observation data (in grib format) are prepared and processed by ncl program. Using ncl program, i separately read in the aboved data and re-write them into a common nc format (referring to a netcdf file in MET_Tutorial:
i.e. sample_obs_2005080703V_03A.nc). it is worth mentioning that all the original data cover a large area, and i have to extract my interested sub-area from the original field in ncl program.

After the above mentioned processes: re-formatting, extracting, two nc files in common format are generated as input files in MODE function. The regridding process is configured in MODEConfigure_default*. No errors are found in runnig MODE, however, the contours of fcst and obs in the output PS figure are really different from the original fields plotted by ncl program (non-native grid, attached). As shown in two figures, the contour is obviously displaced. I check the lat and lon data in the output nc file of MODE, and find that they are not the same as those in the input nc file.

ps: the attached two png figures are plotted by ncl program, and the PS figure is the output of MODE.


----------------------------------------------------------------
  Complete Ticket History
----------------------------------------------------------------

Subject: displacement of objects in the MODE outpur PS figure
From: Julie Prestopnik
Time: Fri Aug 11 11:10:21 2017

I see that you are encountering unexpected behavior from MODE.
Unfortunately, some of our MET experts who are best qualified to
address
these issues are out of the office, but will be back next week.

I think it would be helpful to be able to take a look at your data
files,
both the original files and your modified files.  Could you please
post
data to our anonymous ftp site following these instructions?
   http://www.dtcenter.org/met/
<http://www.dtcenter.org/met/users/support/met_help.php#ftp>
users/support/met_help.php#ftp
<http://www.dtcenter.org/met/users/support/met_help.php#ftp>

Please let me know once you have posted the data, so that I can grab
it
from the ftp server.

Also, can you please tell us the version of MET that you are using?
Please
also tell us the command that you ran for MODE and be sure to upload
any
configuration files that you used as well.

Thanks,
Julie


On Fri, Aug 11, 2017 at 8:01 AM, 沈历都 via RT <met_help at ucar.edu> wrote:

>
> Fri Aug 11 08:01:56 2017: Request 81581 was acted upon.
> Transaction: Ticket created by shenlidu052 at 126.com
>        Queue: met_help
>      Subject: displacement of objects in the MODE outpur PS figure
>        Owner: Nobody
>   Requestors: shenlidu052 at 126.com
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81581 >
>
>
> Hi,
> Upon running MODE, i encounter a weird phenomenon:
> Before running MODE, forcast data (in wrfout* format) and
observation data
> (in grib format) are prepared and processed by ncl program. Using
ncl
> program, i separately read in the aboved data and re-write them into
a
> common nc format (referring to a netcdf file in MET_Tutorial:
> i.e. sample_obs_2005080703V_03A.nc). it is worth mentioning that all
the
> original data cover a large area, and i have to extract my
interested
> sub-area from the original field in ncl program.
>
> After the above mentioned processes: re-formatting, extracting, two
nc
> files in common format are generated as input files in MODE
function. The
> regridding process is configured in MODEConfigure_default*. No
errors are
> found in runnig MODE, however, the contours of fcst and obs in the
output
> PS figure are really different from the original fields plotted by
ncl
> program (non-native grid, attached). As shown in two figures, the
contour
> is obviously displaced. I check the lat and lon data in the output
nc file
> of MODE, and find that they are not the same as those in the input
nc file.
>
> ps: the attached two png figures are plotted by ncl program, and the
PS
> figure is the output of MODE.
>
>

------------------------------------------------
Subject: displacement of objects in the MODE outpur PS figure
From: 沈历都
Time: Mon Aug 14 06:14:01 2017

Hi, Julie
All data requested have uploaded onto the ftp server, directory named
"shen_data", including a configuration file, original fcst/obs data,
my modified NC files (named *_input.nc) and the output of MODE
function (mode_*.nc).
We use MET v6.0, and the command line is "mode  fcst_input.nc
obs_input.nc MODEConfig_shen".


Best Regard !
Lidu Shen
2017-08-14



At 2017-08-12 01:10:22, "Julie Prestopnik via RT" <met_help at ucar.edu>
wrote:
>I see that you are encountering unexpected behavior from MODE.
>Unfortunately, some of our MET experts who are best qualified to
address
>these issues are out of the office, but will be back next week.
>
>I think it would be helpful to be able to take a look at your data
files,
>both the original files and your modified files.  Could you please
post
>data to our anonymous ftp site following these instructions?
>   http://www.dtcenter.org/met/
><http://www.dtcenter.org/met/users/support/met_help.php#ftp>
>users/support/met_help.php#ftp
><http://www.dtcenter.org/met/users/support/met_help.php#ftp>
>
>Please let me know once you have posted the data, so that I can grab
it
>from the ftp server.
>
>Also, can you please tell us the version of MET that you are using?
Please
>also tell us the command that you ran for MODE and be sure to upload
any
>configuration files that you used as well.
>
>Thanks,
>Julie
>


------------------------------------------------
Subject: displacement of objects in the MODE outpur PS figure
From: Julie Prestopnik
Time: Mon Aug 14 09:30:48 2017

Hi.  Unfortunately, there seems to be a problem with the zip file that
you
sent to us.  When running "unzip data.zip", I get:

error [data.zip]:  missing 353 bytes in zipfile
  (attempting to process anyway)
error [data.zip]:  attempt to seek before beginning of zipfile
  (please check that you have transferred or created the zipfile in
the
  appropriate BINARY mode and that you have compiled UnZip properly)
  (attempting to re-compensate)

I see that you have also uploaded the following files in a
shen_data/shen_data directory:

fcst_input.nc
fcst_original_data.nc

however, we still appear to be missing obs_input.nc and
MODEConfig_shen.

Could you please upload those files to the ftp server when you have an
opportunity?

Thank you.

Julie



On Mon, Aug 14, 2017 at 6:14 AM, 沈历都 via RT <met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81581 >
>
> Hi, Julie
> All data requested have uploaded onto the ftp server, directory
named
> "shen_data", including a configuration file, original fcst/obs data,
my
> modified NC files (named *_input.nc) and the output of MODE function
> (mode_*.nc).
> We use MET v6.0, and the command line is "mode  fcst_input.nc
> obs_input.nc MODEConfig_shen".
>
>
> Best Regard !
> Lidu Shen
> 2017-08-14
>
>
>
> At 2017-08-12 01:10:22, "Julie Prestopnik via RT"
<met_help at ucar.edu>
> wrote:
> >I see that you are encountering unexpected behavior from MODE.
> >Unfortunately, some of our MET experts who are best qualified to
address
> >these issues are out of the office, but will be back next week.
> >
> >I think it would be helpful to be able to take a look at your data
files,
> >both the original files and your modified files.  Could you please
post
> >data to our anonymous ftp site following these instructions?
> >   http://www.dtcenter.org/met/
> ><http://www.dtcenter.org/met/users/support/met_help.php#ftp>
> >users/support/met_help.php#ftp
> ><http://www.dtcenter.org/met/users/support/met_help.php#ftp>
> >
> >Please let me know once you have posted the data, so that I can
grab it
> >from the ftp server.
> >
> >Also, can you please tell us the version of MET that you are using?
> Please
> >also tell us the command that you ran for MODE and be sure to
upload any
> >configuration files that you used as well.
> >
> >Thanks,
> >Julie
> >
>
>
>

------------------------------------------------
Subject: displacement of objects in the MODE outpur PS figure
From: 沈历都
Time: Mon Aug 14 20:32:17 2017

Hi, Julie

I have packed and compressed all the data again in *.tar.gz format,
and uploaded (in BINARY mode) it onto the ftp server in "shen_data/"
directory. The other separated files uploaded last day are already
deleted, please take a check and be subject to the version of today's
uploading. Sorry for the inconvinience.

Nice day!

Lidu Shen
2017-08-15

------------------------------------------------
Subject: displacement of objects in the MODE outpur PS figure
From: Julie Prestopnik
Time: Wed Aug 16 14:49:45 2017

Hello.  My apologies for the delay in responding.  I was out of the
office
yesterday.  Thank you for uploading your data. I received all of the
files
with no problems this time.


To take a look at the raw observational wind speed, we used
plot_data_plane:

/usr/local/met-6.0/bin/plot_data_plane \
obs_original_data.GR2 obs_original_data_WIND_Z10.ps 'name="WIND";
level="Z10";' -plot_range 0 15.6


To take a look at the NCL regridded wind speed, we also used
plot_data_plane:

/usr/local/met-6.0/bin/plot_data_plane \
obs_input.nc obs_input_WIND_Z10.ps 'name="WSPD"; level="(*,*)";'
-plot_range 0 15.6


Then, we regridded the NCL wind speed back to the original doman using
regrid_data_plane:

/usr/local/met-6.0/bin/regrid_data_plane \
obs_input.nc obs_original_data.GR2 obs_input_regridded_back_to_or
iginal_domain.nc \
-field 'name="WSPD"; level="(*,*)";' -method NEAREST

and then we plotted this data using plot_data_plane:

/usr/local/met-6.0/bin/plot_data_plane \
obs_input_regridded_back_to_original_domain.nc
obs_input_regridded_back_to_
original_domain.ps 'name="WSPD_10"; level="(*,*)";' -plot_range 0 15.6


We compared the original observation wind speed with the regridded NCL
wind
speed that we regridded back to its original domain and saw an offset
in
the contours.  There appears to be a problem in the regridding logic
used
to regrid the wind speed.

We see with the original forecast data that U10 and V10 are present,
so it
appears that you first used NCL to calculate wind speed and then
regrid
it.

You can proceed by taking a look at the regridding logic and fixing
the
problem with it, or you can let mode do the regridding.  In the MODE
configuration file, there is a "Verification grid" section and
"to_grid"
option under "regrid".  The "regrid" entry specifies regridding logic.
The
"to_grid" entry is described in detail starting on page 49 of our MET
User's Guide for METv6.0:

http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v6.0.pdf

You can let MODE regrid both your forecast and observation file if
you'd
like.  In that case, you'd only need to calculate the wind speed using
NCL
and then give the newly created NetCDF file and the original GRIB2
observation file to MODE and let MODE do the regridding on both data
sets.

I have attached png images that I created from the ps images that were
created using the above tools.  I used "convert" to convert the ps
files to
png files.  For example:

convert -rotate 90 -background white -flatten \
obs_input_regridded_back_to_original_domain.ps \
obs_input_regridded_back_to_original_domain.png

I have attached the png files that we created.

Please let us know how you decide to proceed and how it works out for
you.
I hope this helps!

Kind regards,
Julie

On Mon, Aug 14, 2017 at 8:32 PM, 沈历都 via RT <met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81581 >
>
> Hi, Julie
>
> I have packed and compressed all the data again in *.tar.gz format,
and
> uploaded (in BINARY mode) it onto the ftp server in "shen_data/"
directory.
> The other separated files uploaded last day are already deleted,
please
> take a check and be subject to the version of today's uploading.
Sorry for
> the inconvinience.
>
> Nice day!
>
> Lidu Shen
> 2017-08-15
>

------------------------------------------------


More information about the Met_help mailing list