diff --git a/cfr/proxy.py b/cfr/proxy.py
index 670489b..32aa79e 100644
--- a/cfr/proxy.py
+++ b/cfr/proxy.py
@@ -240,7 +240,6 @@ def to_nc(self, path, verbose=True, **kwargs):
             path (str): the path to save the file.
             verbose (bool, optional): print verbose information. Defaults to False.
         '''
-        p_warning('>>> Warning: this is an experimental feature.')
         da = self.to_da()
 
         try:
@@ -1801,7 +1800,7 @@ def to_df(self):
             
         return df
 
-    def to_ds(self, annualize=False, months=None, verbose=False):
+    def to_ds(self, annualize=False, months=None, verbose=True):
         ''' Convert the proxy database to a `xarray.Dataset`
 
         Args:
@@ -1809,6 +1808,7 @@ def to_ds(self, annualize=False, months=None, verbose=False):
             months (list): months for annulization
             verbose (bool, optional): print verbose information. Defaults to False.
         '''
+        p_warning('>>> Warning: this is an experimental feature.')
         da_dict = {}
         pid_truncated = []
         for pobj in tqdm(self, total=self.nrec):
@@ -1865,6 +1865,7 @@ def to_nc(self, path, annualize=False, months=None, verbose=True, compress_param
             compress_params (dict): the paramters for compression when storing the reconstruction results to netCDF files.
             verbose (bool, optional): print verbose information. Defaults to False.
         '''
+        p_warning('>>> Warning: this is an experimental feature.')
         encoding_dict = {}
         for k in self.records.keys():
             encoding_dict[k] = compress_params