[Met_help] [rt.rap.ucar.edu #85515] History for Diagnostic files diag_conv_ges* and diag_conv_anl*

[John Halley Gotway via RT]

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

Dear MET help desk,

I used MET software to read the diag_conv files output by hybrid 3DEnVAR (
GSI version 3.6).
One output is the ANLY_USE variable:

*what are the differences between the ANLY_USE (1, or -1) for two
diagnostic files?*
*Which ANLY_USE represents the assimilated observations, the diag_conv_ges*
or drag_conv_anl*?*

 I found that the number of rejected observations in *ges files are larger
than that  in *anl files, which is a bit surprising to me.

Thanks for your time in advance.

Regards,
Minghua


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

Subject: Diagnostic files diag_conv_ges* and diag_conv_anl*
From: John Halley Gotway
Time: Wed Jun 13 13:23:33 2018

I see you have a question about the contents of the GSI diagnostic
files.
I asked my colleague about this and have included her response below.

Hope this helps clarify.

Thanks,
John Halley Gotway

The ANLY_USE refers to whether the observation was assimilated or not.
1=
assimilated, -1=not assimilated. The actual assimilated observations
in the
analysis are marked as 1 in the diag_conv_anl file. The reason the
diag_conv_ges (background) has more rejected observations is that this
is
before the first outer loop, where GSI throws out the most
observations.
Observations are brought back in during subsequent outer loops and a
subset
will be rejected again. Therefore, o-a (which is written to
diag_conv_anl)
will have fewer observations flagged as -1 than the o-g
(diag_conv_ges).

If they have more questions on how GSI operates and writes the diag
files,
we can push them over to gsi-help at ucar.edu.

Let me know if you need further clarification.

Thanks,
Kathryn

On Tue, Jun 12, 2018 at 6:11 PM Minghua Zheng via RT
<met_help at ucar.edu>
wrote:

>
> Tue Jun 12 18:11:21 2018: Request 85515 was acted upon.
> Transaction: Ticket created by mzheng at ucsd.edu
>        Queue: met_help
>      Subject: Diagnostic files diag_conv_ges* and diag_conv_anl*
>        Owner: Nobody
>   Requestors: mzheng at ucsd.edu
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=85515 >
>
>
> Dear MET help desk,
>
> I used MET software to read the diag_conv files output by hybrid
3DEnVAR (
> GSI version 3.6).
> One output is the ANLY_USE variable:
>
> *what are the differences between the ANLY_USE (1, or -1) for two
> diagnostic files?*
> *Which ANLY_USE represents the assimilated observations, the
diag_conv_ges*
> or drag_conv_anl*?*
>
>  I found that the number of rejected observations in *ges files are
larger
> than that  in *anl files, which is a bit surprising to me.
>
> Thanks for your time in advance.
>
> Regards,
> Minghua
>
>

------------------------------------------------
Subject: Diagnostic files diag_conv_ges* and diag_conv_anl*
From: Minghua Zheng
Time: Wed Jun 13 13:43:58 2018

Hi John,

Thank you very much for your help. That makes sense to me.
Thank Kathryn for her reply.

Best regards,
Minghua





On Wed, Jun 13, 2018 at 12:23 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:

> I see you have a question about the contents of the GSI diagnostic
files.
> I asked my colleague about this and have included her response
below.
>
> Hope this helps clarify.
>
> Thanks,
> John Halley Gotway
>
> The ANLY_USE refers to whether the observation was assimilated or
not. 1=
> assimilated, -1=not assimilated. The actual assimilated observations
in the
> analysis are marked as 1 in the diag_conv_anl file. The reason the
> diag_conv_ges (background) has more rejected observations is that
this is
> before the first outer loop, where GSI throws out the most
observations.
> Observations are brought back in during subsequent outer loops and a
subset
> will be rejected again. Therefore, o-a (which is written to
diag_conv_anl)
> will have fewer observations flagged as -1 than the o-g
(diag_conv_ges).
>
> If they have more questions on how GSI operates and writes the diag
files,
> we can push them over to gsi-help at ucar.edu.
>
> Let me know if you need further clarification.
>
> Thanks,
> Kathryn
>
> On Tue, Jun 12, 2018 at 6:11 PM Minghua Zheng via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > Tue Jun 12 18:11:21 2018: Request 85515 was acted upon.
> > Transaction: Ticket created by mzheng at ucsd.edu
> >        Queue: met_help
> >      Subject: Diagnostic files diag_conv_ges* and diag_conv_anl*
> >        Owner: Nobody
> >   Requestors: mzheng at ucsd.edu
> >       Status: new
> >  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=85515 >
> >
> >
> > Dear MET help desk,
> >
> > I used MET software to read the diag_conv files output by hybrid
3DEnVAR
> (
> > GSI version 3.6).
> > One output is the ANLY_USE variable:
> >
> > *what are the differences between the ANLY_USE (1, or -1) for two
> > diagnostic files?*
> > *Which ANLY_USE represents the assimilated observations, the
> diag_conv_ges*
> > or drag_conv_anl*?*
> >
> >  I found that the number of rejected observations in *ges files
are
> larger
> > than that  in *anl files, which is a bit surprising to me.
> >
> > Thanks for your time in advance.
> >
> > Regards,
> > Minghua
> >
> >
>
>

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