Opened 3 months ago

Last modified 3 months ago

#223 accepted Support

ASCII output type UNFORMATTED

Reported by: scloux Owned by: pesei
Priority: major Milestone:
Component: FP coding/compilation Version: FLEXPART-WRF
Keywords: ASCII format Cc:

Description

Hello, my name is Sara and I have a question about how to save data in ASCII format. I am running flexpart-wrf 3.3 and the first time I wanted to obtain ASCII outputs I found this error message:

At line 287 of file concoutput_irreg.f90 (unit = 97, file = '/home/usc/fp/scg/lustre/PROYECTOS/FLEXPART/TEST/test_damian/output_back_2/grid_time_19821014170000_001')
Fortran runtime error: Format present for UNFORMATTED data transfer
srun: error: c7026: task 0: Exited with exit code 2
srun: Terminating job step 2165917.0

I looked for that line in the code and, from my point of view, it was difficult the program save ASCII outputs providing that this format was declared as UNFORMATTED. There is a high probability I had been doing something wrong, but I try to solve this incident by changing the concoutput_irreg.f90 file. I would like to share this new file in case this problem persists.

In other terms, I am wondering if it is possible to obtain by any way the partposition file in latlon coordinates if we are working with a Lambert projection and there is not a regular grid for mother domain.

Thank you for your time.
Regards, Sara

Attachments (1)

concoutput_irreg.f90 (28.4 KB) - added by scloux 3 months ago.

Download all attachments as: .zip

Change History (3)

Changed 3 months ago by scloux

comment:1 Changed 3 months ago by pesei

  • Owner set to pesei
  • Priority changed from blocker to major
  • Status changed from new to accepted

comment:2 Changed 3 months ago by pesei

ASCII output is not programmed into FLEXPART because it is less efficient, and as the model may produce huge amounts of data, it would not be practical to produce ASCII output files. It is alwaus possible to postprocess the binary files in order to obtain ASCII output.

If you just want to obtain some x.y.z positions of particles, you can simply add corresponding write statements. You don't need to rewrite all the standard output code.

The line listed in the quoted error message does not correspond to the posted source file as line 287 in this file is "else", has nothing to do with i/o,

Note: See TracTickets for help on using tickets.
hosted by ZAMG