SUBROUTINE LIQC2H6(WAVELN,XNR,XNI) * * * PERS. COMM E. QUIRICO ; eric.quirico@obs.ujf-grenoble.fr * * DIMENSION W(154),XN(154),XK(154) data W/ & 1000. ,13.7832,12.5514,11.7883,11.1533,10.5841,10.0702, & 9.6039, 9.1789, 8.7899, 8.4325, 8.1031, 7.7984, 7.5158, & 7.2530, 7.0079, 6.7789, 6.5643, 6.3629, 6.1383, 5.8776, & 5.6862, 5.3808, 4.4945, 3.8012, 3.6938, 3.6004, 3.5311, & 3.4720, 3.4149, 3.3595, 3.3060, 3.2541, 3.2038, 3.1551, & 3.1078, 3.0620, 3.0167, 2.9724, 2.9286, 2.8869, 2.8469, & 2.8075, 2.7699, 2.7360, 2.7037, 2.6689, 2.6345, 2.5991, & 2.5653, 2.5339, 2.5033, 2.4734, 2.4443, 2.4158, 2.3880, & 2.3608, 2.3343, 2.3083, 2.2829, 2.2580, 2.2337, 2.2099, & 2.1866, 2.1638, 2.1373, 2.1135, 2.0920, 2.0707, 2.0502, & 2.0301, 2.0105, 1.9912, 1.9722, 1.9537, 1.9354, 1.9175, & 1.9000, 1.8827, 1.8658, 1.8492, 1.8328, 1.8168, 1.8010, & 1.7855, 1.7702, 1.7553, 1.7405, 1.7261, 1.7118, 1.6978, & 1.6840, 1.6704, 1.6571, 1.6440, 1.6310, 1.6183, 1.6058, & 1.5934, 1.5813, 1.5693, 1.5575, 1.5457, 1.5342, 1.5229, & 1.5118, 1.5009, 1.4901, 1.4795, 1.4690, 1.4586, 1.4485, & 1.4384, 1.4285, 1.4187, 1.4091, 1.3996, 1.3902, 1.3809, & 1.3718, 1.3628, 1.3539, 1.3451, 1.3365, 1.3279, 1.3195, & 1.3111, 1.3029, 1.2947, 1.2867, 1.2788, 1.2709, 1.2632, & 1.2556, 1.2480, 1.2405, 1.2332, 1.2259, 1.2187, 1.2116, & 1.2045, 1.1976, 1.1907, 1.1839, 1.1772, 1.1705, 1.1640, & 1.1575, 1.1510, 1.1447, 1.1384, 1.1322, 1.1282, 0.0100 &/ data XN/ & 1.3267, 1.3272, 1.3294, 1.2944, 1.3174, 1.3197, 1.3211, & 1.3217, 1.3225, 1.3229, 1.3235, 1.3242, 1.3255, 1.3286, & 1.3266, 1.3408, 1.3029, 1.3101, 1.3149, 1.3171, 1.3188, & 1.3196, 1.3203, 1.3226, 1.3273, 1.3298, 1.3334, 1.3432, & 1.3569, 1.3637, 1.3311, 1.2866, 1.3017, 1.3067, 1.3099, & 1.3120, 1.3133, 1.3139, 1.3148, 1.3153, 1.3159, 1.3162, & 1.3167, 1.3170, 1.3371, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, & 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600, 1.3600 &/ data XK/ &.131E-02,.394E-03,.452E-03,.987E-03,.109E-02,.118E-02,.106E-02, &.107E-02,.102E-02,.992E-03,.997E-03,.855E-03,.849E-03,.751E-03, &.149E-02,.125E-02,.713E-02,.148E-02,.815E-03,.213E-03,.642E-04, &.965E-04,.375E-04,.468E-04,.863E-04,.186E-05,.451E-04,.136E-02, &.428E-02,.237E-01,.162E-01,.273E-02,.144E-02,.416E-03,.155E-03, &.227E-03,.361E-03,.692E-04,.509E-04,.548E-04,.105E-03,.662E-04, &.144E-03,.133E-03,.945E-04,.395E-04,.148E-03,.455E-05,.559E-04, &.187E-04,.929E-04,.348E-04,.896E-04,.480E-03,.471E-03,.285E-03, &.126E-03,.446E-03,.886E-03,.843E-03,.281E-03,.625E-04,.487E-04, &.118E-04,.437E-04,.186E-05,.336E-04,.223E-05,.101E-04,.182E-04, &.436E-04,.171E-04,.127E-04,.186E-04,.181E-04,.243E-04,.252E-04, &.330E-04,.294E-04,.302E-04,.258E-04,.276E-04,.335E-04,.295E-04, &.373E-04,.351E-04,.322E-04,.653E-04,.670E-04,.882E-04,.124E-03, &.261E-04,.973E-05,.144E-04,.356E-05,.228E-05,.475E-05,.414E-05, &.103E-05,.980E-06,.413E-06,.283E-06,.312E-06,.996E-06,.188E-05, &.262E-05,.282E-05,.537E-05,.262E-05,.169E-05,.314E-05,.627E-05, &.147E-04,.975E-05,.866E-05,.161E-04,.225E-04,.315E-04,.206E-04, &.248E-04,.157E-04,.850E-05,.155E-05,.833E-06,.523E-06,.461E-06, &.638E-06,.959E-06,.112E-05,.130E-05,.119E-05,.120E-05,.149E-05, &.180E-05,.183E-05,.157E-05,.190E-05,.264E-05,.178E-05,.339E-05, &.119E-04,.117E-04,.218E-04,.345E-04,.150E-04,.159E-04,.969E-05, &.748E-05,.879E-05,.872E-05,.444E-05,.269E-05,.241E-05,.212E-05 &/ *234567890123456789012345678901234567890123456789012345678901234567890 * * XNR=XN(1) XNI=XK(1) IF (WAVELN .GT. W(1)) RETURN XNR=XN(154) XNI=XK(154) IF (WAVELN .LT. W(154)) RETURN DO 100 I=2,154 IF (WAVELN .GT. W(I) ) GO TO 101 100 CONTINUE 101 CONTINUE C ALL INTERPOLATION IS IN LOG LAMBDA FACTOR= (alog(WAVELN) - alog(W(I)) ) / (alog(W(I-1)) - alog(W(I))) C REAL PART IS LINEARLY INTERPOLATED XNR=XN(I) + FACTOR*(XN(I-1) - XN(I)) C IMAGINARY PART IS LOG INTERPOLATED XNI=alog(XK(I)) + FACTOR*(alog(XK(I-1)) - alog(XK(I))) XNI=exp(XNI) iF(XNI.lt.1.e-7) XNI=1.e-7 RETURN END