diff --git a/wk/fort.90 b/wk/fort.90
index 84051bf6f9444aabf083cef04812c0e70b3c4887..7b37a727913cb19397745149f5a34cfcc012086a 100644
--- a/wk/fort.90
+++ b/wk/fort.90
@@ -1,24 +1,24 @@
 &BASIC
   nrun=100000
-  dt=0.01
-  tmax=50    ! time normalized to 1/omega_pe
+  dt=0.05
+  tmax=100    ! time normalized to 1/omega_pe
 /
 &GRID
-  pmaxe =81    ! number of Hermite moments 
-  jmaxe = 20     ! number of Hermite moments 
-  pmaxi = 81    ! number of Hermite moments 
-  jmaxi = 20     ! number of Hermite moments
+  pmaxe =20    ! Electron Hermite moments 
+  jmaxe = 10     ! Electron Laguerre moments 
+  pmaxi = 20    ! Ion Hermite moments 
+  jmaxi = 10     ! Ion Laguerre moments
   nkr   = 1
   krmin = 0
-  krmax = 0     ! Normalized to cs0/Omega_i
+  krmax = 0
   nkz   = 1
   kzmin = 1
-  kzmax = 1     ! Normalized to cs0/Omega_i
+  kzmax = 1
 /
 &OUTPUT_PAR
   nsave_0d = 0
   nsave_1d = 0
-  nsave_2d = 100
+  nsave_2d = 1
   nsave_5d = 0
   write_Ni00    = .true.
   write_moments = .true.
@@ -28,16 +28,17 @@
 /
 &MODEL_PAR
   ! Collisionality
-  nu      = 0.001       ! Normalized collision frequency normalized to plasma frequency
+  CO      = -1       ! Collision operator (-1:Full Coulomb, 0: Dougherty)
+  nu      = 0       ! Normalized collision frequency normalized to plasma frequency
   tau_e   = 1         ! T_e/T_e
   tau_i   = 1         ! T_i/T_e       temperature ratio
   sigma_e = 0.023338   ! sqrt(m_e/m_i) mass ratio
   sigma_i = 1         ! sqrt(m_i/m_i)
   q_e     = -1         ! Electrons charge
   q_i     = 1         ! Ions charge
-  eta_n   = 1         ! L_perp / L_n Density gradient
-  eta_T   = 0         ! L_perp / L_T Temperature gradient
-  eta_B   = 0.5         ! L_perp / L_B Magnetic gradient and curvature
+  eta_n   = 0         ! L_perp / L_n Density gradient
+  eta_T   = 1         ! L_perp / L_T Temperature gradient
+  eta_B   = 0         ! L_perp / L_B Magnetic gradient and curvature
   lambdaD = 0         ! Debye length
 /
 &INITIAL_CON