f.write(seq_2d_write(seq_2d_xx(instance.md_2d,seq)))
#---v
if instance.type == 'remd' and instance.use_nmr and instance.nmr_format == 'txt':
+ fnmr = open(instance.jobdirname+'/plik3.nmr', 'r')
+ nmr_txt_gamatheta=False
+ for linmr in fnmr:
+ if 'phi' in linmr or 'psi' in linmr :
+ f.write('@include gamatheta.txt\n')
+ nmr_txt_gamatheta=True
+ break
+ fnmr.close()
+
f.write('ndist='+repr(instance.linijek_nmr)+' restr_type=12 scal_peak=')
f.write(repr(instance.nmr_scal)+' fordepth_peak='+repr(instance.nmr_fordepth))
f.write(' slope_peak='+repr(instance.nmr_slope)+'\n')
#---v
if instance.use_nmr and instance.nmr_format == 'txt':
+ if nmr_txt_gamatheta:
+ f.write('@include gamatheta.txt\n')
f.write('ndist='+repr(instance.linijek_nmr)+' restr_type=12 scal_peak=')
f.write(repr(instance.nmr_scal)+' fordepth_peak='+repr(instance.nmr_fordepth))
f.write(' slope_peak='+repr(instance.nmr_slope)+'\n')
#---v
if instance.use_nmr and instance.nmr_format == 'txt':
+ if nmr_txt_gamatheta:
+ f.write('@include gamatheta.txt\n')
f.write('ndist='+repr(instance.linijek_nmr)+' restr_type=12 scal_peak=')
f.write(repr(instance.nmr_scal)+' fordepth_peak='+repr(instance.nmr_fordepth))
f.write(' slope_peak='+repr(instance.nmr_slope)+'\n')