A code for simulating the passive containment heat removal system of nuclear power plants was developed.Code development methodology,subroutine dividing,code structure and physics model were investigated.And related test was done to validate that the code calculation result is reliable.The code development method,physics model and numerical solution method would be used as a support for the future design and analysis and engineering application.