comments etc

EBsep.f90

@@ -82,7 +82,6 @@ program EBsep
         print *, 'smooth fsky=',sum(Mask%TQU)/Mask%npix
     end if
 
-    !/scratch/maja1/planck/repository/exchanges/dx11/maps/hfi/HFI_SkyMap_353_2048_DX11c_full.fits
     if (.not. FileExists(matrix_name)) then
         !First get alm of mask
         !Calculate cut-sky coupling matrix
@@ -105,7 +104,6 @@ program EBsep
                 call Matrix_Mult_SymmLeft(Couplings2%WPAsymm,Proj%M, tmpM)
                 deallocate(Couplings2%WPAsymm)
                 allocate(RotW(Proj%nr,Proj%nr))
-                print *,'mult'
                 call Matrix_Mult_TN(tmpM,tmpM, rotW)
                 deallocate(tmpM)
                 print *,'norm'
@@ -176,6 +174,7 @@ program EBsep
 
     if (in_map /= '') then
         call HealpixMap_Read(LoadMap, in_map)
+        if (LoadMap%nmaps < 3) error stop 'Map has no polarization'
         if (in_map_diff/='') then
             call HealpixMap_Read(InMap, in_map_diff)
             LoadMap%TQU = (LoadMap%TQU - InMap%TQU)/2

EBsep.ini

@@ -1,7 +1,8 @@
 num_threads = 16
-#nside for the input mask and map
+#nside for internal working and output map (will de/upgrade inputs as neccessary)
 nside = 2048
 #max to use for separation, will get all E/B modes at L << lmax
+#lmax=140 takes about 1 hour on a 16core node to calculate coupling matrix (applying it much faster)
 lmax=140
 #Tolerance for separation
 support=0.999

Makefile_EBsep

@@ -20,9 +20,6 @@ OBJFILES= toms760.o inifile.o utils.o spin_alm_tools.o \
 
 OBJFILESEB = $(OBJFILES) Matrix_utils.o AsymmCouplings.o EBsep.o
 
-OBJFILESS = $(OBJFILES)  SimLens.o
-
-OBJFILESW = $(OBJFILES)  Matrix_utils.o PseudoCl.o AsymmCouplings.o SNModes.o TOMS642.o WeightMixer.o
 
 default: EBsep
 all: EBsep
@@ -49,12 +46,6 @@ EBsep.o: HealpixVis.o PseudoCl.o inifile.o AsymmCouplings.o
 	$(F90C) $(F90FLAGS) -c $*.F90
 
 
-simlens: $(OBJFILESS) 	
-	$(F90C) -o simlens $(OBJFILESS) $(F90FLAGS) $(LINKFLAGS)
-
-weightmixer: $(OBJFILESW) 	
-	$(F90C) -o weightmixer $(OBJFILESW) $(F90FLAGS) $(LINKFLAGS)
-
 EBsep: $(OBJFILESEB) 	
 	$(F90C) -o EBsep $(OBJFILESEB) $(F90FLAGS) $(LINKFLAGS)