debug export_to_fsmoshd

Author: arcticsnow

TopoPyScale/topo_export.py

@@ -908,8 +908,8 @@ def write_fsm2oshd_met(ds_pt,
 
         # sampling forest parameters from ds_param
         df_forest[['isfor', 'forcov', 'svf_for', 'CC5', 'CC50', 'CH5', 'LAI5']] = mp.toposub.ds_param.sel(x=xs, y=ys, method='nearest').to_dataframe()[['isfor', 'forcov', 'svf_for', 'CC5', 'CC50', 'CH5', 'LAI5']]
-        pixel_area = np.abs(mp.toposub.ds_param.x.diff(dim='x').mean().values * mp.toposub.ds_param.y.diff(dim='y').mean().values)
-        df_forest['cluster_domain_size'] = pixel_area
+        pixel_size = np.abs(mp.toposub.ds_param.x.diff(dim='x').mean().values) 
+        df_forest['cluster_domain_size'] = pixel_size
 
     else:
         raise ValueError("Sampling method not yet supported. Avail: clusters, points")
@@ -918,7 +918,6 @@ def write_fsm2oshd_met(ds_pt,
     # rename variable columns to match namelist functino varnames
     new_name = {'LAI5':'lai5', 'svf_for':'vfhp', 'CC5':'fveg', 'CC50':'fves', 'CH5':'hcan'}
     df_forest = df_forest.rename(columns=new_name)
-    print(df_forest)
 
     # ----- Loop through all points-------
     # NOTE: eventually this for loop could be parallelized to several cores -----