Commit e8721cb1 authored by Pierre Cazenave's avatar Pierre Cazenave

Merge branch 'master' of gitlab.ecosystem-modelling.pml.ac.uk:pica/fvcom-examples

parents 5759863c 8eac4bf1
......@@ -162,6 +162,7 @@ if strcmpi(inputConf.surface_heat, 'NCEP')
% - Air temperature (air) [celsius]
% - Relative humidity (rhum) [%]
% - Sea level pressure (pres) [Pa]
% - Wind (u and v) (uwnd, vwnd) [ms^{-1}]
%
% The script converts the NCEP data from the OPeNDAP server from
% longitudes in the 0 to 360 range to the -180 to 180 range. It
......@@ -169,7 +170,7 @@ if strcmpi(inputConf.surface_heat, 'NCEP')
% Mobj.lat).
heating = get_NCEP_forcing(Mobj, ...
[inputConf.startDateMJD, inputConf.endDateMJD], ...
'varlist', {'dlwrf', 'dswrf', 'air', 'rhum', 'pres'});
'varlist', {'dlwrf', 'dswrf', 'air', 'rhum', 'pres', 'uwnd', 'vwnd'});
heating.domain_cols = length(heating.lon);
heating.domain_rows = length(heating.lat);
......@@ -286,7 +287,7 @@ if strcmpi(inputConf.surface_heat, 'NCEP')
end
% Interpolate rhum, pres and longwave radiation to hourly values.
todo = {'rhum', 'pres', 'dlwrf', 'air'};
todo = {'rhum', 'pres', 'dlwrf', 'air', 'uwnd', 'vwnd'};
for ff = 1:length(todo)
% Force the data onto the 'global' position arrays. This fixes
% problems with grid2fvcom and the relative humidity grid in
......@@ -334,6 +335,7 @@ if strcmpi(inputConf.surface_heat, 'NCEP')
% [1, 1, size(heating.nshf.data, 3)])) = nan;
interpfields = {'dswrf', 'dlwrf', 'pres', 'air', 'rhum', ...
'uwnd', 'vwnd', ...
'time', 'lon', 'lat', 'x', 'y'};
heating_interp = grid2fvcom(Mobj, interpfields, heating);
if ftbverbose
......
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