+ call reada(controlcard,"D0CM",d0cm,3.78d0)
+ call reada(controlcard,"AKCM",akcm,15.1d0)
+ call reada(controlcard,"AKTH",akth,11.0d0)
+ call reada(controlcard,"AKCT",akct,12.0d0)
+ call reada(controlcard,"V1SS",v1ss,-1.08d0)
+ call reada(controlcard,"V2SS",v2ss,7.61d0)
+ call reada(controlcard,"V3SS",v3ss,13.7d0)
+ call reada(controlcard,"EBR",ebr,-5.50D0)
+ call reada(controlcard,"ATRISS",atriss,0.301D0)
+ call reada(controlcard,"BTRISS",btriss,0.021D0)
+ call reada(controlcard,"CTRISS",ctriss,1.001D0)
+ call reada(controlcard,"DTRISS",dtriss,1.001D0)
+ call reada(controlcard,"SSSCALE",ssscale,1.0D0)
+ dyn_ss=(index(controlcard,'DYN_SS').gt.0)
+