Jinlun, to make sure we will be talking about same code, please look at: http://mitgcm.org/cgi-bin/viewcvs.cgi/MITgcm/pkg/seaice/ http://mitgcm.org/cgi-bin/viewcvs.cgi/MITgcm/pkg/seaice/growth.F http://mitgcm.org/cgi-bin/viewcvs.cgi/*checkout*/MITgcm/pkg/seaice/growth.F?rev=1.32