Hello Everybody, I have writte a new subroutine to calculate nusselt number of the ambient air in KIVA, but when I read it into otape9 it prints all zero values. Can anyone tell me how to proceed with this problem.

Also I wanted to know where the time cycle for KIVA code is mentioned. i.e how to know when the code starts calculating for new cycle.

Thanks in advance.

Sasidhar

Hello Sasidhar,

1. [ I have written a new subroutine to calculate nusselt number of the ambient air in KIVA, but when I read it into otape9 it prints all zero values.]

-> The cause may be the writing format of your nusselt number. Or if you are sure nothing`s wrong with the format, the possible reason can be one of the nominators in the expression for Nu is zero, if not extremely small. Try to check the value of each term in the Nu equation by writing it to a new output file.

2. [ how to know when the code starts calculating for new cycle.]

-> In kiva main program, subroutine 'timstp' is called. This subroutine calculates the timestep required for each cycle. You can check that at the end of this subroutine command lines to write the timestep, time, cycle, and crank angle on the monitor are specified. So if you execute the code, you should easily notice those parameters on your screen. You could also write a new tag like 'new cycle' in the main code (between call timstp and call newcyc):

call timstp

20 write(*,*)'new cycle = ',ncyc

call newcyc

Hope this helps

cheers,

zk