Commit 69a1b05c authored by Pierre Cazenave's avatar Pierre Cazenave

Slight improvement on the way the wind vectors are dealt with

parent 3aa7aff1
...@@ -124,7 +124,7 @@ for i=1:length(suffixes) ...@@ -124,7 +124,7 @@ for i=1:length(suffixes)
netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'title','FVCOM Forcing File') netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'title','FVCOM Forcing File')
netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'institution','Plymouth Marine Laboratory') netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'institution','Plymouth Marine Laboratory')
netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'source','FVCOM grid (unstructured) surface forcing') netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'source','FVCOM grid (unstructured) surface forcing')
netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'history',['File created on ', datestr(now, 'yyyy-mm-dd HH:MM:SS'), ' with write_FVCOM_forcing.m from the fvcom-toolbox (https://github.com/pwcazenave/fvcom-toolbox)']) netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'history',['File created on ', datestr(now, 'yyyy-mm-dd HH:MM:SS'), ' with write_FVCOM_forcing.m from the MATLAB fvcom-toolbox (https://github.com/pwcazenave/fvcom-toolbox)'])
netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'references','http://fvcom.smast.umassd.edu, http://codfish.smast.umassd.edu') netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'references','http://fvcom.smast.umassd.edu, http://codfish.smast.umassd.edu')
netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'Conventions','CF-1.0') netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'Conventions','CF-1.0')
% netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'infos',infos) % netcdf.putAtt(nc,netcdf.getConstant('NC_GLOBAL'),'infos',infos)
...@@ -262,6 +262,11 @@ for i=1:length(suffixes) ...@@ -262,6 +262,11 @@ for i=1:length(suffixes)
used_fnames = [used_fnames, {'uwnd', 'vwnd', 'uwnd', 'vwnd'}]; used_fnames = [used_fnames, {'uwnd', 'vwnd', 'uwnd', 'vwnd'}];
used_dims = [used_dims, {'nElems', 'nElems', 'nElems', 'nElems'}]; used_dims = [used_dims, {'nElems', 'nElems', 'nElems', 'nElems'}];
end end
case {'vwnd', 'v10'}
% We dealt with these in the u component section above, so
% just pass silently.
true;
case 'slp' case 'slp'
if strcmpi(suffixes{i}, '_air_press') || ~multi_out if strcmpi(suffixes{i}, '_air_press') || ~multi_out
...@@ -382,7 +387,7 @@ for i=1:length(suffixes) ...@@ -382,7 +387,7 @@ for i=1:length(suffixes)
otherwise otherwise
if(ftbverbose) if(ftbverbose)
warning('Unknown or unused input data type: %s', fnames{vv}) warning('Unknown or possibly unused input data type: %s', fnames{vv})
end end
end end
end end
...@@ -464,7 +469,9 @@ for i=1:length(suffixes) ...@@ -464,7 +469,9 @@ for i=1:length(suffixes)
fprintf('done.\n') fprintf('done.\n')
end end
end end
if hf_done ~= 4 if hf_done < 4
% hf_done might be higher than four, but unless it is at least
% four, we haven't got everything we need.
warning('Did not have all the required heat flux parameters. Need ''shtfl'', ''lhtfl'', ''nlwrs'' and ''nwsrs''.') warning('Did not have all the required heat flux parameters. Need ''shtfl'', ''lhtfl'', ''nlwrs'' and ''nwsrs''.')
end end
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment