*** Results log of t_fascod_1 (FM FASCOD) *** ./t_fascod NL: state%non_ideal, state%new_ref_op, state%new_bangle_op = F F F Checking FASCOD Scenario: FASCOD_TRO Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_MLS Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_MLW Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_SAS Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_SAW Refractivity passed Bending angle passed **************************** ********** PASS ********** **************************** *** Results log of t_fascod_2 (FM FASCOD -comp) *** ./t_fascod -comp NL: state%non_ideal, state%new_ref_op, state%new_bangle_op = T F F Checking FASCOD Scenario: FASCOD_TRO Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_MLS Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_MLW Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_SAS Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_SAW Refractivity passed Bending angle passed **************************** ********** PASS ********** **************************** *** Results log of t_fascod_3 (FM FASCOD -new_op) *** ./t_fascod -new_op NL: state%non_ideal, state%new_ref_op, state%new_bangle_op = F T T Checking FASCOD Scenario: FASCOD_TRO Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_MLS Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_MLW Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_SAS Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_SAW Refractivity passed Bending angle passed **************************** ********** PASS ********** **************************** *** Results log of t_fascod_4 (FM FASCOD -comp -new_op) *** ./t_fascod -comp -new_op NL: state%non_ideal, state%new_ref_op, state%new_bangle_op = T T T Checking FASCOD Scenario: FASCOD_TRO Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_MLS Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_MLW Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_SAS Refractivity passed Bending angle passed Checking FASCOD Scenario: FASCOD_SAW Refractivity passed Bending angle passed **************************** ********** PASS ********** **************************** *** Results log of t_fascod_ad_1 (FM_AD FASCOD) *** ./t_fascod_ad AD: x%non_ideal, x%new_ref_op, x%new_bangle_op = F F F Checking REFRACTIVITY AD FASCOD_TRO norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 473496185.147 473496185.147 1.0000000 0.00000E+00 Checking BENDING ANGLE AD FASCOD_TRO norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 76.512 76.512 1.0000000 -0.18573E-15 Checking REFRACTIVITY AD FASCOD_MLS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 355823536.657 355823536.657 1.0000000 -0.16751E-15 Checking BENDING ANGLE AD FASCOD_MLS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 144.439 144.439 1.0000000 -0.17710E-14 Checking REFRACTIVITY AD FASCOD_MLW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 474645859.385 474645859.385 1.0000000 -0.50231E-15 Checking BENDING ANGLE AD FASCOD_MLW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 143.528 143.528 1.0000000 -0.23763E-14 Checking REFRACTIVITY AD FASCOD_SAS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 379867322.074 379867322.074 1.0000000 0.00000E+00 Checking BENDING ANGLE AD FASCOD_SAS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 107.062 107.062 1.0000000 0.10619E-14 Checking REFRACTIVITY AD FASCOD_SAW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 475400740.028 475400740.028 1.0000000 -0.25076E-15 Checking BENDING ANGLE AD FASCOD_SAW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 79.056 79.056 1.0000000 0.14381E-14 **************************** ********** PASS ********** **************************** *** Results log of t_fascod_ad_2 (FM_AD FASCOD -comp) *** ./t_fascod_ad -comp AD: x%non_ideal, x%new_ref_op, x%new_bangle_op = T F F Checking REFRACTIVITY AD FASCOD_TRO norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 476145112.147 476145112.147 1.0000000 0.25036E-15 Checking BENDING ANGLE AD FASCOD_TRO norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 76.698 76.698 1.0000000 -0.11117E-14 Checking REFRACTIVITY AD FASCOD_MLS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 358157528.059 358157528.059 1.0000000 0.49926E-15 Checking BENDING ANGLE AD FASCOD_MLS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 143.235 143.235 1.0000000 -0.15874E-14 Checking REFRACTIVITY AD FASCOD_MLW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 476706619.225 476706619.225 1.0000000 0.12503E-15 Checking BENDING ANGLE AD FASCOD_MLW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 141.424 141.424 1.0000000 0.16077E-14 Checking REFRACTIVITY AD FASCOD_SAS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 380318136.194 380318136.194 1.0000000 -0.15672E-15 Checking BENDING ANGLE AD FASCOD_SAS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 105.608 105.608 1.0000000 -0.14802E-14 Checking REFRACTIVITY AD FASCOD_SAW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 479397420.033 479397420.033 1.0000000 -0.62166E-15 Checking BENDING ANGLE AD FASCOD_SAW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 78.015 78.015 1.0000000 0.16394E-14 **************************** ********** PASS ********** **************************** *** Results log of t_fascod_ad_3 (FM_AD FASCOD -new_op) *** ./t_fascod_ad -new_op AD: x%non_ideal, x%new_ref_op, x%new_bangle_op = F T T Checking REFRACTIVITY AD FASCOD_TRO norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 474171710.080 474171710.080 1.0000000 0.25141E-15 Checking BENDING ANGLE AD FASCOD_TRO norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 76.516 76.516 1.0000000 -0.31573E-14 Checking REFRACTIVITY AD FASCOD_MLS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 356516009.418 356516009.418 1.0000000 0.16719E-15 Checking BENDING ANGLE AD FASCOD_MLS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 144.453 144.453 1.0000000 0.00000E+00 Checking REFRACTIVITY AD FASCOD_MLW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 474353041.944 474353041.944 1.0000000 -0.25131E-15 Checking BENDING ANGLE AD FASCOD_MLW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 143.529 143.529 1.0000000 -0.59406E-15 Checking REFRACTIVITY AD FASCOD_SAS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 379820652.360 379820652.360 1.0000000 -0.15693E-15 Checking BENDING ANGLE AD FASCOD_SAS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 107.101 107.101 1.0000000 -0.26537E-15 Checking REFRACTIVITY AD FASCOD_SAW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 476796543.213 476796543.213 1.0000000 0.37503E-15 Checking BENDING ANGLE AD FASCOD_SAW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 79.058 79.058 1.0000000 0.10785E-14 **************************** ********** PASS ********** **************************** *** Results log of t_fascod_ad_4 (FM_AD FASCOD -comp -new_op) *** ./t_fascod_ad -comp -new_op AD: x%non_ideal, x%new_ref_op, x%new_bangle_op = T T T Checking REFRACTIVITY AD FASCOD_TRO norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 476739683.528 476739683.528 1.0000000 -0.12503E-15 Checking BENDING ANGLE AD FASCOD_TRO norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 76.814 76.814 1.0000000 -0.37000E-14 Checking REFRACTIVITY AD FASCOD_MLS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 358898601.347 358898601.347 1.0000000 -0.16608E-15 Checking BENDING ANGLE AD FASCOD_MLS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 143.247 143.247 1.0000000 0.19841E-15 Checking REFRACTIVITY AD FASCOD_MLW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 476429730.446 476429730.446 1.0000000 0.12511E-15 Checking BENDING ANGLE AD FASCOD_MLW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 141.424 141.424 1.0000000 -0.10048E-14 Checking REFRACTIVITY AD FASCOD_SAS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 380293266.393 380293266.393 1.0000000 0.47020E-15 Checking BENDING ANGLE AD FASCOD_SAS norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 105.608 105.608 1.0000000 0.67281E-15 Checking REFRACTIVITY AD FASCOD_SAW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 480702616.199 480702616.199 1.0000000 -0.12399E-15 Checking BENDING ANGLE AD FASCOD_SAW norm1, norm2, norm1/norm2, (norm1-norm2)/norm2 = 78.015 78.015 1.0000000 0.72862E-15 **************************** ********** PASS ********** **************************** *** Results log of t_fascod_tl_1 (FM_TL FASCOD) *** ./t_fascod_tl TL: x%non_ideal, x%new_ref_op, x%new_bangle_op = F F F Checking REFRACTIVITY TL FASCOD_TRO |dy|/|dy_init| cos_angle 100*rel diff 1.0000000000E+00 9.2501426076E-01 8.1817807106E+01 1.0000000000E-01 9.9702326361E-01 1.0402015044E+01 1.0000000000E-02 9.9995571205E-01 1.2995096692E+00 1.0000000000E-03 9.9999951636E-01 1.3724114860E-01 1.0000000000E-04 9.9999999511E-01 1.3813887372E-02 1.0000000000E-05 9.9999999995E-01 1.3823072279E-03 1.0000000000E-06 1.0000000000E+00 1.3823916825E-04 1.0000000000E-07 1.0000000000E+00 1.3818809855E-05 1.0000000000E-08 1.0000000000E+00 1.3318964378E-06 1.0000000000E-09 1.0000000000E+00 3.6207257604E-06 1.0000000000E-10 1.0000000000E+00 3.3375610642E-05 1.0000000000E-11 9.9999999999E-01 3.3200414500E-04 1.0000000000E-12 9.9999999960E-01 2.8283683774E-03 1.0000000000E-13 9.9999994042E-01 3.4527480635E-02 1.0000000000E-14 9.9999504438E-01 3.2317440101E-01 Refractivity TL check passed Checking BENDING ANGLE TL FASCOD_TRO |dy|/|dy_init| cos_angle 100*rel diff 1.0000000000E+00 -1.6726857540E-01 1.0000000012E+02 1.0000000000E-01 -2.0057851944E-01 1.0000000002E+02 1.0000000000E-02 -2.3335030452E-01 1.0000000000E+02 1.0000000000E-03 3.6785347929E-02 1.0000000000E+02 1.0000000000E-04 9.9305762210E-01 1.2328283475E+01 1.0000000000E-05 9.9992887890E-01 1.2400273228E+00 1.0000000000E-06 9.9999927161E-01 1.2534628933E-01 1.0000000000E-07 9.9999999271E-01 1.2540753142E-02 1.0000000000E-08 9.9999999993E-01 1.2519811285E-03 1.0000000000E-09 9.9999999948E-01 3.2218252141E-03 1.0000000000E-10 9.9999993592E-01 3.6724432165E-02 1.0000000000E-11 9.9999443483E-01 3.3365409287E-01 1.0000000000E-12 9.9978246313E-01 2.3084666542E+00 1.0000000000E-13 9.7760754200E-01 2.2496550168E+01 1.0000000000E-14 1.5721749789E-01 1.0170130475E+02 Bending angle TL check passed **************************** ********** PASS ********** **************************** *** Results log of t_fascod_tl_2 (FM_TL FASCOD -comp) *** ./t_fascod_tl -comp TL: x%non_ideal, x%new_ref_op, x%new_bangle_op = T F F Checking REFRACTIVITY TL FASCOD_TRO |dy|/|dy_init| cos_angle 100*rel diff 1.0000000000E+00 9.4108577918E-01 7.1826563522E+01 1.0000000000E-01 9.9735869115E-01 9.5839410436E+00 1.0000000000E-02 9.9995924025E-01 1.2329220313E+00 1.0000000000E-03 9.9999954914E-01 1.3089304601E-01 1.0000000000E-04 9.9999999543E-01 1.3182075104E-02 1.0000000000E-05 9.9999999995E-01 1.3191560638E-03 1.0000000000E-06 1.0000000000E+00 1.3192422068E-04 1.0000000000E-07 1.0000000000E+00 1.3190071552E-05 1.0000000000E-08 1.0000000000E+00 1.4334285105E-06 1.0000000000E-09 1.0000000000E+00 3.8695759740E-06 1.0000000000E-10 1.0000000000E+00 4.3150767195E-05 1.0000000000E-11 9.9999999999E-01 4.1191034763E-04 1.0000000000E-12 9.9999999932E-01 3.6928091270E-03 1.0000000000E-13 9.9999992237E-01 3.9420121192E-02 1.0000000000E-14 9.9999263193E-01 3.9036032933E-01 Refractivity TL check passed Checking BENDING ANGLE TL FASCOD_TRO |dy|/|dy_init| cos_angle 100*rel diff 1.0000000000E+00 -1.6890147313E-01 1.0000000012E+02 1.0000000000E-01 -2.0208592729E-01 1.0000000002E+02 1.0000000000E-02 -2.3532244669E-01 1.0000000000E+02 1.0000000000E-03 3.9599044109E-02 1.0000000000E+02 1.0000000000E-04 9.9205842012E-01 1.3082902837E+01 1.0000000000E-05 9.9992709183E-01 1.2556554769E+00 1.0000000000E-06 9.9999925412E-01 1.2657375460E-01 1.0000000000E-07 9.9999999250E-01 1.2686608847E-02 1.0000000000E-08 9.9999999994E-01 1.1271792322E-03 1.0000000000E-09 9.9999999972E-01 2.3946620638E-03 1.0000000000E-10 9.9999997378E-01 2.4875298384E-02 1.0000000000E-11 9.9999882685E-01 1.6859593577E-01 1.0000000000E-12 9.9980223903E-01 2.1187622575E+00 1.0000000000E-13 9.7079258659E-01 2.4218335304E+01 1.0000000000E-14 5.4376840650E-01 9.6334144114E+01 Bending angle TL check passed **************************** ********** PASS ********** **************************** *** Results log of t_fascod_tl_3 (FM_TL FASCOD -new_op) *** ./t_fascod_tl -new_op TL: x%non_ideal, x%new_ref_op, x%new_bangle_op = F T T Checking REFRACTIVITY TL FASCOD_TRO |dy|/|dy_init| cos_angle 100*rel diff 1.0000000000E+00 9.2405443855E-01 8.1301894831E+01 1.0000000000E-01 9.9381907745E-01 1.4437505940E+01 1.0000000000E-02 9.9893100931E-01 5.5723105450E+00 1.0000000000E-03 9.9954291738E-01 3.2463409581E+00 1.0000000000E-04 9.9968808042E-01 2.5651618482E+00 1.0000000000E-05 9.9987105674E-01 1.6323209430E+00 1.0000000000E-06 9.9999116690E-01 4.2614081892E-01 1.0000000000E-07 9.9999986494E-01 5.2662875942E-02 1.0000000000E-08 9.9999999858E-01 5.4004903967E-03 1.0000000000E-09 9.9999999999E-01 5.4147639888E-04 1.0000000000E-10 1.0000000000E+00 6.1695540125E-05 1.0000000000E-11 9.9999999999E-01 3.6705110535E-04 1.0000000000E-12 9.9999999935E-01 3.6126376984E-03 1.0000000000E-13 9.9999995076E-01 3.1695501416E-02 1.0000000000E-14 9.9999968308E-01 7.9886870900E-02 Refractivity TL check passed Checking BENDING ANGLE TL FASCOD_TRO |dy|/|dy_init| cos_angle 100*rel diff 1.0000000000E+00 -1.6721948670E-01 1.0000000012E+02 1.0000000000E-01 -2.0051372942E-01 1.0000000002E+02 1.0000000000E-02 -2.3327210038E-01 1.0000000000E+02 1.0000000000E-03 3.6787542605E-02 1.0000000000E+02 1.0000000000E-04 9.9305794149E-01 1.2328012995E+01 1.0000000000E-05 9.9992888256E-01 1.2399974506E+00 1.0000000000E-06 9.9999927165E-01 1.2534322822E-01 1.0000000000E-07 9.9999999271E-01 1.2540459655E-02 1.0000000000E-08 9.9999999993E-01 1.2521437558E-03 1.0000000000E-09 9.9999999948E-01 3.2204931606E-03 1.0000000000E-10 9.9999993595E-01 3.6716726836E-02 1.0000000000E-11 9.9999443679E-01 3.3359613532E-01 1.0000000000E-12 9.9978239982E-01 2.3087938687E+00 1.0000000000E-13 9.7759867905E-01 2.2507538844E+01 1.0000000000E-14 1.5680524225E-01 1.0171508220E+02 Bending angle TL check passed **************************** ********** PASS ********** **************************** *** Results log of t_fascod_tl_4 (FM_TL FASCOD -comp -new_op) *** ./t_fascod_tl -comp -new_op TL: x%non_ideal, x%new_ref_op, x%new_bangle_op = T T T Checking REFRACTIVITY TL FASCOD_TRO |dy|/|dy_init| cos_angle 100*rel diff 1.0000000000E+00 9.4034133629E-01 7.1435652079E+01 1.0000000000E-01 9.9446369428E-01 1.3478462557E+01 1.0000000000E-02 9.9896083467E-01 5.4839341670E+00 1.0000000000E-03 9.9955254806E-01 3.2120883888E+00 1.0000000000E-04 9.9969203595E-01 2.5496756872E+00 1.0000000000E-05 9.9987273357E-01 1.6219613878E+00 1.0000000000E-06 9.9999124437E-01 4.2428524674E-01 1.0000000000E-07 9.9999986592E-01 5.2470914172E-02 1.0000000000E-08 9.9999999859E-01 5.3814056346E-03 1.0000000000E-09 9.9999999999E-01 5.3996327686E-04 1.0000000000E-10 1.0000000000E+00 6.4462351602E-05 1.0000000000E-11 9.9999999999E-01 3.5578449872E-04 1.0000000000E-12 9.9999999933E-01 3.6732932851E-03 1.0000000000E-13 9.9999995738E-01 2.9292014373E-02 1.0000000000E-14 9.9999775199E-01 2.1503185437E-01 Refractivity TL check passed Checking BENDING ANGLE TL FASCOD_TRO |dy|/|dy_init| cos_angle 100*rel diff 1.0000000000E+00 -1.6869694382E-01 1.0000000012E+02 1.0000000000E-01 -2.0183344664E-01 1.0000000002E+02 1.0000000000E-02 -2.3502070677E-01 1.0000000000E+02 1.0000000000E-03 3.9572117422E-02 1.0000000000E+02 1.0000000000E-04 9.9206770249E-01 1.3074165278E+01 1.0000000000E-05 9.9992718527E-01 1.2547389005E+00 1.0000000000E-06 9.9999925510E-01 1.2648039600E-01 1.0000000000E-07 9.9999999251E-01 1.2677245884E-02 1.0000000000E-08 9.9999999994E-01 1.1265012195E-03 1.0000000000E-09 9.9999999972E-01 2.3943977120E-03 1.0000000000E-10 9.9999997377E-01 2.4891678795E-02 1.0000000000E-11 9.9999882935E-01 1.6851825365E-01 1.0000000000E-12 9.9980255479E-01 2.1175361555E+00 1.0000000000E-13 9.7081929322E-01 2.4208730436E+01 1.0000000000E-14 5.4459222770E-01 9.6272676643E+01 Bending angle TL check passed **************************** ********** PASS ********** **************************** *** Results log of t_fm_1D_1 (FM 1D; default options) *** ./../tools/ropp_fm_bg2ro_1d -f ../data/bgr*_XXXX.nc -o ropp_fm_bgr20090401.nc -d ----------------------------------------------------------------------- ROPP Forward Model ----------------------------------------------------------------------- INFO (from ropp_fm_bg2ro_1d): Processing profile 1 of 1 INFO (from ropp_fm_bg2ro_1d): (BG_20090401000329_UNKN_U999_DMI_) WARNING (from set_obs_levels_bangle): RoC missing from data structure ... setting equal to effective radius INFO (from ropp_fm_bg2ro_1d): Processing profile 1 of 1 INFO (from ropp_fm_bg2ro_1d): (BG_20090401000601_UNKN_U999_DMI_) WARNING (from set_obs_levels_bangle): RoC missing from data structure ... setting equal to effective radius INFO (from ropp_fm_bg2ro_1d): Processing profile 1 of 1 INFO (from ropp_fm_bg2ro_1d): (BG_20090401000630_UNKN_U999_DMI_) WARNING (from set_obs_levels_bangle): RoC missing from data structure ... setting equal to effective radius INFO (from ropp_fm_bg2ro_1d): Processing profile 1 of 1 INFO (from ropp_fm_bg2ro_1d): (BG_20090401000758_UNKN_U999_DMI_) WARNING (from set_obs_levels_bangle): RoC missing from data structure ... setting equal to effective radius INFO (from ropp_fm_bg2ro_1d): Processing profile 1 of 1 INFO (from ropp_fm_bg2ro_1d): (BG_20090401000849_UNKN_U999_DMI_) WARNING (from set_obs_levels_bangle): RoC missing from data structure ... setting equal to effective radius *** Running IDL to plot results IDL Version 8.2 (linux x86_64 m64). (c) 2012, Exelis Visual Information Solutions, Inc. Installation number: 404878-1. Licensed for use by: Met Office Automatic journal file: $HOME/idl_journal_files/idlidculvFriAug1711:57:402018.jou % Compiled module: PLOT_FM. PLOT_FM: Plot forward modelled observables for file: ropp_fm_bgr20090401.nc % Loaded DLM: NCDF. % Compiled module: NCDF_GETVAR. Generating output image file ropp_fm_bgr20090401.jpg % Compiled module: OPEN_JPGFILE. % Compiled module: CLOSE_JPGFILE. % Loaded DLM: JPEG. *** Finished IDL *** ropp_fm_bgr20090401.jpg can be compared against ../data/ropp_fm_bgr20090401_reference.jpg *** Results log of t_fm_1D_2 (FM 1D; compress factors) *** ./../tools/ropp_fm_bg2ro_1d -f -comp ../data/bgr*_XXXX.nc -o ropp_fm_bgr20090401_comp.nc -d ----------------------------------------------------------------------- ROPP Forward Model ----------------------------------------------------------------------- INFO (from ropp_fm_bg2ro_1d): Processing profile 1 of 1 INFO (from ropp_fm_bg2ro_1d): (BG_20090401000329_UNKN_U999_DMI_) WARNING (from set_obs_levels_bangle): RoC missing from data structure ... setting equal to effective radius INFO (from ropp_fm_bg2ro_1d): Processing profile 1 of 1 INFO (from ropp_fm_bg2ro_1d): (BG_20090401000601_UNKN_U999_DMI_) WARNING (from set_obs_levels_bangle): RoC missing from data structure ... setting equal to effective radius INFO (from ropp_fm_bg2ro_1d): Processing profile 1 of 1 INFO (from ropp_fm_bg2ro_1d): (BG_20090401000630_UNKN_U999_DMI_) WARNING (from set_obs_levels_bangle): RoC missing from data structure ... setting equal to effective radius INFO (from ropp_fm_bg2ro_1d): Processing profile 1 of 1 INFO (from ropp_fm_bg2ro_1d): (BG_20090401000758_UNKN_U999_DMI_) WARNING (from set_obs_levels_bangle): RoC missing from data structure ... setting equal to effective radius INFO (from ropp_fm_bg2ro_1d): Processing profile 1 of 1 INFO (from ropp_fm_bg2ro_1d): (BG_20090401000849_UNKN_U999_DMI_) WARNING (from set_obs_levels_bangle): RoC missing from data structure ... setting equal to effective radius *** Running IDL to plot results IDL Version 8.2 (linux x86_64 m64). (c) 2012, Exelis Visual Information Solutions, Inc. Installation number: 404878-1. Licensed for use by: Met Office Automatic journal file: $HOME/idl_journal_files/idlidculvFriAug1711:57:412018.jou % Compiled module: PLOT_FM. PLOT_FM: Plot forward modelled observables for file: ropp_fm_bgr20090401_comp.nc % Loaded DLM: NCDF. % Compiled module: NCDF_GETVAR. Generating output image file ropp_fm_bgr20090401_comp.jpg % Compiled module: OPEN_JPGFILE. % Compiled module: CLOSE_JPGFILE. % Loaded DLM: JPEG. *** Finished IDL *** ropp_fm_bgr20090401_comp.jpg can be compared against ../data/ropp_fm_bgr20090401_comp_reference.jpg *** Results log of t_fm_2D_1 (FM 2D; default options) *** ./../tools/ropp_fm_bg2ro_2d ../data/ECMWF_10_OCCS.nc -o ECMWF_10_OCCS_out.nc -d ----------------------------------------------------------------------- ROPP Forward Model (2D operator) ----------------------------------------------------------------------- INFO: Processing profile 1 of 10 INFO: (OC_20081031215329_744__U999_UNKN) INFO: Processing profile 2 of 10 INFO: (OC_20081101033345_745__U999_UNKN) INFO: Processing profile 3 of 10 INFO: (OC_20081101034701_745__U999_UNKN) INFO: Processing profile 4 of 10 INFO: (OC_20081101042122_745__U999_UNKN) INFO: Processing profile 5 of 10 INFO: (OC_20081101080844_745__U999_UNKN) INFO: Processing profile 6 of 10 INFO: (OC_20081101085533_745__U999_UNKN) INFO: Processing profile 7 of 10 INFO: (OC_20081101014405_741__U999_UNKN) INFO: Processing profile 8 of 10 INFO: (OC_20081101014418_741__U999_UNKN) INFO: Processing profile 9 of 10 INFO: (OC_20081101022155_743__U999_UNKN) INFO: Processing profile 10 of 10 INFO: (OC_20081101031244_743__U999_UNKN) *** Running IDL to plot results IDL Version 8.2 (linux x86_64 m64). (c) 2012, Exelis Visual Information Solutions, Inc. Installation number: 404878-1. Licensed for use by: Met Office Automatic journal file: $HOME/idl_journal_files/idlidculvFriAug1711:57:412018.jou % Compiled module: PLOT_FM_TWOD. PLOT_FM: Plot forward modelled observables for file: ECMWF_10_OCCS_out.nc % Loaded DLM: NCDF. % Compiled module: NCDF_GETVAR. Generating output image file ECMWF_10_OCCS_out.jpg % Compiled module: OPEN_JPGFILE. % Compiled module: CLOSE_JPGFILE. % Loaded DLM: JPEG. *** Finished IDL *** ECMWF_10_OCCS_out.jpg can be compared against ../data/ECMWF_10_OCCS_out_reference.jpg *** Results log of t_fm_2D_2 (FM 2D; compress factors) *** ./../tools/ropp_fm_bg2ro_2d -comp ../data/ECMWF_10_OCCS.nc -o ECMWF_10_OCCS_out_comp.nc -d ----------------------------------------------------------------------- ROPP Forward Model (2D operator) ----------------------------------------------------------------------- INFO: Processing profile 1 of 10 INFO: (OC_20081031215329_744__U999_UNKN) INFO: Processing profile 2 of 10 INFO: (OC_20081101033345_745__U999_UNKN) INFO: Processing profile 3 of 10 INFO: (OC_20081101034701_745__U999_UNKN) INFO: Processing profile 4 of 10 INFO: (OC_20081101042122_745__U999_UNKN) INFO: Processing profile 5 of 10 INFO: (OC_20081101080844_745__U999_UNKN) INFO: Processing profile 6 of 10 INFO: (OC_20081101085533_745__U999_UNKN) INFO: Processing profile 7 of 10 INFO: (OC_20081101014405_741__U999_UNKN) INFO: Processing profile 8 of 10 INFO: (OC_20081101014418_741__U999_UNKN) INFO: Processing profile 9 of 10 INFO: (OC_20081101022155_743__U999_UNKN) INFO: Processing profile 10 of 10 INFO: (OC_20081101031244_743__U999_UNKN) *** Running IDL to plot results IDL Version 8.2 (linux x86_64 m64). (c) 2012, Exelis Visual Information Solutions, Inc. Installation number: 404878-1. Licensed for use by: Met Office Automatic journal file: $HOME/idl_journal_files/idlidculvFriAug1711:57:422018.jou % Compiled module: PLOT_FM_TWOD. PLOT_FM: Plot forward modelled observables for file: ECMWF_10_OCCS_out_comp.nc % Loaded DLM: NCDF. % Compiled module: NCDF_GETVAR. Generating output image file ECMWF_10_OCCS_out_comp.jpg % Compiled module: OPEN_JPGFILE. % Compiled module: CLOSE_JPGFILE. % Loaded DLM: JPEG. *** Finished IDL *** ECMWF_10_OCCS_out_comp.jpg can be compared against ../data/ECMWF_10_OCCS_out_comp_reference.jpg *** Results log of t_fm_iono_1 (FM iono; L_neutral) *** ./../tools/ropp_fm_bg2ro_1d -f ../data/bgr20090401_000329_M02_2030337800_N0007_YYYY.nc -o bgr20090401_000329_M02_2030337800_N0007_YYYY_neut.nc -d ----------------------------------------------------------------------- ROPP Forward Model ----------------------------------------------------------------------- INFO (from ropp_fm_bg2ro_1d): Processing profile 1 of 1 INFO (from ropp_fm_bg2ro_1d): (BG_20090401000329_UNKN_U999_DMI_) *** Results log of t_fm_iono_2 (FM iono; L1 and L2) *** ./../tools/ropp_fm_bg2ro_1d -f -direct_ion ../data/bgr20090401_000329_M02_2030337800_N0007_YYYY.nc -o bgr20090401_000329_M02_2030337800_N0007_YYYY_iono.nc -d ----------------------------------------------------------------------- ROPP Forward Model ----------------------------------------------------------------------- INFO (from ropp_fm_bg2ro_1d): Processing profile 1 of 1 INFO (from ropp_fm_bg2ro_1d): (BG_20090401000329_UNKN_U999_DMI_) INFO (from ropp_fm_iono_set_default): -direct_ion requested but sigma(Ne_max) not set ... defaulting to 200.0e9 m-3. INFO (from ropp_fm_iono_set_default): -direct_ion requested but sigma(H_peak) not set ... defaulting to 150.0 km. INFO (from ropp_fm_iono_set_default): -direct_ion requested but sigma(H_width) not set ... defaulting to 25.0 km. *** Results log of t_iono_1 (FM iono2) *** ./t_iono Ne_max, H_peak, H_width = 3.00000E+11 3.00000E+05 7.50000E+04 RoC, und, lat = 6.35003E+06 1.29449E+01 -3.65360E+01 ii impact bangle_L1(test) bangle_L1(ref) bangle_L2(test) bangle_L2(ref) 1 2365.81 0.28312904E-01 0.28312904E-01 0.28322297E-01 0.28322297E-01 11 3764.98 0.24506589E-01 0.24506589E-01 0.24516059E-01 0.24516059E-01 21 5409.42 0.11645305E-01 0.11645305E-01 0.11654867E-01 0.11654867E-01 31 7178.12 0.11537644E-01 0.11537644E-01 0.11547307E-01 0.11547307E-01 41 8972.60 0.86063310E-02 0.86063310E-02 0.86160992E-02 0.86160992E-02 51 10824.44 0.68233446E-02 0.68233446E-02 0.68332236E-02 0.68332236E-02 61 12702.58 0.56243477E-02 0.56243477E-02 0.56343416E-02 0.56343416E-02 71 14598.90 0.42547598E-02 0.42547598E-02 0.42648727E-02 0.42648727E-02 81 16523.92 0.33219036E-02 0.33219036E-02 0.33321400E-02 0.33321400E-02 91 18463.34 0.24479817E-02 0.24479817E-02 0.24583458E-02 0.24583458E-02 101 20422.41 0.16899809E-02 0.16899809E-02 0.17004771E-02 0.17004771E-02 111 22400.60 0.11921099E-02 0.11921099E-02 0.12027430E-02 0.12027430E-02 121 24390.83 0.84512897E-03 0.84512897E-03 0.85590349E-03 0.85590349E-03 131 26390.37 0.61557073E-03 0.61557073E-03 0.62649114E-03 0.62649114E-03 141 28396.00 0.44076546E-03 0.44076546E-03 0.45183620E-03 0.45183620E-03 151 30407.02 0.32362821E-03 0.32362821E-03 0.33485386E-03 0.33485386E-03 161 32421.82 0.24237560E-03 0.24237560E-03 0.25376081E-03 0.25376081E-03 171 34439.65 0.18108515E-03 0.18108515E-03 0.19263471E-03 0.19263471E-03 181 36460.11 0.13572330E-03 0.13572330E-03 0.14744216E-03 0.14744216E-03 191 38482.89 0.10409341E-03 0.10409341E-03 0.11598673E-03 0.11598673E-03 201 40507.67 0.82296391E-04 0.82296391E-04 0.94369473E-04 0.94369473E-04 211 42534.18 0.66367951E-04 0.66367951E-04 0.78626329E-04 0.78626329E-04 221 44562.34 0.53847123E-04 0.53847123E-04 0.66296544E-04 0.66296544E-04 231 46592.07 0.45112649E-04 0.45112649E-04 0.57759096E-04 0.57759096E-04 241 48623.30 0.39406818E-04 0.39406818E-04 0.52256517E-04 0.52256517E-04 251 50655.94 0.35230215E-04 0.35230215E-04 0.48289641E-04 0.48289641E-04 261 52689.96 0.32089115E-04 0.32089115E-04 0.45365016E-04 0.45365016E-04 271 54725.35 0.29725204E-04 0.29725204E-04 0.43224602E-04 0.43224602E-04 281 56762.09 0.28120508E-04 0.28120508E-04 0.41850718E-04 0.41850718E-04 291 58800.17 0.26963013E-04 0.26963013E-04 0.40931650E-04 0.40931650E-04 Comparing against bending angles in ../data/bgr20090401_000329_M02_2030337800_N 0007_YYYY_iono_ref.nc **************************** ********** PASS ********** **************************** *** Results log of t_iono_ad_1 (FM_AD iono2) *** ./t_iono_ad Ne_max, H_peak, H_width = 3.00000E+11 3.00000E+05 7.50000E+04 RoC, und, lat = 6.35003E+06 1.29449E+01 -3.65360E+01 |dx|/|dx|_init norm1=(Kdx|Kdx)_y norm2=(K^TKdx|dx)_x (norm1-norm2)/norm2 1.0000000000E+00 1.9403327627E-06 1.9403327627E-06 -4.3654004279E-16 1.0000000000E-01 1.7655001277E-08 1.7655001277E-08 1.8740992415E-16 1.0000000000E-02 1.7495993808E-10 1.7495993808E-10 0.0000000000E+00 1.0000000000E-03 1.7480235624E-12 1.7480235624E-12 -2.3105911852E-16 1.0000000000E-04 1.7478661217E-14 1.7478661217E-14 -1.8053119639E-16 1.0000000000E-05 1.7478503791E-16 1.7478503791E-16 4.2312380253E-16 1.0000000000E-06 1.7478488048E-18 1.7478488048E-18 4.4075435795E-16 1.0000000000E-07 1.7478486474E-20 1.7478486474E-20 5.1650905974E-16 1.0000000000E-08 1.7478486316E-22 1.7478486316E-22 -6.7253784426E-16 1.0000000000E-09 1.7478486301E-24 1.7478486301E-24 -2.1016807652E-16 **************************** ********** PASS ********** **************************** *** Results log of t_iono_tl_1 (FM_TL iono2) *** ./t_iono_tl Ne_max, H_peak, H_width = 3.00000E+11 3.00000E+05 7.50000E+04 RoC, und, lat = 6.35003E+06 1.29449E+01 -3.65360E+01 |dx|/|dx|_init |dy=Kdx| (Kdx|H(x+dx)-H(x))/|.||.| |H(x+dx)-H(x)-Kdx|/|H(x+dx)-H(x)| 1.0000000000E+00 3.0402137299E-03 -9.9176755758E-01 2.9294560229E+01 1.0000000000E-01 1.0735125085E-03 -9.9697502043E-01 9.0126511363E+00 1.0000000000E-02 7.9490714808E-04 9.7527854993E-01 2.4560144272E-01 1.0000000000E-03 9.3864010470E-05 9.9999805741E-01 8.3315521206E-02 1.0000000000E-04 8.6807772836E-06 9.9999993471E-01 8.7863610557E-03 1.0000000000E-05 8.6121593486E-07 9.9999999928E-01 8.8218798585E-04 1.0000000000E-06 8.6053275443E-08 9.9999999999E-01 8.8253272338E-05 1.0000000000E-07 8.6046446733E-09 1.0000000000E+00 8.8257196731E-06 1.0000000000E-08 8.6045763880E-10 1.0000000000E+00 8.8249087300E-07 1.0000000000E-09 8.6045694862E-11 1.0000000000E+00 1.2441719449E-07 1.0000000000E-10 8.6045688317E-12 1.0000000000E+00 7.9425641862E-07 1.0000000000E-11 8.6045612046E-13 9.9999999997E-01 7.6992770969E-06 1.0000000000E-12 8.6045454125E-14 9.9999999596E-01 8.9911954737E-05 1.0000000000E-13 8.6054424252E-15 9.9999957831E-01 9.2390587184E-04 1.0000000000E-14 8.6046718036E-16 9.9995429958E-01 9.5603286279E-03 **************************** ********** PASS ********** **************************** *** Results log of t_twodad_1 (FM_AD TWOD) *** ./t_twodad norm1: 5.57977E+02 norm2: 5.57977E+02 norm1/norm2 - 1: 3.48206E-12 **************************** ********** PASS ********** **************************** *** Results log of t_twodad_2 (FM_AD TWOD -comp) *** ./t_twodad -comp norm1: 6.56316E+02 norm2: 6.56316E+02 norm1/norm2 - 1: -2.90749E-12 **************************** ********** PASS ********** **************************** *** Results log of t_twodop_1 (FM TWOD) *** ./t_twodop All fractional differences in bending angle < 1.00000E-04 **************************** ********** PASS ********** **************************** *** Results log of t_twodop_2 (FM TWOD -comp) *** ./t_twodop -comp All fractional differences in bending angle < 1.00000E-04 **************************** ********** PASS ********** **************************** *** Results log of t_twodtl_1 (FM_TL TWOD) *** ./t_twodtl I COS(alpha) SIN(alpha) REL ERR ---------------------------------------------------------------------- 1 -0.154911786950484 0.987928306236747 1.00000000367772 2 -0.171487744573838 0.985186253183111 1.00000000045867 3 -0.457299692900489 0.889312650800110 1.00000000040829 4 -0.503856805871227 0.863787195539066 1.00000000005323 5 0.393631005679664 0.919268530608775 5.68496150900049 6 0.999994288744531 0.003379715715790 0.115913367628674E-01 7 0.999999940388039 0.000345288166711 0.114691994060953E-02 8 0.999999999400268 0.000034633276790 0.114784737498998E-03 9 0.999999999994801 0.000003224681518 0.926060197152841E-05 10 0.999999999988387 0.000004819276044 0.137568214955166E-04 11 0.999999999422730 0.000033978526172 0.207467529756747E-03 12 0.999999952082007 0.000309573875164 0.880552915655245E-03 13 0.999977403878036 0.006722479701958 0.362091586341675E-01 14 0.998393643098385 0.056658039330127 0.264728815636092 15 -0.257266829037393 0.966340405176688 3.70291315432597 **************************** ********** PASS ********** **************************** *** Results log of t_twodtl_2 (FM_TL TWOD -comp) *** ./t_twodtl -comp I COS(alpha) SIN(alpha) REL ERR ---------------------------------------------------------------------- 1 -0.151595428390545 0.988442626605656 1.00000000392314 2 -0.167048985545289 0.985948597254588 1.00000000048457 3 -0.451615203454038 0.892212815425315 1.00000000043730 4 -0.501967461464587 0.864886505635738 1.00000000005751 5 -0.927021770695764 0.375007515466143 2.48952846332201 6 0.999995619171922 0.002960006243847 0.122925432687321E-01 7 0.999999953443850 0.000305143076338 0.125539041934801E-02 8 0.999999999523139 0.000030882389643 0.127974346639637E-03 9 0.999999999994722 0.000003248999926 0.131700206537189E-04 10 0.999999999980238 0.000006286842114 0.330890417155269E-04 11 0.999999999513112 0.000031205397557 0.100885837538163E-03 12 0.999999959720376 0.000283829606328 0.283848434120066E-03 13 0.999985562222865 0.005373578492963 0.303742787489081E-01 14 0.999129619365075 0.041713351668253 0.190715480770831 15 0.374597825705778 0.927187396903400 2.55507624886755 **************************** ********** PASS ********** ****************************