This paper presents a numerical calculation method for the static performance of gas journal bearing based on finite element analysis(FEM) by python programming. The dimensionless form of compressible Reynolds lubrication equation of gas journal bearing is derived and is transformed into standard elliptic partial differential equation(PDE) by substitution method. Based on Galerkin method, the finite element method program is developed to calculation the pressure distribution and static performance of the gas bearing through python code. The FEM program of bearing calculation includes three parts: preprocessing, calculation and post-processing. The class of FEM model some useful class methods are built by object-oriented programming to keep and deal with the nodes and element information through out the three parts of the programming, which makes the program modularization and turns the calculation process in several procedures of model operation. After the solving of Reynolds equation, we can go further to calculate the other static performance such as static bearing capacity and friction torque.