[Met_help] [rt.rap.ucar.edu #94517] History for AERONET AOD processing
Howard Soh via RT
met_help at ucar.edu
Tue Mar 17 12:26:05 MDT 2020
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello,
Is a formatting code available for version 3 AERONET? (https://aeronet.gsfc.nasa.gov/new_web/data.html <https://aeronet.gsfc.nasa.gov/new_web/data.html>)
(If yes) Does the code allow for selecting chosen wavelengths? Is capable of obtaining selected channels (e.g. 550 nm) from other wavelengths?
Thanks,
Mariusz
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: AERONET AOD processing
From: Minna Win
Time: Thu Mar 12 13:57:48 2020
Hello Mariusz,
I checked the MET v8.1 release notes and for ASCII2NC, support to the
"-format" command line was added to handle both aeronetv2 and
aeronetv3.
Is this what you were hoping for?
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
wrote:
>
> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
> Queue: met_help
> Subject: AERONET AOD processing
> Owner: Nobody
> Requestors: Mariusz.Pagowski at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>
>
> Hello,
> Is a formatting code available for version 3 AERONET? (
> https://aeronet.gsfc.nasa.gov/new_web/data.html <
> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
> (If yes) Does the code allow for selecting chosen wavelengths? Is
capable
> of obtaining selected channels (e.g. 550 nm) from other wavelengths?
> Thanks,
> Mariusz
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #94517] AERONET AOD processing
From: Mariusz Pagowski
Time: Thu Mar 12 14:01:41 2020
Yes, and also the info what wavelengths it extracts and if it
interpolates.
Thanks,
Mariusz
> On Mar 12, 2020, at 1:57 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>
> Hello Mariusz,
>
> I checked the MET v8.1 release notes and for ASCII2NC, support to
the
> "-format" command line was added to handle both aeronetv2 and
aeronetv3.
> Is this what you were hoping for?
>
>
> Regards,
> Minna
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
>
>
>
> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
> wrote:
>
>>
>> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
>> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
>> Queue: met_help
>> Subject: AERONET AOD processing
>> Owner: Nobody
>> Requestors: Mariusz.Pagowski at noaa.gov
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>
>>
>> Hello,
>> Is a formatting code available for version 3 AERONET? (
>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
>> (If yes) Does the code allow for selecting chosen wavelengths? Is
capable
>> of obtaining selected channels (e.g. 550 nm) from other
wavelengths?
>> Thanks,
>> Mariusz
>>
>>
>>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Howard Soh
Time: Thu Mar 12 14:02:43 2020
Internal only
ascii2nc reads a text format AERONET 3 data and produces a NetCDF
output (MET's point observation data file). ascii2nc interpolates
550nm wavelength from 440 nm and 675nm if not exists.
Cheers,
Howard
On Thu Mar 12 13:27:57 2020, Mariusz.Pagowski at noaa.gov wrote:
> Hello,
> Is a formatting code available for version 3 AERONET?
> (https://aeronet.gsfc.nasa.gov/new_web/data.html
> <https://aeronet.gsfc.nasa.gov/new_web/data.html>)
> (If yes) Does the code allow for selecting chosen wavelengths? Is
> capable of obtaining selected channels (e.g. 550 nm) from other
> wavelengths?
> Thanks,
> Mariusz
------------------------------------------------
Subject: AERONET AOD processing
From: Minna Win
Time: Thu Mar 12 14:06:56 2020
Hi Mariusz,
According to our resident expert Howard Soh:
"ascii2nc reads a text format AERONET 3 data and produces a NetCDF
output
(MET's point observation data file). ascii2nc interpolates 550nm
wavelength
from 440 nm and 675nm if not exists."
I hope this helps.
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
wrote:
>
> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
> Queue: met_help
> Subject: AERONET AOD processing
> Owner: Nobody
> Requestors: Mariusz.Pagowski at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>
>
> Hello,
> Is a formatting code available for version 3 AERONET? (
> https://aeronet.gsfc.nasa.gov/new_web/data.html <
> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
> (If yes) Does the code allow for selecting chosen wavelengths? Is
capable
> of obtaining selected channels (e.g. 550 nm) from other wavelengths?
> Thanks,
> Mariusz
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #94517] AERONET AOD processing
From: Mariusz Pagowski
Time: Thu Mar 12 14:10:44 2020
Minna,
thanks very much - just one follow-up: how is the interpolation done
(linear/log)?
Mariusz
> On Mar 12, 2020, at 2:06 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>
> Hi Mariusz,
>
> According to our resident expert Howard Soh:
>
> "ascii2nc reads a text format AERONET 3 data and produces a NetCDF
output
> (MET's point observation data file). ascii2nc interpolates 550nm
wavelength
> from 440 nm and 675nm if not exists."
>
> I hope this helps.
>
> Regards,
> Minna
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
>
>
>
> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
> wrote:
>
>>
>> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
>> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
>> Queue: met_help
>> Subject: AERONET AOD processing
>> Owner: Nobody
>> Requestors: Mariusz.Pagowski at noaa.gov
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>
>>
>> Hello,
>> Is a formatting code available for version 3 AERONET? (
>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
>> (If yes) Does the code allow for selecting chosen wavelengths? Is
capable
>> of obtaining selected channels (e.g. 550 nm) from other
wavelengths?
>> Thanks,
>> Mariusz
>>
>>
>>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Minna Win
Time: Thu Mar 12 14:37:28 2020
Hi Mariusz,
According to Howard:
It's a log interpolation: angstrom_power_interpolation
double angstrom_power_interplation(double value_1, double value_2,
double level_1, double level_2, double target_level) {
double angstrom_log =
-log10(value_1/value_2)/log10(level_1/level_2);
double angstrom_value = value_2 *
pow((target_level/level_2),-angstrom_log);
return angstrom_value;
}
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>
> Minna,
> thanks very much - just one follow-up: how is the interpolation done
> (linear/log)?
> Mariusz
>
>
> > On Mar 12, 2020, at 2:06 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
> >
> > Hi Mariusz,
> >
> > According to our resident expert Howard Soh:
> >
> > "ascii2nc reads a text format AERONET 3 data and produces a NetCDF
output
> > (MET's point observation data file). ascii2nc interpolates 550nm
> wavelength
> > from 440 nm and 675nm if not exists."
> >
> > I hope this helps.
> >
> > Regards,
> > Minna
> > ---------------
> > Minna Win
> > National Center for Atmospheric Research
> > Developmental Testbed Center
> > Phone: 303-497-8423
> > Fax: 303-497-8401
> >
> >
> >
> > On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT <
> met_help at ucar.edu>
> > wrote:
> >
> >>
> >> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
> >> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
> >> Queue: met_help
> >> Subject: AERONET AOD processing
> >> Owner: Nobody
> >> Requestors: Mariusz.Pagowski at noaa.gov
> >> Status: new
> >> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>
> >>
> >> Hello,
> >> Is a formatting code available for version 3 AERONET? (
> >> https://aeronet.gsfc.nasa.gov/new_web/data.html <
> >> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
> >> (If yes) Does the code allow for selecting chosen wavelengths? Is
> capable
> >> of obtaining selected channels (e.g. 550 nm) from other
wavelengths?
> >> Thanks,
> >> Mariusz
> >>
> >>
> >>
> >
>
>
>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Minna Win
Time: Thu Mar 12 15:56:02 2020
Hi Mariusz,
Upon further checking from Howard, the implementation is actually:
The Fortran code I referenced
(https://github.com/NCAR/MET/issues/1121) was
> implemented that way:
> a_675_440 = -log(aerdata(i)%aot_675nm/ &
> aerdata(i)%aot_440nm)/log(675./440.)
> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* & ((550./440.)**(-
a_675_440))
Just in case you need to see the correct implementation code. The
first code snippet (in previous email) he cited was replaced with the
above code snippet, which is associated with Github issue 1121...
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister
<minnawin at ucar.edu>
wrote:
> Hi Mariusz,
>
> According to Howard:
>
> It's a log interpolation: angstrom_power_interpolation
>
> double angstrom_power_interplation(double value_1, double value_2,
> double level_1, double level_2, double target_level) {
> double angstrom_log =
-log10(value_1/value_2)/log10(level_1/level_2);
> double angstrom_value = value_2 *
> pow((target_level/level_2),-angstrom_log);
> return angstrom_value;
> }
>
> Regards,
> Minna
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
>
>
>
> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>
>> Minna,
>> thanks very much - just one follow-up: how is the interpolation
done
>> (linear/log)?
>> Mariusz
>>
>>
>> > On Mar 12, 2020, at 2:06 PM, Minna Win via RT <met_help at ucar.edu>
>> wrote:
>> >
>> > Hi Mariusz,
>> >
>> > According to our resident expert Howard Soh:
>> >
>> > "ascii2nc reads a text format AERONET 3 data and produces a
NetCDF
>> output
>> > (MET's point observation data file). ascii2nc interpolates 550nm
>> wavelength
>> > from 440 nm and 675nm if not exists."
>> >
>> > I hope this helps.
>> >
>> > Regards,
>> > Minna
>> > ---------------
>> > Minna Win
>> > National Center for Atmospheric Research
>> > Developmental Testbed Center
>> > Phone: 303-497-8423
>> > Fax: 303-497-8401
>> >
>> >
>> >
>> > On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT <
>> met_help at ucar.edu>
>> > wrote:
>> >
>> >>
>> >> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
>> >> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
>> >> Queue: met_help
>> >> Subject: AERONET AOD processing
>> >> Owner: Nobody
>> >> Requestors: Mariusz.Pagowski at noaa.gov
>> >> Status: new
>> >> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>> >>
>> >>
>> >> Hello,
>> >> Is a formatting code available for version 3 AERONET? (
>> >> https://aeronet.gsfc.nasa.gov/new_web/data.html <
>> >> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
>> >> (If yes) Does the code allow for selecting chosen wavelengths?
Is
>> capable
>> >> of obtaining selected channels (e.g. 550 nm) from other
wavelengths?
>> >> Thanks,
>> >> Mariusz
>> >>
>> >>
>> >>
>> >
>>
>>
>>
>>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #94517] AERONET AOD processing
From: Mariusz Pagowski
Time: Thu Mar 12 16:02:37 2020
Thanks Mina,
just wanted to be sure that it is not a linear interpolation. Still I
wonder what happens when 675 and 440 are not available.
So if my understanding is correct MET only outputs 550nm.
Mariusz
> On Mar 12, 2020, at 3:56 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>
> Hi Mariusz,
>
> Upon further checking from Howard, the implementation is actually:
>
>
> The Fortran code I referenced
(https://github.com/NCAR/MET/issues/1121) was
>> implemented that way:
>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
>> aerdata(i)%aot_440nm)/log(675./440.)
>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* & ((550./440.)**(-
a_675_440))
>
>
> Just in case you need to see the correct implementation code. The
> first code snippet (in previous email) he cited was replaced with
the
> above code snippet, which is associated with Github issue 1121...
>
>
> Regards,
>
> Minna
>
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
>
>
>
> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister
<minnawin at ucar.edu>
> wrote:
>
>> Hi Mariusz,
>>
>> According to Howard:
>>
>> It's a log interpolation: angstrom_power_interpolation
>>
>> double angstrom_power_interplation(double value_1, double value_2,
>> double level_1, double level_2, double target_level) {
>> double angstrom_log =
-log10(value_1/value_2)/log10(level_1/level_2);
>> double angstrom_value = value_2 *
>> pow((target_level/level_2),-angstrom_log);
>> return angstrom_value;
>> }
>>
>> Regards,
>> Minna
>> ---------------
>> Minna Win
>> National Center for Atmospheric Research
>> Developmental Testbed Center
>> Phone: 303-497-8423
>> Fax: 303-497-8401
>>
>>
>>
>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
>> wrote:
>>
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>
>>> Minna,
>>> thanks very much - just one follow-up: how is the interpolation
done
>>> (linear/log)?
>>> Mariusz
>>>
>>>
>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT <met_help at ucar.edu>
>>> wrote:
>>>>
>>>> Hi Mariusz,
>>>>
>>>> According to our resident expert Howard Soh:
>>>>
>>>> "ascii2nc reads a text format AERONET 3 data and produces a
NetCDF
>>> output
>>>> (MET's point observation data file). ascii2nc interpolates 550nm
>>> wavelength
>>>> from 440 nm and 675nm if not exists."
>>>>
>>>> I hope this helps.
>>>>
>>>> Regards,
>>>> Minna
>>>> ---------------
>>>> Minna Win
>>>> National Center for Atmospheric Research
>>>> Developmental Testbed Center
>>>> Phone: 303-497-8423
>>>> Fax: 303-497-8401
>>>>
>>>>
>>>>
>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT <
>>> met_help at ucar.edu>
>>>> wrote:
>>>>
>>>>>
>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
>>>>> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
>>>>> Queue: met_help
>>>>> Subject: AERONET AOD processing
>>>>> Owner: Nobody
>>>>> Requestors: Mariusz.Pagowski at noaa.gov
>>>>> Status: new
>>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>
>>>>>
>>>>> Hello,
>>>>> Is a formatting code available for version 3 AERONET? (
>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
>>>>> (If yes) Does the code allow for selecting chosen wavelengths?
Is
>>> capable
>>>>> of obtaining selected channels (e.g. 550 nm) from other
wavelengths?
>>>>> Thanks,
>>>>> Mariusz
>>>>>
>>>>>
>>>>>
>>>>
>>>
>>>
>>>
>>>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Minna Win
Time: Fri Mar 13 08:03:24 2020
Hi Mariusz,
Here's what happens when 675 and 440 nm are unavailable:
ascii2nc does not save 550 nm wavelength into NetCDF output if 675 and
440
> are not available.
> ascii2nc figures out the variable list and associated wavelengths
from the
> header record and knows where AOD at 440 and 650 nm values are.
> ascii2nc saves AOD of all available
> wavelengths. And ascii2nc checks AOD at 440, 550 and
> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550 does not
exist
> from data record and AOD values at 675 and at 440 are valid.
I hope this is useful.
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>
> Thanks Mina,
> just wanted to be sure that it is not a linear interpolation. Still
I
> wonder what happens when 675 and 440 are not available.
> So if my understanding is correct MET only outputs 550nm.
> Mariusz
>
>
>
> > On Mar 12, 2020, at 3:56 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
> >
> > Hi Mariusz,
> >
> > Upon further checking from Howard, the implementation is actually:
> >
> >
> > The Fortran code I referenced
(https://github.com/NCAR/MET/issues/1121)
> was
> >> implemented that way:
> >> a_675_440 = -log(aerdata(i)%aot_675nm/ &
> >> aerdata(i)%aot_440nm)/log(675./440.)
> >> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
> ((550./440.)**(-a_675_440))
> >
> >
> > Just in case you need to see the correct implementation code. The
> > first code snippet (in previous email) he cited was replaced with
the
> > above code snippet, which is associated with Github issue 1121...
> >
> >
> > Regards,
> >
> > Minna
> >
> > ---------------
> > Minna Win
> > National Center for Atmospheric Research
> > Developmental Testbed Center
> > Phone: 303-497-8423
> > Fax: 303-497-8401
> >
> >
> >
> > On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
> minnawin at ucar.edu>
> > wrote:
> >
> >> Hi Mariusz,
> >>
> >> According to Howard:
> >>
> >> It's a log interpolation: angstrom_power_interpolation
> >>
> >> double angstrom_power_interplation(double value_1, double
value_2,
> >> double level_1, double level_2, double target_level) {
> >> double angstrom_log =
-log10(value_1/value_2)/log10(level_1/level_2);
> >> double angstrom_value = value_2 *
> >> pow((target_level/level_2),-angstrom_log);
> >> return angstrom_value;
> >> }
> >>
> >> Regards,
> >> Minna
> >> ---------------
> >> Minna Win
> >> National Center for Atmospheric Research
> >> Developmental Testbed Center
> >> Phone: 303-497-8423
> >> Fax: 303-497-8401
> >>
> >>
> >>
> >> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT <
> met_help at ucar.edu>
> >> wrote:
> >>
> >>>
> >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>>
> >>> Minna,
> >>> thanks very much - just one follow-up: how is the interpolation
done
> >>> (linear/log)?
> >>> Mariusz
> >>>
> >>>
> >>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT
<met_help at ucar.edu>
> >>> wrote:
> >>>>
> >>>> Hi Mariusz,
> >>>>
> >>>> According to our resident expert Howard Soh:
> >>>>
> >>>> "ascii2nc reads a text format AERONET 3 data and produces a
NetCDF
> >>> output
> >>>> (MET's point observation data file). ascii2nc interpolates
550nm
> >>> wavelength
> >>>> from 440 nm and 675nm if not exists."
> >>>>
> >>>> I hope this helps.
> >>>>
> >>>> Regards,
> >>>> Minna
> >>>> ---------------
> >>>> Minna Win
> >>>> National Center for Atmospheric Research
> >>>> Developmental Testbed Center
> >>>> Phone: 303-497-8423
> >>>> Fax: 303-497-8401
> >>>>
> >>>>
> >>>>
> >>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT <
> >>> met_help at ucar.edu>
> >>>> wrote:
> >>>>
> >>>>>
> >>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
> >>>>> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
> >>>>> Queue: met_help
> >>>>> Subject: AERONET AOD processing
> >>>>> Owner: Nobody
> >>>>> Requestors: Mariusz.Pagowski at noaa.gov
> >>>>> Status: new
> >>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
> >
> >>>>>
> >>>>>
> >>>>> Hello,
> >>>>> Is a formatting code available for version 3 AERONET? (
> >>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
> >>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
> >>>>> (If yes) Does the code allow for selecting chosen wavelengths?
Is
> >>> capable
> >>>>> of obtaining selected channels (e.g. 550 nm) from other
wavelengths?
> >>>>> Thanks,
> >>>>> Mariusz
> >>>>>
> >>>>>
> >>>>>
> >>>>
> >>>
> >>>
> >>>
> >>>
> >
>
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #94517] AERONET AOD processing
From: Mariusz Pagowski
Time: Fri Mar 13 11:22:57 2020
Minna,
thanks for your dedication to answering these questions.
We are also interested in outputting 440 and 675 nm from AERONET, and
even possibly 870 nm.
We would be ready to work with someone from MET to help implement
these as options.
Pls let me know if that is possible,
Mariusz
> On Mar 13, 2020, at 8:03 AM, Minna Win via RT <met_help at ucar.edu>
wrote:
>
> Hi Mariusz,
>
> Here's what happens when 675 and 440 nm are unavailable:
>
> ascii2nc does not save 550 nm wavelength into NetCDF output if 675
and 440
>> are not available.
>> ascii2nc figures out the variable list and associated wavelengths
from the
>> header record and knows where AOD at 440 and 650 nm values are.
>> ascii2nc saves AOD of all available
>> wavelengths. And ascii2nc checks AOD at 440, 550 and
>> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550 does not
exist
>> from data record and AOD values at 675 and at 440 are valid.
>
>
> I hope this is useful.
>
> Regards,
> Minna
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
>
>
>
> On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>
>> Thanks Mina,
>> just wanted to be sure that it is not a linear interpolation. Still
I
>> wonder what happens when 675 and 440 are not available.
>> So if my understanding is correct MET only outputs 550nm.
>> Mariusz
>>
>>
>>
>>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>>>
>>> Hi Mariusz,
>>>
>>> Upon further checking from Howard, the implementation is actually:
>>>
>>>
>>> The Fortran code I referenced
(https://github.com/NCAR/MET/issues/1121)
>> was
>>>> implemented that way:
>>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
>>>> aerdata(i)%aot_440nm)/log(675./440.)
>>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
>> ((550./440.)**(-a_675_440))
>>>
>>>
>>> Just in case you need to see the correct implementation code. The
>>> first code snippet (in previous email) he cited was replaced with
the
>>> above code snippet, which is associated with Github issue 1121...
>>>
>>>
>>> Regards,
>>>
>>> Minna
>>>
>>> ---------------
>>> Minna Win
>>> National Center for Atmospheric Research
>>> Developmental Testbed Center
>>> Phone: 303-497-8423
>>> Fax: 303-497-8401
>>>
>>>
>>>
>>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
>> minnawin at ucar.edu>
>>> wrote:
>>>
>>>> Hi Mariusz,
>>>>
>>>> According to Howard:
>>>>
>>>> It's a log interpolation: angstrom_power_interpolation
>>>>
>>>> double angstrom_power_interplation(double value_1, double
value_2,
>>>> double level_1, double level_2, double target_level) {
>>>> double angstrom_log =
-log10(value_1/value_2)/log10(level_1/level_2);
>>>> double angstrom_value = value_2 *
>>>> pow((target_level/level_2),-angstrom_log);
>>>> return angstrom_value;
>>>> }
>>>>
>>>> Regards,
>>>> Minna
>>>> ---------------
>>>> Minna Win
>>>> National Center for Atmospheric Research
>>>> Developmental Testbed Center
>>>> Phone: 303-497-8423
>>>> Fax: 303-497-8401
>>>>
>>>>
>>>>
>>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT <
>> met_help at ucar.edu>
>>>> wrote:
>>>>
>>>>>
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>
>>>>> Minna,
>>>>> thanks very much - just one follow-up: how is the interpolation
done
>>>>> (linear/log)?
>>>>> Mariusz
>>>>>
>>>>>
>>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT
<met_help at ucar.edu>
>>>>> wrote:
>>>>>>
>>>>>> Hi Mariusz,
>>>>>>
>>>>>> According to our resident expert Howard Soh:
>>>>>>
>>>>>> "ascii2nc reads a text format AERONET 3 data and produces a
NetCDF
>>>>> output
>>>>>> (MET's point observation data file). ascii2nc interpolates
550nm
>>>>> wavelength
>>>>>> from 440 nm and 675nm if not exists."
>>>>>>
>>>>>> I hope this helps.
>>>>>>
>>>>>> Regards,
>>>>>> Minna
>>>>>> ---------------
>>>>>> Minna Win
>>>>>> National Center for Atmospheric Research
>>>>>> Developmental Testbed Center
>>>>>> Phone: 303-497-8423
>>>>>> Fax: 303-497-8401
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT <
>>>>> met_help at ucar.edu>
>>>>>> wrote:
>>>>>>
>>>>>>>
>>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
>>>>>>> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
>>>>>>> Queue: met_help
>>>>>>> Subject: AERONET AOD processing
>>>>>>> Owner: Nobody
>>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
>>>>>>> Status: new
>>>>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>>>
>>>>>>>
>>>>>>>
>>>>>>> Hello,
>>>>>>> Is a formatting code available for version 3 AERONET? (
>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
>>>>>>> (If yes) Does the code allow for selecting chosen wavelengths?
Is
>>>>> capable
>>>>>>> of obtaining selected channels (e.g. 550 nm) from other
wavelengths?
>>>>>>> Thanks,
>>>>>>> Mariusz
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>
>>
>>
>>
>>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Minna Win
Time: Fri Mar 13 11:50:18 2020
Hi Mariusz,
Let me see if this request can be submitted in the form of a GitHub
issue.
We're currently ready to release MET v9.0, so if accepted, it would
make
its way into a future release. I'll check with John to see if this
can be
done within the current budget. What time-frame constraints do you
have?
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
On Fri, Mar 13, 2020 at 11:23 AM Mariusz Pagowski via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>
> Minna,
> thanks for your dedication to answering these questions.
> We are also interested in outputting 440 and 675 nm from AERONET,
and even
> possibly 870 nm.
> We would be ready to work with someone from MET to help implement
these as
> options.
> Pls let me know if that is possible,
> Mariusz
>
>
> > On Mar 13, 2020, at 8:03 AM, Minna Win via RT <met_help at ucar.edu>
wrote:
> >
> > Hi Mariusz,
> >
> > Here's what happens when 675 and 440 nm are unavailable:
> >
> > ascii2nc does not save 550 nm wavelength into NetCDF output if 675
and
> 440
> >> are not available.
> >> ascii2nc figures out the variable list and associated wavelengths
from
> the
> >> header record and knows where AOD at 440 and 650 nm values are.
> >> ascii2nc saves AOD of all available
> >> wavelengths. And ascii2nc checks AOD at 440, 550 and
> >> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550 does
not
> exist
> >> from data record and AOD values at 675 and at 440 are valid.
> >
> >
> > I hope this is useful.
> >
> > Regards,
> > Minna
> > ---------------
> > Minna Win
> > National Center for Atmospheric Research
> > Developmental Testbed Center
> > Phone: 303-497-8423
> > Fax: 303-497-8401
> >
> >
> >
> > On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT <
> met_help at ucar.edu>
> > wrote:
> >
> >>
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>
> >> Thanks Mina,
> >> just wanted to be sure that it is not a linear interpolation.
Still I
> >> wonder what happens when 675 and 440 are not available.
> >> So if my understanding is correct MET only outputs 550nm.
> >> Mariusz
> >>
> >>
> >>
> >>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT
<met_help at ucar.edu>
> wrote:
> >>>
> >>> Hi Mariusz,
> >>>
> >>> Upon further checking from Howard, the implementation is
actually:
> >>>
> >>>
> >>> The Fortran code I referenced
(https://github.com/NCAR/MET/issues/1121
> )
> >> was
> >>>> implemented that way:
> >>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
> >>>> aerdata(i)%aot_440nm)/log(675./440.)
> >>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
> >> ((550./440.)**(-a_675_440))
> >>>
> >>>
> >>> Just in case you need to see the correct implementation code.
The
> >>> first code snippet (in previous email) he cited was replaced
with the
> >>> above code snippet, which is associated with Github issue
1121...
> >>>
> >>>
> >>> Regards,
> >>>
> >>> Minna
> >>>
> >>> ---------------
> >>> Minna Win
> >>> National Center for Atmospheric Research
> >>> Developmental Testbed Center
> >>> Phone: 303-497-8423
> >>> Fax: 303-497-8401
> >>>
> >>>
> >>>
> >>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
> >> minnawin at ucar.edu>
> >>> wrote:
> >>>
> >>>> Hi Mariusz,
> >>>>
> >>>> According to Howard:
> >>>>
> >>>> It's a log interpolation: angstrom_power_interpolation
> >>>>
> >>>> double angstrom_power_interplation(double value_1, double
value_2,
> >>>> double level_1, double level_2, double target_level) {
> >>>> double angstrom_log =
-log10(value_1/value_2)/log10(level_1/level_2);
> >>>> double angstrom_value = value_2 *
> >>>> pow((target_level/level_2),-angstrom_log);
> >>>> return angstrom_value;
> >>>> }
> >>>>
> >>>> Regards,
> >>>> Minna
> >>>> ---------------
> >>>> Minna Win
> >>>> National Center for Atmospheric Research
> >>>> Developmental Testbed Center
> >>>> Phone: 303-497-8423
> >>>> Fax: 303-497-8401
> >>>>
> >>>>
> >>>>
> >>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT <
> >> met_help at ucar.edu>
> >>>> wrote:
> >>>>
> >>>>>
> >>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>
> >>>>>
> >>>>> Minna,
> >>>>> thanks very much - just one follow-up: how is the
interpolation done
> >>>>> (linear/log)?
> >>>>> Mariusz
> >>>>>
> >>>>>
> >>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT
<met_help at ucar.edu>
> >>>>> wrote:
> >>>>>>
> >>>>>> Hi Mariusz,
> >>>>>>
> >>>>>> According to our resident expert Howard Soh:
> >>>>>>
> >>>>>> "ascii2nc reads a text format AERONET 3 data and produces a
NetCDF
> >>>>> output
> >>>>>> (MET's point observation data file). ascii2nc interpolates
550nm
> >>>>> wavelength
> >>>>>> from 440 nm and 675nm if not exists."
> >>>>>>
> >>>>>> I hope this helps.
> >>>>>>
> >>>>>> Regards,
> >>>>>> Minna
> >>>>>> ---------------
> >>>>>> Minna Win
> >>>>>> National Center for Atmospheric Research
> >>>>>> Developmental Testbed Center
> >>>>>> Phone: 303-497-8423
> >>>>>> Fax: 303-497-8401
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT <
> >>>>> met_help at ucar.edu>
> >>>>>> wrote:
> >>>>>>
> >>>>>>>
> >>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
> >>>>>>> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
> >>>>>>> Queue: met_help
> >>>>>>> Subject: AERONET AOD processing
> >>>>>>> Owner: Nobody
> >>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
> >>>>>>> Status: new
> >>>>>>> Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
> >>>
> >>>>>>>
> >>>>>>>
> >>>>>>> Hello,
> >>>>>>> Is a formatting code available for version 3 AERONET? (
> >>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
> >>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
> >>>>>>> (If yes) Does the code allow for selecting chosen
wavelengths? Is
> >>>>> capable
> >>>>>>> of obtaining selected channels (e.g. 550 nm) from other
> wavelengths?
> >>>>>>> Thanks,
> >>>>>>> Mariusz
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>
> >>
> >>
> >>
> >>
> >
>
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #94517] AERONET AOD processing
From: Mariusz Pagowski
Time: Fri Mar 13 11:54:38 2020
Minna,
our time-frame is asap and would be happy to help with coding if
someone provided guidance.
Pls let me know what your and John’s thinking is,
Mariusz
> On Mar 13, 2020, at 11:50 AM, Minna Win via RT <met_help at ucar.edu>
wrote:
>
> Hi Mariusz,
>
> Let me see if this request can be submitted in the form of a GitHub
issue.
> We're currently ready to release MET v9.0, so if accepted, it would
make
> its way into a future release. I'll check with John to see if this
can be
> done within the current budget. What time-frame constraints do you
have?
>
> Regards,
> Minna
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
>
>
>
> On Fri, Mar 13, 2020 at 11:23 AM Mariusz Pagowski via RT
<met_help at ucar.edu>
> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>
>> Minna,
>> thanks for your dedication to answering these questions.
>> We are also interested in outputting 440 and 675 nm from AERONET,
and even
>> possibly 870 nm.
>> We would be ready to work with someone from MET to help implement
these as
>> options.
>> Pls let me know if that is possible,
>> Mariusz
>>
>>
>>> On Mar 13, 2020, at 8:03 AM, Minna Win via RT <met_help at ucar.edu>
wrote:
>>>
>>> Hi Mariusz,
>>>
>>> Here's what happens when 675 and 440 nm are unavailable:
>>>
>>> ascii2nc does not save 550 nm wavelength into NetCDF output if 675
and
>> 440
>>>> are not available.
>>>> ascii2nc figures out the variable list and associated wavelengths
from
>> the
>>>> header record and knows where AOD at 440 and 650 nm values are.
>>>> ascii2nc saves AOD of all available
>>>> wavelengths. And ascii2nc checks AOD at 440, 550 and
>>>> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550 does
not
>> exist
>>>> from data record and AOD values at 675 and at 440 are valid.
>>>
>>>
>>> I hope this is useful.
>>>
>>> Regards,
>>> Minna
>>> ---------------
>>> Minna Win
>>> National Center for Atmospheric Research
>>> Developmental Testbed Center
>>> Phone: 303-497-8423
>>> Fax: 303-497-8401
>>>
>>>
>>>
>>> On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT <
>> met_help at ucar.edu>
>>> wrote:
>>>
>>>>
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>
>>>> Thanks Mina,
>>>> just wanted to be sure that it is not a linear interpolation.
Still I
>>>> wonder what happens when 675 and 440 are not available.
>>>> So if my understanding is correct MET only outputs 550nm.
>>>> Mariusz
>>>>
>>>>
>>>>
>>>>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT
<met_help at ucar.edu>
>> wrote:
>>>>>
>>>>> Hi Mariusz,
>>>>>
>>>>> Upon further checking from Howard, the implementation is
actually:
>>>>>
>>>>>
>>>>> The Fortran code I referenced
(https://github.com/NCAR/MET/issues/1121
>> )
>>>> was
>>>>>> implemented that way:
>>>>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
>>>>>> aerdata(i)%aot_440nm)/log(675./440.)
>>>>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
>>>> ((550./440.)**(-a_675_440))
>>>>>
>>>>>
>>>>> Just in case you need to see the correct implementation code.
The
>>>>> first code snippet (in previous email) he cited was replaced
with the
>>>>> above code snippet, which is associated with Github issue
1121...
>>>>>
>>>>>
>>>>> Regards,
>>>>>
>>>>> Minna
>>>>>
>>>>> ---------------
>>>>> Minna Win
>>>>> National Center for Atmospheric Research
>>>>> Developmental Testbed Center
>>>>> Phone: 303-497-8423
>>>>> Fax: 303-497-8401
>>>>>
>>>>>
>>>>>
>>>>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
>>>> minnawin at ucar.edu>
>>>>> wrote:
>>>>>
>>>>>> Hi Mariusz,
>>>>>>
>>>>>> According to Howard:
>>>>>>
>>>>>> It's a log interpolation: angstrom_power_interpolation
>>>>>>
>>>>>> double angstrom_power_interplation(double value_1, double
value_2,
>>>>>> double level_1, double level_2, double target_level) {
>>>>>> double angstrom_log =
-log10(value_1/value_2)/log10(level_1/level_2);
>>>>>> double angstrom_value = value_2 *
>>>>>> pow((target_level/level_2),-angstrom_log);
>>>>>> return angstrom_value;
>>>>>> }
>>>>>>
>>>>>> Regards,
>>>>>> Minna
>>>>>> ---------------
>>>>>> Minna Win
>>>>>> National Center for Atmospheric Research
>>>>>> Developmental Testbed Center
>>>>>> Phone: 303-497-8423
>>>>>> Fax: 303-497-8401
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT <
>>>> met_help at ucar.edu>
>>>>>> wrote:
>>>>>>
>>>>>>>
>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>
>>>>>>>
>>>>>>> Minna,
>>>>>>> thanks very much - just one follow-up: how is the
interpolation done
>>>>>>> (linear/log)?
>>>>>>> Mariusz
>>>>>>>
>>>>>>>
>>>>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT
<met_help at ucar.edu>
>>>>>>> wrote:
>>>>>>>>
>>>>>>>> Hi Mariusz,
>>>>>>>>
>>>>>>>> According to our resident expert Howard Soh:
>>>>>>>>
>>>>>>>> "ascii2nc reads a text format AERONET 3 data and produces a
NetCDF
>>>>>>> output
>>>>>>>> (MET's point observation data file). ascii2nc interpolates
550nm
>>>>>>> wavelength
>>>>>>>> from 440 nm and 675nm if not exists."
>>>>>>>>
>>>>>>>> I hope this helps.
>>>>>>>>
>>>>>>>> Regards,
>>>>>>>> Minna
>>>>>>>> ---------------
>>>>>>>> Minna Win
>>>>>>>> National Center for Atmospheric Research
>>>>>>>> Developmental Testbed Center
>>>>>>>> Phone: 303-497-8423
>>>>>>>> Fax: 303-497-8401
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT <
>>>>>>> met_help at ucar.edu>
>>>>>>>> wrote:
>>>>>>>>
>>>>>>>>>
>>>>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
>>>>>>>>> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
>>>>>>>>> Queue: met_help
>>>>>>>>> Subject: AERONET AOD processing
>>>>>>>>> Owner: Nobody
>>>>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
>>>>>>>>> Status: new
>>>>>>>>> Ticket <URL:
>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Hello,
>>>>>>>>> Is a formatting code available for version 3 AERONET? (
>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
>>>>>>>>> (If yes) Does the code allow for selecting chosen
wavelengths? Is
>>>>>>> capable
>>>>>>>>> of obtaining selected channels (e.g. 550 nm) from other
>> wavelengths?
>>>>>>>>> Thanks,
>>>>>>>>> Mariusz
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>
>>>>
>>>>
>>>>
>>>>
>>>
>>
>>
>>
>>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Minna Win
Time: Fri Mar 13 13:50:58 2020
Hi Mariusz,
I created a GitHub issue for MET:
https://github.com/NCAR/MET/issues/1278
Please feel free to add any additional comments to the GitHub issue.
I've
assigned the issue to Howard Soh and assigned this with a high
priority
(Blocker).
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
On Fri, Mar 13, 2020 at 12:05 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>
> Minna,
> our time-frame is asap and would be happy to help with coding if
someone
> provided guidance.
> Pls let me know what your and John’s thinking is,
> Mariusz
>
>
> > On Mar 13, 2020, at 11:50 AM, Minna Win via RT <met_help at ucar.edu>
> wrote:
> >
> > Hi Mariusz,
> >
> > Let me see if this request can be submitted in the form of a
GitHub
> issue.
> > We're currently ready to release MET v9.0, so if accepted, it
would make
> > its way into a future release. I'll check with John to see if
this can
> be
> > done within the current budget. What time-frame constraints do
you have?
> >
> > Regards,
> > Minna
> > ---------------
> > Minna Win
> > National Center for Atmospheric Research
> > Developmental Testbed Center
> > Phone: 303-497-8423
> > Fax: 303-497-8401
> >
> >
> >
> > On Fri, Mar 13, 2020 at 11:23 AM Mariusz Pagowski via RT <
> met_help at ucar.edu>
> > wrote:
> >
> >>
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>
> >> Minna,
> >> thanks for your dedication to answering these questions.
> >> We are also interested in outputting 440 and 675 nm from AERONET,
and
> even
> >> possibly 870 nm.
> >> We would be ready to work with someone from MET to help implement
these
> as
> >> options.
> >> Pls let me know if that is possible,
> >> Mariusz
> >>
> >>
> >>> On Mar 13, 2020, at 8:03 AM, Minna Win via RT
<met_help at ucar.edu>
> wrote:
> >>>
> >>> Hi Mariusz,
> >>>
> >>> Here's what happens when 675 and 440 nm are unavailable:
> >>>
> >>> ascii2nc does not save 550 nm wavelength into NetCDF output if
675 and
> >> 440
> >>>> are not available.
> >>>> ascii2nc figures out the variable list and associated
wavelengths from
> >> the
> >>>> header record and knows where AOD at 440 and 650 nm values are.
> >>>> ascii2nc saves AOD of all available
> >>>> wavelengths. And ascii2nc checks AOD at 440, 550 and
> >>>> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550 does
not
> >> exist
> >>>> from data record and AOD values at 675 and at 440 are valid.
> >>>
> >>>
> >>> I hope this is useful.
> >>>
> >>> Regards,
> >>> Minna
> >>> ---------------
> >>> Minna Win
> >>> National Center for Atmospheric Research
> >>> Developmental Testbed Center
> >>> Phone: 303-497-8423
> >>> Fax: 303-497-8401
> >>>
> >>>
> >>>
> >>> On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT <
> >> met_help at ucar.edu>
> >>> wrote:
> >>>
> >>>>
> >>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>>>
> >>>> Thanks Mina,
> >>>> just wanted to be sure that it is not a linear interpolation.
Still I
> >>>> wonder what happens when 675 and 440 are not available.
> >>>> So if my understanding is correct MET only outputs 550nm.
> >>>> Mariusz
> >>>>
> >>>>
> >>>>
> >>>>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT
<met_help at ucar.edu>
> >> wrote:
> >>>>>
> >>>>> Hi Mariusz,
> >>>>>
> >>>>> Upon further checking from Howard, the implementation is
actually:
> >>>>>
> >>>>>
> >>>>> The Fortran code I referenced (
> https://github.com/NCAR/MET/issues/1121
> >> )
> >>>> was
> >>>>>> implemented that way:
> >>>>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
> >>>>>> aerdata(i)%aot_440nm)/log(675./440.)
> >>>>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
> >>>> ((550./440.)**(-a_675_440))
> >>>>>
> >>>>>
> >>>>> Just in case you need to see the correct implementation code.
The
> >>>>> first code snippet (in previous email) he cited was replaced
with the
> >>>>> above code snippet, which is associated with Github issue
1121...
> >>>>>
> >>>>>
> >>>>> Regards,
> >>>>>
> >>>>> Minna
> >>>>>
> >>>>> ---------------
> >>>>> Minna Win
> >>>>> National Center for Atmospheric Research
> >>>>> Developmental Testbed Center
> >>>>> Phone: 303-497-8423
> >>>>> Fax: 303-497-8401
> >>>>>
> >>>>>
> >>>>>
> >>>>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
> >>>> minnawin at ucar.edu>
> >>>>> wrote:
> >>>>>
> >>>>>> Hi Mariusz,
> >>>>>>
> >>>>>> According to Howard:
> >>>>>>
> >>>>>> It's a log interpolation: angstrom_power_interpolation
> >>>>>>
> >>>>>> double angstrom_power_interplation(double value_1, double
value_2,
> >>>>>> double level_1, double level_2, double target_level) {
> >>>>>> double angstrom_log =
> -log10(value_1/value_2)/log10(level_1/level_2);
> >>>>>> double angstrom_value = value_2 *
> >>>>>> pow((target_level/level_2),-angstrom_log);
> >>>>>> return angstrom_value;
> >>>>>> }
> >>>>>>
> >>>>>> Regards,
> >>>>>> Minna
> >>>>>> ---------------
> >>>>>> Minna Win
> >>>>>> National Center for Atmospheric Research
> >>>>>> Developmental Testbed Center
> >>>>>> Phone: 303-497-8423
> >>>>>> Fax: 303-497-8401
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT <
> >>>> met_help at ucar.edu>
> >>>>>> wrote:
> >>>>>>
> >>>>>>>
> >>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>>>>>>
> >>>>>>> Minna,
> >>>>>>> thanks very much - just one follow-up: how is the
interpolation
> done
> >>>>>>> (linear/log)?
> >>>>>>> Mariusz
> >>>>>>>
> >>>>>>>
> >>>>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT
<met_help at ucar.edu>
> >>>>>>> wrote:
> >>>>>>>>
> >>>>>>>> Hi Mariusz,
> >>>>>>>>
> >>>>>>>> According to our resident expert Howard Soh:
> >>>>>>>>
> >>>>>>>> "ascii2nc reads a text format AERONET 3 data and produces a
NetCDF
> >>>>>>> output
> >>>>>>>> (MET's point observation data file). ascii2nc interpolates
550nm
> >>>>>>> wavelength
> >>>>>>>> from 440 nm and 675nm if not exists."
> >>>>>>>>
> >>>>>>>> I hope this helps.
> >>>>>>>>
> >>>>>>>> Regards,
> >>>>>>>> Minna
> >>>>>>>> ---------------
> >>>>>>>> Minna Win
> >>>>>>>> National Center for Atmospheric Research
> >>>>>>>> Developmental Testbed Center
> >>>>>>>> Phone: 303-497-8423
> >>>>>>>> Fax: 303-497-8401
> >>>>>>>>
> >>>>>>>>
> >>>>>>>>
> >>>>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT <
> >>>>>>> met_help at ucar.edu>
> >>>>>>>> wrote:
> >>>>>>>>
> >>>>>>>>>
> >>>>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
> >>>>>>>>> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
> >>>>>>>>> Queue: met_help
> >>>>>>>>> Subject: AERONET AOD processing
> >>>>>>>>> Owner: Nobody
> >>>>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
> >>>>>>>>> Status: new
> >>>>>>>>> Ticket <URL:
> >> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
> >>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>> Hello,
> >>>>>>>>> Is a formatting code available for version 3 AERONET? (
> >>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
> >>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
> >>>>>>>>> (If yes) Does the code allow for selecting chosen
wavelengths? Is
> >>>>>>> capable
> >>>>>>>>> of obtaining selected channels (e.g. 550 nm) from other
> >> wavelengths?
> >>>>>>>>> Thanks,
> >>>>>>>>> Mariusz
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>
> >>>>
> >>>>
> >>>>
> >>>>
> >>>
> >>
> >>
> >>
> >>
> >
>
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #94517] AERONET AOD processing
From: Mariusz Pagowski
Time: Fri Mar 13 13:53:04 2020
Thanks a lot, pls let me know if we can be helpful in any way,
Mariusz
> On Mar 13, 2020, at 1:50 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>
> Hi Mariusz,
>
> I created a GitHub issue for MET:
> https://github.com/NCAR/MET/issues/1278
>
> Please feel free to add any additional comments to the GitHub issue.
I've
> assigned the issue to Howard Soh and assigned this with a high
priority
> (Blocker).
>
> Regards,
> Minna
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
>
>
>
> On Fri, Mar 13, 2020 at 12:05 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>
>> Minna,
>> our time-frame is asap and would be happy to help with coding if
someone
>> provided guidance.
>> Pls let me know what your and John’s thinking is,
>> Mariusz
>>
>>
>>> On Mar 13, 2020, at 11:50 AM, Minna Win via RT <met_help at ucar.edu>
>> wrote:
>>>
>>> Hi Mariusz,
>>>
>>> Let me see if this request can be submitted in the form of a
GitHub
>> issue.
>>> We're currently ready to release MET v9.0, so if accepted, it
would make
>>> its way into a future release. I'll check with John to see if
this can
>> be
>>> done within the current budget. What time-frame constraints do
you have?
>>>
>>> Regards,
>>> Minna
>>> ---------------
>>> Minna Win
>>> National Center for Atmospheric Research
>>> Developmental Testbed Center
>>> Phone: 303-497-8423
>>> Fax: 303-497-8401
>>>
>>>
>>>
>>> On Fri, Mar 13, 2020 at 11:23 AM Mariusz Pagowski via RT <
>> met_help at ucar.edu>
>>> wrote:
>>>
>>>>
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>
>>>> Minna,
>>>> thanks for your dedication to answering these questions.
>>>> We are also interested in outputting 440 and 675 nm from AERONET,
and
>> even
>>>> possibly 870 nm.
>>>> We would be ready to work with someone from MET to help implement
these
>> as
>>>> options.
>>>> Pls let me know if that is possible,
>>>> Mariusz
>>>>
>>>>
>>>>> On Mar 13, 2020, at 8:03 AM, Minna Win via RT
<met_help at ucar.edu>
>> wrote:
>>>>>
>>>>> Hi Mariusz,
>>>>>
>>>>> Here's what happens when 675 and 440 nm are unavailable:
>>>>>
>>>>> ascii2nc does not save 550 nm wavelength into NetCDF output if
675 and
>>>> 440
>>>>>> are not available.
>>>>>> ascii2nc figures out the variable list and associated
wavelengths from
>>>> the
>>>>>> header record and knows where AOD at 440 and 650 nm values are.
>>>>>> ascii2nc saves AOD of all available
>>>>>> wavelengths. And ascii2nc checks AOD at 440, 550 and
>>>>>> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550 does
not
>>>> exist
>>>>>> from data record and AOD values at 675 and at 440 are valid.
>>>>>
>>>>>
>>>>> I hope this is useful.
>>>>>
>>>>> Regards,
>>>>> Minna
>>>>> ---------------
>>>>> Minna Win
>>>>> National Center for Atmospheric Research
>>>>> Developmental Testbed Center
>>>>> Phone: 303-497-8423
>>>>> Fax: 303-497-8401
>>>>>
>>>>>
>>>>>
>>>>> On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT <
>>>> met_help at ucar.edu>
>>>>> wrote:
>>>>>
>>>>>>
>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>>
>>>>>> Thanks Mina,
>>>>>> just wanted to be sure that it is not a linear interpolation.
Still I
>>>>>> wonder what happens when 675 and 440 are not available.
>>>>>> So if my understanding is correct MET only outputs 550nm.
>>>>>> Mariusz
>>>>>>
>>>>>>
>>>>>>
>>>>>>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT
<met_help at ucar.edu>
>>>> wrote:
>>>>>>>
>>>>>>> Hi Mariusz,
>>>>>>>
>>>>>>> Upon further checking from Howard, the implementation is
actually:
>>>>>>>
>>>>>>>
>>>>>>> The Fortran code I referenced (
>> https://github.com/NCAR/MET/issues/1121
>>>> )
>>>>>> was
>>>>>>>> implemented that way:
>>>>>>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
>>>>>>>> aerdata(i)%aot_440nm)/log(675./440.)
>>>>>>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
>>>>>> ((550./440.)**(-a_675_440))
>>>>>>>
>>>>>>>
>>>>>>> Just in case you need to see the correct implementation code.
The
>>>>>>> first code snippet (in previous email) he cited was replaced
with the
>>>>>>> above code snippet, which is associated with Github issue
1121...
>>>>>>>
>>>>>>>
>>>>>>> Regards,
>>>>>>>
>>>>>>> Minna
>>>>>>>
>>>>>>> ---------------
>>>>>>> Minna Win
>>>>>>> National Center for Atmospheric Research
>>>>>>> Developmental Testbed Center
>>>>>>> Phone: 303-497-8423
>>>>>>> Fax: 303-497-8401
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
>>>>>> minnawin at ucar.edu>
>>>>>>> wrote:
>>>>>>>
>>>>>>>> Hi Mariusz,
>>>>>>>>
>>>>>>>> According to Howard:
>>>>>>>>
>>>>>>>> It's a log interpolation: angstrom_power_interpolation
>>>>>>>>
>>>>>>>> double angstrom_power_interplation(double value_1, double
value_2,
>>>>>>>> double level_1, double level_2, double target_level) {
>>>>>>>> double angstrom_log =
>> -log10(value_1/value_2)/log10(level_1/level_2);
>>>>>>>> double angstrom_value = value_2 *
>>>>>>>> pow((target_level/level_2),-angstrom_log);
>>>>>>>> return angstrom_value;
>>>>>>>> }
>>>>>>>>
>>>>>>>> Regards,
>>>>>>>> Minna
>>>>>>>> ---------------
>>>>>>>> Minna Win
>>>>>>>> National Center for Atmospheric Research
>>>>>>>> Developmental Testbed Center
>>>>>>>> Phone: 303-497-8423
>>>>>>>> Fax: 303-497-8401
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT <
>>>>>> met_help at ucar.edu>
>>>>>>>> wrote:
>>>>>>>>
>>>>>>>>>
>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>>>>>
>>>>>>>>> Minna,
>>>>>>>>> thanks very much - just one follow-up: how is the
interpolation
>> done
>>>>>>>>> (linear/log)?
>>>>>>>>> Mariusz
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT
<met_help at ucar.edu>
>>>>>>>>> wrote:
>>>>>>>>>>
>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>
>>>>>>>>>> According to our resident expert Howard Soh:
>>>>>>>>>>
>>>>>>>>>> "ascii2nc reads a text format AERONET 3 data and produces a
NetCDF
>>>>>>>>> output
>>>>>>>>>> (MET's point observation data file). ascii2nc interpolates
550nm
>>>>>>>>> wavelength
>>>>>>>>>> from 440 nm and 675nm if not exists."
>>>>>>>>>>
>>>>>>>>>> I hope this helps.
>>>>>>>>>>
>>>>>>>>>> Regards,
>>>>>>>>>> Minna
>>>>>>>>>> ---------------
>>>>>>>>>> Minna Win
>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT <
>>>>>>>>> met_help at ucar.edu>
>>>>>>>>>> wrote:
>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
>>>>>>>>>>> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
>>>>>>>>>>> Queue: met_help
>>>>>>>>>>> Subject: AERONET AOD processing
>>>>>>>>>>> Owner: Nobody
>>>>>>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
>>>>>>>>>>> Status: new
>>>>>>>>>>> Ticket <URL:
>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> Hello,
>>>>>>>>>>> Is a formatting code available for version 3 AERONET? (
>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
>>>>>>>>>>> (If yes) Does the code allow for selecting chosen
wavelengths? Is
>>>>>>>>> capable
>>>>>>>>>>> of obtaining selected channels (e.g. 550 nm) from other
>>>> wavelengths?
>>>>>>>>>>> Thanks,
>>>>>>>>>>> Mariusz
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>
>>>>
>>>>
>>>>
>>>
>>
>>
>>
>>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Minna Win
Time: Fri Mar 13 13:53:47 2020
Hi Mariusz,
I put down your affiliation as NOAA/GSD on the GitHub issue. Is that
correct?
Thanks,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
On Fri, Mar 13, 2020 at 1:50 PM Minna Win-Gildenmeister
<minnawin at ucar.edu>
wrote:
> Hi Mariusz,
>
> I created a GitHub issue for MET:
> https://github.com/NCAR/MET/issues/1278
>
> Please feel free to add any additional comments to the GitHub issue.
I've
> assigned the issue to Howard Soh and assigned this with a high
priority
> (Blocker).
>
> Regards,
> Minna
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
>
>
>
> On Fri, Mar 13, 2020 at 12:05 PM Mariusz Pagowski via RT <
> met_help at ucar.edu> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>
>> Minna,
>> our time-frame is asap and would be happy to help with coding if
someone
>> provided guidance.
>> Pls let me know what your and John’s thinking is,
>> Mariusz
>>
>>
>> > On Mar 13, 2020, at 11:50 AM, Minna Win via RT
<met_help at ucar.edu>
>> wrote:
>> >
>> > Hi Mariusz,
>> >
>> > Let me see if this request can be submitted in the form of a
GitHub
>> issue.
>> > We're currently ready to release MET v9.0, so if accepted, it
would make
>> > its way into a future release. I'll check with John to see if
this can
>> be
>> > done within the current budget. What time-frame constraints do
you
>> have?
>> >
>> > Regards,
>> > Minna
>> > ---------------
>> > Minna Win
>> > National Center for Atmospheric Research
>> > Developmental Testbed Center
>> > Phone: 303-497-8423
>> > Fax: 303-497-8401
>> >
>> >
>> >
>> > On Fri, Mar 13, 2020 at 11:23 AM Mariusz Pagowski via RT <
>> met_help at ucar.edu>
>> > wrote:
>> >
>> >>
>> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>> >>
>> >> Minna,
>> >> thanks for your dedication to answering these questions.
>> >> We are also interested in outputting 440 and 675 nm from
AERONET, and
>> even
>> >> possibly 870 nm.
>> >> We would be ready to work with someone from MET to help
implement
>> these as
>> >> options.
>> >> Pls let me know if that is possible,
>> >> Mariusz
>> >>
>> >>
>> >>> On Mar 13, 2020, at 8:03 AM, Minna Win via RT
<met_help at ucar.edu>
>> wrote:
>> >>>
>> >>> Hi Mariusz,
>> >>>
>> >>> Here's what happens when 675 and 440 nm are unavailable:
>> >>>
>> >>> ascii2nc does not save 550 nm wavelength into NetCDF output if
675 and
>> >> 440
>> >>>> are not available.
>> >>>> ascii2nc figures out the variable list and associated
wavelengths
>> from
>> >> the
>> >>>> header record and knows where AOD at 440 and 650 nm values
are.
>> >>>> ascii2nc saves AOD of all available
>> >>>> wavelengths. And ascii2nc checks AOD at 440, 550 and
>> >>>> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550
does not
>> >> exist
>> >>>> from data record and AOD values at 675 and at 440 are valid.
>> >>>
>> >>>
>> >>> I hope this is useful.
>> >>>
>> >>> Regards,
>> >>> Minna
>> >>> ---------------
>> >>> Minna Win
>> >>> National Center for Atmospheric Research
>> >>> Developmental Testbed Center
>> >>> Phone: 303-497-8423
>> >>> Fax: 303-497-8401
>> >>>
>> >>>
>> >>>
>> >>> On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT <
>> >> met_help at ucar.edu>
>> >>> wrote:
>> >>>
>> >>>>
>> >>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>
>> >>>>
>> >>>> Thanks Mina,
>> >>>> just wanted to be sure that it is not a linear interpolation.
Still I
>> >>>> wonder what happens when 675 and 440 are not available.
>> >>>> So if my understanding is correct MET only outputs 550nm.
>> >>>> Mariusz
>> >>>>
>> >>>>
>> >>>>
>> >>>>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT
<met_help at ucar.edu>
>> >> wrote:
>> >>>>>
>> >>>>> Hi Mariusz,
>> >>>>>
>> >>>>> Upon further checking from Howard, the implementation is
actually:
>> >>>>>
>> >>>>>
>> >>>>> The Fortran code I referenced (
>> https://github.com/NCAR/MET/issues/1121
>> >> )
>> >>>> was
>> >>>>>> implemented that way:
>> >>>>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
>> >>>>>> aerdata(i)%aot_440nm)/log(675./440.)
>> >>>>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
>> >>>> ((550./440.)**(-a_675_440))
>> >>>>>
>> >>>>>
>> >>>>> Just in case you need to see the correct implementation code.
The
>> >>>>> first code snippet (in previous email) he cited was replaced
with
>> the
>> >>>>> above code snippet, which is associated with Github issue
1121...
>> >>>>>
>> >>>>>
>> >>>>> Regards,
>> >>>>>
>> >>>>> Minna
>> >>>>>
>> >>>>> ---------------
>> >>>>> Minna Win
>> >>>>> National Center for Atmospheric Research
>> >>>>> Developmental Testbed Center
>> >>>>> Phone: 303-497-8423
>> >>>>> Fax: 303-497-8401
>> >>>>>
>> >>>>>
>> >>>>>
>> >>>>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
>> >>>> minnawin at ucar.edu>
>> >>>>> wrote:
>> >>>>>
>> >>>>>> Hi Mariusz,
>> >>>>>>
>> >>>>>> According to Howard:
>> >>>>>>
>> >>>>>> It's a log interpolation: angstrom_power_interpolation
>> >>>>>>
>> >>>>>> double angstrom_power_interplation(double value_1, double
value_2,
>> >>>>>> double level_1, double level_2, double target_level) {
>> >>>>>> double angstrom_log =
>> -log10(value_1/value_2)/log10(level_1/level_2);
>> >>>>>> double angstrom_value = value_2 *
>> >>>>>> pow((target_level/level_2),-angstrom_log);
>> >>>>>> return angstrom_value;
>> >>>>>> }
>> >>>>>>
>> >>>>>> Regards,
>> >>>>>> Minna
>> >>>>>> ---------------
>> >>>>>> Minna Win
>> >>>>>> National Center for Atmospheric Research
>> >>>>>> Developmental Testbed Center
>> >>>>>> Phone: 303-497-8423
>> >>>>>> Fax: 303-497-8401
>> >>>>>>
>> >>>>>>
>> >>>>>>
>> >>>>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT <
>> >>>> met_help at ucar.edu>
>> >>>>>> wrote:
>> >>>>>>
>> >>>>>>>
>> >>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>> >>>>>>>
>> >>>>>>> Minna,
>> >>>>>>> thanks very much - just one follow-up: how is the
interpolation
>> done
>> >>>>>>> (linear/log)?
>> >>>>>>> Mariusz
>> >>>>>>>
>> >>>>>>>
>> >>>>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT
<met_help at ucar.edu
>> >
>> >>>>>>> wrote:
>> >>>>>>>>
>> >>>>>>>> Hi Mariusz,
>> >>>>>>>>
>> >>>>>>>> According to our resident expert Howard Soh:
>> >>>>>>>>
>> >>>>>>>> "ascii2nc reads a text format AERONET 3 data and produces
a
>> NetCDF
>> >>>>>>> output
>> >>>>>>>> (MET's point observation data file). ascii2nc interpolates
550nm
>> >>>>>>> wavelength
>> >>>>>>>> from 440 nm and 675nm if not exists."
>> >>>>>>>>
>> >>>>>>>> I hope this helps.
>> >>>>>>>>
>> >>>>>>>> Regards,
>> >>>>>>>> Minna
>> >>>>>>>> ---------------
>> >>>>>>>> Minna Win
>> >>>>>>>> National Center for Atmospheric Research
>> >>>>>>>> Developmental Testbed Center
>> >>>>>>>> Phone: 303-497-8423
>> >>>>>>>> Fax: 303-497-8401
>> >>>>>>>>
>> >>>>>>>>
>> >>>>>>>>
>> >>>>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT <
>> >>>>>>> met_help at ucar.edu>
>> >>>>>>>> wrote:
>> >>>>>>>>
>> >>>>>>>>>
>> >>>>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
>> >>>>>>>>> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
>> >>>>>>>>> Queue: met_help
>> >>>>>>>>> Subject: AERONET AOD processing
>> >>>>>>>>> Owner: Nobody
>> >>>>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
>> >>>>>>>>> Status: new
>> >>>>>>>>> Ticket <URL:
>> >> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>> >>>>>
>> >>>>>>>>>
>> >>>>>>>>>
>> >>>>>>>>> Hello,
>> >>>>>>>>> Is a formatting code available for version 3 AERONET? (
>> >>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
>> >>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
>> >>>>>>>>> (If yes) Does the code allow for selecting chosen
wavelengths?
>> Is
>> >>>>>>> capable
>> >>>>>>>>> of obtaining selected channels (e.g. 550 nm) from other
>> >> wavelengths?
>> >>>>>>>>> Thanks,
>> >>>>>>>>> Mariusz
>> >>>>>>>>>
>> >>>>>>>>>
>> >>>>>>>>>
>> >>>>>>>>
>> >>>>>>>
>> >>>>>>>
>> >>>>>>>
>> >>>>>>>
>> >>>>>
>> >>>>
>> >>>>
>> >>>>
>> >>>>
>> >>>
>> >>
>> >>
>> >>
>> >>
>> >
>>
>>
>>
>>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #94517] AERONET AOD processing
From: Mariusz Pagowski
Time: Fri Mar 13 13:56:18 2020
Yes, thanks,
Mariusz
> On Mar 13, 2020, at 1:53 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>
> Hi Mariusz,
>
> I put down your affiliation as NOAA/GSD on the GitHub issue. Is
that
> correct?
>
> Thanks,
> Minna
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
>
>
>
> On Fri, Mar 13, 2020 at 1:50 PM Minna Win-Gildenmeister
<minnawin at ucar.edu>
> wrote:
>
>> Hi Mariusz,
>>
>> I created a GitHub issue for MET:
>> https://github.com/NCAR/MET/issues/1278
>>
>> Please feel free to add any additional comments to the GitHub
issue. I've
>> assigned the issue to Howard Soh and assigned this with a high
priority
>> (Blocker).
>>
>> Regards,
>> Minna
>> ---------------
>> Minna Win
>> National Center for Atmospheric Research
>> Developmental Testbed Center
>> Phone: 303-497-8423
>> Fax: 303-497-8401
>>
>>
>>
>> On Fri, Mar 13, 2020 at 12:05 PM Mariusz Pagowski via RT <
>> met_help at ucar.edu> wrote:
>>
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>
>>> Minna,
>>> our time-frame is asap and would be happy to help with coding if
someone
>>> provided guidance.
>>> Pls let me know what your and John’s thinking is,
>>> Mariusz
>>>
>>>
>>>> On Mar 13, 2020, at 11:50 AM, Minna Win via RT
<met_help at ucar.edu>
>>> wrote:
>>>>
>>>> Hi Mariusz,
>>>>
>>>> Let me see if this request can be submitted in the form of a
GitHub
>>> issue.
>>>> We're currently ready to release MET v9.0, so if accepted, it
would make
>>>> its way into a future release. I'll check with John to see if
this can
>>> be
>>>> done within the current budget. What time-frame constraints do
you
>>> have?
>>>>
>>>> Regards,
>>>> Minna
>>>> ---------------
>>>> Minna Win
>>>> National Center for Atmospheric Research
>>>> Developmental Testbed Center
>>>> Phone: 303-497-8423
>>>> Fax: 303-497-8401
>>>>
>>>>
>>>>
>>>> On Fri, Mar 13, 2020 at 11:23 AM Mariusz Pagowski via RT <
>>> met_help at ucar.edu>
>>>> wrote:
>>>>
>>>>>
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>
>>>>> Minna,
>>>>> thanks for your dedication to answering these questions.
>>>>> We are also interested in outputting 440 and 675 nm from
AERONET, and
>>> even
>>>>> possibly 870 nm.
>>>>> We would be ready to work with someone from MET to help
implement
>>> these as
>>>>> options.
>>>>> Pls let me know if that is possible,
>>>>> Mariusz
>>>>>
>>>>>
>>>>>> On Mar 13, 2020, at 8:03 AM, Minna Win via RT
<met_help at ucar.edu>
>>> wrote:
>>>>>>
>>>>>> Hi Mariusz,
>>>>>>
>>>>>> Here's what happens when 675 and 440 nm are unavailable:
>>>>>>
>>>>>> ascii2nc does not save 550 nm wavelength into NetCDF output if
675 and
>>>>> 440
>>>>>>> are not available.
>>>>>>> ascii2nc figures out the variable list and associated
wavelengths
>>> from
>>>>> the
>>>>>>> header record and knows where AOD at 440 and 650 nm values
are.
>>>>>>> ascii2nc saves AOD of all available
>>>>>>> wavelengths. And ascii2nc checks AOD at 440, 550 and
>>>>>>> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550
does not
>>>>> exist
>>>>>>> from data record and AOD values at 675 and at 440 are valid.
>>>>>>
>>>>>>
>>>>>> I hope this is useful.
>>>>>>
>>>>>> Regards,
>>>>>> Minna
>>>>>> ---------------
>>>>>> Minna Win
>>>>>> National Center for Atmospheric Research
>>>>>> Developmental Testbed Center
>>>>>> Phone: 303-497-8423
>>>>>> Fax: 303-497-8401
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT <
>>>>> met_help at ucar.edu>
>>>>>> wrote:
>>>>>>
>>>>>>>
>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>
>>>>>>>
>>>>>>> Thanks Mina,
>>>>>>> just wanted to be sure that it is not a linear interpolation.
Still I
>>>>>>> wonder what happens when 675 and 440 are not available.
>>>>>>> So if my understanding is correct MET only outputs 550nm.
>>>>>>> Mariusz
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT
<met_help at ucar.edu>
>>>>> wrote:
>>>>>>>>
>>>>>>>> Hi Mariusz,
>>>>>>>>
>>>>>>>> Upon further checking from Howard, the implementation is
actually:
>>>>>>>>
>>>>>>>>
>>>>>>>> The Fortran code I referenced (
>>> https://github.com/NCAR/MET/issues/1121
>>>>> )
>>>>>>> was
>>>>>>>>> implemented that way:
>>>>>>>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
>>>>>>>>> aerdata(i)%aot_440nm)/log(675./440.)
>>>>>>>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
>>>>>>> ((550./440.)**(-a_675_440))
>>>>>>>>
>>>>>>>>
>>>>>>>> Just in case you need to see the correct implementation code.
The
>>>>>>>> first code snippet (in previous email) he cited was replaced
with
>>> the
>>>>>>>> above code snippet, which is associated with Github issue
1121...
>>>>>>>>
>>>>>>>>
>>>>>>>> Regards,
>>>>>>>>
>>>>>>>> Minna
>>>>>>>>
>>>>>>>> ---------------
>>>>>>>> Minna Win
>>>>>>>> National Center for Atmospheric Research
>>>>>>>> Developmental Testbed Center
>>>>>>>> Phone: 303-497-8423
>>>>>>>> Fax: 303-497-8401
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
>>>>>>> minnawin at ucar.edu>
>>>>>>>> wrote:
>>>>>>>>
>>>>>>>>> Hi Mariusz,
>>>>>>>>>
>>>>>>>>> According to Howard:
>>>>>>>>>
>>>>>>>>> It's a log interpolation: angstrom_power_interpolation
>>>>>>>>>
>>>>>>>>> double angstrom_power_interplation(double value_1, double
value_2,
>>>>>>>>> double level_1, double level_2, double target_level) {
>>>>>>>>> double angstrom_log =
>>> -log10(value_1/value_2)/log10(level_1/level_2);
>>>>>>>>> double angstrom_value = value_2 *
>>>>>>>>> pow((target_level/level_2),-angstrom_log);
>>>>>>>>> return angstrom_value;
>>>>>>>>> }
>>>>>>>>>
>>>>>>>>> Regards,
>>>>>>>>> Minna
>>>>>>>>> ---------------
>>>>>>>>> Minna Win
>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>> Developmental Testbed Center
>>>>>>>>> Phone: 303-497-8423
>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT <
>>>>>>> met_help at ucar.edu>
>>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>>>>>>
>>>>>>>>>> Minna,
>>>>>>>>>> thanks very much - just one follow-up: how is the
interpolation
>>> done
>>>>>>>>>> (linear/log)?
>>>>>>>>>> Mariusz
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT
<met_help at ucar.edu
>>>>
>>>>>>>>>> wrote:
>>>>>>>>>>>
>>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>>
>>>>>>>>>>> According to our resident expert Howard Soh:
>>>>>>>>>>>
>>>>>>>>>>> "ascii2nc reads a text format AERONET 3 data and produces
a
>>> NetCDF
>>>>>>>>>> output
>>>>>>>>>>> (MET's point observation data file). ascii2nc interpolates
550nm
>>>>>>>>>> wavelength
>>>>>>>>>>> from 440 nm and 675nm if not exists."
>>>>>>>>>>>
>>>>>>>>>>> I hope this helps.
>>>>>>>>>>>
>>>>>>>>>>> Regards,
>>>>>>>>>>> Minna
>>>>>>>>>>> ---------------
>>>>>>>>>>> Minna Win
>>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT <
>>>>>>>>>> met_help at ucar.edu>
>>>>>>>>>>> wrote:
>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
>>>>>>>>>>>> Transaction: Ticket created by Mariusz.Pagowski at noaa.gov
>>>>>>>>>>>> Queue: met_help
>>>>>>>>>>>> Subject: AERONET AOD processing
>>>>>>>>>>>> Owner: Nobody
>>>>>>>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
>>>>>>>>>>>> Status: new
>>>>>>>>>>>> Ticket <URL:
>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> Hello,
>>>>>>>>>>>> Is a formatting code available for version 3 AERONET? (
>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
>>>>>>>>>>>> (If yes) Does the code allow for selecting chosen
wavelengths?
>>> Is
>>>>>>>>>> capable
>>>>>>>>>>>> of obtaining selected channels (e.g. 550 nm) from other
>>>>> wavelengths?
>>>>>>>>>>>> Thanks,
>>>>>>>>>>>> Mariusz
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>
>>>
>>>
>>>
>>>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Minna Win
Time: Fri Mar 13 14:08:19 2020
Hi Mariusz,
I updated the GitHub Issue. Howard reminded me that MET 9.0 has
support
for 440nm and 675nm AERONET3. It looks like 870nm interpolation
support
will need to be added.
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
On Fri, Mar 13, 2020 at 2:05 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>
> Yes, thanks,
> Mariusz
>
> > On Mar 13, 2020, at 1:53 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
> >
> > Hi Mariusz,
> >
> > I put down your affiliation as NOAA/GSD on the GitHub issue. Is
that
> > correct?
> >
> > Thanks,
> > Minna
> > ---------------
> > Minna Win
> > National Center for Atmospheric Research
> > Developmental Testbed Center
> > Phone: 303-497-8423
> > Fax: 303-497-8401
> >
> >
> >
> > On Fri, Mar 13, 2020 at 1:50 PM Minna Win-Gildenmeister <
> minnawin at ucar.edu>
> > wrote:
> >
> >> Hi Mariusz,
> >>
> >> I created a GitHub issue for MET:
> >> https://github.com/NCAR/MET/issues/1278
> >>
> >> Please feel free to add any additional comments to the GitHub
issue.
> I've
> >> assigned the issue to Howard Soh and assigned this with a high
priority
> >> (Blocker).
> >>
> >> Regards,
> >> Minna
> >> ---------------
> >> Minna Win
> >> National Center for Atmospheric Research
> >> Developmental Testbed Center
> >> Phone: 303-497-8423
> >> Fax: 303-497-8401
> >>
> >>
> >>
> >> On Fri, Mar 13, 2020 at 12:05 PM Mariusz Pagowski via RT <
> >> met_help at ucar.edu> wrote:
> >>
> >>>
> >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>>
> >>> Minna,
> >>> our time-frame is asap and would be happy to help with coding if
> someone
> >>> provided guidance.
> >>> Pls let me know what your and John’s thinking is,
> >>> Mariusz
> >>>
> >>>
> >>>> On Mar 13, 2020, at 11:50 AM, Minna Win via RT
<met_help at ucar.edu>
> >>> wrote:
> >>>>
> >>>> Hi Mariusz,
> >>>>
> >>>> Let me see if this request can be submitted in the form of a
GitHub
> >>> issue.
> >>>> We're currently ready to release MET v9.0, so if accepted, it
would
> make
> >>>> its way into a future release. I'll check with John to see if
this
> can
> >>> be
> >>>> done within the current budget. What time-frame constraints do
you
> >>> have?
> >>>>
> >>>> Regards,
> >>>> Minna
> >>>> ---------------
> >>>> Minna Win
> >>>> National Center for Atmospheric Research
> >>>> Developmental Testbed Center
> >>>> Phone: 303-497-8423
> >>>> Fax: 303-497-8401
> >>>>
> >>>>
> >>>>
> >>>> On Fri, Mar 13, 2020 at 11:23 AM Mariusz Pagowski via RT <
> >>> met_help at ucar.edu>
> >>>> wrote:
> >>>>
> >>>>>
> >>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>
> >>>>>
> >>>>> Minna,
> >>>>> thanks for your dedication to answering these questions.
> >>>>> We are also interested in outputting 440 and 675 nm from
AERONET, and
> >>> even
> >>>>> possibly 870 nm.
> >>>>> We would be ready to work with someone from MET to help
implement
> >>> these as
> >>>>> options.
> >>>>> Pls let me know if that is possible,
> >>>>> Mariusz
> >>>>>
> >>>>>
> >>>>>> On Mar 13, 2020, at 8:03 AM, Minna Win via RT
<met_help at ucar.edu>
> >>> wrote:
> >>>>>>
> >>>>>> Hi Mariusz,
> >>>>>>
> >>>>>> Here's what happens when 675 and 440 nm are unavailable:
> >>>>>>
> >>>>>> ascii2nc does not save 550 nm wavelength into NetCDF output
if 675
> and
> >>>>> 440
> >>>>>>> are not available.
> >>>>>>> ascii2nc figures out the variable list and associated
wavelengths
> >>> from
> >>>>> the
> >>>>>>> header record and knows where AOD at 440 and 650 nm values
are.
> >>>>>>> ascii2nc saves AOD of all available
> >>>>>>> wavelengths. And ascii2nc checks AOD at 440, 550 and
> >>>>>>> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550
does not
> >>>>> exist
> >>>>>>> from data record and AOD values at 675 and at 440 are valid.
> >>>>>>
> >>>>>>
> >>>>>> I hope this is useful.
> >>>>>>
> >>>>>> Regards,
> >>>>>> Minna
> >>>>>> ---------------
> >>>>>> Minna Win
> >>>>>> National Center for Atmospheric Research
> >>>>>> Developmental Testbed Center
> >>>>>> Phone: 303-497-8423
> >>>>>> Fax: 303-497-8401
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>> On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT <
> >>>>> met_help at ucar.edu>
> >>>>>> wrote:
> >>>>>>
> >>>>>>>
> >>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>>>>>>
> >>>>>>> Thanks Mina,
> >>>>>>> just wanted to be sure that it is not a linear
interpolation.
> Still I
> >>>>>>> wonder what happens when 675 and 440 are not available.
> >>>>>>> So if my understanding is correct MET only outputs 550nm.
> >>>>>>> Mariusz
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT
<met_help at ucar.edu>
> >>>>> wrote:
> >>>>>>>>
> >>>>>>>> Hi Mariusz,
> >>>>>>>>
> >>>>>>>> Upon further checking from Howard, the implementation is
actually:
> >>>>>>>>
> >>>>>>>>
> >>>>>>>> The Fortran code I referenced (
> >>> https://github.com/NCAR/MET/issues/1121
> >>>>> )
> >>>>>>> was
> >>>>>>>>> implemented that way:
> >>>>>>>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
> >>>>>>>>> aerdata(i)%aot_440nm)/log(675./440.)
> >>>>>>>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
> >>>>>>> ((550./440.)**(-a_675_440))
> >>>>>>>>
> >>>>>>>>
> >>>>>>>> Just in case you need to see the correct implementation
code. The
> >>>>>>>> first code snippet (in previous email) he cited was
replaced with
> >>> the
> >>>>>>>> above code snippet, which is associated with Github issue
1121...
> >>>>>>>>
> >>>>>>>>
> >>>>>>>> Regards,
> >>>>>>>>
> >>>>>>>> Minna
> >>>>>>>>
> >>>>>>>> ---------------
> >>>>>>>> Minna Win
> >>>>>>>> National Center for Atmospheric Research
> >>>>>>>> Developmental Testbed Center
> >>>>>>>> Phone: 303-497-8423
> >>>>>>>> Fax: 303-497-8401
> >>>>>>>>
> >>>>>>>>
> >>>>>>>>
> >>>>>>>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
> >>>>>>> minnawin at ucar.edu>
> >>>>>>>> wrote:
> >>>>>>>>
> >>>>>>>>> Hi Mariusz,
> >>>>>>>>>
> >>>>>>>>> According to Howard:
> >>>>>>>>>
> >>>>>>>>> It's a log interpolation: angstrom_power_interpolation
> >>>>>>>>>
> >>>>>>>>> double angstrom_power_interplation(double value_1, double
> value_2,
> >>>>>>>>> double level_1, double level_2, double target_level) {
> >>>>>>>>> double angstrom_log =
> >>> -log10(value_1/value_2)/log10(level_1/level_2);
> >>>>>>>>> double angstrom_value = value_2 *
> >>>>>>>>> pow((target_level/level_2),-angstrom_log);
> >>>>>>>>> return angstrom_value;
> >>>>>>>>> }
> >>>>>>>>>
> >>>>>>>>> Regards,
> >>>>>>>>> Minna
> >>>>>>>>> ---------------
> >>>>>>>>> Minna Win
> >>>>>>>>> National Center for Atmospheric Research
> >>>>>>>>> Developmental Testbed Center
> >>>>>>>>> Phone: 303-497-8423
> >>>>>>>>> Fax: 303-497-8401
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT <
> >>>>>>> met_help at ucar.edu>
> >>>>>>>>> wrote:
> >>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>>>>>>>>>
> >>>>>>>>>> Minna,
> >>>>>>>>>> thanks very much - just one follow-up: how is the
interpolation
> >>> done
> >>>>>>>>>> (linear/log)?
> >>>>>>>>>> Mariusz
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT <
> met_help at ucar.edu
> >>>>
> >>>>>>>>>> wrote:
> >>>>>>>>>>>
> >>>>>>>>>>> Hi Mariusz,
> >>>>>>>>>>>
> >>>>>>>>>>> According to our resident expert Howard Soh:
> >>>>>>>>>>>
> >>>>>>>>>>> "ascii2nc reads a text format AERONET 3 data and
produces a
> >>> NetCDF
> >>>>>>>>>> output
> >>>>>>>>>>> (MET's point observation data file). ascii2nc
interpolates
> 550nm
> >>>>>>>>>> wavelength
> >>>>>>>>>>> from 440 nm and 675nm if not exists."
> >>>>>>>>>>>
> >>>>>>>>>>> I hope this helps.
> >>>>>>>>>>>
> >>>>>>>>>>> Regards,
> >>>>>>>>>>> Minna
> >>>>>>>>>>> ---------------
> >>>>>>>>>>> Minna Win
> >>>>>>>>>>> National Center for Atmospheric Research
> >>>>>>>>>>> Developmental Testbed Center
> >>>>>>>>>>> Phone: 303-497-8423
> >>>>>>>>>>> Fax: 303-497-8401
> >>>>>>>>>>>
> >>>>>>>>>>>
> >>>>>>>>>>>
> >>>>>>>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT
<
> >>>>>>>>>> met_help at ucar.edu>
> >>>>>>>>>>> wrote:
> >>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
> >>>>>>>>>>>> Transaction: Ticket created by
Mariusz.Pagowski at noaa.gov
> >>>>>>>>>>>> Queue: met_help
> >>>>>>>>>>>> Subject: AERONET AOD processing
> >>>>>>>>>>>> Owner: Nobody
> >>>>>>>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
> >>>>>>>>>>>> Status: new
> >>>>>>>>>>>> Ticket <URL:
> >>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
> >>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>> Hello,
> >>>>>>>>>>>> Is a formatting code available for version 3 AERONET? (
> >>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
> >>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
> >>>>>>>>>>>> (If yes) Does the code allow for selecting chosen
wavelengths?
> >>> Is
> >>>>>>>>>> capable
> >>>>>>>>>>>> of obtaining selected channels (e.g. 550 nm) from other
> >>>>> wavelengths?
> >>>>>>>>>>>> Thanks,
> >>>>>>>>>>>> Mariusz
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>
> >>>
> >>>
> >>>
> >>>
> >
>
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #94517] AERONET AOD processing
From: Mariusz Pagowski
Time: Fri Mar 13 14:12:22 2020
Minna,
To clarify:
1. MET 9.0 can output 440 nm and 675 nm? (along with interpolated
550nm)
2. 870nm would be an addition - it would be fine only to output if
exists
Mariusz
> On Mar 13, 2020, at 2:08 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>
> Hi Mariusz,
>
> I updated the GitHub Issue. Howard reminded me that MET 9.0 has
support
> for 440nm and 675nm AERONET3. It looks like 870nm interpolation
support
> will need to be added.
>
> Regards,
> Minna
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
>
>
>
> On Fri, Mar 13, 2020 at 2:05 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>
>> Yes, thanks,
>> Mariusz
>>
>>> On Mar 13, 2020, at 1:53 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>>>
>>> Hi Mariusz,
>>>
>>> I put down your affiliation as NOAA/GSD on the GitHub issue. Is
that
>>> correct?
>>>
>>> Thanks,
>>> Minna
>>> ---------------
>>> Minna Win
>>> National Center for Atmospheric Research
>>> Developmental Testbed Center
>>> Phone: 303-497-8423
>>> Fax: 303-497-8401
>>>
>>>
>>>
>>> On Fri, Mar 13, 2020 at 1:50 PM Minna Win-Gildenmeister <
>> minnawin at ucar.edu>
>>> wrote:
>>>
>>>> Hi Mariusz,
>>>>
>>>> I created a GitHub issue for MET:
>>>> https://github.com/NCAR/MET/issues/1278
>>>>
>>>> Please feel free to add any additional comments to the GitHub
issue.
>> I've
>>>> assigned the issue to Howard Soh and assigned this with a high
priority
>>>> (Blocker).
>>>>
>>>> Regards,
>>>> Minna
>>>> ---------------
>>>> Minna Win
>>>> National Center for Atmospheric Research
>>>> Developmental Testbed Center
>>>> Phone: 303-497-8423
>>>> Fax: 303-497-8401
>>>>
>>>>
>>>>
>>>> On Fri, Mar 13, 2020 at 12:05 PM Mariusz Pagowski via RT <
>>>> met_help at ucar.edu> wrote:
>>>>
>>>>>
>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>
>>>>> Minna,
>>>>> our time-frame is asap and would be happy to help with coding if
>> someone
>>>>> provided guidance.
>>>>> Pls let me know what your and John’s thinking is,
>>>>> Mariusz
>>>>>
>>>>>
>>>>>> On Mar 13, 2020, at 11:50 AM, Minna Win via RT
<met_help at ucar.edu>
>>>>> wrote:
>>>>>>
>>>>>> Hi Mariusz,
>>>>>>
>>>>>> Let me see if this request can be submitted in the form of a
GitHub
>>>>> issue.
>>>>>> We're currently ready to release MET v9.0, so if accepted, it
would
>> make
>>>>>> its way into a future release. I'll check with John to see if
this
>> can
>>>>> be
>>>>>> done within the current budget. What time-frame constraints do
you
>>>>> have?
>>>>>>
>>>>>> Regards,
>>>>>> Minna
>>>>>> ---------------
>>>>>> Minna Win
>>>>>> National Center for Atmospheric Research
>>>>>> Developmental Testbed Center
>>>>>> Phone: 303-497-8423
>>>>>> Fax: 303-497-8401
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Fri, Mar 13, 2020 at 11:23 AM Mariusz Pagowski via RT <
>>>>> met_help at ucar.edu>
>>>>>> wrote:
>>>>>>
>>>>>>>
>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>
>>>>>>>
>>>>>>> Minna,
>>>>>>> thanks for your dedication to answering these questions.
>>>>>>> We are also interested in outputting 440 and 675 nm from
AERONET, and
>>>>> even
>>>>>>> possibly 870 nm.
>>>>>>> We would be ready to work with someone from MET to help
implement
>>>>> these as
>>>>>>> options.
>>>>>>> Pls let me know if that is possible,
>>>>>>> Mariusz
>>>>>>>
>>>>>>>
>>>>>>>> On Mar 13, 2020, at 8:03 AM, Minna Win via RT
<met_help at ucar.edu>
>>>>> wrote:
>>>>>>>>
>>>>>>>> Hi Mariusz,
>>>>>>>>
>>>>>>>> Here's what happens when 675 and 440 nm are unavailable:
>>>>>>>>
>>>>>>>> ascii2nc does not save 550 nm wavelength into NetCDF output
if 675
>> and
>>>>>>> 440
>>>>>>>>> are not available.
>>>>>>>>> ascii2nc figures out the variable list and associated
wavelengths
>>>>> from
>>>>>>> the
>>>>>>>>> header record and knows where AOD at 440 and 650 nm values
are.
>>>>>>>>> ascii2nc saves AOD of all available
>>>>>>>>> wavelengths. And ascii2nc checks AOD at 440, 550 and
>>>>>>>>> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550
does not
>>>>>>> exist
>>>>>>>>> from data record and AOD values at 675 and at 440 are valid.
>>>>>>>>
>>>>>>>>
>>>>>>>> I hope this is useful.
>>>>>>>>
>>>>>>>> Regards,
>>>>>>>> Minna
>>>>>>>> ---------------
>>>>>>>> Minna Win
>>>>>>>> National Center for Atmospheric Research
>>>>>>>> Developmental Testbed Center
>>>>>>>> Phone: 303-497-8423
>>>>>>>> Fax: 303-497-8401
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT <
>>>>>>> met_help at ucar.edu>
>>>>>>>> wrote:
>>>>>>>>
>>>>>>>>>
>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>>>>>
>>>>>>>>> Thanks Mina,
>>>>>>>>> just wanted to be sure that it is not a linear
interpolation.
>> Still I
>>>>>>>>> wonder what happens when 675 and 440 are not available.
>>>>>>>>> So if my understanding is correct MET only outputs 550nm.
>>>>>>>>> Mariusz
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT
<met_help at ucar.edu>
>>>>>>> wrote:
>>>>>>>>>>
>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>
>>>>>>>>>> Upon further checking from Howard, the implementation is
actually:
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> The Fortran code I referenced (
>>>>> https://github.com/NCAR/MET/issues/1121
>>>>>>> )
>>>>>>>>> was
>>>>>>>>>>> implemented that way:
>>>>>>>>>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
>>>>>>>>>>> aerdata(i)%aot_440nm)/log(675./440.)
>>>>>>>>>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
>>>>>>>>> ((550./440.)**(-a_675_440))
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> Just in case you need to see the correct implementation
code. The
>>>>>>>>>> first code snippet (in previous email) he cited was
replaced with
>>>>> the
>>>>>>>>>> above code snippet, which is associated with Github issue
1121...
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> Regards,
>>>>>>>>>>
>>>>>>>>>> Minna
>>>>>>>>>>
>>>>>>>>>> ---------------
>>>>>>>>>> Minna Win
>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
>>>>>>>>> minnawin at ucar.edu>
>>>>>>>>>> wrote:
>>>>>>>>>>
>>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>>
>>>>>>>>>>> According to Howard:
>>>>>>>>>>>
>>>>>>>>>>> It's a log interpolation: angstrom_power_interpolation
>>>>>>>>>>>
>>>>>>>>>>> double angstrom_power_interplation(double value_1, double
>> value_2,
>>>>>>>>>>> double level_1, double level_2, double target_level) {
>>>>>>>>>>> double angstrom_log =
>>>>> -log10(value_1/value_2)/log10(level_1/level_2);
>>>>>>>>>>> double angstrom_value = value_2 *
>>>>>>>>>>> pow((target_level/level_2),-angstrom_log);
>>>>>>>>>>> return angstrom_value;
>>>>>>>>>>> }
>>>>>>>>>>>
>>>>>>>>>>> Regards,
>>>>>>>>>>> Minna
>>>>>>>>>>> ---------------
>>>>>>>>>>> Minna Win
>>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT <
>>>>>>>>> met_help at ucar.edu>
>>>>>>>>>>> wrote:
>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>>>>>>>>
>>>>>>>>>>>> Minna,
>>>>>>>>>>>> thanks very much - just one follow-up: how is the
interpolation
>>>>> done
>>>>>>>>>>>> (linear/log)?
>>>>>>>>>>>> Mariusz
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT <
>> met_help at ucar.edu
>>>>>>
>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>
>>>>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>>>>
>>>>>>>>>>>>> According to our resident expert Howard Soh:
>>>>>>>>>>>>>
>>>>>>>>>>>>> "ascii2nc reads a text format AERONET 3 data and
produces a
>>>>> NetCDF
>>>>>>>>>>>> output
>>>>>>>>>>>>> (MET's point observation data file). ascii2nc
interpolates
>> 550nm
>>>>>>>>>>>> wavelength
>>>>>>>>>>>>> from 440 nm and 675nm if not exists."
>>>>>>>>>>>>>
>>>>>>>>>>>>> I hope this helps.
>>>>>>>>>>>>>
>>>>>>>>>>>>> Regards,
>>>>>>>>>>>>> Minna
>>>>>>>>>>>>> ---------------
>>>>>>>>>>>>> Minna Win
>>>>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT
<
>>>>>>>>>>>> met_help at ucar.edu>
>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted upon.
>>>>>>>>>>>>>> Transaction: Ticket created by
Mariusz.Pagowski at noaa.gov
>>>>>>>>>>>>>> Queue: met_help
>>>>>>>>>>>>>> Subject: AERONET AOD processing
>>>>>>>>>>>>>> Owner: Nobody
>>>>>>>>>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
>>>>>>>>>>>>>> Status: new
>>>>>>>>>>>>>> Ticket <URL:
>>>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Hello,
>>>>>>>>>>>>>> Is a formatting code available for version 3 AERONET? (
>>>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
>>>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
>>>>>>>>>>>>>> (If yes) Does the code allow for selecting chosen
wavelengths?
>>>>> Is
>>>>>>>>>>>> capable
>>>>>>>>>>>>>> of obtaining selected channels (e.g. 550 nm) from other
>>>>>>> wavelengths?
>>>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>>> Mariusz
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>
>>
>>
>>
>>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Minna Win
Time: Fri Mar 13 14:21:53 2020
Hi Mariusz,
Can you take a look at the comments for GitHub issue #1121:
https://github.com/NCAR/MET/issues/1121
It looks like it is *using* 440nm and 675 nm to *output* 550nm, but
you
also want to *output* 440nm and 675 nm in addition to 550nm-is that
correct
? My apologies for the confusion, I'm not familiar with ascii2nc so I
need
Howard to verify if I've interpreted the request correctly.
Thanks,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
On Fri, Mar 13, 2020 at 2:12 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>
> Minna,
> To clarify:
> 1. MET 9.0 can output 440 nm and 675 nm? (along with interpolated
550nm)
> 2. 870nm would be an addition - it would be fine only to output if
exists
> Mariusz
>
>
> > On Mar 13, 2020, at 2:08 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
> >
> > Hi Mariusz,
> >
> > I updated the GitHub Issue. Howard reminded me that MET 9.0 has
support
> > for 440nm and 675nm AERONET3. It looks like 870nm interpolation
support
> > will need to be added.
> >
> > Regards,
> > Minna
> > ---------------
> > Minna Win
> > National Center for Atmospheric Research
> > Developmental Testbed Center
> > Phone: 303-497-8423
> > Fax: 303-497-8401
> >
> >
> >
> > On Fri, Mar 13, 2020 at 2:05 PM Mariusz Pagowski via RT <
> met_help at ucar.edu>
> > wrote:
> >
> >>
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>
> >> Yes, thanks,
> >> Mariusz
> >>
> >>> On Mar 13, 2020, at 1:53 PM, Minna Win via RT
<met_help at ucar.edu>
> wrote:
> >>>
> >>> Hi Mariusz,
> >>>
> >>> I put down your affiliation as NOAA/GSD on the GitHub issue. Is
that
> >>> correct?
> >>>
> >>> Thanks,
> >>> Minna
> >>> ---------------
> >>> Minna Win
> >>> National Center for Atmospheric Research
> >>> Developmental Testbed Center
> >>> Phone: 303-497-8423
> >>> Fax: 303-497-8401
> >>>
> >>>
> >>>
> >>> On Fri, Mar 13, 2020 at 1:50 PM Minna Win-Gildenmeister <
> >> minnawin at ucar.edu>
> >>> wrote:
> >>>
> >>>> Hi Mariusz,
> >>>>
> >>>> I created a GitHub issue for MET:
> >>>> https://github.com/NCAR/MET/issues/1278
> >>>>
> >>>> Please feel free to add any additional comments to the GitHub
issue.
> >> I've
> >>>> assigned the issue to Howard Soh and assigned this with a high
> priority
> >>>> (Blocker).
> >>>>
> >>>> Regards,
> >>>> Minna
> >>>> ---------------
> >>>> Minna Win
> >>>> National Center for Atmospheric Research
> >>>> Developmental Testbed Center
> >>>> Phone: 303-497-8423
> >>>> Fax: 303-497-8401
> >>>>
> >>>>
> >>>>
> >>>> On Fri, Mar 13, 2020 at 12:05 PM Mariusz Pagowski via RT <
> >>>> met_help at ucar.edu> wrote:
> >>>>
> >>>>>
> >>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>
> >>>>>
> >>>>> Minna,
> >>>>> our time-frame is asap and would be happy to help with coding
if
> >> someone
> >>>>> provided guidance.
> >>>>> Pls let me know what your and John’s thinking is,
> >>>>> Mariusz
> >>>>>
> >>>>>
> >>>>>> On Mar 13, 2020, at 11:50 AM, Minna Win via RT
<met_help at ucar.edu>
> >>>>> wrote:
> >>>>>>
> >>>>>> Hi Mariusz,
> >>>>>>
> >>>>>> Let me see if this request can be submitted in the form of a
GitHub
> >>>>> issue.
> >>>>>> We're currently ready to release MET v9.0, so if accepted, it
would
> >> make
> >>>>>> its way into a future release. I'll check with John to see
if this
> >> can
> >>>>> be
> >>>>>> done within the current budget. What time-frame constraints
do you
> >>>>> have?
> >>>>>>
> >>>>>> Regards,
> >>>>>> Minna
> >>>>>> ---------------
> >>>>>> Minna Win
> >>>>>> National Center for Atmospheric Research
> >>>>>> Developmental Testbed Center
> >>>>>> Phone: 303-497-8423
> >>>>>> Fax: 303-497-8401
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>> On Fri, Mar 13, 2020 at 11:23 AM Mariusz Pagowski via RT <
> >>>>> met_help at ucar.edu>
> >>>>>> wrote:
> >>>>>>
> >>>>>>>
> >>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>>>>>>
> >>>>>>> Minna,
> >>>>>>> thanks for your dedication to answering these questions.
> >>>>>>> We are also interested in outputting 440 and 675 nm from
AERONET,
> and
> >>>>> even
> >>>>>>> possibly 870 nm.
> >>>>>>> We would be ready to work with someone from MET to help
implement
> >>>>> these as
> >>>>>>> options.
> >>>>>>> Pls let me know if that is possible,
> >>>>>>> Mariusz
> >>>>>>>
> >>>>>>>
> >>>>>>>> On Mar 13, 2020, at 8:03 AM, Minna Win via RT
<met_help at ucar.edu>
> >>>>> wrote:
> >>>>>>>>
> >>>>>>>> Hi Mariusz,
> >>>>>>>>
> >>>>>>>> Here's what happens when 675 and 440 nm are unavailable:
> >>>>>>>>
> >>>>>>>> ascii2nc does not save 550 nm wavelength into NetCDF output
if 675
> >> and
> >>>>>>> 440
> >>>>>>>>> are not available.
> >>>>>>>>> ascii2nc figures out the variable list and associated
wavelengths
> >>>>> from
> >>>>>>> the
> >>>>>>>>> header record and knows where AOD at 440 and 650 nm values
are.
> >>>>>>>>> ascii2nc saves AOD of all available
> >>>>>>>>> wavelengths. And ascii2nc checks AOD at 440, 550 and
> >>>>>>>>> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550
does
> not
> >>>>>>> exist
> >>>>>>>>> from data record and AOD values at 675 and at 440 are
valid.
> >>>>>>>>
> >>>>>>>>
> >>>>>>>> I hope this is useful.
> >>>>>>>>
> >>>>>>>> Regards,
> >>>>>>>> Minna
> >>>>>>>> ---------------
> >>>>>>>> Minna Win
> >>>>>>>> National Center for Atmospheric Research
> >>>>>>>> Developmental Testbed Center
> >>>>>>>> Phone: 303-497-8423
> >>>>>>>> Fax: 303-497-8401
> >>>>>>>>
> >>>>>>>>
> >>>>>>>>
> >>>>>>>> On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT <
> >>>>>>> met_help at ucar.edu>
> >>>>>>>> wrote:
> >>>>>>>>
> >>>>>>>>>
> >>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>>>>>>>>
> >>>>>>>>> Thanks Mina,
> >>>>>>>>> just wanted to be sure that it is not a linear
interpolation.
> >> Still I
> >>>>>>>>> wonder what happens when 675 and 440 are not available.
> >>>>>>>>> So if my understanding is correct MET only outputs 550nm.
> >>>>>>>>> Mariusz
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT <
> met_help at ucar.edu>
> >>>>>>> wrote:
> >>>>>>>>>>
> >>>>>>>>>> Hi Mariusz,
> >>>>>>>>>>
> >>>>>>>>>> Upon further checking from Howard, the implementation is
> actually:
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>> The Fortran code I referenced (
> >>>>> https://github.com/NCAR/MET/issues/1121
> >>>>>>> )
> >>>>>>>>> was
> >>>>>>>>>>> implemented that way:
> >>>>>>>>>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
> >>>>>>>>>>> aerdata(i)%aot_440nm)/log(675./440.)
> >>>>>>>>>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
> >>>>>>>>> ((550./440.)**(-a_675_440))
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>> Just in case you need to see the correct implementation
code.
> The
> >>>>>>>>>> first code snippet (in previous email) he cited was
replaced
> with
> >>>>> the
> >>>>>>>>>> above code snippet, which is associated with Github issue
> 1121...
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>> Regards,
> >>>>>>>>>>
> >>>>>>>>>> Minna
> >>>>>>>>>>
> >>>>>>>>>> ---------------
> >>>>>>>>>> Minna Win
> >>>>>>>>>> National Center for Atmospheric Research
> >>>>>>>>>> Developmental Testbed Center
> >>>>>>>>>> Phone: 303-497-8423
> >>>>>>>>>> Fax: 303-497-8401
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
> >>>>>>>>> minnawin at ucar.edu>
> >>>>>>>>>> wrote:
> >>>>>>>>>>
> >>>>>>>>>>> Hi Mariusz,
> >>>>>>>>>>>
> >>>>>>>>>>> According to Howard:
> >>>>>>>>>>>
> >>>>>>>>>>> It's a log interpolation: angstrom_power_interpolation
> >>>>>>>>>>>
> >>>>>>>>>>> double angstrom_power_interplation(double value_1,
double
> >> value_2,
> >>>>>>>>>>> double level_1, double level_2, double target_level) {
> >>>>>>>>>>> double angstrom_log =
> >>>>> -log10(value_1/value_2)/log10(level_1/level_2);
> >>>>>>>>>>> double angstrom_value = value_2 *
> >>>>>>>>>>> pow((target_level/level_2),-angstrom_log);
> >>>>>>>>>>> return angstrom_value;
> >>>>>>>>>>> }
> >>>>>>>>>>>
> >>>>>>>>>>> Regards,
> >>>>>>>>>>> Minna
> >>>>>>>>>>> ---------------
> >>>>>>>>>>> Minna Win
> >>>>>>>>>>> National Center for Atmospheric Research
> >>>>>>>>>>> Developmental Testbed Center
> >>>>>>>>>>> Phone: 303-497-8423
> >>>>>>>>>>> Fax: 303-497-8401
> >>>>>>>>>>>
> >>>>>>>>>>>
> >>>>>>>>>>>
> >>>>>>>>>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT
<
> >>>>>>>>> met_help at ucar.edu>
> >>>>>>>>>>> wrote:
> >>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
> >
> >>>>>>>>>>>>
> >>>>>>>>>>>> Minna,
> >>>>>>>>>>>> thanks very much - just one follow-up: how is the
> interpolation
> >>>>> done
> >>>>>>>>>>>> (linear/log)?
> >>>>>>>>>>>> Mariusz
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT <
> >> met_help at ucar.edu
> >>>>>>
> >>>>>>>>>>>> wrote:
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> Hi Mariusz,
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> According to our resident expert Howard Soh:
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> "ascii2nc reads a text format AERONET 3 data and
produces a
> >>>>> NetCDF
> >>>>>>>>>>>> output
> >>>>>>>>>>>>> (MET's point observation data file). ascii2nc
interpolates
> >> 550nm
> >>>>>>>>>>>> wavelength
> >>>>>>>>>>>>> from 440 nm and 675nm if not exists."
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> I hope this helps.
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> Regards,
> >>>>>>>>>>>>> Minna
> >>>>>>>>>>>>> ---------------
> >>>>>>>>>>>>> Minna Win
> >>>>>>>>>>>>> National Center for Atmospheric Research
> >>>>>>>>>>>>> Developmental Testbed Center
> >>>>>>>>>>>>> Phone: 303-497-8423
> >>>>>>>>>>>>> Fax: 303-497-8401
> >>>>>>>>>>>>>
> >>>>>>>>>>>>>
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via
RT <
> >>>>>>>>>>>> met_help at ucar.edu>
> >>>>>>>>>>>>> wrote:
> >>>>>>>>>>>>>
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted
upon.
> >>>>>>>>>>>>>> Transaction: Ticket created by
Mariusz.Pagowski at noaa.gov
> >>>>>>>>>>>>>> Queue: met_help
> >>>>>>>>>>>>>> Subject: AERONET AOD processing
> >>>>>>>>>>>>>> Owner: Nobody
> >>>>>>>>>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
> >>>>>>>>>>>>>> Status: new
> >>>>>>>>>>>>>> Ticket <URL:
> >>>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
> >>>>>>>>>>
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> Hello,
> >>>>>>>>>>>>>> Is a formatting code available for version 3 AERONET?
(
> >>>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
> >>>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
> >>>>>>>>>>>>>> (If yes) Does the code allow for selecting chosen
> wavelengths?
> >>>>> Is
> >>>>>>>>>>>> capable
> >>>>>>>>>>>>>> of obtaining selected channels (e.g. 550 nm) from
other
> >>>>>>> wavelengths?
> >>>>>>>>>>>>>> Thanks,
> >>>>>>>>>>>>>> Mariusz
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>
> >>
> >>
> >>
> >>
> >
>
>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #94517] AERONET AOD processing
From: Mariusz Pagowski
Time: Fri Mar 13 14:25:49 2020
Yes, from GitHub it looks like MET only outputs 550 nm interpolated
from 440 and 675 nms.
It would be good for us to output AOD at 440, 550, and 675 nms (and
870 nm in the future),
Mariusz
> On Mar 13, 2020, at 2:21 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>
> Hi Mariusz,
>
> Can you take a look at the comments for GitHub issue #1121:
>
> https://github.com/NCAR/MET/issues/1121
>
> It looks like it is *using* 440nm and 675 nm to *output* 550nm, but
you
> also want to *output* 440nm and 675 nm in addition to 550nm-is that
correct
> ? My apologies for the confusion, I'm not familiar with ascii2nc so
I need
> Howard to verify if I've interpreted the request correctly.
>
> Thanks,
> Minna
> ---------------
> Minna Win
> National Center for Atmospheric Research
> Developmental Testbed Center
> Phone: 303-497-8423
> Fax: 303-497-8401
>
>
>
> On Fri, Mar 13, 2020 at 2:12 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>
>> Minna,
>> To clarify:
>> 1. MET 9.0 can output 440 nm and 675 nm? (along with interpolated
550nm)
>> 2. 870nm would be an addition - it would be fine only to output if
exists
>> Mariusz
>>
>>
>>> On Mar 13, 2020, at 2:08 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>>>
>>> Hi Mariusz,
>>>
>>> I updated the GitHub Issue. Howard reminded me that MET 9.0 has
support
>>> for 440nm and 675nm AERONET3. It looks like 870nm interpolation
support
>>> will need to be added.
>>>
>>> Regards,
>>> Minna
>>> ---------------
>>> Minna Win
>>> National Center for Atmospheric Research
>>> Developmental Testbed Center
>>> Phone: 303-497-8423
>>> Fax: 303-497-8401
>>>
>>>
>>>
>>> On Fri, Mar 13, 2020 at 2:05 PM Mariusz Pagowski via RT <
>> met_help at ucar.edu>
>>> wrote:
>>>
>>>>
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>
>>>> Yes, thanks,
>>>> Mariusz
>>>>
>>>>> On Mar 13, 2020, at 1:53 PM, Minna Win via RT
<met_help at ucar.edu>
>> wrote:
>>>>>
>>>>> Hi Mariusz,
>>>>>
>>>>> I put down your affiliation as NOAA/GSD on the GitHub issue. Is
that
>>>>> correct?
>>>>>
>>>>> Thanks,
>>>>> Minna
>>>>> ---------------
>>>>> Minna Win
>>>>> National Center for Atmospheric Research
>>>>> Developmental Testbed Center
>>>>> Phone: 303-497-8423
>>>>> Fax: 303-497-8401
>>>>>
>>>>>
>>>>>
>>>>> On Fri, Mar 13, 2020 at 1:50 PM Minna Win-Gildenmeister <
>>>> minnawin at ucar.edu>
>>>>> wrote:
>>>>>
>>>>>> Hi Mariusz,
>>>>>>
>>>>>> I created a GitHub issue for MET:
>>>>>> https://github.com/NCAR/MET/issues/1278
>>>>>>
>>>>>> Please feel free to add any additional comments to the GitHub
issue.
>>>> I've
>>>>>> assigned the issue to Howard Soh and assigned this with a high
>> priority
>>>>>> (Blocker).
>>>>>>
>>>>>> Regards,
>>>>>> Minna
>>>>>> ---------------
>>>>>> Minna Win
>>>>>> National Center for Atmospheric Research
>>>>>> Developmental Testbed Center
>>>>>> Phone: 303-497-8423
>>>>>> Fax: 303-497-8401
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Fri, Mar 13, 2020 at 12:05 PM Mariusz Pagowski via RT <
>>>>>> met_help at ucar.edu> wrote:
>>>>>>
>>>>>>>
>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>
>>>>>>>
>>>>>>> Minna,
>>>>>>> our time-frame is asap and would be happy to help with coding
if
>>>> someone
>>>>>>> provided guidance.
>>>>>>> Pls let me know what your and John’s thinking is,
>>>>>>> Mariusz
>>>>>>>
>>>>>>>
>>>>>>>> On Mar 13, 2020, at 11:50 AM, Minna Win via RT
<met_help at ucar.edu>
>>>>>>> wrote:
>>>>>>>>
>>>>>>>> Hi Mariusz,
>>>>>>>>
>>>>>>>> Let me see if this request can be submitted in the form of a
GitHub
>>>>>>> issue.
>>>>>>>> We're currently ready to release MET v9.0, so if accepted, it
would
>>>> make
>>>>>>>> its way into a future release. I'll check with John to see
if this
>>>> can
>>>>>>> be
>>>>>>>> done within the current budget. What time-frame constraints
do you
>>>>>>> have?
>>>>>>>>
>>>>>>>> Regards,
>>>>>>>> Minna
>>>>>>>> ---------------
>>>>>>>> Minna Win
>>>>>>>> National Center for Atmospheric Research
>>>>>>>> Developmental Testbed Center
>>>>>>>> Phone: 303-497-8423
>>>>>>>> Fax: 303-497-8401
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> On Fri, Mar 13, 2020 at 11:23 AM Mariusz Pagowski via RT <
>>>>>>> met_help at ucar.edu>
>>>>>>>> wrote:
>>>>>>>>
>>>>>>>>>
>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>>>>>
>>>>>>>>> Minna,
>>>>>>>>> thanks for your dedication to answering these questions.
>>>>>>>>> We are also interested in outputting 440 and 675 nm from
AERONET,
>> and
>>>>>>> even
>>>>>>>>> possibly 870 nm.
>>>>>>>>> We would be ready to work with someone from MET to help
implement
>>>>>>> these as
>>>>>>>>> options.
>>>>>>>>> Pls let me know if that is possible,
>>>>>>>>> Mariusz
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>> On Mar 13, 2020, at 8:03 AM, Minna Win via RT
<met_help at ucar.edu>
>>>>>>> wrote:
>>>>>>>>>>
>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>
>>>>>>>>>> Here's what happens when 675 and 440 nm are unavailable:
>>>>>>>>>>
>>>>>>>>>> ascii2nc does not save 550 nm wavelength into NetCDF output
if 675
>>>> and
>>>>>>>>> 440
>>>>>>>>>>> are not available.
>>>>>>>>>>> ascii2nc figures out the variable list and associated
wavelengths
>>>>>>> from
>>>>>>>>> the
>>>>>>>>>>> header record and knows where AOD at 440 and 650 nm values
are.
>>>>>>>>>>> ascii2nc saves AOD of all available
>>>>>>>>>>> wavelengths. And ascii2nc checks AOD at 440, 550 and
>>>>>>>>>>> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550
does
>> not
>>>>>>>>> exist
>>>>>>>>>>> from data record and AOD values at 675 and at 440 are
valid.
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> I hope this is useful.
>>>>>>>>>>
>>>>>>>>>> Regards,
>>>>>>>>>> Minna
>>>>>>>>>> ---------------
>>>>>>>>>> Minna Win
>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT <
>>>>>>>>> met_help at ucar.edu>
>>>>>>>>>> wrote:
>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>>>>>>>
>>>>>>>>>>> Thanks Mina,
>>>>>>>>>>> just wanted to be sure that it is not a linear
interpolation.
>>>> Still I
>>>>>>>>>>> wonder what happens when 675 and 440 are not available.
>>>>>>>>>>> So if my understanding is correct MET only outputs 550nm.
>>>>>>>>>>> Mariusz
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT <
>> met_help at ucar.edu>
>>>>>>>>> wrote:
>>>>>>>>>>>>
>>>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>>>
>>>>>>>>>>>> Upon further checking from Howard, the implementation is
>> actually:
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> The Fortran code I referenced (
>>>>>>> https://github.com/NCAR/MET/issues/1121
>>>>>>>>> )
>>>>>>>>>>> was
>>>>>>>>>>>>> implemented that way:
>>>>>>>>>>>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
>>>>>>>>>>>>> aerdata(i)%aot_440nm)/log(675./440.)
>>>>>>>>>>>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
>>>>>>>>>>> ((550./440.)**(-a_675_440))
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> Just in case you need to see the correct implementation
code.
>> The
>>>>>>>>>>>> first code snippet (in previous email) he cited was
replaced
>> with
>>>>>>> the
>>>>>>>>>>>> above code snippet, which is associated with Github issue
>> 1121...
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> Regards,
>>>>>>>>>>>>
>>>>>>>>>>>> Minna
>>>>>>>>>>>>
>>>>>>>>>>>> ---------------
>>>>>>>>>>>> Minna Win
>>>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
>>>>>>>>>>> minnawin at ucar.edu>
>>>>>>>>>>>> wrote:
>>>>>>>>>>>>
>>>>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>>>>
>>>>>>>>>>>>> According to Howard:
>>>>>>>>>>>>>
>>>>>>>>>>>>> It's a log interpolation: angstrom_power_interpolation
>>>>>>>>>>>>>
>>>>>>>>>>>>> double angstrom_power_interplation(double value_1,
double
>>>> value_2,
>>>>>>>>>>>>> double level_1, double level_2, double target_level) {
>>>>>>>>>>>>> double angstrom_log =
>>>>>>> -log10(value_1/value_2)/log10(level_1/level_2);
>>>>>>>>>>>>> double angstrom_value = value_2 *
>>>>>>>>>>>>> pow((target_level/level_2),-angstrom_log);
>>>>>>>>>>>>> return angstrom_value;
>>>>>>>>>>>>> }
>>>>>>>>>>>>>
>>>>>>>>>>>>> Regards,
>>>>>>>>>>>>> Minna
>>>>>>>>>>>>> ---------------
>>>>>>>>>>>>> Minna Win
>>>>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT
<
>>>>>>>>>>> met_help at ucar.edu>
>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Minna,
>>>>>>>>>>>>>> thanks very much - just one follow-up: how is the
>> interpolation
>>>>>>> done
>>>>>>>>>>>>>> (linear/log)?
>>>>>>>>>>>>>> Mariusz
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT <
>>>> met_help at ucar.edu
>>>>>>>>
>>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> According to our resident expert Howard Soh:
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> "ascii2nc reads a text format AERONET 3 data and
produces a
>>>>>>> NetCDF
>>>>>>>>>>>>>> output
>>>>>>>>>>>>>>> (MET's point observation data file). ascii2nc
interpolates
>>>> 550nm
>>>>>>>>>>>>>> wavelength
>>>>>>>>>>>>>>> from 440 nm and 675nm if not exists."
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> I hope this helps.
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Regards,
>>>>>>>>>>>>>>> Minna
>>>>>>>>>>>>>>> ---------------
>>>>>>>>>>>>>>> Minna Win
>>>>>>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via
RT <
>>>>>>>>>>>>>> met_help at ucar.edu>
>>>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted
upon.
>>>>>>>>>>>>>>>> Transaction: Ticket created by
Mariusz.Pagowski at noaa.gov
>>>>>>>>>>>>>>>> Queue: met_help
>>>>>>>>>>>>>>>> Subject: AERONET AOD processing
>>>>>>>>>>>>>>>> Owner: Nobody
>>>>>>>>>>>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
>>>>>>>>>>>>>>>> Status: new
>>>>>>>>>>>>>>>> Ticket <URL:
>>>>>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Hello,
>>>>>>>>>>>>>>>> Is a formatting code available for version 3 AERONET?
(
>>>>>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
>>>>>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
>>>>>>>>>>>>>>>> (If yes) Does the code allow for selecting chosen
>> wavelengths?
>>>>>>> Is
>>>>>>>>>>>>>> capable
>>>>>>>>>>>>>>>> of obtaining selected channels (e.g. 550 nm) from
other
>>>>>>>>> wavelengths?
>>>>>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>>>>> Mariusz
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>
>>>>
>>>>
>>>>
>>>>
>>>
>>
>>
>>
>>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Minna Win
Time: Fri Mar 13 14:42:27 2020
Hi Mariusz,
Thanks, I'll update the GitHub issue to reflect what you stated.
Regards,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
On Fri, Mar 13, 2020 at 2:35 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>
> Yes, from GitHub it looks like MET only outputs 550 nm interpolated
from
> 440 and 675 nms.
> It would be good for us to output AOD at 440, 550, and 675 nms (and
870 nm
> in the future),
> Mariusz
>
>
> > On Mar 13, 2020, at 2:21 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
> >
> > Hi Mariusz,
> >
> > Can you take a look at the comments for GitHub issue #1121:
> >
> > https://github.com/NCAR/MET/issues/1121
> >
> > It looks like it is *using* 440nm and 675 nm to *output* 550nm,
but you
> > also want to *output* 440nm and 675 nm in addition to 550nm-is
that
> correct
> > ? My apologies for the confusion, I'm not familiar with ascii2nc
so I
> need
> > Howard to verify if I've interpreted the request correctly.
> >
> > Thanks,
> > Minna
> > ---------------
> > Minna Win
> > National Center for Atmospheric Research
> > Developmental Testbed Center
> > Phone: 303-497-8423
> > Fax: 303-497-8401
> >
> >
> >
> > On Fri, Mar 13, 2020 at 2:12 PM Mariusz Pagowski via RT <
> met_help at ucar.edu>
> > wrote:
> >
> >>
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>
> >> Minna,
> >> To clarify:
> >> 1. MET 9.0 can output 440 nm and 675 nm? (along with interpolated
550nm)
> >> 2. 870nm would be an addition - it would be fine only to output
if
> exists
> >> Mariusz
> >>
> >>
> >>> On Mar 13, 2020, at 2:08 PM, Minna Win via RT
<met_help at ucar.edu>
> wrote:
> >>>
> >>> Hi Mariusz,
> >>>
> >>> I updated the GitHub Issue. Howard reminded me that MET 9.0 has
> support
> >>> for 440nm and 675nm AERONET3. It looks like 870nm interpolation
> support
> >>> will need to be added.
> >>>
> >>> Regards,
> >>> Minna
> >>> ---------------
> >>> Minna Win
> >>> National Center for Atmospheric Research
> >>> Developmental Testbed Center
> >>> Phone: 303-497-8423
> >>> Fax: 303-497-8401
> >>>
> >>>
> >>>
> >>> On Fri, Mar 13, 2020 at 2:05 PM Mariusz Pagowski via RT <
> >> met_help at ucar.edu>
> >>> wrote:
> >>>
> >>>>
> >>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>>>
> >>>> Yes, thanks,
> >>>> Mariusz
> >>>>
> >>>>> On Mar 13, 2020, at 1:53 PM, Minna Win via RT
<met_help at ucar.edu>
> >> wrote:
> >>>>>
> >>>>> Hi Mariusz,
> >>>>>
> >>>>> I put down your affiliation as NOAA/GSD on the GitHub issue.
Is that
> >>>>> correct?
> >>>>>
> >>>>> Thanks,
> >>>>> Minna
> >>>>> ---------------
> >>>>> Minna Win
> >>>>> National Center for Atmospheric Research
> >>>>> Developmental Testbed Center
> >>>>> Phone: 303-497-8423
> >>>>> Fax: 303-497-8401
> >>>>>
> >>>>>
> >>>>>
> >>>>> On Fri, Mar 13, 2020 at 1:50 PM Minna Win-Gildenmeister <
> >>>> minnawin at ucar.edu>
> >>>>> wrote:
> >>>>>
> >>>>>> Hi Mariusz,
> >>>>>>
> >>>>>> I created a GitHub issue for MET:
> >>>>>> https://github.com/NCAR/MET/issues/1278
> >>>>>>
> >>>>>> Please feel free to add any additional comments to the GitHub
issue.
> >>>> I've
> >>>>>> assigned the issue to Howard Soh and assigned this with a
high
> >> priority
> >>>>>> (Blocker).
> >>>>>>
> >>>>>> Regards,
> >>>>>> Minna
> >>>>>> ---------------
> >>>>>> Minna Win
> >>>>>> National Center for Atmospheric Research
> >>>>>> Developmental Testbed Center
> >>>>>> Phone: 303-497-8423
> >>>>>> Fax: 303-497-8401
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>> On Fri, Mar 13, 2020 at 12:05 PM Mariusz Pagowski via RT <
> >>>>>> met_help at ucar.edu> wrote:
> >>>>>>
> >>>>>>>
> >>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>>>>>>
> >>>>>>> Minna,
> >>>>>>> our time-frame is asap and would be happy to help with
coding if
> >>>> someone
> >>>>>>> provided guidance.
> >>>>>>> Pls let me know what your and John’s thinking is,
> >>>>>>> Mariusz
> >>>>>>>
> >>>>>>>
> >>>>>>>> On Mar 13, 2020, at 11:50 AM, Minna Win via RT
<met_help at ucar.edu
> >
> >>>>>>> wrote:
> >>>>>>>>
> >>>>>>>> Hi Mariusz,
> >>>>>>>>
> >>>>>>>> Let me see if this request can be submitted in the form of
a
> GitHub
> >>>>>>> issue.
> >>>>>>>> We're currently ready to release MET v9.0, so if accepted,
it
> would
> >>>> make
> >>>>>>>> its way into a future release. I'll check with John to see
if
> this
> >>>> can
> >>>>>>> be
> >>>>>>>> done within the current budget. What time-frame
constraints do
> you
> >>>>>>> have?
> >>>>>>>>
> >>>>>>>> Regards,
> >>>>>>>> Minna
> >>>>>>>> ---------------
> >>>>>>>> Minna Win
> >>>>>>>> National Center for Atmospheric Research
> >>>>>>>> Developmental Testbed Center
> >>>>>>>> Phone: 303-497-8423
> >>>>>>>> Fax: 303-497-8401
> >>>>>>>>
> >>>>>>>>
> >>>>>>>>
> >>>>>>>> On Fri, Mar 13, 2020 at 11:23 AM Mariusz Pagowski via RT <
> >>>>>>> met_help at ucar.edu>
> >>>>>>>> wrote:
> >>>>>>>>
> >>>>>>>>>
> >>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
> >>>>>>>>>
> >>>>>>>>> Minna,
> >>>>>>>>> thanks for your dedication to answering these questions.
> >>>>>>>>> We are also interested in outputting 440 and 675 nm from
AERONET,
> >> and
> >>>>>>> even
> >>>>>>>>> possibly 870 nm.
> >>>>>>>>> We would be ready to work with someone from MET to help
implement
> >>>>>>> these as
> >>>>>>>>> options.
> >>>>>>>>> Pls let me know if that is possible,
> >>>>>>>>> Mariusz
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>> On Mar 13, 2020, at 8:03 AM, Minna Win via RT <
> met_help at ucar.edu>
> >>>>>>> wrote:
> >>>>>>>>>>
> >>>>>>>>>> Hi Mariusz,
> >>>>>>>>>>
> >>>>>>>>>> Here's what happens when 675 and 440 nm are unavailable:
> >>>>>>>>>>
> >>>>>>>>>> ascii2nc does not save 550 nm wavelength into NetCDF
output if
> 675
> >>>> and
> >>>>>>>>> 440
> >>>>>>>>>>> are not available.
> >>>>>>>>>>> ascii2nc figures out the variable list and associated
> wavelengths
> >>>>>>> from
> >>>>>>>>> the
> >>>>>>>>>>> header record and knows where AOD at 440 and 650 nm
values are.
> >>>>>>>>>>> ascii2nc saves AOD of all available
> >>>>>>>>>>> wavelengths. And ascii2nc checks AOD at 440, 550 and
> >>>>>>>>>>> 675. ascii2nc computes & saves AOD at 550 if If AOD at
550 does
> >> not
> >>>>>>>>> exist
> >>>>>>>>>>> from data record and AOD values at 675 and at 440 are
valid.
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>> I hope this is useful.
> >>>>>>>>>>
> >>>>>>>>>> Regards,
> >>>>>>>>>> Minna
> >>>>>>>>>> ---------------
> >>>>>>>>>> Minna Win
> >>>>>>>>>> National Center for Atmospheric Research
> >>>>>>>>>> Developmental Testbed Center
> >>>>>>>>>> Phone: 303-497-8423
> >>>>>>>>>> Fax: 303-497-8401
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>> On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT <
> >>>>>>>>> met_help at ucar.edu>
> >>>>>>>>>> wrote:
> >>>>>>>>>>
> >>>>>>>>>>>
> >>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
> >
> >>>>>>>>>>>
> >>>>>>>>>>> Thanks Mina,
> >>>>>>>>>>> just wanted to be sure that it is not a linear
interpolation.
> >>>> Still I
> >>>>>>>>>>> wonder what happens when 675 and 440 are not available.
> >>>>>>>>>>> So if my understanding is correct MET only outputs
550nm.
> >>>>>>>>>>> Mariusz
> >>>>>>>>>>>
> >>>>>>>>>>>
> >>>>>>>>>>>
> >>>>>>>>>>>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT <
> >> met_help at ucar.edu>
> >>>>>>>>> wrote:
> >>>>>>>>>>>>
> >>>>>>>>>>>> Hi Mariusz,
> >>>>>>>>>>>>
> >>>>>>>>>>>> Upon further checking from Howard, the implementation
is
> >> actually:
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>> The Fortran code I referenced (
> >>>>>>> https://github.com/NCAR/MET/issues/1121
> >>>>>>>>> )
> >>>>>>>>>>> was
> >>>>>>>>>>>>> implemented that way:
> >>>>>>>>>>>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
> >>>>>>>>>>>>> aerdata(i)%aot_440nm)/log(675./440.)
> >>>>>>>>>>>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
> >>>>>>>>>>> ((550./440.)**(-a_675_440))
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>> Just in case you need to see the correct implementation
code.
> >> The
> >>>>>>>>>>>> first code snippet (in previous email) he cited was
replaced
> >> with
> >>>>>>> the
> >>>>>>>>>>>> above code snippet, which is associated with Github
issue
> >> 1121...
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>> Regards,
> >>>>>>>>>>>>
> >>>>>>>>>>>> Minna
> >>>>>>>>>>>>
> >>>>>>>>>>>> ---------------
> >>>>>>>>>>>> Minna Win
> >>>>>>>>>>>> National Center for Atmospheric Research
> >>>>>>>>>>>> Developmental Testbed Center
> >>>>>>>>>>>> Phone: 303-497-8423
> >>>>>>>>>>>> Fax: 303-497-8401
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister
<
> >>>>>>>>>>> minnawin at ucar.edu>
> >>>>>>>>>>>> wrote:
> >>>>>>>>>>>>
> >>>>>>>>>>>>> Hi Mariusz,
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> According to Howard:
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> It's a log interpolation: angstrom_power_interpolation
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> double angstrom_power_interplation(double value_1,
double
> >>>> value_2,
> >>>>>>>>>>>>> double level_1, double level_2, double target_level) {
> >>>>>>>>>>>>> double angstrom_log =
> >>>>>>> -log10(value_1/value_2)/log10(level_1/level_2);
> >>>>>>>>>>>>> double angstrom_value = value_2 *
> >>>>>>>>>>>>> pow((target_level/level_2),-angstrom_log);
> >>>>>>>>>>>>> return angstrom_value;
> >>>>>>>>>>>>> }
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> Regards,
> >>>>>>>>>>>>> Minna
> >>>>>>>>>>>>> ---------------
> >>>>>>>>>>>>> Minna Win
> >>>>>>>>>>>>> National Center for Atmospheric Research
> >>>>>>>>>>>>> Developmental Testbed Center
> >>>>>>>>>>>>> Phone: 303-497-8423
> >>>>>>>>>>>>> Fax: 303-497-8401
> >>>>>>>>>>>>>
> >>>>>>>>>>>>>
> >>>>>>>>>>>>>
> >>>>>>>>>>>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via
RT <
> >>>>>>>>>>> met_help at ucar.edu>
> >>>>>>>>>>>>> wrote:
> >>>>>>>>>>>>>
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
> >>>
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>> Minna,
> >>>>>>>>>>>>>> thanks very much - just one follow-up: how is the
> >> interpolation
> >>>>>>> done
> >>>>>>>>>>>>>> (linear/log)?
> >>>>>>>>>>>>>> Mariusz
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT <
> >>>> met_help at ucar.edu
> >>>>>>>>
> >>>>>>>>>>>>>> wrote:
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> Hi Mariusz,
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> According to our resident expert Howard Soh:
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> "ascii2nc reads a text format AERONET 3 data and
produces a
> >>>>>>> NetCDF
> >>>>>>>>>>>>>> output
> >>>>>>>>>>>>>>> (MET's point observation data file). ascii2nc
interpolates
> >>>> 550nm
> >>>>>>>>>>>>>> wavelength
> >>>>>>>>>>>>>>> from 440 nm and 675nm if not exists."
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> I hope this helps.
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> Regards,
> >>>>>>>>>>>>>>> Minna
> >>>>>>>>>>>>>>> ---------------
> >>>>>>>>>>>>>>> Minna Win
> >>>>>>>>>>>>>>> National Center for Atmospheric Research
> >>>>>>>>>>>>>>> Developmental Testbed Center
> >>>>>>>>>>>>>>> Phone: 303-497-8423
> >>>>>>>>>>>>>>> Fax: 303-497-8401
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via
RT <
> >>>>>>>>>>>>>> met_help at ucar.edu>
> >>>>>>>>>>>>>>> wrote:
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted
upon.
> >>>>>>>>>>>>>>>> Transaction: Ticket created by
Mariusz.Pagowski at noaa.gov
> >>>>>>>>>>>>>>>> Queue: met_help
> >>>>>>>>>>>>>>>> Subject: AERONET AOD processing
> >>>>>>>>>>>>>>>> Owner: Nobody
> >>>>>>>>>>>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
> >>>>>>>>>>>>>>>> Status: new
> >>>>>>>>>>>>>>>> Ticket <URL:
> >>>>>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
> >>>>>>>>>>>>
> >>>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>> Hello,
> >>>>>>>>>>>>>>>> Is a formatting code available for version 3
AERONET? (
> >>>>>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
> >>>>>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
> >>>>>>>>>>>>>>>> (If yes) Does the code allow for selecting chosen
> >> wavelengths?
> >>>>>>> Is
> >>>>>>>>>>>>>> capable
> >>>>>>>>>>>>>>>> of obtaining selected channels (e.g. 550 nm) from
other
> >>>>>>>>> wavelengths?
> >>>>>>>>>>>>>>>> Thanks,
> >>>>>>>>>>>>>>>> Mariusz
> >>>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>>
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>>>
> >>>>>>>>>>>>
> >>>>>>>>>>>
> >>>>>>>>>>>
> >>>>>>>>>>>
> >>>>>>>>>>>
> >>>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>>
> >>>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>>>
> >>>>>
> >>>>
> >>>>
> >>>>
> >>>>
> >>>
> >>
> >>
> >>
> >>
> >
>
>
>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Howard Soh
Time: Fri Mar 13 15:05:07 2020
We need which wavelengths for interpolating for AOD_870nm.
The sample AERONET3 file has valid values for following wavelengths:
Date(dd:mm:yyyy) Time(hh:mm:ss) Day_of_Year AOD_1640nm AOD_1020nm
AOD_870nm AOD_675nm AOD_500nm AOD_440nm AOD_380nm AOD_340nm
0.779572 12:00:00 261 0.022949 0.031322 0.038461 0.049653
0.079107 0.097381 0.128962 0.139122
0.821238 12:00:00 262 0.01301 0.019142 0.022565 0.02795
0.041882
0.051215 0.063218 0.068337
0.862905 12:00:00 263 0.011324 0.016164 0.0199 0.023205
0.032876
0.039101 0.051348 0.054325
0.904572 12:00:00 264 0.010447 0.014793 0.018324 0.019833
0.026401 0.030903 0.043605 0.041629
0.987905 12:00:00 266 0.034807 0.046698 0.053709 0.06755
0.100127
0.118879 0.148089 0.164597
24:09:2015 12:00:00 267 0.02107 0.02921 0.034945 0.043239
0.065506 0.079168 0.104625 0.112846
25:09:2015 12:00:00 268 0.021043 0.028205 0.03332 0.040212
0.060867 0.074257 0.095984 0.105028
26:09:2015 12:00:00 269 0.015853 0.022548 0.027219 0.032956
0.051153 0.063412 0.084131 0.089649
They are saved into ascii2nc output as is (AOD for 8 wavelengths).
If there are missing values at AOD_675nm column, interpolated value is
saved.
Cheers,
Howard
On 3/13/2020 2:12 PM, Mariusz Pagowski via RT wrote:
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>
> Minna,
> To clarify:
> 1. MET 9.0 can output 440 nm and 675 nm? (along with interpolated
550nm)
> 2. 870nm would be an addition - it would be fine only to output if
exists
> Mariusz
>
>
>> On Mar 13, 2020, at 2:08 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>>
>> Hi Mariusz,
>>
>> I updated the GitHub Issue. Howard reminded me that MET 9.0 has
support
>> for 440nm and 675nm AERONET3. It looks like 870nm interpolation
support
>> will need to be added.
>>
>> Regards,
>> Minna
>> ---------------
>> Minna Win
>> National Center for Atmospheric Research
>> Developmental Testbed Center
>> Phone: 303-497-8423
>> Fax: 303-497-8401
>>
>>
>>
>> On Fri, Mar 13, 2020 at 2:05 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
>> wrote:
>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>
>>> Yes, thanks,
>>> Mariusz
>>>
>>>> On Mar 13, 2020, at 1:53 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>>>>
>>>> Hi Mariusz,
>>>>
>>>> I put down your affiliation as NOAA/GSD on the GitHub issue. Is
that
>>>> correct?
>>>>
>>>> Thanks,
>>>> Minna
>>>> ---------------
>>>> Minna Win
>>>> National Center for Atmospheric Research
>>>> Developmental Testbed Center
>>>> Phone: 303-497-8423
>>>> Fax: 303-497-8401
>>>>
>>>>
>>>>
>>>> On Fri, Mar 13, 2020 at 1:50 PM Minna Win-Gildenmeister <
>>> minnawin at ucar.edu>
>>>> wrote:
>>>>
>>>>> Hi Mariusz,
>>>>>
>>>>> I created a GitHub issue for MET:
>>>>> https://github.com/NCAR/MET/issues/1278
>>>>>
>>>>> Please feel free to add any additional comments to the GitHub
issue.
>>> I've
>>>>> assigned the issue to Howard Soh and assigned this with a high
priority
>>>>> (Blocker).
>>>>>
>>>>> Regards,
>>>>> Minna
>>>>> ---------------
>>>>> Minna Win
>>>>> National Center for Atmospheric Research
>>>>> Developmental Testbed Center
>>>>> Phone: 303-497-8423
>>>>> Fax: 303-497-8401
>>>>>
>>>>>
>>>>>
>>>>> On Fri, Mar 13, 2020 at 12:05 PM Mariusz Pagowski via RT <
>>>>> met_help at ucar.edu> wrote:
>>>>>
>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>>
>>>>>> Minna,
>>>>>> our time-frame is asap and would be happy to help with coding
if
>>> someone
>>>>>> provided guidance.
>>>>>> Pls let me know what your and John’s thinking is,
>>>>>> Mariusz
>>>>>>
>>>>>>
>>>>>>> On Mar 13, 2020, at 11:50 AM, Minna Win via RT
<met_help at ucar.edu>
>>>>>> wrote:
>>>>>>> Hi Mariusz,
>>>>>>>
>>>>>>> Let me see if this request can be submitted in the form of a
GitHub
>>>>>> issue.
>>>>>>> We're currently ready to release MET v9.0, so if accepted, it
would
>>> make
>>>>>>> its way into a future release. I'll check with John to see if
this
>>> can
>>>>>> be
>>>>>>> done within the current budget. What time-frame constraints
do you
>>>>>> have?
>>>>>>> Regards,
>>>>>>> Minna
>>>>>>> ---------------
>>>>>>> Minna Win
>>>>>>> National Center for Atmospheric Research
>>>>>>> Developmental Testbed Center
>>>>>>> Phone: 303-497-8423
>>>>>>> Fax: 303-497-8401
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> On Fri, Mar 13, 2020 at 11:23 AM Mariusz Pagowski via RT <
>>>>>> met_help at ucar.edu>
>>>>>>> wrote:
>>>>>>>
>>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>
>>>>>>>>
>>>>>>>> Minna,
>>>>>>>> thanks for your dedication to answering these questions.
>>>>>>>> We are also interested in outputting 440 and 675 nm from
AERONET, and
>>>>>> even
>>>>>>>> possibly 870 nm.
>>>>>>>> We would be ready to work with someone from MET to help
implement
>>>>>> these as
>>>>>>>> options.
>>>>>>>> Pls let me know if that is possible,
>>>>>>>> Mariusz
>>>>>>>>
>>>>>>>>
>>>>>>>>> On Mar 13, 2020, at 8:03 AM, Minna Win via RT
<met_help at ucar.edu>
>>>>>> wrote:
>>>>>>>>> Hi Mariusz,
>>>>>>>>>
>>>>>>>>> Here's what happens when 675 and 440 nm are unavailable:
>>>>>>>>>
>>>>>>>>> ascii2nc does not save 550 nm wavelength into NetCDF output
if 675
>>> and
>>>>>>>> 440
>>>>>>>>>> are not available.
>>>>>>>>>> ascii2nc figures out the variable list and associated
wavelengths
>>>>>> from
>>>>>>>> the
>>>>>>>>>> header record and knows where AOD at 440 and 650 nm values
are.
>>>>>>>>>> ascii2nc saves AOD of all available
>>>>>>>>>> wavelengths. And ascii2nc checks AOD at 440, 550 and
>>>>>>>>>> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550
does not
>>>>>>>> exist
>>>>>>>>>> from data record and AOD values at 675 and at 440 are
valid.
>>>>>>>>>
>>>>>>>>> I hope this is useful.
>>>>>>>>>
>>>>>>>>> Regards,
>>>>>>>>> Minna
>>>>>>>>> ---------------
>>>>>>>>> Minna Win
>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>> Developmental Testbed Center
>>>>>>>>> Phone: 303-497-8423
>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT <
>>>>>>>> met_help at ucar.edu>
>>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>>>>>>
>>>>>>>>>> Thanks Mina,
>>>>>>>>>> just wanted to be sure that it is not a linear
interpolation.
>>> Still I
>>>>>>>>>> wonder what happens when 675 and 440 are not available.
>>>>>>>>>> So if my understanding is correct MET only outputs 550nm.
>>>>>>>>>> Mariusz
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT
<met_help at ucar.edu>
>>>>>>>> wrote:
>>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>>
>>>>>>>>>>> Upon further checking from Howard, the implementation is
actually:
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> The Fortran code I referenced (
>>>>>> https://github.com/NCAR/MET/issues/1121
>>>>>>>> )
>>>>>>>>>> was
>>>>>>>>>>>> implemented that way:
>>>>>>>>>>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
>>>>>>>>>>>> aerdata(i)%aot_440nm)/log(675./440.)
>>>>>>>>>>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
>>>>>>>>>> ((550./440.)**(-a_675_440))
>>>>>>>>>>>
>>>>>>>>>>> Just in case you need to see the correct implementation
code. The
>>>>>>>>>>> first code snippet (in previous email) he cited was
replaced with
>>>>>> the
>>>>>>>>>>> above code snippet, which is associated with Github issue
1121...
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> Regards,
>>>>>>>>>>>
>>>>>>>>>>> Minna
>>>>>>>>>>>
>>>>>>>>>>> ---------------
>>>>>>>>>>> Minna Win
>>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
>>>>>>>>>> minnawin at ucar.edu>
>>>>>>>>>>> wrote:
>>>>>>>>>>>
>>>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>>>
>>>>>>>>>>>> According to Howard:
>>>>>>>>>>>>
>>>>>>>>>>>> It's a log interpolation: angstrom_power_interpolation
>>>>>>>>>>>>
>>>>>>>>>>>> double angstrom_power_interplation(double value_1, double
>>> value_2,
>>>>>>>>>>>> double level_1, double level_2, double target_level) {
>>>>>>>>>>>> double angstrom_log =
>>>>>> -log10(value_1/value_2)/log10(level_1/level_2);
>>>>>>>>>>>> double angstrom_value = value_2 *
>>>>>>>>>>>> pow((target_level/level_2),-angstrom_log);
>>>>>>>>>>>> return angstrom_value;
>>>>>>>>>>>> }
>>>>>>>>>>>>
>>>>>>>>>>>> Regards,
>>>>>>>>>>>> Minna
>>>>>>>>>>>> ---------------
>>>>>>>>>>>> Minna Win
>>>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT <
>>>>>>>>>> met_help at ucar.edu>
>>>>>>>>>>>> wrote:
>>>>>>>>>>>>
>>>>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>>>>>>>>>
>>>>>>>>>>>>> Minna,
>>>>>>>>>>>>> thanks very much - just one follow-up: how is the
interpolation
>>>>>> done
>>>>>>>>>>>>> (linear/log)?
>>>>>>>>>>>>> Mariusz
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT <
>>> met_help at ucar.edu
>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> According to our resident expert Howard Soh:
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> "ascii2nc reads a text format AERONET 3 data and
produces a
>>>>>> NetCDF
>>>>>>>>>>>>> output
>>>>>>>>>>>>>> (MET's point observation data file). ascii2nc
interpolates
>>> 550nm
>>>>>>>>>>>>> wavelength
>>>>>>>>>>>>>> from 440 nm and 675nm if not exists."
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> I hope this helps.
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Regards,
>>>>>>>>>>>>>> Minna
>>>>>>>>>>>>>> ---------------
>>>>>>>>>>>>>> Minna Win
>>>>>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via RT
<
>>>>>>>>>>>>> met_help at ucar.edu>
>>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted
upon.
>>>>>>>>>>>>>>> Transaction: Ticket created by
Mariusz.Pagowski at noaa.gov
>>>>>>>>>>>>>>> Queue: met_help
>>>>>>>>>>>>>>> Subject: AERONET AOD processing
>>>>>>>>>>>>>>> Owner: Nobody
>>>>>>>>>>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
>>>>>>>>>>>>>>> Status: new
>>>>>>>>>>>>>>> Ticket <URL:
>>>>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Hello,
>>>>>>>>>>>>>>> Is a formatting code available for version 3 AERONET?
(
>>>>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
>>>>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
>>>>>>>>>>>>>>> (If yes) Does the code allow for selecting chosen
wavelengths?
>>>>>> Is
>>>>>>>>>>>>> capable
>>>>>>>>>>>>>>> of obtaining selected channels (e.g. 550 nm) from
other
>>>>>>>> wavelengths?
>>>>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>>>> Mariusz
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>
>>>
>>>
>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Mariusz Pagowski
Time: Fri Mar 13 15:20:42 2020
If all values of AOD from AERONET are saved to ascii2nc output we are
good.
We will use 870nm whenever it is available i.e. MET will skip the
missing values in calculating statistics.
I will try running ascii2nc soon.
I have a cutout of AERONET 3 file that I recently downloaded - is this
the format ascii2nc can handle?
Thanks,
Mariusz
------------------------------------------------
Subject: AERONET AOD processing
From: Mariusz Pagowski
Time: Fri Mar 13 15:20:42 2020
> On Mar 13, 2020, at 3:05 PM, Howard Soh via RT <met_help at ucar.edu>
wrote:
>
> We need which wavelengths for interpolating for AOD_870nm.
>
> The sample AERONET3 file has valid values for following wavelengths:
>
> Date(dd:mm:yyyy) Time(hh:mm:ss) Day_of_Year AOD_1640nm
AOD_1020nm
> AOD_870nm AOD_675nm AOD_500nm AOD_440nm AOD_380nm AOD_340nm
> 0.779572 12:00:00 261 0.022949 0.031322 0.038461 0.049653
> 0.079107 0.097381 0.128962 0.139122
> 0.821238 12:00:00 262 0.01301 0.019142 0.022565 0.02795
0.041882
> 0.051215 0.063218 0.068337
> 0.862905 12:00:00 263 0.011324 0.016164 0.0199 0.023205
0.032876
> 0.039101 0.051348 0.054325
> 0.904572 12:00:00 264 0.010447 0.014793 0.018324 0.019833
> 0.026401 0.030903 0.043605 0.041629
> 0.987905 12:00:00 266 0.034807 0.046698 0.053709 0.06755
0.100127
> 0.118879 0.148089 0.164597
> 24:09:2015 12:00:00 267 0.02107 0.02921 0.034945 0.043239
> 0.065506 0.079168 0.104625 0.112846
> 25:09:2015 12:00:00 268 0.021043 0.028205 0.03332 0.040212
> 0.060867 0.074257 0.095984 0.105028
> 26:09:2015 12:00:00 269 0.015853 0.022548 0.027219 0.032956
> 0.051153 0.063412 0.084131 0.089649
>
>
> They are saved into ascii2nc output as is (AOD for 8 wavelengths).
> If there are missing values at AOD_675nm column, interpolated value
is
> saved.
>
> Cheers,
> Howard
>
>
> On 3/13/2020 2:12 PM, Mariusz Pagowski via RT wrote:
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>
>> Minna,
>> To clarify:
>> 1. MET 9.0 can output 440 nm and 675 nm? (along with interpolated
550nm)
>> 2. 870nm would be an addition - it would be fine only to output if
exists
>> Mariusz
>>
>>
>>> On Mar 13, 2020, at 2:08 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>>>
>>> Hi Mariusz,
>>>
>>> I updated the GitHub Issue. Howard reminded me that MET 9.0 has
support
>>> for 440nm and 675nm AERONET3. It looks like 870nm interpolation
support
>>> will need to be added.
>>>
>>> Regards,
>>> Minna
>>> ---------------
>>> Minna Win
>>> National Center for Atmospheric Research
>>> Developmental Testbed Center
>>> Phone: 303-497-8423
>>> Fax: 303-497-8401
>>>
>>>
>>>
>>> On Fri, Mar 13, 2020 at 2:05 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
>>> wrote:
>>>
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>
>>>> Yes, thanks,
>>>> Mariusz
>>>>
>>>>> On Mar 13, 2020, at 1:53 PM, Minna Win via RT
<met_help at ucar.edu> wrote:
>>>>>
>>>>> Hi Mariusz,
>>>>>
>>>>> I put down your affiliation as NOAA/GSD on the GitHub issue. Is
that
>>>>> correct?
>>>>>
>>>>> Thanks,
>>>>> Minna
>>>>> ---------------
>>>>> Minna Win
>>>>> National Center for Atmospheric Research
>>>>> Developmental Testbed Center
>>>>> Phone: 303-497-8423
>>>>> Fax: 303-497-8401
>>>>>
>>>>>
>>>>>
>>>>> On Fri, Mar 13, 2020 at 1:50 PM Minna Win-Gildenmeister <
>>>> minnawin at ucar.edu>
>>>>> wrote:
>>>>>
>>>>>> Hi Mariusz,
>>>>>>
>>>>>> I created a GitHub issue for MET:
>>>>>> https://github.com/NCAR/MET/issues/1278
>>>>>>
>>>>>> Please feel free to add any additional comments to the GitHub
issue.
>>>> I've
>>>>>> assigned the issue to Howard Soh and assigned this with a high
priority
>>>>>> (Blocker).
>>>>>>
>>>>>> Regards,
>>>>>> Minna
>>>>>> ---------------
>>>>>> Minna Win
>>>>>> National Center for Atmospheric Research
>>>>>> Developmental Testbed Center
>>>>>> Phone: 303-497-8423
>>>>>> Fax: 303-497-8401
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Fri, Mar 13, 2020 at 12:05 PM Mariusz Pagowski via RT <
>>>>>> met_help at ucar.edu> wrote:
>>>>>>
>>>>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>
>>>>>>>
>>>>>>> Minna,
>>>>>>> our time-frame is asap and would be happy to help with coding
if
>>>> someone
>>>>>>> provided guidance.
>>>>>>> Pls let me know what your and John’s thinking is,
>>>>>>> Mariusz
>>>>>>>
>>>>>>>
>>>>>>>> On Mar 13, 2020, at 11:50 AM, Minna Win via RT
<met_help at ucar.edu>
>>>>>>> wrote:
>>>>>>>> Hi Mariusz,
>>>>>>>>
>>>>>>>> Let me see if this request can be submitted in the form of a
GitHub
>>>>>>> issue.
>>>>>>>> We're currently ready to release MET v9.0, so if accepted, it
would
>>>> make
>>>>>>>> its way into a future release. I'll check with John to see
if this
>>>> can
>>>>>>> be
>>>>>>>> done within the current budget. What time-frame constraints
do you
>>>>>>> have?
>>>>>>>> Regards,
>>>>>>>> Minna
>>>>>>>> ---------------
>>>>>>>> Minna Win
>>>>>>>> National Center for Atmospheric Research
>>>>>>>> Developmental Testbed Center
>>>>>>>> Phone: 303-497-8423
>>>>>>>> Fax: 303-497-8401
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> On Fri, Mar 13, 2020 at 11:23 AM Mariusz Pagowski via RT <
>>>>>>> met_help at ucar.edu>
>>>>>>>> wrote:
>>>>>>>>
>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>>>>>
>>>>>>>>> Minna,
>>>>>>>>> thanks for your dedication to answering these questions.
>>>>>>>>> We are also interested in outputting 440 and 675 nm from
AERONET, and
>>>>>>> even
>>>>>>>>> possibly 870 nm.
>>>>>>>>> We would be ready to work with someone from MET to help
implement
>>>>>>> these as
>>>>>>>>> options.
>>>>>>>>> Pls let me know if that is possible,
>>>>>>>>> Mariusz
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>> On Mar 13, 2020, at 8:03 AM, Minna Win via RT
<met_help at ucar.edu>
>>>>>>> wrote:
>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>
>>>>>>>>>> Here's what happens when 675 and 440 nm are unavailable:
>>>>>>>>>>
>>>>>>>>>> ascii2nc does not save 550 nm wavelength into NetCDF output
if 675
>>>> and
>>>>>>>>> 440
>>>>>>>>>>> are not available.
>>>>>>>>>>> ascii2nc figures out the variable list and associated
wavelengths
>>>>>>> from
>>>>>>>>> the
>>>>>>>>>>> header record and knows where AOD at 440 and 650 nm values
are.
>>>>>>>>>>> ascii2nc saves AOD of all available
>>>>>>>>>>> wavelengths. And ascii2nc checks AOD at 440, 550 and
>>>>>>>>>>> 675. ascii2nc computes & saves AOD at 550 if If AOD at 550
does not
>>>>>>>>> exist
>>>>>>>>>>> from data record and AOD values at 675 and at 440 are
valid.
>>>>>>>>>>
>>>>>>>>>> I hope this is useful.
>>>>>>>>>>
>>>>>>>>>> Regards,
>>>>>>>>>> Minna
>>>>>>>>>> ---------------
>>>>>>>>>> Minna Win
>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> On Thu, Mar 12, 2020 at 4:03 PM Mariusz Pagowski via RT <
>>>>>>>>> met_help at ucar.edu>
>>>>>>>>>> wrote:
>>>>>>>>>>
>>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>>>>>>>
>>>>>>>>>>> Thanks Mina,
>>>>>>>>>>> just wanted to be sure that it is not a linear
interpolation.
>>>> Still I
>>>>>>>>>>> wonder what happens when 675 and 440 are not available.
>>>>>>>>>>> So if my understanding is correct MET only outputs 550nm.
>>>>>>>>>>> Mariusz
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>> On Mar 12, 2020, at 3:56 PM, Minna Win via RT
<met_help at ucar.edu>
>>>>>>>>> wrote:
>>>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>>>
>>>>>>>>>>>> Upon further checking from Howard, the implementation is
actually:
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> The Fortran code I referenced (
>>>>>>> https://github.com/NCAR/MET/issues/1121
>>>>>>>>> )
>>>>>>>>>>> was
>>>>>>>>>>>>> implemented that way:
>>>>>>>>>>>>> a_675_440 = -log(aerdata(i)%aot_675nm/ &
>>>>>>>>>>>>> aerdata(i)%aot_440nm)/log(675./440.)
>>>>>>>>>>>>> aerdata(i)%aot_550nm = aerdata(i)%aot_440nm* &
>>>>>>>>>>> ((550./440.)**(-a_675_440))
>>>>>>>>>>>>
>>>>>>>>>>>> Just in case you need to see the correct implementation
code. The
>>>>>>>>>>>> first code snippet (in previous email) he cited was
replaced with
>>>>>>> the
>>>>>>>>>>>> above code snippet, which is associated with Github issue
1121...
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> Regards,
>>>>>>>>>>>>
>>>>>>>>>>>> Minna
>>>>>>>>>>>>
>>>>>>>>>>>> ---------------
>>>>>>>>>>>> Minna Win
>>>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>>
>>>>>>>>>>>> On Thu, Mar 12, 2020 at 2:36 PM Minna Win-Gildenmeister <
>>>>>>>>>>> minnawin at ucar.edu>
>>>>>>>>>>>> wrote:
>>>>>>>>>>>>
>>>>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>>>>
>>>>>>>>>>>>> According to Howard:
>>>>>>>>>>>>>
>>>>>>>>>>>>> It's a log interpolation: angstrom_power_interpolation
>>>>>>>>>>>>>
>>>>>>>>>>>>> double angstrom_power_interplation(double value_1,
double
>>>> value_2,
>>>>>>>>>>>>> double level_1, double level_2, double target_level) {
>>>>>>>>>>>>> double angstrom_log =
>>>>>>> -log10(value_1/value_2)/log10(level_1/level_2);
>>>>>>>>>>>>> double angstrom_value = value_2 *
>>>>>>>>>>>>> pow((target_level/level_2),-angstrom_log);
>>>>>>>>>>>>> return angstrom_value;
>>>>>>>>>>>>> }
>>>>>>>>>>>>>
>>>>>>>>>>>>> Regards,
>>>>>>>>>>>>> Minna
>>>>>>>>>>>>> ---------------
>>>>>>>>>>>>> Minna Win
>>>>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>>
>>>>>>>>>>>>> On Thu, Mar 12, 2020 at 2:11 PM Mariusz Pagowski via RT
<
>>>>>>>>>>> met_help at ucar.edu>
>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>
>>>>>>>>>>>>>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>>>>>>>>>>>>>>
>>>>>>>>>>>>>> Minna,
>>>>>>>>>>>>>> thanks very much - just one follow-up: how is the
interpolation
>>>>>>> done
>>>>>>>>>>>>>> (linear/log)?
>>>>>>>>>>>>>> Mariusz
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> On Mar 12, 2020, at 2:06 PM, Minna Win via RT <
>>>> met_help at ucar.edu
>>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>>> Hi Mariusz,
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> According to our resident expert Howard Soh:
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> "ascii2nc reads a text format AERONET 3 data and
produces a
>>>>>>> NetCDF
>>>>>>>>>>>>>> output
>>>>>>>>>>>>>>> (MET's point observation data file). ascii2nc
interpolates
>>>> 550nm
>>>>>>>>>>>>>> wavelength
>>>>>>>>>>>>>>> from 440 nm and 675nm if not exists."
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> I hope this helps.
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> Regards,
>>>>>>>>>>>>>>> Minna
>>>>>>>>>>>>>>> ---------------
>>>>>>>>>>>>>>> Minna Win
>>>>>>>>>>>>>>> National Center for Atmospheric Research
>>>>>>>>>>>>>>> Developmental Testbed Center
>>>>>>>>>>>>>>> Phone: 303-497-8423
>>>>>>>>>>>>>>> Fax: 303-497-8401
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>> On Thu, Mar 12, 2020 at 1:27 PM Mariusz Pagowski via
RT <
>>>>>>>>>>>>>> met_help at ucar.edu>
>>>>>>>>>>>>>>> wrote:
>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Thu Mar 12 13:27:57 2020: Request 94517 was acted
upon.
>>>>>>>>>>>>>>>> Transaction: Ticket created by
Mariusz.Pagowski at noaa.gov
>>>>>>>>>>>>>>>> Queue: met_help
>>>>>>>>>>>>>>>> Subject: AERONET AOD processing
>>>>>>>>>>>>>>>> Owner: Nobody
>>>>>>>>>>>>>>>> Requestors: Mariusz.Pagowski at noaa.gov
>>>>>>>>>>>>>>>> Status: new
>>>>>>>>>>>>>>>> Ticket <URL:
>>>>>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>> Hello,
>>>>>>>>>>>>>>>> Is a formatting code available for version 3 AERONET?
(
>>>>>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html <
>>>>>>>>>>>>>>>> https://aeronet.gsfc.nasa.gov/new_web/data.html>)
>>>>>>>>>>>>>>>> (If yes) Does the code allow for selecting chosen
wavelengths?
>>>>>>> Is
>>>>>>>>>>>>>> capable
>>>>>>>>>>>>>>>> of obtaining selected channels (e.g. 550 nm) from
other
>>>>>>>>> wavelengths?
>>>>>>>>>>>>>>>> Thanks,
>>>>>>>>>>>>>>>> Mariusz
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>
>>>>
>>>>
>>
>>
>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Howard Soh
Time: Fri Mar 13 16:13:10 2020
The first column of the attached file is "Site_location(m)".
ascii2nc only supports AERONET V3 data file with first three columns
for date info: "Date(dd:mm:yyyy)", "Time(hh:mm:ss)", and "Day_of_Year"
columns.
ascii2nc does not save any observations because the valid time from
the first three columns is 0 which MET considers it as an invalid
record.
Cheers,
Howard
On Fri Mar 13 15:20:42 2020, Mariusz.Pagowski at noaa.gov wrote:
> If all values of AOD from AERONET are saved to ascii2nc output we
are
> good.
> We will use 870nm whenever it is available i.e. MET will skip the
> missing values in calculating statistics.
> I will try running ascii2nc soon.
> I have a cutout of AERONET 3 file that I recently downloaded - is
this
> the format ascii2nc can handle?
> Thanks,
> Mariusz
------------------------------------------------
Subject: AERONET AOD processing
From: Mariusz Pagowski
Time: Fri Mar 13 16:21:55 2020
Ok, this file was downloaded from NASA website and is formatted for
AERONET 3 with the following command
wget --no-check-certificate -q -O AERONET_NA_2015.lev20 "
https://aeronet.gsfc.nasa.gov/cgi-
bin/print_web_data_v3?year=2015&month=1&day=1&year2=2015&month2=12&day2=13&AOD20=1&AVG=20&lat1=12.&lon1=-
175.&lat2=71.&lon2=-60
."
Where/how to download data to concur with the ascii2nc format?
Mariusz
On Fri, Mar 13, 2020 at 4:13 PM Howard Soh via RT <met_help at ucar.edu>
wrote:
> The first column of the attached file is "Site_location(m)".
> ascii2nc only supports AERONET V3 data file with first three columns
for
> date info: "Date(dd:mm:yyyy)", "Time(hh:mm:ss)", and "Day_of_Year"
columns.
> ascii2nc does not save any observations because the valid time from
the
> first three columns is 0 which MET considers it as an invalid
record.
>
> Cheers,
> Howard
>
> On Fri Mar 13 15:20:42 2020, Mariusz.Pagowski at noaa.gov wrote:
> > If all values of AOD from AERONET are saved to ascii2nc output we
are
> > good.
> > We will use 870nm whenever it is available i.e. MET will skip the
> > missing values in calculating statistics.
> > I will try running ascii2nc soon.
> > I have a cutout of AERONET 3 file that I recently downloaded - is
this
> > the format ascii2nc can handle?
> > Thanks,
> > Mariusz
>
>
>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Howard Soh
Time: Fri Mar 13 16:57:35 2020
I downloaded the sample data files by using a web browser from
https://aeronet.gsfc.nasa.gov/cgi-bin/draw_map_display_aod_v3?level=3
Please consult it with AERONET 3.
Cheers,
Howard
On Fri Mar 13 16:21:55 2020, Mariusz.Pagowski at noaa.gov wrote:
> Ok, this file was downloaded from NASA website and is formatted for
> AERONET 3 with the following command
> wget --no-check-certificate -q -O AERONET_NA_2015.lev20 "
> https://aeronet.gsfc.nasa.gov/cgi-
>
bin/print_web_data_v3?year=2015&month=1&day=1&year2=2015&month2=12&day2=13&AOD20=1&AVG=20&lat1=12.&lon1=-
> 175.&lat2=71.&lon2=-60
> ."
> Where/how to download data to concur with the ascii2nc format?
> Mariusz
>
>
> On Fri, Mar 13, 2020 at 4:13 PM Howard Soh via RT
<met_help at ucar.edu>
> wrote:
>
> > The first column of the attached file is "Site_location(m)".
> > ascii2nc only supports AERONET V3 data file with first three
columns
> > for
> > date info: "Date(dd:mm:yyyy)", "Time(hh:mm:ss)", and "Day_of_Year"
> > columns.
> > ascii2nc does not save any observations because the valid time
from
> > the
> > first three columns is 0 which MET considers it as an invalid
record.
> >
> > Cheers,
> > Howard
> >
> > On Fri Mar 13 15:20:42 2020, Mariusz.Pagowski at noaa.gov wrote:
> > > If all values of AOD from AERONET are saved to ascii2nc output
we
> > > are
> > > good.
> > > We will use 870nm whenever it is available i.e. MET will skip
the
> > > missing values in calculating statistics.
> > > I will try running ascii2nc soon.
> > > I have a cutout of AERONET 3 file that I recently downloaded -
is
> > > this
> > > the format ascii2nc can handle?
> > > Thanks,
> > > Mariusz
> >
> >
> >
> >
------------------------------------------------
Subject: AERONET AOD processing
From: Mariusz Pagowski
Time: Fri Mar 13 17:01:28 2020
Ok, we need to sort it out because I received this command from the
official AERONET distribution office.
Mariusz
On Fri, Mar 13, 2020 at 4:57 PM Howard Soh via RT <met_help at ucar.edu>
wrote:
> I downloaded the sample data files by using a web browser from
> https://aeronet.gsfc.nasa.gov/cgi-
bin/draw_map_display_aod_v3?level=3
>
> Please consult it with AERONET 3.
>
> Cheers,
> Howard
>
> On Fri Mar 13 16:21:55 2020, Mariusz.Pagowski at noaa.gov wrote:
> > Ok, this file was downloaded from NASA website and is formatted
for
> > AERONET 3 with the following command
> > wget --no-check-certificate -q -O AERONET_NA_2015.lev20 "
> > https://aeronet.gsfc.nasa.gov/cgi-
> >
>
bin/print_web_data_v3?year=2015&month=1&day=1&year2=2015&month2=12&day2=13&AOD20=1&AVG=20&lat1=12.&lon1=-
> > 175.&lat2=71.&lon2=-60
> > ."
> > Where/how to download data to concur with the ascii2nc format?
> > Mariusz
> >
> >
> > On Fri, Mar 13, 2020 at 4:13 PM Howard Soh via RT
<met_help at ucar.edu>
> > wrote:
> >
> > > The first column of the attached file is "Site_location(m)".
> > > ascii2nc only supports AERONET V3 data file with first three
columns
> > > for
> > > date info: "Date(dd:mm:yyyy)", "Time(hh:mm:ss)", and
"Day_of_Year"
> > > columns.
> > > ascii2nc does not save any observations because the valid time
from
> > > the
> > > first three columns is 0 which MET considers it as an invalid
record.
> > >
> > > Cheers,
> > > Howard
> > >
> > > On Fri Mar 13 15:20:42 2020, Mariusz.Pagowski at noaa.gov wrote:
> > > > If all values of AOD from AERONET are saved to ascii2nc output
we
> > > > are
> > > > good.
> > > > We will use 870nm whenever it is available i.e. MET will skip
the
> > > > missing values in calculating statistics.
> > > > I will try running ascii2nc soon.
> > > > I have a cutout of AERONET 3 file that I recently downloaded -
is
> > > > this
> > > > the format ascii2nc can handle?
> > > > Thanks,
> > > > Mariusz
> > >
> > >
> > >
> > >
>
>
>
>
------------------------------------------------
Subject: AERONET AOD processing
From: Mariusz Pagowski
Time: Tue Mar 17 12:10:36 2020
Howard:
Barry has python scripts for obtaining AERONET data and can format
them
appropriately for MET.
Can you provide input format that you can process with your converter
and
its output for MET in netCDF so that it can be used
as template in Barry's code?
Thanks,
Mariusz
On Fri, Mar 13, 2020 at 5:01 PM Mariusz Pagowski - NOAA Affiliate <
mariusz.pagowski at noaa.gov> wrote:
> Ok, we need to sort it out because I received this command from the
> official AERONET distribution office.
> Mariusz
>
> On Fri, Mar 13, 2020 at 4:57 PM Howard Soh via RT
<met_help at ucar.edu>
> wrote:
>
>> I downloaded the sample data files by using a web browser from
>> https://aeronet.gsfc.nasa.gov/cgi-
bin/draw_map_display_aod_v3?level=3
>>
>> Please consult it with AERONET 3.
>>
>> Cheers,
>> Howard
>>
>> On Fri Mar 13 16:21:55 2020, Mariusz.Pagowski at noaa.gov wrote:
>> > Ok, this file was downloaded from NASA website and is formatted
for
>> > AERONET 3 with the following command
>> > wget --no-check-certificate -q -O AERONET_NA_2015.lev20 "
>> > https://aeronet.gsfc.nasa.gov/cgi-
>> >
>>
bin/print_web_data_v3?year=2015&month=1&day=1&year2=2015&month2=12&day2=13&AOD20=1&AVG=20&lat1=12.&lon1=-
>> > 175.&lat2=71.&lon2=-60
>> > ."
>> > Where/how to download data to concur with the ascii2nc format?
>> > Mariusz
>> >
>> >
>> > On Fri, Mar 13, 2020 at 4:13 PM Howard Soh via RT
<met_help at ucar.edu>
>> > wrote:
>> >
>> > > The first column of the attached file is "Site_location(m)".
>> > > ascii2nc only supports AERONET V3 data file with first three
columns
>> > > for
>> > > date info: "Date(dd:mm:yyyy)", "Time(hh:mm:ss)", and
"Day_of_Year"
>> > > columns.
>> > > ascii2nc does not save any observations because the valid time
from
>> > > the
>> > > first three columns is 0 which MET considers it as an invalid
record.
>> > >
>> > > Cheers,
>> > > Howard
>> > >
>> > > On Fri Mar 13 15:20:42 2020, Mariusz.Pagowski at noaa.gov wrote:
>> > > > If all values of AOD from AERONET are saved to ascii2nc
output we
>> > > > are
>> > > > good.
>> > > > We will use 870nm whenever it is available i.e. MET will skip
the
>> > > > missing values in calculating statistics.
>> > > > I will try running ascii2nc soon.
>> > > > I have a cutout of AERONET 3 file that I recently downloaded
- is
>> > > > this
>> > > > the format ascii2nc can handle?
>> > > > Thanks,
>> > > > Mariusz
>> > >
>> > >
>> > >
>> > >
>>
>>
>>
>>
------------------------------------------------
Subject: AERONET AOD processing
From: Tara Jensen
Time: Tue Mar 17 12:18:00 2020
Mariusz,
We have spoken with Barry, as recently as last month, about making
connections with METplus. Now that the METplus 3.0 release is
officially
out, he can start working with us to call MONET scripts via our tools.
This seems like a really good use-case to work on.
Howard, please close this ticket. Mariusz, I'm going to take this off-
line
and open a line of communication between you, Barry, and myself.
Sound good?
Thanks so much, Tara
On Tue, Mar 17, 2020 at 12:10 PM Mariusz Pagowski via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94517 >
>
> Howard:
> Barry has python scripts for obtaining AERONET data and can format
them
> appropriately for MET.
> Can you provide input format that you can process with your
converter and
> its output for MET in netCDF so that it can be used
> as template in Barry's code?
> Thanks,
> Mariusz
>
>
> On Fri, Mar 13, 2020 at 5:01 PM Mariusz Pagowski - NOAA Affiliate <
> mariusz.pagowski at noaa.gov> wrote:
>
> > Ok, we need to sort it out because I received this command from
the
> > official AERONET distribution office.
> > Mariusz
> >
> > On Fri, Mar 13, 2020 at 4:57 PM Howard Soh via RT
<met_help at ucar.edu>
> > wrote:
> >
> >> I downloaded the sample data files by using a web browser from
> >> https://aeronet.gsfc.nasa.gov/cgi-
bin/draw_map_display_aod_v3?level=3
> >>
> >> Please consult it with AERONET 3.
> >>
> >> Cheers,
> >> Howard
> >>
> >> On Fri Mar 13 16:21:55 2020, Mariusz.Pagowski at noaa.gov wrote:
> >> > Ok, this file was downloaded from NASA website and is formatted
for
> >> > AERONET 3 with the following command
> >> > wget --no-check-certificate -q -O AERONET_NA_2015.lev20 "
> >> > https://aeronet.gsfc.nasa.gov/cgi-
> >> >
> >>
>
bin/print_web_data_v3?year=2015&month=1&day=1&year2=2015&month2=12&day2=13&AOD20=1&AVG=20&lat1=12.&lon1=-
> >> > 175.&lat2=71.&lon2=-60
> >> > ."
> >> > Where/how to download data to concur with the ascii2nc format?
> >> > Mariusz
> >> >
> >> >
> >> > On Fri, Mar 13, 2020 at 4:13 PM Howard Soh via RT
<met_help at ucar.edu>
> >> > wrote:
> >> >
> >> > > The first column of the attached file is "Site_location(m)".
> >> > > ascii2nc only supports AERONET V3 data file with first three
columns
> >> > > for
> >> > > date info: "Date(dd:mm:yyyy)", "Time(hh:mm:ss)", and
"Day_of_Year"
> >> > > columns.
> >> > > ascii2nc does not save any observations because the valid
time from
> >> > > the
> >> > > first three columns is 0 which MET considers it as an invalid
> record.
> >> > >
> >> > > Cheers,
> >> > > Howard
> >> > >
> >> > > On Fri Mar 13 15:20:42 2020, Mariusz.Pagowski at noaa.gov wrote:
> >> > > > If all values of AOD from AERONET are saved to ascii2nc
output we
> >> > > > are
> >> > > > good.
> >> > > > We will use 870nm whenever it is available i.e. MET will
skip the
> >> > > > missing values in calculating statistics.
> >> > > > I will try running ascii2nc soon.
> >> > > > I have a cutout of AERONET 3 file that I recently
downloaded - is
> >> > > > this
> >> > > > the format ascii2nc can handle?
> >> > > > Thanks,
> >> > > > Mariusz
> >> > >
> >> > >
> >> > >
> >> > >
> >>
> >>
> >>
> >>
>
>
--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Tara Jensen
Project Manager II
NCAR RAL and DTC
PO Box 3000, Boulder, Colorado 80307 USA
+1 303-497-8479 jensen at ucar.edu
------------------------------------------------
Subject: AERONET AOD processing
From: Howard Soh
Time: Tue Mar 17 12:21:52 2020
The attached file was modified from the file you sent.
The file you sent us did not work because
- It was a HTML file (not valid input because of extra HTML tags)
- The second line did not have the site name
Cheers,
Howard
On Tue Mar 17 12:10:36 2020, Mariusz.Pagowski at noaa.gov wrote:
> Howard:
> Barry has python scripts for obtaining AERONET data and can format
> them
> appropriately for MET.
> Can you provide input format that you can process with your
converter
> and
> its output for MET in netCDF so that it can be used
> as template in Barry's code?
> Thanks,
> Mariusz
>
>
> On Fri, Mar 13, 2020 at 5:01 PM Mariusz Pagowski - NOAA Affiliate <
> mariusz.pagowski at noaa.gov> wrote:
>
> > Ok, we need to sort it out because I received this command from
the
> > official AERONET distribution office.
> > Mariusz
> >
> > On Fri, Mar 13, 2020 at 4:57 PM Howard Soh via RT
<met_help at ucar.edu>
> > wrote:
> >
> >> I downloaded the sample data files by using a web browser from
> >> https://aeronet.gsfc.nasa.gov/cgi-
> >> bin/draw_map_display_aod_v3?level=3
> >>
> >> Please consult it with AERONET 3.
> >>
> >> Cheers,
> >> Howard
> >>
> >> On Fri Mar 13 16:21:55 2020, Mariusz.Pagowski at noaa.gov wrote:
> >> > Ok, this file was downloaded from NASA website and is formatted
> >> > for
> >> > AERONET 3 with the following command
> >> > wget --no-check-certificate -q -O AERONET_NA_2015.lev20 "
> >> > https://aeronet.gsfc.nasa.gov/cgi-
> >> >
> >>
bin/print_web_data_v3?year=2015&month=1&day=1&year2=2015&month2=12&day2=13&AOD20=1&AVG=20&lat1=12.&lon1=-
> >> > 175.&lat2=71.&lon2=-60
> >> > ."
> >> > Where/how to download data to concur with the ascii2nc format?
> >> > Mariusz
> >> >
> >> >
> >> > On Fri, Mar 13, 2020 at 4:13 PM Howard Soh via RT
> >> > <met_help at ucar.edu>
> >> > wrote:
> >> >
> >> > > The first column of the attached file is "Site_location(m)".
> >> > > ascii2nc only supports AERONET V3 data file with first three
> >> > > columns
> >> > > for
> >> > > date info: "Date(dd:mm:yyyy)", "Time(hh:mm:ss)", and
> >> > > "Day_of_Year"
> >> > > columns.
> >> > > ascii2nc does not save any observations because the valid
time
> >> > > from
> >> > > the
> >> > > first three columns is 0 which MET considers it as an invalid
> >> > > record.
> >> > >
> >> > > Cheers,
> >> > > Howard
> >> > >
> >> > > On Fri Mar 13 15:20:42 2020, Mariusz.Pagowski at noaa.gov wrote:
> >> > > > If all values of AOD from AERONET are saved to ascii2nc
output
> >> > > > we
> >> > > > are
> >> > > > good.
> >> > > > We will use 870nm whenever it is available i.e. MET will
skip
> >> > > > the
> >> > > > missing values in calculating statistics.
> >> > > > I will try running ascii2nc soon.
> >> > > > I have a cutout of AERONET 3 file that I recently
downloaded -
> >> > > > is
> >> > > > this
> >> > > > the format ascii2nc can handle?
> >> > > > Thanks,
> >> > > > Mariusz
> >> > >
> >> > >
> >> > >
> >> > >
> >>
> >>
> >>
> >>
------------------------------------------------
Subject: AERONET AOD processing
From: Howard Soh
Time: Tue Mar 17 12:26:04 2020
Closed this ticket based on Tara's comment.
------------------------------------------------
More information about the Met_help
mailing list