Add h and land_mask to ncml file

app/write_ncml/file/write_ncml_file.py

@@ -7,6 +7,9 @@ fvcom_file_loc = sys.argv[1]
 output_file = sys.argv[2]
 write_vars = sys.argv[3].split(',')
 
+write_vars.append('h')
+write_vars.append('mask')
+
 all_attributes = {'lon':{'var_name':'longitude','standard_name':'longitude', 'units':'degrees_east', 'long_name':'longitude', 'shape':'longitude'},
             'lat':{'var_name':'latitude','standard_name':'latitude', 'units':'degrees_north', 'long_name':'latitude', 'shape':'latitude'},
             'depth':{'standard_name':'depth', 'units':'m', 'long_name':'depth, measured downwards from free surface', 'axis':'Z', 'positive':'down', 'shape':'depth'},
@@ -16,7 +19,10 @@ all_attributes = {'lon':{'var_name':'longitude','standard_name':'longitude', 'un
             'salinity':{'var_name':'so','standard_name':'sea_water_salinity','units':'psu', 'missing_value':-32768, 'long_name':'Sea Water Salinity', 'shape':'time depth latitude longitude'},
             'u':{'var_name':'uo', 'standard_name':'eastward_sea_water_velocity','units':'m s-1', 'missing_value':-32768, 'long_name':'Eastward Current Velocity', 'shape':'time depth latitude longitude'},
             'v':{'var_name':'vo','standard_name':'northward_sea_water_velocity','units':'m s-1', 'missing_value':-32768, 'long_name':'Northward Current Velocity', 'shape':'time depth latitude longitude'},
-            'zeta':{'var_name':'zos','standard_name':'sea_surface_height_above_geoid','units':'m', 'missing_value':-32768,'long_name':'Sea surface height above geoid', 'shape':'time latitude longitude'}}
+            'zeta':{'var_name':'zos','standard_name':'sea_surface_height_above_geoid','units':'m', 'missing_value':-32768,'long_name':'Sea surface height above geoid', 'shape':'time latitude longitude'},
+            'mask':{'standard_name':'sea_binary_mask','units':'1', 'long_name':'Land-sea mask: sea = 1 ; land = 0', 'shape':'depth latitude longitude'},
+            'h':{'standard_name':'h', 'units':'m', 'long_name':'model bathymetry depth, measured downwards from geoid', 'axis':'Z', 'positive':'down', 'shape':'depth latitude longitude' }
+            }
 
 all_type = {'add_offset':'float', 'scale_factor':'float', '_FillValue':'double'}