C INCLUDE FILE FOR XIAM H.Hartwig Feb. 1994 - Jul. 1996 C C Important array dimensions, to be set before compiling C C maximal J quantum number integer DIMJ c parameter (DIMJ=30) c parameter (DIMJ=50) zk XIAM3 parameter (DIMJ=70) zk XIAM4 C maximal number of v-blocks in the last diagonalisation. C the size of the matrix is (2*DIMJ+1)*DIMVV integer DIMVV C parameter (DIMVV=3) zk parameter (DIMVV=8) C maximal number of internal rotors (max. 3, see MAXTOP) integer DIMTOP C parameter (DIMTOP=3) parameter (DIMTOP=1) zk XIAM3 C maximal internal rotation symmetry number sigma for one top integer DIMSIG parameter (DIMSIG=1) C maximal number of basis funcions for the pure internal rotation C Hamiltonian integer DIMM C parameter (DIMM=31) c parameter (DIMM=51) zk XIAM3 parameter (DIMM=71) zk XIAM4 C maximum number of v basis funktions for one top integer DIMV c parameter (DIMV=3) parameter (DIMV=8) zk C maximum number of transitions integer DIMLIN c parameter (DIMLIN=500) parameter (DIMLIN=3000) zk XIAM3 C maximum number of different sets of paramaters integer DIMVB parameter(DIMVB=5) C ------------------------------------------------------------------- C please change the following parameters only if you really know what you C are doing ! C ------------------------------------------------------------------- C max. allowed number of tops integer MAXTOP C parameter (MAXTOP=3) parameter (MAXTOP=3) zk integer DIMTOT,DIM2J1 parameter (DIM2J1=2*DIMJ+1) parameter (DIMTOT=DIM2J1 * DIMVV) integer DIMGAM parameter (DIMGAM=(3**DIMTOP)/2+1) integer gamma(DIMGAM,0:DIMTOP) C array dimensions for data integer DIMPRR, DIMPIR, DIMPAR, DIMPM, MAXPAR parameter (DIMPRR=32) parameter (DIMPIR=13) parameter (DIMPAR=DIMPRR + DIMPIR*DIMTOP) parameter (MAXPAR=DIMPRR + DIMPIR*MAXTOP) C see also DIMLMF (=DIMFIT) in 'iamfit.f' integer DIMFIT, DIMPLC parameter (DIMFIT=60, DIMPLC=DIMTOP*DIMVB) C Quantum Number Array for transitions integer Q_UP, Q_LO parameter (Q_UP=1) parameter (Q_LO=2) integer Q_K,Q_T,Q_V1,Q_V2,Q_B,Q_K2,Q_GK,Q_J,Q_S,Q_STAT,Q_F integer Q_REF,Q_AVG,Q_V,Q_TJ parameter (Q_J=1) parameter (Q_K=2) parameter (Q_V1=3) parameter (Q_V=Q_V1-1) parameter (Q_V2=4) parameter (Q_B=5) parameter (Q_S=6) parameter (Q_F=7) parameter (Q_T=8) parameter (Q_TJ=9) parameter (Q_STAT=10) parameter (Q_K2=11) parameter (Q_GK=12) parameter (Q_REF=13) parameter (Q_AVG=14) integer DIMQLP, DIMTDO parameter (DIMQLP=14) C parameter (DIMTDO=200) C parameter (DIMTDO=500) c parameter (DIMTDO=1000) zk parameter (DIMTDO=3000) zk XIAM4 integer qlin(DIMLIN, DIMQLP, Q_UP:Q_LO) integer todo(DIMTDO, DIMQLP) integer ncalc integer LN_FREQ, LN_ERR, LN_PSI, LN_INT parameter (LN_FREQ=1) parameter (LN_ERR=2) parameter (LN_PSI=3) parameter (LN_INT=4) ! Intensity integer DIMDLN parameter (DIMDLN=LN_INT) real*8 dln(DIMLIN, DIMDLN) integer NV_ENG, NV_DEF, NV_DEL parameter (NV_ENG=1) ! energyparameter of one niveau parameter (NV_DEF=1) ! first derivativ of energy in array dnv parameter (NV_DEL=DIMFIT+NV_ENG) ! last derivativ of energy integer DIMDNV parameter (DIMDNV=NV_DEL) real*8 dnv(DIMLIN, DIMDNV, Q_UP:Q_LO) integer NaQN,BINUNIT,UI,GU parameter (NaQN=-12345,BINUNIT=11,UI=5,GU=1) real*8 NOFIT parameter (NOFIT=1.0d12) character*3 BINSTR parameter (BINSTR='iam') integer S_H, S_K, S_VV integer S_G integer S_V, S_S, S_MINV, S_MAXV integer S_MAXK,S_NB integer S_FIRV,S_MAXM parameter (S_H =1) parameter (S_K =2) parameter (S_VV =3) parameter (S_MAXK=4) parameter (S_NB =5) parameter (S_G =6) parameter (S_V =6) ! S_V+itop= 7, 8, 9 parameter (S_S =S_V+1*DIMTOP) ! S_S+itop=10,11,12 parameter (S_MINV=S_V+2*DIMTOP) ! 13,14,15 parameter (S_MAXV=S_V+3*DIMTOP) ! 16,17,18 parameter (S_FIRV=S_V+4*DIMTOP) ! 19,20,21 parameter (S_MAXM=S_V+5*DIMTOP) ! 22,23,24 integer DIMSIZ parameter (DIMSIZ=S_MAXM+6*DIMTOP) integer size(DIMSIZ) C array used in IAMM parameter (DIMPM =4) C for DIMOVV, for DIMPM integer DIMOVV parameter (DIMOVV=08) integer PM_F, PM_VN1, PM_VN2, PM_RHO C integer PM_DPI4, PM_DPIC integer PM_PI, PM_COS, PM_SIN C integer PM_M, PM_RK integer PM_PI2 parameter (PM_F =1) parameter (PM_VN1 =2) parameter (PM_VN2 =3) parameter (PM_RHO =4) c parameter (PM_DPI4=5) c parameter (PM_DPIC=6) C DIMENSION DIMPM parameter (PM_PI =5) parameter (PM_COS =6) parameter (PM_SIN =7) parameter (PM_PI2 =8) C DIMENSION DIMOVV integer PI_VN1, PI_VN2, PI_F, PI_RHO C integer PI_DPI4,PI_DPIC integer PI_BETA,PI_GAMA integer PI_DPIJ,PI_DPIK,PI_DPID integer PI_DC3J integer PI_F0, PI_ANGX,PI_ANGZ parameter (PI_F = 1) parameter (PI_VN1 = 2) parameter (PI_VN2 = 3) parameter (PI_RHO = 4) c parameter (PI_DPI4= 5) c parameter (PI_DPIC= 6) parameter (PI_BETA= 5) parameter (PI_GAMA= 6) parameter (PI_DC3J= 7) parameter (PI_DPIJ= 8) parameter (PI_DPIK= 9) parameter (PI_DPID=10) parameter (PI_F0 =11) parameter (PI_ANGX=12) parameter (PI_ANGZ=13) C ******** DIMENSION: DIMPIR integer P_BJ, P_BK, P_BD, P_DJ, P_DJK, P_DK integer P_DJD, P_DKD, P_R6, P_FF integer P_QZ , P_QD, P_QXY, P_QXZ, P_QYZ integer P_VCC, P_VSS integer P_HJ, P_HJK, P_HKJ, P_HK, P_HJD, P_HJKD,P_HKD integer P_CP, P_CZ, P_CD integer P_MUX, P_MUY, P_MUZ integer P_PX, P_PY, P_PZ integer P1_VN1, P1_VN2, P1_F, P1_RHO, P1_BETA,P1_GAMA integer P1_DC3J,P1_DPIJ,P1_DPIK,P1_DPID integer P1_F0, P1_ANGX,P1_ANGZ integer P2_VN1, P2_VN2, P2_F, P2_RHO, P2_BETA,P2_GAMA integer P2_DC3J,P2_DPIJ,P2_DPIK,P2_DPID integer P2_F0, P2_ANGX,P2_ANGZ integer P3_VN1, P3_VN2, P3_F, P3_RHO, P3_BETA,P3_GAMA integer P3_DC3J,P3_DPIJ,P3_DPIK,P3_DPID integer P3_F0, P3_ANGX,P3_ANGZ parameter (P_BJ = 1) parameter (P_BK = 2) parameter (P_BD = 3) parameter (P_DJ = 4) parameter (P_DJK = 5) parameter (P_DK = 6) parameter (P_DJD = 7) parameter (P_DKD = 8) C parameter (P_R6 = 9) parameter (P_FF =09) parameter (P_QZ =10) parameter (P_QD =11) parameter (P_QXY =12) parameter (P_QXZ =13) parameter (P_QYZ =14) parameter (P_VSS =15) parameter (P_VCC =16) parameter (P_HJ =17) parameter (P_HJK =18) parameter (P_HKJ =19) parameter (P_HK =20) parameter (P_HJD =21) parameter (P_HJKD=22) parameter (P_HKD =23) parameter (P_CP =24) parameter (P_CZ =25) parameter (P_CD =26) parameter (P_MUX =27) parameter (P_MUY =28) parameter (P_MUZ =29) parameter (P_PX =30) parameter (P_PY =31) parameter (P_PZ =32) C ******** DIMENSION: DIMPRR parameter (P1_VN1 = DIMPRR+PI_VN1 ) parameter (P1_VN2 = DIMPRR+PI_VN2 ) parameter (P1_F = DIMPRR+PI_F ) parameter (P1_RHO = DIMPRR+PI_RHO ) c parameter (P1_DPI4= DIMPRR+PI_DPI4) c parameter (P1_DPIC= DIMPRR+PI_DPIC) parameter (P1_BETA= DIMPRR+PI_BETA) parameter (P1_GAMA= DIMPRR+PI_GAMA) parameter (P1_DPIJ= DIMPRR+PI_DPIJ) parameter (P1_DPIK= DIMPRR+PI_DPIK) parameter (P1_DPID= DIMPRR+PI_DPID) parameter (P1_DC3J= DIMPRR+PI_DC3J) parameter (P1_F0 = DIMPRR+PI_F0 ) parameter (P1_ANGX= DIMPRR+PI_ANGX) parameter (P1_ANGZ= DIMPRR+PI_ANGZ) parameter (P2_VN1 = DIMPRR+DIMPIR+PI_VN1 ) parameter (P2_VN2 = DIMPRR+DIMPIR+PI_VN2 ) parameter (P2_F = DIMPRR+DIMPIR+PI_F ) parameter (P2_RHO = DIMPRR+DIMPIR+PI_RHO ) c parameter (P2_DPI4= DIMPRR+DIMPIR+PI_DPI4) c parameter (P2_DPIC= DIMPRR+DIMPIR+PI_DPIC) parameter (P2_BETA= DIMPRR+DIMPIR+PI_BETA) parameter (P2_GAMA= DIMPRR+DIMPIR+PI_GAMA) parameter (P2_DPIJ= DIMPRR+DIMPIR+PI_DPIJ) parameter (P2_DPIK= DIMPRR+DIMPIR+PI_DPIK) parameter (P2_DPID= DIMPRR+DIMPIR+PI_DPID) parameter (P2_DC3J= DIMPRR+DIMPIR+PI_DC3J) parameter (P2_F0 = DIMPRR+DIMPIR+PI_F0 ) parameter (P2_ANGX= DIMPRR+DIMPIR+PI_ANGX) parameter (P2_ANGZ= DIMPRR+DIMPIR+PI_ANGZ) parameter (P3_VN1 = DIMPRR+2*DIMPIR+PI_VN1 ) parameter (P3_VN2 = DIMPRR+2*DIMPIR+PI_VN2 ) parameter (P3_F = DIMPRR+2*DIMPIR+PI_F ) parameter (P3_RHO = DIMPRR+2*DIMPIR+PI_RHO ) c parameter (P3_DPI4= DIMPRR+2*DIMPIR+PI_DPI4) c parameter (P3_DPIC= DIMPRR+2*DIMPIR+PI_DPIC) parameter (P3_BETA= DIMPRR+2*DIMPIR+PI_BETA) parameter (P3_GAMA= DIMPRR+2*DIMPIR+PI_GAMA) parameter (P3_DPIJ= DIMPRR+2*DIMPIR+PI_DPIJ) parameter (P3_DPIK= DIMPRR+2*DIMPIR+PI_DPIK) parameter (P3_DPID= DIMPRR+2*DIMPIR+PI_DPID) parameter (P3_DC3J= DIMPRR+2*DIMPIR+PI_DC3J) parameter (P3_F0 = DIMPRR+2*DIMPIR+PI_F0 ) parameter (P3_ANGX= DIMPRR+2*DIMPIR+PI_ANGX) parameter (P3_ANGZ= DIMPRR+2*DIMPIR+PI_ANGZ) C control integers integer C_NZYK, C_PRINT,C_EVAL, C_DFRQ, C_ORGER, C_INTS integer C_APR, C_PRI ,C_XPR ,C_NCYCL,C_RED integer C_WOODS,C_ADJF integer C_NDATA,C_NFOLD,C_SPIN, C_NTOP integer C_SVDER,C_LAST, C_FITSC integer C_MAXM, C_MAXV C integer C_MAXM1,C_MAXM2,C_MAXM3,C_MAXM C integer C_MAXV1,C_MAXV2,C_MAXV3,C_MAXV C integer C_WOOD1,C_WOOD2,C_WOOD3,C_WOOD4,C_WOODS C integer C_ADJ1 ,C_ADJ2 ,C_ADJ3 ,C_ADJ4 ,C_ADJF C parameter (C_MAXM1= 7, C_MAXM2= 8, C_MAXM3= 9) C parameter (C_WOOD1=10, C_WOOD2=11, C_WOOD3=12, C_WOOD4=13 ) C parameter (C_MAXV1=16, C_MAXV2=17, C_MAXV3=18) C parameter (C_ADJ1 =28, C_ADJ2 =29, C_ADJ3 =30, C_ADJ4 =31) C parameter (C_MAXM=C_MAXM1-1) C parameter (C_MAXV=C_MAXV1-1) parameter (C_NZYK = 1, C_PRINT= 2, C_EVAL = 3) parameter (C_DFRQ = 4, C_ORGER= 5, C_INTS = 6) parameter (C_MAXM =07) parameter (C_WOODS=08, C_NDATA=09) parameter (C_NFOLD=10, C_SPIN =11, C_NTOP =12) parameter (C_ADJF =13) parameter (C_MAXV =14) parameter (C_PRI =15, C_XPR=16, C_NCYCL=17) parameter (C_SVDER=18, C_FITSC=19) parameter (C_RED =20) parameter (C_LAST =20) integer C_ROFIT, C_EPS, C_DEFER, C_WEIGF, C_CNVG integer C_LMBDA, C_FRQLO,C_FRQUP,C_INTLM, C_TEMP parameter (C_ROFIT=C_LAST+1, C_EPS =C_LAST+2) parameter (C_DEFER=C_LAST+3, C_WEIGF=C_LAST+4) parameter (C_CNVG =C_LAST+5, C_LMBDA=C_LAST+6) parameter (C_FRQLO=C_LAST+7) parameter (C_FRQUP=C_LAST+8, C_INTLM=C_LAST+9) parameter (C_TEMP=C_LAST+10) integer DIMCINT, DIMCPAR parameter (DIMCINT=C_LAST) parameter (DIMCPAR=10) integer CTLINT(1:DIMCINT) real*8 CTLPAR(DIMCINT+1:DIMCINT+DIMCPAR) integer qvv(DIMVV,DIMTOP,DIMVB) C woods(1-DIMVB) adjust(1-DIMVB) integer CB_WDS, CB_ADJ parameter (CB_WDS=1, CB_ADJ=2) integer ctlnb(2,DIMVB) integer AP_PL ! Parameter List (0..3) integer AP_TF ! List of Transitions zero cycle integer AP_TL ! List of Transitions integer AP_TE ! List of Transitions extended integer AP_PC ! additional Parameter information integer AP_IO ! Input control integer AP_LT ! Latex Output (at the end only) integer AP_SV ! SVD-Information integer AP_ST ! status integer AP_TI ! Torsional Integrals < v K | v' K' > integer AP_RM ! Rotation Matrix D integer AP_EO ! Eigenenergies of one top operator integer AP_MO ! Matrixelements of one top operator integer AP_EH ! Eigenenergies / vectors integer AP_MH ! Matrixelements integer XP_FI ! first cycle integer XP_LA ! last cycle integer XP_CC ! every converged cycle integer XP_EC,xcy ! every cycle integer XP_DE,xde ! for every dqu derivativ parameter(AP_PL= 7 ) parameter(AP_TF= 8, AP_TL= 16, AP_TE= 32) parameter(AP_PC= 64, AP_IO= 128, AP_LT= 256) parameter(AP_SV= 512, AP_ST= 1024, AP_TI= 2048) parameter(AP_RM= 4096, AP_EO= 8192, AP_MO= 16384) parameter(AP_EH= 32768, AP_MH= 65536) parameter(XP_FI= 1, XP_LA= 2, XP_CC= 4) parameter(XP_EC= 8, XP_DE= 16) common /data/qlin,dln,dnv,ctlpar,ctlint $ ,ctlnb $ ,size,gamma $ ,todo,ncalc,qvv,xcy,xde integer DIMQC,MAXQC parameter (MAXQC=30) parameter (DIMQC=MAXQC+DIMGAM)