@@ -795,7 +795,7 @@ def pcol_2dxy(self, x, y):
795795    #date_str, date_end = 19980107, 19991110 #  
796796    #date_str, date_end = 20081107, 20100110 #  
797797    #date_str, date_end = 19921014, 20120718 #  
798-     date_str , date_end  =  20020101 , 20021231  #  
798+     date_str , date_end  =  20020101 , 20020630  #  
799799
800800    # Choose type of diagnostic: either q-parameter ('Q') or sea level anomaly ('sla') 
801801    #diag_type = 'Q' #<<< not implemented in 1.2.0 
@@ -807,11 +807,12 @@ def pcol_2dxy(self, x, y):
807807    #savedir = '/marula/emason/aviso_eddy_tracking/pablo_exp/' 
808808    #savedir = '/marula/emason/aviso_eddy_tracking/new_AVISO_test/' 
809809    #savedir = '/marula/emason/aviso_eddy_tracking/new_AVISO_SUBSAMP-3days/' 
810-     savedir  =  '/marula/emason/aviso_eddy_tracking/junk/' 
810+     # savedir = '/marula/emason/aviso_eddy_tracking/junk/'
811811    #savedir = '/marula/emason/aviso_eddy_tracking/new_AVISO_test/BlackSea/' 
812812    #savedir = '/marula/emason/aviso_eddy_tracking/Corridor_V3_Dec2014/' 
813813    #savedir = '/marula/emason/aviso_eddy_tracking/CLS_test_1_acd66ba338a9/' 
814814    #savedir = '/marula/emason/aviso_eddy_tracking/CLS_test_2/' 
815+     savedir  =  '/marula/emason/aviso_eddy_tracking/CLS_test_3/' 
815816    #savedir = '/path/to/save/your/outputs/' 
816817
817818
@@ -908,10 +909,10 @@ def pcol_2dxy(self, x, y):
908909        #latmin = -47. 
909910        #latmax = -24. 
910911
911-         #lonmin = -32 .     # small test 
912+         #lonmin = -38 .     # small test 
912913        #lonmax = -25. 
913-         #latmin = 20 . 
914-         #latmax = 28 . 
914+         #latmin = 30 . 
915+         #latmax = 38 . 
915916
916917
917918    elif  the_domain  in  'MedSea' :
@@ -1005,8 +1006,8 @@ def pcol_2dxy(self, x, y):
10051006        hanning_passes  =  5 
10061007
10071008    elif  'sla'  in  diag_type :
1008-         # shape_err = 55. * np.ones(slaparameter.size)
1009-         shape_err  =  1000.  *  np .ones (slaparameter .size )
1009+         shape_err  =  55.  *  np .ones (slaparameter .size )
1010+         # shape_err = 1000. * np.ones(slaparameter.size)
10101011        #shape_err = np.power(np.linspace(85., 40,  slaparameter.size), 2) / 100. 
10111012        #shape_err[shape_err < 50.] = 50. 
10121013
0 commit comments