Changeset 28065

Timestamp:

01/24/24 10:39:30 (14 months ago)

Author:

bgetraer

Message:

CHG: interpISMIP6 function inputs define model with single string, simplify nc file identification

Location:

issm/trunk-jpl/src/m/parameterization

Files:

• interpISMIP6AntarcticaOcn.m (modified) (2 diffs)
• interpISMIP6AntarcticaSMB.m (modified) (2 diffs)

issm/trunk-jpl/src/m/parameterization/interpISMIP6AntarcticaOcn.m

@@ -1 @@
-function basalforcings = interpISMIP6AntarcticaOcn(md,model_name,scenario)
+function basalforcings = interpISMIP6AntarcticaOcn(md,model_name)
 %interpISMIP6AntarcticaOcn - interpolate chosen ISMIP6 atmospheric forcing to model
 %
 %   Input:
 %     - md (model object)
-%     - model_name  (string): name of the climate model 
-%                             Examples: cnrm-esm2-1 ccsm4cesm2 cnrm-cm6-1 csiro-mk3-6-0
-%                                       hadgem2-es ipsl-cm5a-mr miroc-esm-chem noresm1-m ukesm1-0-ll
+%     - model_name (string): name of the climate model and scenario
+%       - suppported options from /totten_1/ModelData/ISMIP6/Projections/AIS/Ocean_Forcing/
+%             2.6 scenario             8.5 scenario
+%             ---------------------------------------------
+%                                      ccsm4_rcp8.5
+%                                      cesm2_ssp585
+%             cnrm-cm6-1_ssp126        cnrm-cm6-1_ssp585
+%                                      cnrm-esm2-1_ssp585
+%                                      csiro-mk3-6-0_rcp8.5
+%                                      hadgem2-es_rcp8.5
+%             ipsl-cm5a-mr_rcp2.6      ipsl-cm5a-mr_rcp8.5
+%                                      miroc-esm-chem_rcp8.5
+%             noresm1-m_rcp2.6         noresm1-m_rcp8.5
+%                                      ukesm1-0-ll_ssp585
 %
-%     - scenario    (string): name of the climate scenario
-%                             Examples: rcp2.6, rcp8.5, ssp126, ssp585
 %   Output:
-%     - basalforcings: prepared to be input directly into md.smb
+%     - basalforcings: prepared to be input directly into md.basalforcings
+%                      time series from 1995-2100
 %
 %   Examples:
-%      md.basalforcings = interpISMIP6AntarcticaOcn(md,'miroc-esm-chem','rcp8.5');
+%      md.basalforcings = interpISMIP6AntarcticaOcn(md,'miroc-esm-chem_rcp8.5');
 
 % Find appropriate directory
@@ -24 +34 @@
 end
 
-%Find forcing file
-rootname = [path model_name '_' scenario '/1995-2100'];
-list = dir(rootname);
-found = false;
-for i=1:numel(list)
-        if strcmpi(list(i).name, [model_name '_thermal_forcing_8km_x_60m.nc'])
-                disp(['Found tf file: ' list(i).name]);
-                tfnc = [rootname '/' list(i).name];
-                found = true;
-        end
+% search for thermal forcing file in the ISMIP climate model directory
+rootname=[path model_name '/1995-2100/']; % root directory for the climate model files
+tffile=dir([rootname '*_thermal_forcing_8km_x_60m.nc']); % thermal forcing file if found
+
+% throw error if file not found, or if the file search is not unique
+if length(tffile)~=1
+   error(['this path does not exist or is not unique under ' rootname]);
 end
-if ~found
-   error(['this path does not exist or the ' model_name ' and ' scenario ' are not available in this combination.']);
-end
+
+% save the full path of the found files
+tfnc=[rootname tffile.name];
 
 %load TF data

issm/trunk-jpl/src/m/parameterization/interpISMIP6AntarcticaSMB.m

@@ -1 @@
-function smb = interpISMIP6AntarcticaSMB(md,model_name,scenario)
+function smb = interpISMIP6AntarcticaSMB(md,model_name)
 %interpISMIP6AntarcticaSMB - interpolate chosen ISMIP6 atmospheric forcing to model
 %
 %   Input:
 %     - md (model object)
-%     - model_name  (string): name of the climate model 
-%                             Examples: CESM2 CNRM_CM6 CNRM_ESM2 CSIRO-Mk3-6-0 HadGEM2-ES IPSL-CM5A-MR 
-%                                       ccsm4 miroc-esm-chem noresm1-m
-%     - scenario    (string): name of the climate scenario
-%                             Examples: rcp26, rcp85, ssp126, ssp585
+%     - model_name (string): name of the climate model and scenario
+%       - suppported options from /totten_1/ModelData/ISMIP6/Projections/AIS/Atmosphere_Forcing/
+%             2.6 scenario             8.5 scenario
+%             ---------------------------------------------
+%             ccsm4_rcp2.6             ccsm4_rcp8.5
+%                                      CESM2_ssp585 
+%             CNRM_CM6_ssp126          CNRM_CM6_ssp585 
+%                                      CNRM_ESM2_ssp585
+%                                      CSIRO-Mk3-6-0_rcp85
+%                                      HadGEM2-ES_rcp85
+%             IPSL-CM5A-MR_rcp26       IPSL-CM5A-MR_rcp85
+%             miroc-esm-chem_rcp2.6    miroc-esm-chem_rcp8.5
+%             noresm1-m_rcp2.6         noresm1-m_rcp8.5
+%
 %   Output:
 %     - smb: prepared to be input directly into md.smb
+%            time series from 1995-2100
 %
 %   Examples:
-%      md.smb = interpISMIP6AntarcticaSMB(md,'miroc-esm-chem','rcp8.5');
+%      md.smb = interpISMIP6AntarcticaSMB(md,'miroc-esm-chem_rcp8.5');
 
 % Find appropriate directory
@@ -23 +33 @@
 end
 
-%find which files are available and do a case insensitive test
-rootname = [path model_name '_' scenario '/'];
-list = dir([rootname '/Regridded_2km/']);
-found1 = false; found2 = false;
-for i=1:numel(list)
-        if strcmpi(list(i).name, [model_name '_2km_clim_1995-2014.nc'])
-                disp(['Found clim file: ' list(i).name]);
-                smbclimnc = [rootname '/Regridded_2km/' list(i).name];
-                found1 = true;
-        end
-        if strcmpi(list(i).name, [model_name '_2km_anomaly_1995-2100.nc'])
-                disp(['Found anom file: ' list(i).name]);
-                smbanomnc = [rootname '/Regridded_2km/' list(i).name];
-                found2 = true;
-        end
-end
-if ~found1 || ~found2
-   error(['this path does not exist or the ' model_name ' and ' scenario ' are not available in this combination.']);
+% search for smb files in the ISMIP climate model directory
+rootname=[path model_name '/Regridded_2km/']; % root directory for the climate model files
+smbclimfile=dir([rootname '*_2km_clim_1995-2014.nc']); % climatology file if found
+smbanomfile=dir([rootname '*_2km_anomaly_1995-2100.nc']); % anomaly file if found
+
+% throw error if files are not found, or if the file search is not unique
+if length(smbclimfile)~=1 || length(smbanomfile)~=1
+   error(['this path does not exist or is not unique under ' rootname]);
 end
 
+% save the full path of the found files
+smbclimnc=[rootname smbclimfile.name];
+smbanomnc=[rootname smbanomfile.name];
+
+% load data from files
 disp('   == loading TS and SMB climatology data');
 lat                 = double(ncread(smbclimnc,'lat'));