Merge pull request #244 from ZPYin/dev

Bugfixes of wrong assiged aerosol backscatter in in-complete overlap zone

.gitignore

@@ -13,6 +13,7 @@ todo_filelist
 test
 netcdf4
 log
+analysis
 done_filelist
 
 # Sphinx build

lib/interface/picassoProcV3.m

@@ -3327,7 +3327,7 @@
 
             % optical thickness (OT)
             aExt355 = aerExt355_klett(iGrp, :);
-            aExt355(1:hIndBase) = aerExt355_klett(hIndBase);
+            aExt355(1:hIndBase) = aerExt355_klett(iGrp, hIndBase);
             aBsc355 = aerBsc355_klett(iGrp, :);
             aOT355 = nancumsum(aExt355 .* [data.distance0(1), diff(data.distance0)]);
             mOT355 = nancumsum(mExt355 .* [data.distance0(1), diff(data.distance0)]);
@@ -3366,7 +3366,7 @@
 
             % optical thickness (OT)
             aExt532 = aerExt532_klett(iGrp, :);
-            aExt532(1:hIndBase) = aerExt532_klett(hIndBase);
+            aExt532(1:hIndBase) = aerExt532_klett(iGrp, hIndBase);
             aBsc532 = aerBsc532_klett(iGrp, :);
             aOT532 = nancumsum(aExt532 .* [data.distance0(1), diff(data.distance0)]);
             mOT532 = nancumsum(mExt532 .* [data.distance0(1), diff(data.distance0)]);
@@ -3405,7 +3405,7 @@
 
             % optical thickness (OT)
             aExt1064 = aerExt1064_klett(iGrp, :);
-            aExt1064(1:hIndBase) = aerExt1064_klett(hIndBase);
+            aExt1064(1:hIndBase) = aerExt1064_klett(iGrp, hIndBase);
             aBsc1064 = aerBsc1064_klett(iGrp, :);
             aOT1064 = nancumsum(aExt1064 .* [data.distance0(1), diff(data.distance0)]);
             mOT1064 = nancumsum(mExt1064 .* [data.distance0(1), diff(data.distance0)]);
@@ -3632,7 +3632,7 @@
 
             % optical thickness (OT)
             aExt355 = aerExt355_aeronet(iGrp, :);
-            aExt355(1:hIndBase) = aerExt355_aeronet(hIndBase);
+            aExt355(1:hIndBase) = aerExt355_aeronet(iGrp, hIndBase);
             aBsc355 = aerBsc355_aeronet(iGrp, :);
             aOT355 = nancumsum(aExt355 .* [data.distance0(1), diff(data.distance0)]);
             mOT355 = nancumsum(mExt355 .* [data.distance0(1), diff(data.distance0)]);
@@ -3671,7 +3671,7 @@
 
             % optical thickness (OT)
             aExt532 = aerExt532_aeronet(iGrp, :);
-            aExt532(1:hIndBase) = aerExt532_aeronet(hIndBase);
+            aExt532(1:hIndBase) = aerExt532_aeronet(iGrp, hIndBase);
             aBsc532 = aerBsc532_aeronet(iGrp, :);
             aOT532 = nancumsum(aExt532 .* [data.distance0(1), diff(data.distance0)]);
             mOT532 = nancumsum(mExt532 .* [data.distance0(1), diff(data.distance0)]);
@@ -3710,7 +3710,7 @@
 
             % optical thickness (OT)
             aExt1064 = aerExt1064_aeronet(iGrp, :);
-            aExt1064(1:hIndBase) = aerExt1064_aeronet(hIndBase);
+            aExt1064(1:hIndBase) = aerExt1064_aeronet(iGrp, hIndBase);
             aBsc1064 = aerBsc1064_aeronet(iGrp, :);
             aOT1064 = nancumsum(aExt1064 .* [data.distance0(1), diff(data.distance0)]);
             mOT1064 = nancumsum(mExt1064 .* [data.distance0(1), diff(data.distance0)]);