POLIPHON 2 - first part

[AthenaAugousta]

Here comes the separation part of the 2nd POLIPHON step. Have a look and let Jonas know if anything needs to be changed!

[HolgerPollyNet]

obviously you have not used the latest dev version, there are things in the dev branch which are not in your code

[HolgerPollyNet]

what is this cryptic line for?

[AthenaAugousta]

obviously you have not used the latest dev version, there are things in the dev branch which are not in your code

I think this happened because we copied the old main file instead of adding the respective POLIPHON 2 part. @jn76lewy could you please just add the POLIPHON 2 lines in the current picassoProcV3.m file and then push this file only again?

[AthenaAugousta]

@HolgerPollyNet could you please check again? @jn76lewy has just pushed the correct version of the main script!

[HolgerPollyNet]

Looks good. I will switch to his branch on rsd2 for some testing. Any objections @ulysses78 ?

[HolgerPollyNet]

I get the following error on rsd2:

Index in position 2 exceeds array bounds. Index must not exceed 4000.

Error in poliphon_two (line 154)
if pdr355_klett(i,n)<dndfd_e(1,n)

Error in picassoProcV3 (line 3292)
[data.POLIPHON2] = poliphon_two ...

Error in picassoProcTodolist (line 62)
reportTmp = picassoProcV3(pollyDataFile, pollyDataTasks.pollyType{iTask}, ...

Any idea?

[AthenaAugousta]

I get the following error on rsd2:

Index in position 2 exceeds array bounds. Index must not exceed 4000.

Error in poliphon_two (line 154) if pdr355_klett(i,n)<dndfd_e(1,n)

Error in picassoProcV3 (line 3292) [data.POLIPHON2] = poliphon_two ...

Error in picassoProcTodolist (line 62) reportTmp = picassoProcV3(pollyDataFile, pollyDataTasks.pollyType{iTask}, ...

Any idea?

n is the height index, which indeed should not exceed 4000. I do not know why it does in this case- which station/date did you try? Maybe @jn76lewy can run locally the same case to examine more closely why it happens

[AthenaAugousta]

I get the following error on rsd2:
Index in position 2 exceeds array bounds. Index must not exceed 4000.
Error in poliphon_two (line 154) if pdr355_klett(i,n)<dndfd_e(1,n)
Error in picassoProcV3 (line 3292) [data.POLIPHON2] = poliphon_two ...
Error in picassoProcTodolist (line 62) reportTmp = picassoProcV3(pollyDataFile, pollyDataTasks.pollyType{iTask}, ...
Any idea?

n is the height index, which indeed should not exceed 4000. I do not know why it does in this case- which station/date did you try? Maybe @jn76lewy can run locally the same case to examine more closely why it happens

@jn76lewy, the issue is on how we currently create the mr variable (line 129: mr = (0:2.500625156289073e-04:1);)
As it is, the size of mr is always 1x4000 and sometimes the height might exceed index 4000. Instead, it should always follow the size of the respective Bsc profile (e.g., aerBsc532_raman).

I quickly tested the following lines, they should work:
step = size(aerBsc355_klett,2)
mr = linspace(0, 1, step);

You could try it locally & if it runs through properly, then push the updated poliphon_2 function only.

[AthenaAugousta]

@HolgerPollyNet Try once more please

[jn76lewy]

i change the code and it worked for me locally

[HolgerPollyNet]

its currently running on rsd2, so far without problems ;-)

[AthenaAugousta]

its currently running on rsd2, so far without problems ;-)

Nice!! You can have a look at the new nc file and tell us if the variables names etc are ok or if something needs to be changed!

[HolgerPollyNet]

I think the updated global config with the new keyword for saving is missing...but I added it manually....

[AthenaAugousta]

Thanks! We will also take care of the visualization together with @ulysses78 soon

[jn76lewy]

I made a first approach with plotting the poliphon 2 data. This is TJK 20:00 - 21:30 (13.09.2025) at 532 nm (355 and 1064 have the same style but in blue/ red). Let me know what you think! :)