write_FVCOM_bath

 Write bathymetry to FVCOM format bathymetry file

 function write_FVCOM_bath(Mobj,filename)

 DESCRIPTION:
    Generate an ascii FVCOM 3.x format bathymetry from Mesh object

 INPUT 
   Mobj     = Mesh object
   filename = FVCOM bathymetry file name

 OUTPUT:
    FVCOM bathymetry file: filename

 EXAMPLE USAGE
    write_FVCOM_bath(Mobj,'tst_cor.dat')   

 Author(s):  
    Geoff Cowles (University of Massachusetts Dartmouth)

 Revision history
   
==============================================================================

0001 function write_FVCOM_bath(Mobj,filename) 
0002 
0003 % Write bathymetry to FVCOM format bathymetry file
0004 %
0005 % function write_FVCOM_bath(Mobj,filename)
0006 %
0007 % DESCRIPTION:
0008 %    Generate an ascii FVCOM 3.x format bathymetry from Mesh object
0009 %
0010 % INPUT
0011 %   Mobj     = Mesh object
0012 %   filename = FVCOM bathymetry file name
0013 %
0014 % OUTPUT:
0015 %    FVCOM bathymetry file: filename
0016 %
0017 % EXAMPLE USAGE
0018 %    write_FVCOM_bath(Mobj,'tst_cor.dat')
0019 %
0020 % Author(s):
0021 %    Geoff Cowles (University of Massachusetts Dartmouth)
0022 %
0023 % Revision history
0024 %
0025 %==============================================================================
0026 subname = 'write_FVCOM_bath';
0027 global ftbverbose
0028 if(ftbverbose)
0029   fprintf('\n'); fprintf(['begin : ' subname '\n']);
0030 end;
0031 
0032 %------------------------------------------------------------------------------
0033 % Parse input arguments
0034 %------------------------------------------------------------------------------
0035 if(exist('Mobj')*exist('filename')==0)
0036     error('arguments to write_FVCOM_cor are incorrect')
0037 end;
0038 
0039 %------------------------------------------------------------------------------
0040 % Dump the file
0041 %------------------------------------------------------------------------------
0042 if(lower(Mobj.nativeCoords(1:3)) == 'car')
0043     x = Mobj.x;
0044     y = Mobj.y;
0045 else
0046     x = Mobj.lon;
0047     y = Mobj.lat;
0048 end;
0049 if(Mobj.have_bath)
0050     if(ftbverbose); fprintf('writing FVCOM bathymetry file %s\n',filename); end;
0051     fid = fopen(filename,'w');
0052     fprintf(fid,'Node Number = %d\n',Mobj.nVerts);
0053     for i=1:Mobj.nVerts
0054       fprintf(fid,'%f %f %f\n',x(i),y(i),Mobj.h(i));
0055     end;
0056     fclose(fid);
0057 else
0058     error('can''t write bathymetry to file, mesh object has no bathymetry')
0059 end;
0060 
0061 if(ftbverbose)
0062   fprintf(['end   : ' subname '\n'])
0063 end;
0064 
0065