[ncl-talk] Error when writing out a text file with write_table: "fatal:process_list: can not print list in list yet."

[Maxime Colin]

Dear NCL-talk list,

I'm trying to output an ASCII .txt file with some basic stats on my data. The idea is to print a list of variables (vertically), for various time periods (horizontally).

I've already done this successfully in the past, so I don't know what I'm doing wrong this time.

The header and footer print out well, but not the data, with should have 14 rows and 4 columns. To try to fix the problem I also made a test with meaningless data that I have commented out, such as [/11,12,13,14/]. But this test fails as well.

I'm receiving the following error, and I don't see what's wrong with my list:

"
fatal:process_list: can not print list in list yet.

fatal:process_list: can not print list in list yet.
"

The relevant bit of code looks like this:

  variable_string_list = [/"PW","CAPE","CIN","LCL","LFC","temp_500hPa","rh_500hPa","mse_500hPa","qvapor_500hPa","theta_get_subcloud","rh_subcloud","mse_subcloud","qvapor_subcloud","temp_subcloud"/]
  ;variable_string_list = [/"s1","s2","s3","s4"/]

  all_day_average      = [/pw_manual_members_ave,cape_0D_members_ave,cin_0D_members_ave,LCL_members_ave,LFC_members_ave,temp_500hPa_members_ave,rh_500hPa_members_ave,mse_500hPa_members_ave,qvapor_500hPa_members_ave,theta_get_subcloud_layer_ave_members_ave,rh_subcloud_layer_ave_members_ave,mse_subcloud_layer_ave_members_ave,qvapor_subcloud_layer_ave_members_ave,temp_subcloud_layer_ave_members_ave/]
  ;all_day_average      = [/11,12,13,14/]

  exp_growth_days_ave  = [/avg(pw_manual_members(0:last_growth_day-1)),avg(cape_0D_members(0:last_growth_day-1)),avg(cin_0D_members(0:last_growth_day-1)),avg(LCL_members(0:last_growth_day-1)),avg(LFC_members(0:last_growth_day-1)),avg(temp_500hPa_members(0:last_growth_day-1)),avg(rh_500hPa_members(0:last_growth_day-1)),avg(mse_500hPa_members(0:last_growth_day-1)),avg(qvapor_500hPa_members(0:last_growth_day-1)),avg(theta_get_subcloud_layer_ave_members(0:last_growth_day-1)),avg(rh_subcloud_layer_ave_members(0:last_growth_day-1)),avg(mse_subcloud_layer_ave_members(0:last_growth_day-1)),avg(qvapor_subcloud_layer_ave_members(0:last_growth_day-1)),avg(temp_subcloud_layer_ave_members(0:last_growth_day-1))/]
  ;exp_growth_days_ave  = [/5,6,7,8/]

  lin = num_day_max-1
  exp_decay_days_ave   = [/avg(pw_manual_members(last_growth_day:lin)),avg(cape_0D_members(last_growth_day:lin)),avg(cin_0D_members(last_growth_day:lin)),avg(LCL_members(last_growth_day:lin)),avg(LFC_members(last_growth_day:lin)),avg(temp_500hPa_members(last_growth_day:lin)),avg(rh_500hPa_members(last_growth_day:lin)),avg(mse_500hPa_members(last_growth_day:lin)),avg(qvapor_500hPa_members(last_growth_day:lin)),avg(theta_get_subcloud_layer_ave_members(last_growth_day:lin)),avg(rh_subcloud_layer_ave_members(last_growth_day:lin)),avg(mse_subcloud_layer_ave_members(last_growth_day:lin)),avg(qvapor_subcloud_layer_ave_members(last_growth_day:lin)),avg(temp_subcloud_layer_ave_members(last_growth_day:lin))/]
  print ("pw_manual_members(last_growth_day:lin) : " + pw_manual_members(last_growth_day:lin))
  ;exp_decay_days_ave   = [/1,2,3,4/]



  alist = [/variable_string_list, all_day_average, exp_growth_days_ave, exp_decay_days_ave/]

  dir_for_plots_idea = data_dir_for_plots(0)
  header = (/"--------------------------------", \
            "Stats about the target profiles, depending on exp growth or decay", \
            alist_string, \
            "for the targets computed with " + dir_for_plots_idea, \
            "--------------------------------"/)
  footer = (/"--------------------------------", \
            "This is a file footer", \
            "--------------------------------"/)

  hlist = [/header/]
  flist = [/footer/]



  write_table("./C02_target_profile_analysis_without_getvar/stats_for_all_day_and_exp_growth_days_and_exp_decay_days.txt", "w", hlist, "%s")
  write_table("./C02_target_profile_analysis_without_getvar/stats_for_all_day_and_exp_growth_days_and_exp_decay_days.txt", "a", alist, "%s%10.10f%10.10f%10.10f")
  write_table("./C02_target_profile_analysis_without_getvar/stats_for_all_day_and_exp_growth_days_and_exp_decay_days.txt", "a", flist, "%s")



Does anyone have a clue on where the problem comes from?

What else can I test to identify the culprit?

Thank you very much for any suggestions.

Maxime.




Maxime Colin
---------------------------------------------
<https://www.normalesup.org/phare/squirrelmail/src/compose.php?send_to=maxime.colin%40normalesup.org>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200915/6e3d2a6e/attachment.html>