search

NOAA-SWPC / GSMWAM-IPE

Coupled NEMS app for WAM-IPE
5 stars 11 forks source link

CCMC neutal density product validation #75

Closed ZhuxiaoLi66 closed 1 year ago

ZhuxiaoLi66 commented 1 year ago

The neutral density product data files for Feb.3 2022 storm delivered by CCMC have been checked. The comparison plots with our original product files are posted for further discussion and record.

ZhuxiaoLi66 commented 1 year ago

DEN_CCMC_250km_4p_20220203 DEN_product_250km_4p_20220203 DEN_CCMC_400km_4p_20220203 DEN_product_400km_4p_20220203 DEN_CCMC_800km_4p_20220203 DEN_product_800km_4p_20220203 DEN_CCMC_1000km_4p_20220203 DEN_product_1000km_4p_20220203 DEN_CCMC_800km_24p_20220203 DEN_product_800km_24p_20220203 DEN_CCMC_400km_24p_20220203 DEN_product_400km_24p_20220203 DEN_CCMC_600km_4p_20220203 DEN_product_600km_4p_20220203

ZhuxiaoLi66 commented 1 year ago

Hi All, would like to tell you some good news, after I revised my Fortran code and let it do the neutral density expolation since WAM Level 149, which means excluding Level 150 for the interpolation, the results become very close to CCMC data with Adam's package. In conclusion, now I think we shouldn't have any concerns about the neutral density product from CCMC. Would like to post the plots for the new result at 800km here, also attached the CCMC plots for comparison.

DEN_product_L149_800km_4p_20220203 DEN_CCMC_800km_4p_20220203 DEN_product_L149_800km_24p_20220203 DEN_CCMC_800km_24p_20220203

ZhuxiaoLi66 commented 1 year ago

The following plots are the density difference (CCMC-product) ratio to the neutral product [(CCMC-product)/product]. larger along the altitude, I think they are small enough even for 800km. DEN_diff_perc_CCMC_product_250km_4p_20220203 DEN_diff_perc_CCMC_product_400km_4p_20220203 DEN_diff_perc_CCMC_product_800km_4p_20220203

ZhuxiaoLi66 commented 1 year ago

The other thing I found in CCMC files is that the unit of the neutral density is "m^-3", which should be "kg*m^-3". since the formal name of the variable is "thermosphere_mass_density".

twfang commented 1 year ago

Thanks for catching that Zhuxiao.

akubaryk commented 1 year ago

The other thing I found in CCMC files is that the unit of the neutral density is "m^-3", which should be "kg*m^-3". since the formal name of the variable is "thermosphere_mass_density".

A fix has been pushed to address this.

ZhuxiaoLi66 commented 1 year ago

@akubaryk nice to know, thanks.