================================================================== lmdif_s.f:L817 FDJAC2() ------------------------------------------------------------------ double precision aaa FJAC(I,J) = (WA(I) - FVEC(I))/H aaa = (dble(WA(I)) - dble(FVEC(I)))/dble(H) print *, j, i, wa(i), fvec(i), h, FJAC(I,J) print *, j, i, dble(wa(i)), dble(fvec(i)), dble(h), aaa ------------------------------------------------------------------ PGI 2 3 -0.1004647 -0.1004647 3.4526672E-04 8.6316810E-05 2 3 -0.1004647016525269 -0.1004647314548492 3.4526671515777707E-004 8.6316812711229632E-005 ------------------------------------------------------------------ INTEL 2 3 -0.1004647 -0.1004647 3.4526672E-04 4.3158405E-05 2 3 -0.100464716553688 -0.100464731454849 3.452667151577771E-004 4.315840635561482E-005 ================================================================== ================================================================== lmdif_s.f:L813 FDJAC2() ------------------------------------------------------------------ CALL FCN(M,N,X,WA,IFLAG) print *, j,WA print *, j,dble(WA) ------------------------------------------------------------------ PGI 2 -8.6558852E-03 4.5458078E-03 -0.1004647 2 -8.6558852344751358E-003 4.5458078384399414E-003 -0.1004647016525269 ------------------------------------------------------------------ INTEL ================================================================== 2 -8.6558852E-03 4.5458078E-03 -0.1004647 2 -8.655885234475136E-003 4.545807838439941E-003 -0.100464716553688 ================================================================== ================================================================== aerosol_SGLI_3ch.3param.QA.201807.f:L1463 fn_land() ------------------------------------------------------------------ F_Vec(3) = si380 - Obs_R380 print *, si380, Obs_R380, F_Vec(3) print *, dble(si380), dble(Obs_R380), dble(F_Vec(3)) ------------------------------------------------------------------ PGI 0.2059973 0.3064620 -0.1004647 0.2059973180294037 0.3064620196819305 -0.1004647016525269 ------------------------------------------------------------------ INTEL 0.2059973 0.3064620 -0.1004647 0.205997332930565 0.306462049484253 -0.100464716553688 ================================================================== ================================================================== aerosol_SGLI_3ch.3param.QA.201807.f:L1404 fn_land() ------------------------------------------------------------------ si380 = LIP_3param ( $ LUT_TC_CAI( 1, 1 , 1, IX38), $ Tau_A, CMF, Ref_R, dTAi, dCMFi, dRefi_R, $ E_Tau_a, E_CMF, E_Ref_R, ix,iy,iz, n_Tau_A, n_CMF, n_M) print *, si380, ATMOS_WITH_BRDF3 print *, dble(si380), dble(ATMOS_WITH_BRDF3) si380 = si380 + ATMOS_WITH_BRDF3 print *, si380 print *, dble(si380) ------------------------------------------------------------------ PGI 0.1912796 1.4717661E-02 0.1912796497344971 1.4717660844326019E-002 0.2059973 0.2059973180294037 ----------------------------------------------------------------- INTEL 0.1912797 1.4717662E-02 0.191279664635658 1.471766177564859E-002 0.2059973 0.205997332930565 ================================================================== ================================================================== globals_cai_parasol.f:L1360 main() ------------------------------------------------------------------ do i=1,21 write(*,*) 'OBS_CAI = ', i, OBS_CAI(i), dble(OBS_CAI(i)) enddo do i=2,20 write(*,*) 'dir_info= ', i, dir_info(Num_of_dir,i), $ dble(dir_info(Num_of_dir,i)) enddo write(*,*) 'xTau_Aero = ', xTau_Aero write(*,*) 'xCMF= ', xCMF write(*,*) 'xRef_R= ', xRef_R write(*,*) 'xQA_Flag= ', xQA_Flag ------------------------------------------------------------------ PGI OBS_CAI = 1 28.83000 28.82999992370605 OBS_CAI = 2 27.70000 27.69999885559082 OBS_CAI = 3 131.6400 131.6399993896484 OBS_CAI = 4 98.81000 98.80999755859375 OBS_CAI = 5 0.3064620 0.3064620196819305 OBS_CAI = 6 0.1609124 0.1609123796224594 OBS_CAI = 7 0.3931367 0.3931366801261902 OBS_CAI = 8 0.1934716 0.1934715658426285 OBS_CAI = 9 1.6325001E-02 1.6325000673532486E-002 OBS_CAI = 10 1.9993749E-02 1.9993748515844345E-002 OBS_CAI = 11 0.3675187 0.3675187230110168 OBS_CAI = 12 0.000000 0.000000000000000 OBS_CAI = 13 0.000000 0.000000000000000 OBS_CAI = 14 0.000000 0.000000000000000 OBS_CAI = 15 0.000000 0.000000000000000 OBS_CAI = 16 0.000000 0.000000000000000 OBS_CAI = 17 0.000000 0.000000000000000 OBS_CAI = 18 0.000000 0.000000000000000 OBS_CAI = 19 1.000000 1.000000000000000 OBS_CAI = 20 396.0000 396.0000000000000 OBS_CAI = 21 20.00000 20.00000000000000 dir_info= 2 28.83000 28.82999992370605 dir_info= 3 56.52000 56.52000045776367 dir_info= 4 81.62000 81.61999511718750 dir_info= 5 0.5516459 0.5516458749771118 dir_info= 6 0.8760543 0.8760543465614319 dir_info= 7 115.1289 115.1288681030273 dir_info= 8 57.18775 57.18774795532227 dir_info= 9 0.2436072 0.2436071783304214 dir_info= 10 0.3842940 0.3842939734458923 dir_info= 11 2.2175968E-02 2.2175967693328857E-002 dir_info= 12 4.0993690E-03 4.0993690490722656E-003 dir_info= 13 -7.4311495E-03 -7.4311494827270508E-003 dir_info= 14 -1.6677380E-04 -1.6677379608154297E-004 dir_info= 15 9.600676 9.600675582885742 dir_info= 16 1.067610 1.067609786987305 dir_info= 17 2.3387935E-02 2.3387935012578964E-002 dir_info= 18 4.1027600E-03 4.1027599945664406E-003 dir_info= 19 0.000000 0.000000000000000 dir_info= 20 0.000000 0.000000000000000 xTau_Aero = 0.7377689 xCMF= 1.000000 xRef_R= 1.400000 xQA_Flag= 0.5035146 ------------------------------------------------------------------ INTEL OBS_CAI = 1 28.83000 28.8299999237061 OBS_CAI = 2 27.70000 27.6999988555908 OBS_CAI = 3 131.6400 131.639999389648 OBS_CAI = 4 98.81000 98.8099975585938 OBS_CAI = 5 0.3064620 0.306462049484253 OBS_CAI = 6 0.1609124 0.160912394523621 OBS_CAI = 7 0.3931367 0.393136709928513 OBS_CAI = 8 0.1934716 0.193471580743790 OBS_CAI = 9 1.6325001E-02 1.632500067353249E-002 OBS_CAI = 10 1.9993749E-02 1.999374851584435E-002 OBS_CAI = 11 0.3675187 0.367518723011017 OBS_CAI = 12 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 13 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 14 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 15 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 16 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 17 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 18 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 19 1.000000 1.00000000000000 OBS_CAI = 20 396.0000 396.000000000000 OBS_CAI = 21 20.00000 20.0000000000000 dir_info= 2 28.83000 28.8299999237061 dir_info= 3 56.52000 56.5200004577637 dir_info= 4 81.62000 81.6199951171875 dir_info= 5 0.5516459 0.551645874977112 dir_info= 6 0.8760543 0.876054346561432 dir_info= 7 115.1289 115.128868103027 dir_info= 8 57.18775 57.1877479553223 dir_info= 9 0.2436072 0.243607178330421 dir_info= 10 0.3842940 0.384293973445892 dir_info= 11 2.2175968E-02 2.217596769332886E-002 dir_info= 12 4.0993690E-03 4.099369049072266E-003 dir_info= 13 -7.4311495E-03 -7.431149482727051E-003 dir_info= 14 -1.6677380E-04 -1.667737960815430E-004 dir_info= 15 9.600676 9.60067558288574 dir_info= 16 1.067610 1.06760978698730 dir_info= 17 2.3387935E-02 2.338793501257896E-002 dir_info= 18 4.1027600E-03 4.102759994566441E-003 dir_info= 19 0.0000000E+00 0.000000000000000E+000 dir_info= 20 0.0000000E+00 0.000000000000000E+000 xTau_Aero = 0.4521422 xCMF= 1.000000 xRef_R= 1.400000 xQA_Flag= 0.3174434 ================================================================== ================================================================== globals_cai_parasol.f:L988 main() ------------------------------------------------------------------ OBS_CAI(i)=HDF_ref(k,m,i) / cos( OBS_CAI(1) * rad ) print *, i,OBS_CAI(i),HDF_ref(k,m,i), $ rad, OBS_CAI(1) *rad, $ cos( OBS_CAI(1) *rad ) print *, i,dble(OBS_CAI(i)),dble(HDF_ref(k,m,i)), $ dble(rad), dble(OBS_CAI(1) * (rad)), $ dble(cos( OBS_CAI(1) *rad )) ------------------------------------------------------------------ PGI 5 0.3064620 0.2684774 1.7453292E-02 0.5031784 0.8760543 5 0.3064620196819305 0.2684773802757263 1.7453292384743690E-002 0.5031784176826477 0.8760543465614319 6 0.1609124 0.1409680 1.7453292E-02 0.5031784 0.8760543 6 0.1609123796224594 0.1409679949283600 1.7453292384743690E-002 0.5031784176826477 0.8760543465614319 7 0.3931367 0.3444091 1.7453292E-02 0.5031784 0.8760543 7 0.3931366801261902 0.3444091081619263 1.7453292384743690E-002 0.5031784176826477 0.8760543465614319 8 0.1934716 0.1694916 1.7453292E-02 0.5031784 0.8760543 8 0.1934715658426285 0.1694916039705276 1.7453292384743690E-002 0.5031784176826477 0.8760543465614319 ----------------------------------------------------------------- INTEL 5 0.3064620 0.2684774 1.7453292E-02 0.5031784 0.8760543 5 0.306462049484253 0.268477380275726 1.745329238474369E-002 0.503178417682648 0.876054286956787 6 0.1609124 0.1409680 1.7453292E-02 0.5031784 0.8760543 6 0.160912394523621 0.140967994928360 1.745329238474369E-002 0.503178417682648 0.876054286956787 7 0.3931367 0.3444091 1.7453292E-02 0.5031784 0.8760543 7 0.393136709928513 0.344409108161926 1.745329238474369E-002 0.503178417682648 0.876054286956787 8 0.1934716 0.1694916 1.7453292E-02 0.5031784 0.8760543 8 0.193471580743790 0.169491603970528 1.745329238474369E-002 0.503178417682648 0.876054286956787 ================================================================== ================================================================== globals_cai_parasol.f:L988 main() ------------------------------------------------------------------ OBS_CAI(i)=HDF_ref(k,m,i) / cos( OBS_CAI(1) * rad ) print *, i,OBS_CAI(i),HDF_ref(k,m,i), $ cos( OBS_CAI(1) *rad ) OBS_CAI(i)=HDF_ref(k,m,i) / dcos( dble(OBS_CAI(1) * rad )) print *, i,dble(OBS_CAI(i)),dble(HDF_ref(k,m,i)), $ dcos( dble(OBS_CAI(1) *rad )) ------------------------------------------------------------------ PGI 5 0.3064620 0.2684774 0.8760543 5 0.3064620196819305 0.2684773802757263 0.8760543170321176 6 0.1609124 0.1409680 0.8760543 6 0.1609123945236206 0.1409679949283600 0.8760543170321176 7 0.3931367 0.3444091 0.8760543 7 0.3931367099285126 0.3444091081619263 0.8760543170321176 8 0.1934716 0.1694916 0.8760543 8 0.1934715658426285 0.1694916039705276 0.8760543170321176 ----------------------------------------------------------------- INTEL 5 0.3064620 0.2684774 0.8760543 5 0.306462019681931 0.268477380275726 0.876054317032118 6 0.1609124 0.1409680 0.8760543 6 0.160912394523621 0.140967994928360 0.876054317032118 7 0.3931367 0.3444091 0.8760543 7 0.393136709928513 0.344409108161926 0.876054317032118 8 0.1934716 0.1694916 0.8760543 8 0.193471565842628 0.169491603970528 0.876054317032118 ================================================================== ================================================================== globals_cai_parasol.f:L988 main() ------------------------------------------------------------------ OBS_CAI(i)=HDF_ref(k,m,i) / cos( dble(OBS_CAI(1) * rad)) ... do i=1,21 write(*,*) 'OBS_CAI = ', i, OBS_CAI(i), dble(OBS_CAI(i)) enddo do i=2,20 write(*,*) 'dir_info= ', i, dir_info(Num_of_dir,i), $ dble(dir_info(Num_of_dir,i)) enddo write(*,*) 'xTau_Aero = ', xTau_Aero write(*,*) 'xCMF= ', xCMF write(*,*) 'xRef_R= ', xRef_R ------------------------------------------------------------------ PGI OBS_CAI = 1 28.83000 28.82999992370605 OBS_CAI = 2 27.70000 27.69999885559082 OBS_CAI = 3 131.6400 131.6399993896484 OBS_CAI = 4 98.81000 98.80999755859375 OBS_CAI = 5 0.3064620 0.3064620196819305 OBS_CAI = 6 0.1609124 0.1609123945236206 OBS_CAI = 7 0.3931367 0.3931367099285126 OBS_CAI = 8 0.1934716 0.1934715658426285 OBS_CAI = 9 1.6325001E-02 1.6325000673532486E-002 OBS_CAI = 10 1.9993749E-02 1.9993748515844345E-002 OBS_CAI = 11 0.3675187 0.3675187230110168 OBS_CAI = 12 0.000000 0.000000000000000 OBS_CAI = 13 0.000000 0.000000000000000 OBS_CAI = 14 0.000000 0.000000000000000 OBS_CAI = 15 0.000000 0.000000000000000 OBS_CAI = 16 0.000000 0.000000000000000 OBS_CAI = 17 0.000000 0.000000000000000 OBS_CAI = 18 0.000000 0.000000000000000 OBS_CAI = 19 1.000000 1.000000000000000 OBS_CAI = 20 396.0000 396.0000000000000 OBS_CAI = 21 20.00000 20.00000000000000 dir_info= 2 28.83000 28.82999992370605 dir_info= 3 56.52000 56.52000045776367 dir_info= 4 81.62000 81.61999511718750 dir_info= 5 0.5516459 0.5516458749771118 dir_info= 6 0.8760543 0.8760543465614319 dir_info= 7 115.1289 115.1288681030273 dir_info= 8 57.18775 57.18774795532227 dir_info= 9 0.2436072 0.2436071783304214 dir_info= 10 0.3842940 0.3842939734458923 dir_info= 11 2.2175968E-02 2.2175967693328857E-002 dir_info= 12 4.0993690E-03 4.0993690490722656E-003 dir_info= 13 -7.4311495E-03 -7.4311494827270508E-003 dir_info= 14 -1.6677380E-04 -1.6677379608154297E-004 dir_info= 15 9.600676 9.600675582885742 dir_info= 16 1.067610 1.067609786987305 dir_info= 17 2.3387935E-02 2.3387935012578964E-002 dir_info= 18 4.1027600E-03 4.1027599945664406E-003 dir_info= 19 0.000000 0.000000000000000 dir_info= 20 0.000000 0.000000000000000 xTau_Aero = 0.7377689 xCMF= 1.000000 xRef_R= 1.400000 xQA_Flag= 0.5035146 ----------------------------------------------------------------- INTEL OBS_CAI = 1 28.83000 28.8299999237061 OBS_CAI = 2 27.70000 27.6999988555908 OBS_CAI = 3 131.6400 131.639999389648 OBS_CAI = 4 98.81000 98.8099975585938 OBS_CAI = 5 0.3064620 0.306462019681931 OBS_CAI = 6 0.1609124 0.160912394523621 OBS_CAI = 7 0.3931367 0.393136709928513 OBS_CAI = 8 0.1934716 0.193471565842628 OBS_CAI = 9 1.6325001E-02 1.632500067353249E-002 OBS_CAI = 10 1.9993749E-02 1.999374851584435E-002 OBS_CAI = 11 0.3675187 0.367518723011017 OBS_CAI = 12 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 13 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 14 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 15 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 16 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 17 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 18 0.0000000E+00 0.000000000000000E+000 OBS_CAI = 19 1.000000 1.00000000000000 OBS_CAI = 20 396.0000 396.000000000000 OBS_CAI = 21 20.00000 20.0000000000000 dir_info= 2 28.83000 28.8299999237061 dir_info= 3 56.52000 56.5200004577637 dir_info= 4 81.62000 81.6199951171875 dir_info= 5 0.5516459 0.551645874977112 dir_info= 6 0.8760543 0.876054346561432 dir_info= 7 115.1289 115.128868103027 dir_info= 8 57.18775 57.1877479553223 dir_info= 9 0.2436072 0.243607178330421 dir_info= 10 0.3842940 0.384293973445892 dir_info= 11 2.2175968E-02 2.217596769332886E-002 dir_info= 12 4.0993690E-03 4.099369049072266E-003 dir_info= 13 -7.4311495E-03 -7.431149482727051E-003 dir_info= 14 -1.6677380E-04 -1.667737960815430E-004 dir_info= 15 9.600676 9.60067558288574 dir_info= 16 1.067610 1.06760978698730 dir_info= 17 2.3387935E-02 2.338793501257896E-002 dir_info= 18 4.1027600E-03 4.102759994566441E-003 dir_info= 19 0.0000000E+00 0.000000000000000E+000 dir_info= 20 0.0000000E+00 0.000000000000000E+000 xTau_Aero = 0.4519970 xCMF= 1.000000 xRef_R= 1.400000 xQA_Flag= 0.3173314 ================================================================== ================================================================== globals_cai_parasol.f:L288 main() ------------------------------------------------------------------ pi = acos(-1.0) rad = pi / 180.0 deg = 180.0 / pi print *, pi, rad print *, dble(pi), dble(rad) ------------------------------------------------------------------ PGI 3.141593 1.7453292E-02 57.29578 3.141592741012573 1.7453292384743690E-002 57.29577636718750 ----------------------------------------------------------------- INTEL 3.141593 1.7453292E-02 57.29578 3.14159274101257 1.745329238474369E-002 57.2957763671875 ================================================================== ================================================================== aerosol_SGLI_3ch.3param.QA.201807.f:L1864 LIP_3param() ------------------------------------------------------------------ print *, Ref_R_IP(2), Ref_R_IP(1), dRef_R print *, dble(Ref_R_IP(2)), dble(Ref_R_IP(1)), dble(dRef_R) LIP_3param = (Ref_R_IP(2)-Ref_R_IP(1)) * dRef_R + Ref_R_IP(1) print *, LIP_3param print *, dble(LIP_3param) ------------------------------------------------------------------ PGI 0.1912796 0.1912796 0.000000 0.1912796497344971 0.1912796497344971 0.000000000000000 0.1912796 0.1912796497344971 ----------------------------------------------------------------- INTEL 0.1912797 0.1912797 0.0000000E+00 0.191279664635658 0.191279664635658 0.000000000000000E+000 0.1912797 0.191279664635658 ================================================================== ================================================================== aerosol_SGLI_3ch.3param.QA.201807.f:L1856 LIP_3param() ------------------------------------------------------------------ print *, D(1,1,iz ) print *, dble(D(1,1,iz )) print *, z1 print *, dble(z1) print *, z2 print *, dble(z2) print *, dz1i print *, dble(dz1i) print *, dz2i print *, dble(dz2i) print *, x, y print *, dble(x), dble(y) Ref_R_IP(1) = bi_dim_IP2( D(1,1,iz ), z1, z2, dz1i, dz2i, $ x, y, ix, iy, n, n2) ------------------------------------------------------------------ PGI 0.1912796 0.1912796497344971 0.000000 5.0000001E-02 0.1000000 0.2000000 0.4000000 0.8000000 1.600000 2.400000 3.200000 4.800000 6.400000 0.000000000000000 5.0000000745058060E-002 0.1000000014901161 0.2000000029802322 0.4000000059604645 0.8000000119209290 1.600000023841858 2.400000095367432 3.200000047683716 4.800000190734863 6.400000095367432 --- 0.000000 0.2500000 0.5000000 0.7000000 0.9000000 0.9600000 1.000000 0.000000000000000 0.2500000000000000 0.5000000000000000 0.6999999880790710 0.8999999761581421 0.9599999785423279 1.000000000000000 --- 20.00000 20.00000 10.00000 5.000000 2.500000 1.250000 1.250000 1.250000 0.6249999 0.6250001 20.00000000000000 20.00000000000000 10.00000000000000 5.000000000000000 2.500000000000000 1.250000000000000 1.249999880790710 1.250000119209290 0.6249999403953552 0.6250000596046448 --- 4.000000 4.000000 5.000000 5.000000 16.66667 24.99999 4.000000000000000 4.000000000000000 5.000000476837158 5.000000476837158 16.66666603088379 24.99998664855957 --- 0.000000 0.000000 0.000000000000000 0.000000000000000 ----------------------------------------------------------------- INTEL 0.1912797 0.191279664635658 --- 0.0000000E+00 5.0000001E-02 0.1000000 0.2000000 0.4000000 0.8000000 1.600000 2.400000 3.200000 4.800000 6.400000 0.000000000000000E+000 5.000000074505806E-002 0.100000001490116 0.200000002980232 0.400000005960464 0.800000011920929 1.60000002384186 2.40000009536743 3.20000004768372 4.80000019073486 6.40000009536743 --- 0.0000000E+00 0.2500000 0.5000000 0.7000000 0.9000000 0.9600000 1.000000 0.000000000000000E+000 0.250000000000000 0.500000000000000 0.699999988079071 0.899999976158142 0.959999978542328 1.00000000000000 --- 20.00000 20.00000 10.00000 5.000000 2.500000 1.250000 1.250000 1.250000 0.6249999 0.6250001 20.0000000000000 20.0000000000000 10.0000000000000 5.00000000000000 2.50000000000000 1.25000000000000 1.24999988079071 1.25000011920929 0.624999940395355 0.625000059604645 --- 4.000000 4.000000 5.000000 5.000000 16.66667 24.99999 4.00000000000000 4.00000000000000 5.00000047683716 5.00000047683716 16.6666660308838 24.9999866485596 --- 0.0000000E+00 0.0000000E+00 0.000000000000000E+000 0.000000000000000E+000 ================================================================== ================================================================== aerosol_SGLI_3ch.3param.QA.201807.f:L1404 fn_land() ------------------------------------------------------------------ print *, LUT_TC_CAI( 1, 1 , 1, IX38), IX38 print *, dble(LUT_TC_CAI( 1, 1 , 1, IX38)), IX38 si380 = LIP_3param ( $ LUT_TC_CAI( 1, 1 , 1, IX38), $ Tau_A, CMF, Ref_R, dTAi, dCMFi, dRefi_R, $ E_Tau_a, E_CMF, E_Ref_R, ix,iy,iz, n_Tau_A, n_CMF, n_M) ------------------------------------------------------------------ PGI 0.1912796 1 0.1912796497344971 1 ----------------------------------------------------------------- INTEL 0.1912797 1 0.191279664635658 1 ================================================================== ================================================================== lmdif_s.f:L365 LMDIF() ------------------------------------------------------------------ CALL FCN(M,N,X,FVEC,IFLAG) print *, FVEC print *, dble(FVEC) ------------------------------------------------------------------ PGI -8.6558852E-03 4.5458078E-03 -0.1004647 -8.6558852344751358E-003 4.5458078384399414E-003 -0.1004647314548492 ------------------------------------------------------------------ INTEL -8.6558852E-03 4.5458078E-03 -0.1004647 -8.655885234475136E-003 4.545807838439941E-003 -0.100464731454849 ==================================================================