Changeset 27532
- Timestamp:
- 01/19/23 07:36:24 (2 years ago)
- Location:
- issm/trunk-jpl/src/m
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/clusters/pace.m
r27531 r27532 9 9 properties (SetAccess=public) 10 10 % {{{ 11 name = 'login-phoenix- slurm.pace.gatech.edu' %Phoenix cluster name11 name = 'login-phoenix-4.pace.gatech.edu' %Phoenix cluster name 12 12 login = ''; %personal login 13 13 numnodes = 1; %number of nodes requested … … 17 17 queue = 'inferno'; %queue 18 18 time = 60; %time requested per run [minutes] 19 accountname = ' gts-arobel3-atlas'; %group account name19 accountname = 'GT-arobel3-atlas'; %group account name 20 20 codepath = ''; %path to issm binaries 21 21 executionpath = ''; %path for execution folder … … 62 62 fid=fopen([modelname '.queue'],'w'); 63 63 fprintf(fid,'#!/bin/sh\n'); 64 fprintf(fid,'#PBS -l walltime=%i\n',cluster.time*60); %conversion of walltime from minutes to seconds. 65 fprintf(fid,'#PBS -N %s\n',modelname); 66 fprintf(fid,'#PBS -l nodes=1:ppn=%i\n',cluster.np); 67 fprintf(fid,'#PBS -l pmem=%igb\n',cluster.mem); 68 fprintf(fid,'#PBS -q %s\n',cluster.queue); 69 fprintf(fid,'#PBS -A %s\n',cluster.accountname); 70 71 fprintf(fid,'#PBS -o %s/%s/%s.outlog \n',cluster.executionpath,dirname,modelname); 72 fprintf(fid,'#PBS -e %s/%s/%s.errlog \n\n',cluster.executionpath,dirname,modelname); 64 73 65 fprintf(fid,'#SBATCH -t%i\n',cluster.time); 66 fprintf(fid,'#SBATCH -J%s\n',modelname); 67 fprintf(fid,'#SBATCH -N 1 --ntasks-per-node=%i\n',cluster.np); 68 %fprintf(fid,'#SBATCH -N %i\n',cluster.numnodes); 69 %fprintf(fid,'#SBATCH --ntasks=1\n'); 70 %fprintf(fid,'#SBATCH --cpus-per-task=%i\n',cluster.np); 71 fprintf(fid,'#SBATCH --mem-per-cpu=%iG\n',cluster.mem); 72 fprintf(fid,'#SBATCH -p%s\n',cluster.queue); 73 fprintf(fid,'#SBATCH -A %s\n',cluster.accountname); 74 fprintf(fid,'#SBATCH -o%s/%s/%s.outlog \n',cluster.executionpath,dirname,modelname); 75 fprintf(fid,'#SBATCH -e%s/%s/%s.errlog \n\n',cluster.executionpath,dirname,modelname); 76 fprintf(fid,'export SLURM_SUBMIT_DIR=%s\n',[cluster.executionpath '/' dirname]); 77 fprintf(fid,'cd $SLURM_SUBMIT_DIR\n'); 78 fprintf(fid,'export LD_LIBRARY_PATH=/opt/slurm/current/lib:/opt/pmix/current/lib:$LD_LIBRARY_PATH \n'); 79 fprintf(fid,'srun --mpi=pmi2 -n %i %s/%s %s %s %s \n',cluster.np,cluster.codepath,executable,solution,[cluster.executionpath '/' dirname],modelname); 74 fprintf(fid,'export PBS_O_WORKDIR=%s\n',[cluster.executionpath '/' dirname]); 75 fprintf(fid,'cd $PBS_O_WORKDIR\n'); 76 fprintf(fid,'mpiexec -np %i %s/%s %s %s %s \n',cluster.np,cluster.codepath,executable,solution,[cluster.executionpath '/' dirname],modelname); 80 77 81 78 fclose(fid); … … 100 97 disp('launching solution sequence on remote cluster'); 101 98 if ~isempty(restart) 102 launchcommand=['cd ' cluster.executionpath ' && cd ' dirname ' && sbatch' modelname '.queue '];99 launchcommand=['cd ' cluster.executionpath ' && cd ' dirname ' && qsub ' modelname '.queue ']; 103 100 else 104 101 launchcommand=['cd ' cluster.executionpath ' && rm -rf ./' dirname ' && mkdir ' dirname ... 105 ' && cd ' dirname ' && mv ../' dirname '.tar.gz ./ && tar -zxf ' dirname '.tar.gz && sbatch' modelname '.queue '];102 ' && cd ' dirname ' && mv ../' dirname '.tar.gz ./ && tar -zxf ' dirname '.tar.gz && qsub ' modelname '.queue ']; 106 103 end 107 104 issmssh(cluster.name,cluster.login,cluster.port,launchcommand); -
issm/trunk-jpl/src/m/classes/hydrologyarmapw.m
r27531 r27532 18 18 malag_coefs = NaN; 19 19 datebreaks = NaN; 20 requested_outputs = {};21 20 end 22 21 methods … … 185 184 WriteData(fid,prefix,'data',dbreaks,'name','md.hydrology.datebreaks','format','DoubleMat','scale',yts); 186 185 WriteData(fid,prefix,'data',tempmonthlyfactors,'name','md.hydrology.monthlyfactors','format','DoubleMat'); 187 188 %process requested outputs 189 outputs = self.requested_outputs; 190 pos = find(ismember(outputs,'default')); 191 if ~isempty(pos), 192 outputs(pos) = []; %remove 'default' from outputs 193 outputs = [outputs defaultoutputs(self,md)]; %add defaults 194 end 195 WriteData(fid,prefix,'data',outputs,'name','md.hydrology.requested_outputs','format','StringArray'); 186 WriteData(fid,prefix,'data',{'FrictionWaterPressure'},'name','md.hydrology.requested_outputs','format','StringArray'); 196 187 end % }}} 197 188 end -
issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpBedmachineGreenland.m
r27531 r27532 45 45 ['/Users/larour/ModelData/BedMachine/' basename '-' ncdate '.nc'],... 46 46 ['./' basename '-' ncdate '.nc'],... 47 '/media/vincent/TOSH4TB/GeorgiaTech/DataSearch/BedMachine/BedMachineGreenland-2021-04-20.nc',...48 47 }; 49 48 … … 90 89 if strcmp(string,'mask') | strcmp(string,'source'), 91 90 %Need nearest neighbor to avoid interpolation between 0 and 2 92 %output = InterpFromGrid(xdata,ydata,data,double(X),double(Y),'nearest'); 93 output = InterpFromGridToMesh(xdata,flipud(ydata),flipud(data),double(X),double(Y),NaN); %VV 91 output = InterpFromGrid(xdata,ydata,data,double(X),double(Y),'nearest'); 94 92 else 95 %output = InterpFromGrid(xdata,ydata,data,double(X),double(Y)); 96 output = InterpFromGridToMesh(xdata,flipud(ydata),flipud(data),double(X),double(Y),NaN); %VV 93 output = InterpFromGrid(xdata,ydata,data,double(X),double(Y)); 97 94 end 98 95 -
issm/trunk-jpl/src/m/contrib/morlighem/modeldata/interpSeaRISE.m
r27531 r27532 58 58 searisenc='/totten_1/ModelData/SeaRISE/Antarctica_5km_dev1.0.nc'; 59 59 end 60 case {'vincentxps159500'}61 if hemisphere==162 searisenc = '/media/vincent/TOSH4TB/GeorgiaTech/ISSM/fromSeroussi/GrIS04Feb2022/seariseData/Greenland_5km_dev1.2.nc';63 end64 60 otherwise 65 61 error('hostname not supported yet'); … … 84 80 if verbose, disp([' -- SeaRISE: interpolating ' string]); end 85 81 if strcmpi(string,'LandMask'); 86 %dataout = InterpFromGrid(xdata,ydata,data,xproj,yproj,'nearest'); 87 dataout = InterpFromGridToMesh(xdata,ydata,data,xproj,yproj,NaN); %VV 82 dataout = InterpFromGrid(xdata,ydata,data,xproj,yproj,'nearest'); 88 83 else 89 %dataout = InterpFromGrid(xdata,ydata,data,xproj,yproj); 90 dataout = InterpFromGridToMesh(xdata,ydata,data,xproj,yproj,NaN); %VV 84 dataout = InterpFromGrid(xdata,ydata,data,xproj,yproj); 91 85 end 92 86 dataout = reshape(dataout,size(X,1),size(X,2));
Note:
See TracChangeset
for help on using the changeset viewer.