set output frequency using usnpst.f90
Posted: Thu May 27, 2010 7:40 pm
Dear all,
I would like to set my output frequency without the GUI using usnpst.f90. I am using a variable in time and uniform in space time step so I would like to control my output frequency using ttcabs (in order to have an output every second between 300 and 500s for example). Consequently, I put the following in usnpst.f90:
call pstntc(ntcabs) ! indwri = 0 indact = 0 ! time = 0.d0 ! do i = 2, 252, 50 time = dble(i) if (ttcabs.ge.time.and.ttpabs.lt.time) then indwri = 0 indact = 1 endif enddo ! do i = 300, 500 time = dble(i) if (ttcabs.ge.time.and.ttpabs.lt.time) then indwri = 0 indact = 1 endif enddo ! call pstact(indwri, indact)
But the result I get is that the output files are written when the physical time step is almost equal to 2 (that's what I want) but then, I got an output at every time step and not when the physical time is almost equal to 52.
Please let me know if you have any suggestions, my being to have an output every 50 s between 2s and 252s and then an output every second between 300s and 500s.
Regards.
Nicolas
I would like to set my output frequency without the GUI using usnpst.f90. I am using a variable in time and uniform in space time step so I would like to control my output frequency using ttcabs (in order to have an output every second between 300 and 500s for example). Consequently, I put the following in usnpst.f90:
call pstntc(ntcabs) ! indwri = 0 indact = 0 ! time = 0.d0 ! do i = 2, 252, 50 time = dble(i) if (ttcabs.ge.time.and.ttpabs.lt.time) then indwri = 0 indact = 1 endif enddo ! do i = 300, 500 time = dble(i) if (ttcabs.ge.time.and.ttpabs.lt.time) then indwri = 0 indact = 1 endif enddo ! call pstact(indwri, indact)
But the result I get is that the output files are written when the physical time step is almost equal to 2 (that's what I want) but then, I got an output at every time step and not when the physical time is almost equal to 52.
Please let me know if you have any suggestions, my being to have an output every 50 s between 2s and 252s and then an output every second between 300s and 500s.
Regards.
Nicolas