calculate current by saved transmission spectrum report

[jaco]

Dear Professor:
     i have got the transmission spectrum report from a .nc document, and named it as .nc document. the origin .nc document can not be used. so i hope i can use this  transmission spectrum report to calculate the current. here i make a attempt to do that, but it does not work. the wrong message is
     Traceback (most recent call last):
  File "c:\docume~1\admini~1\locals~1\temp\0991215239625338.py", line 5, in <module>
    transmission_spectrum = nlread('F:/cnb/cnb/0000.nc', TransmissionSpectrum)

  File ".\zipdir\NL\IO\NLSaveUtilities.py", line 225, in nlread
  File ".\zipdir\NL\IO\NetCDFFile.py", line 156, in __init__
  File ".\zipdir\NL\IO\NetCDFFile.py", line 362, in _read
AssertionError: Error: F:/cnb/cnb/0000.nc is not a valid NetCDF 3 file
     so what can i do next?
your sincerely,
jaco

[Nordland]

I think I know what is the problem, but can you send me the 0000.nc, so I can confirm the problem?

[jaco]

yes, of course, the 0000.nc was saved by ever transmission spectrum in the big .nc document. thanks a lot. this .nc report can hardly be open. I means their are a save button for every transmission spectrum in the big .nc document. so i save it and named it as .nc document.  besides, i doubt that if i get the message from the .log document(as we know, the .log document can save the transmission spectrum), how can i calculate the current from these.
thanks  a lot
Great regards,
jaco

[Nordland]

This is not a netcdf file. This is a txt file saved with the extension .nc

The netcdf file is obtained from the calculation, and is not something that you can obtain from renaming a log file to a nc file.

[jaco]

Dear Norland
I think i got your ideas. but the question become that if i know the Transmission spectrum in a form of .log or .txt. how can i calculate the current?

[kstokbro]

In this case you will need to write you own routine which can integrate the transmission spectrum for the bias window

[Anders Blom]

But there is no need to do that - you obviously have the original NC file with the TransmissionSpectrum. All you need is to do is use your original script but read the real NC file instead. Use

Code: python

transmission_spectrum = nlread('F:/cnb/cnb/calculation.nc', object_id="gID002")[0]
current = transmission_spectrum.current()
print 'current = ', current

where "gID002" should be replaced with the ID of the TransmissionSpectrum that you clicked in VNL to save the data from in the first place.

[jaco]

Thanks all of you,  you kind reply make me feel the sunshine in the cold. maybe i did not express my idea exactly. the origin nc document can not work, it is too big. so i have to use that to calculate the current.
I get Professor kstokbro's ideas, but it is hard for me to do that. So i ask if any script that can be used in ATK and can used to calculate this.Or how can i change the format of saved transmission spectrum into the nc format?
I am appreciated your answer as soon as possible.
Best Regards,
jaco
 

[Anders Blom]

Please specify the symptom of "origin nc document can not work". The file may be big, but what happens when you run the script I suggested? The point of pursuing this is that the manual approach is cumbersome and it's in your own interest to make the proper way work. So it's worthwhile investing time in that.