*.s3d files

61 views
Skip to first unread message

Leobardo Valera

unread,
Oct 19, 2021, 2:58:55 PM10/19/21
to fds...@googlegroups.com
Hello everyone,

I would like to duplicate the information stored on the *.s3d files but on text files. Could any of you let me know where that writing is done to add a new line to write the same information but in ASCII format?

Thanks in advance,

 Leobardo Valera  

dr_jfloyd

unread,
Oct 19, 2021, 3:12:28 PM10/19/21
to FDS and Smokeview Discussions
Outputs are processed and written in dump.f90.

Leobardo Valera

unread,
Oct 19, 2021, 3:17:18 PM10/19/21
to fds...@googlegroups.com
Thanks for your reply Dr. Floyd. I understand that the information is written in thos line
WRITE(FN_SMOKE3D(N,NM),  '(A,A,I4.4,A,I2.2,A)') TRIM(CHID),'_',NM,'_',N,'.s3d'
How can I modify the line to write the same information on an ASCII file?

Thanks again


--
You received this message because you are subscribed to the Google Groups "FDS and Smokeview Discussions" group.
To unsubscribe from this group and stop receiving emails from it, send an email to fds-smv+u...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/fds-smv/18cc2587-3b2d-4418-b29d-16672aa9be7fn%40googlegroups.com.

Kevin McGrattan

unread,
Oct 19, 2021, 4:43:32 PM10/19/21
to fds...@googlegroups.com
The data you want is found in the array QQ in these lines from dump.f90

   ! Pack the data into a 1-D array and call routine that writes the file

    ALLOCATE(QQ_PACK(IBP1*JBP1*KBP1))
   QQ_PACK = PACK(QQ(0:IBAR,0:JBAR,0:KBAR,1),MASK=.TRUE.)
   CALL SMOKE3D_TO_FILE(LU_SMOKE3D(N,NM),LU_SMOKE3D(N+4,NM),STIME,DXX,MASS_EXT_COEF,N,QQ_PACK,IBP1,JBP1,KBP1,&
                        HRRPUV_MAX_SMV,TEMP_MAX_SMV)

Leobardo Valera

unread,
Oct 19, 2021, 4:45:20 PM10/19/21
to fds...@googlegroups.com
Thanks a lot Dr. MacGrattan

--
You received this message because you are subscribed to the Google Groups "FDS and Smokeview Discussions" group.
To unsubscribe from this group and stop receiving emails from it, send an email to fds-smv+u...@googlegroups.com.

Leobardo Valera

unread,
Oct 21, 2021, 2:53:40 AM10/21/21
to fds...@googlegroups.com
Hello,
I tried to replicate these lines
CALL SMOKE3D_TO_FILE(LU_SMOKE3D(N,NM),LU_SMOKE3D(N+4,NM),STIME,DXX,MASS_EXT_COEF,N,QQ_PACK,IBP1,JBP1,KBP1,&
                        HRRPUV_MAX_SMV,TEMP_MAX_SMV)
but using a formatted file to have the information written in an ascii file. I did not get any compilation error, but when I checked the files
they were empty. 

After some experiments, I realized two things:
1. When I do have only one mesh, everything is ok
2. When I have more than one mesh the number is not consistent.

Any suggestions? Thanks in advance,

Leobardo Valera




On Tue, Oct 19, 2021 at 2:43 PM Kevin McGrattan <mcgr...@gmail.com> wrote:
--
You received this message because you are subscribed to the Google Groups "FDS and Smokeview Discussions" group.
To unsubscribe from this group and stop receiving emails from it, send an email to fds-smv+u...@googlegroups.com.

Kevin McGrattan

unread,
Oct 21, 2021, 9:19:57 AM10/21/21
to fds...@googlegroups.com
All you need to do is add a line before the allocate statement like this

WRITE(999,*) QQ

which will dump the entire array into a text file. You could also do this

DO K=0,KBP1
   DO J=0,JBP1
      DO I=0,IBP1
         WRITE(999,*) I,J,K,QQ(I,J,K,1)
      ENDDO
   ENDDO
ENDDO

The resulting file will be called fort.999. If you have multiple meshes, adjust the 999 to account for different mesh numbers.

Leobardo Valera

unread,
Oct 21, 2021, 9:25:15 AM10/21/21
to fds...@googlegroups.com
Thanks a lot. I will try it and let you know.

Best regards

--
You received this message because you are subscribed to the Google Groups "FDS and Smokeview Discussions" group.
To unsubscribe from this group and stop receiving emails from it, send an email to fds-smv+u...@googlegroups.com.

Leobardo Valera

unread,
Oct 21, 2021, 1:27:04 PM10/21/21
to fds...@googlegroups.com
Thanks again for the help. It is working well, so I have just one more question: when I check a file corresponding to a mesh, I noticed that I have (KBP1+1)*(JBP1+1)*(IBP1+1)*(a_number) lines. May I assume that (a_number) is the number of t_times?
Best regards

Leobardo Valera



On Thu, Oct 21, 2021 at 7:19 AM Kevin McGrattan <mcgr...@gmail.com> wrote:
--
You received this message because you are subscribed to the Google Groups "FDS and Smokeview Discussions" group.
To unsubscribe from this group and stop receiving emails from it, send an email to fds-smv+u...@googlegroups.com.

Kevin McGrattan

unread,
Oct 21, 2021, 1:50:06 PM10/21/21
to fds...@googlegroups.com
Yes
Reply all
Reply to author
Forward
0 new messages