_RL area_reg_sq, hice_reg_sq
area_reg_sq = SEAICE_area_reg * SEAICE_area_reg
& + area_reg_sq)
tmpscal1 = SQRT(AREApreTH(I,J)* AREApreTH(I,J) + area_reg_sq)
& /SQRT( AREAITD(i,j,k,bi,bj)**2 + area_reg_sq )