BISMG:SarahS/gauss newton method

SEt up a matrix $$ J_{i,j} = [d smb(1000m) /d lapse(1000m) ]$$