write_FVCOM_grid

 Write grid and connectivity to FVCOM format grid file

 function write_FVCOM_grid(Mobj,filename)

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

 INPUT 
   Mobj     = Mesh object
   filename = FVCOM grid file name

 OUTPUT:
    FVCOM grid file: filename

 EXAMPLE USAGE
    write_FVCOM_grid(Mobj,'tst_grd.dat')   

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

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

0001 function write_FVCOM_grid(Mobj,filename) 
0002 
0003 % Write grid and connectivity to FVCOM format grid file
0004 %
0005 % function write_FVCOM_grid(Mobj,filename)
0006 %
0007 % DESCRIPTION:
0008 %    Generate an ascii FVCOM 3.x format gridfile from Mesh object
0009 %
0010 % INPUT
0011 %   Mobj     = Mesh object
0012 %   filename = FVCOM grid file name
0013 %
0014 % OUTPUT:
0015 %    FVCOM grid file: filename
0016 %
0017 % EXAMPLE USAGE
0018 %    write_FVCOM_grid(Mobj,'tst_grd.dat')
0019 %
0020 % Author(s):
0021 %    Geoff Cowles (University of Massachusetts Dartmouth)
0022 %
0023 % Revision history
0024 %
0025 %==============================================================================
0026 subname = 'write_FVCOM_grid';
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_grid 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(ftbverbose);  fprintf('writing FVCOM gridfile %s\n',filename); end;
0050 fid = fopen(filename,'w');
0051 fprintf(fid,'Node Number = %d\n',Mobj.nVerts);
0052 fprintf(fid,'Cell Number = %d\n',Mobj.nElems);
0053 for i=1:Mobj.nElems
0054   fprintf(fid,'%d %d %d %d %d\n',i,Mobj.tri(i,1:3),i);
0055 end;
0056 for i=1:Mobj.nVerts
0057   fprintf(fid,'%d %f %f %f\n',i,x(i),y(i),Mobj.h(i));
0058 end;
0059 fclose(fid);
0060 
0061 if(ftbverbose);
0062   fprintf(['end   : ' subname '\n'])
0063 end;
0064