@@ -1601,7 +1601,7 @@ EXTERN_MSC int GMT_grdseamount (void *V_API, int mode, void *args) {
1601
1601
else /* Circular features */
1602
1602
major = minor = S [smt ].radius ; /* Radius in m */
1603
1603
/* Compute area, volume, mean amplitude */
1604
- shape_func [S [smt ].build_mode ] (major , minor , S [smt ].height , Ctrl -> L .value , S [smt ].f , & area , & volume , & height );
1604
+ shape_func [S [smt ].build_mode ](major , minor , S [smt ].height , Ctrl -> L .value , S [smt ].f , & area , & volume , & height );
1605
1605
V [smt ] = volume ;
1606
1606
h [smt ] = amplitude ;
1607
1607
if (map ) { /* Report values in km^2, km^3, and m */
@@ -1613,13 +1613,14 @@ EXTERN_MSC int GMT_grdseamount (void *V_API, int mode, void *args) {
1613
1613
out [GMT_X ] = area ;
1614
1614
out [GMT_Y ] = volume ;
1615
1615
out [GMT_Z ] = height ;
1616
- GMT_Put_Record (API , GMT_WRITE_DATA , Out ); /* Write this to output */
1616
+ API -> GMT -> current .io .col_type [1 ][0 ] = API -> GMT -> current .io .col_type [1 ][1 ] = GMT_IS_FLOAT ; /* If set to geog Put_Record would error. */
1617
+ GMT_Put_Record (API , GMT_WRITE_DATA , Out ); /* Write this to output */
1617
1618
}
1618
1619
GMT_Report (API , GMT_MSG_INFORMATION , "Seamount %" PRIu64 " [%c] area, volume, mean height: %g %g %g\n" , n_smts , S [smt ].code , area , volume , height );
1619
1620
}
1620
1621
if (Ctrl -> L .active ) { /* OK, that was all we wanted */
1621
- gmt_M_free (GMT , Out );
1622
- if (GMT_End_IO (API , GMT_OUT , 0 ) != GMT_NOERROR ) /* Disables further data output */
1622
+ gmt_M_free (GMT , Out );
1623
+ if (GMT_End_IO (API , GMT_OUT , 0 ) != GMT_NOERROR ) /* Disables further data output */
1623
1624
goto wrap_up ;
1624
1625
goto wrap_up ;
1625
1626
}
0 commit comments