#!/usr/bin/env looper
# Calculation of the universal conductance curve G(T/T_K, B/T_K)
#PRELUDE: $Nz=8;
#AUTOLOOP: nrginit ; nrgrun

[extra]
Nz=8

B=1.05e-4

[param]
symtype=ISOSZ
model=SIAM
variant=MAGFIELD

U=1.0
Gamma=0.05
delta=0.0

discretization=Z
Lambda=2.0
@$z = 1/$Nz; $z <= 1.00001; $z += 1/$Nz
z=$z
!print "$cnt $z\n";

Tmin=1e-11

# High-accuracy calculation
keepmin=500
keepenergy=12
keep=10000

ops=A_d n_d SZd
specgt=A_d-A_d-u A_d-A_d-d

strategy=kept
safeguard=1e-4

gtp=0.7

calc0=true

log=c

tri=cpp
prec=100
preccpp=2000