Changeset 14587
- Timestamp:
- 04/15/13 15:21:20 (13 years ago)
- Location:
- issm/trunk-jpl/src/m/classes
- Files:
-
- 2 edited
-
clusters/greenplanet.m (modified) (1 diff)
-
hydrologydc.m (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/clusters/greenplanet.m
r14501 r14587 62 62 if isempty(cluster.executionpath), md = checkmessage(md,'executionpath empty'); end 63 63 64 end 65 %}}} 66 function BuildKrigingQueueScript(cluster,modelname,solution,io_gather,isvalgrind,isgprof) % {{{ 67 68 if(isvalgrind), disp('valgrind not supported by cluster, ignoring...'); end 69 if(isgprof), disp('gprof not supported by cluster, ignoring...'); end 70 71 %compute number of processors 72 cluster.np=cluster.numnodes*cluster.cpuspernode; 73 74 %write queuing script 75 fid=fopen([modelname '.queue'],'w'); 76 fprintf(fid,'#PBS -S /bin/bash\n'); 77 fprintf(fid,'#PBS -N %s\n',modelname); 78 fprintf(fid,'#PBS -q %s \n',cluster.queue); 79 fprintf(fid,'#PBS -l nodes=%i:ppn=%i\n',cluster.numnodes,cluster.cpuspernode); 80 fprintf(fid,'#PBS -m bea\n'); 81 fprintf(fid,'#PBS -M mmorligh@uci.edu\n'); 82 fprintf(fid,'#PBS -o %s.outlog \n',modelname); 83 fprintf(fid,'#PBS -e %s.errlog \n\n',modelname); 84 %fprintf(fid,'cd %s/%s\n\n',cluster.executionpath,dirname); 85 fprintf(fid,'mpiexec -np %i %s/kriging.exe %s %s\n',cluster.np,cluster.codepath,[cluster.executionpath '/' modelname],modelname); 86 if ~io_gather, %concatenate the output files: 87 fprintf(fid,'cat %s.outbin.* > %s.outbin',modelname,modelname); 88 end 89 fclose(fid); 64 90 end 65 91 %}}} -
issm/trunk-jpl/src/m/classes/hydrologydc.m
r14579 r14587 5 5 6 6 classdef hydrologydc 7 properties (SetAccess=public)8 spcsediment_head = NaN;9 sediment_compressibility = 0;10 sediment_porosity = 0;11 sediment_thickness = 0;12 sediment_transmitivity = 0;13 water_compressibility = 0;14 water_density = 0;7 properties (SetAccess=public) 8 spcsediment_head = NaN; 9 sediment_compressibility = 0; 10 sediment_porosity = 0; 11 sediment_thickness = 0; 12 sediment_transmitivity = 0; 13 water_compressibility = 0; 14 water_density = 0; 15 15 end 16 16 methods 17 function obj = hydrologydc(varargin) % {{{18 switch nargin19 case 020 obj=setdefaultparameters(obj);21 otherwise22 error('constructor not supported');23 end24 end % }}}25 function obj = setdefaultparameters(obj) % {{{26 27 %Parameters from de Fleurian 201328 obj.sediment_compressibility = 1.0e-08;29 obj.sediment_porosity = .4;30 obj.sediment_thickness = 20.0;31 obj.sediment_transmitivity = 8.0e-04;32 obj.water_compressibility = 5.04e-10;33 obj.water_density = 1.0e+03;34 35 end % }}}36 function md = checkconsistency(obj,md,solution,analyses) % {{{37 38 %Early return39 if ~ismember(HydrologyAnalysisEnum(),analyses), return; end40 41 md = checkfield(md,'hydrology.spcsediment_head','forcing',1);42 md = checkfield(md,'hydrology.sediment_compressibility','>',0,'numel',1);43 md = checkfield(md,'hydrology.sediment_porosity','>',0,'numel',1);44 md = checkfield(md,'hydrology.sediment_thickness','>',0,'numel',1);45 md = checkfield(md,'hydrology.sediment_transmitivity','>',0,'numel',1);46 md = checkfield(md,'hydrology.water_compressibility','>',0,'numel',1);47 md = checkfield(md,'hydrology.water_density','>',0,'numel',1);48 end % }}}49 function disp(obj) % {{{50 disp(sprintf(' hydrology Dual Porous Continuum Equivalent parameters:'));51 52 fielddisplay(obj,'spcsediment_head','sediment water head constraints (NaN means no constraint, meters above MSL)');53 fielddisplay(obj,'sediment_compressibility','sediment compressibility (Pa -1)');54 fielddisplay(obj,'sediment_porosity','sediment(dimensionless)');55 fielddisplay(obj,'sediment_thickness','sediment thickness (meters)');56 fielddisplay(obj,'sediment_transmitivity','sediment transmitivity (m2 s-1)');57 fielddisplay(obj,'water_compressibility','water compressibility (Pa -1)');58 fielddisplay(obj,'water_density','water density (kg m 3)');59 60 end % }}}61 function marshall(obj,fid) % {{{62 WriteData(fid,'enum',HydrologyEnum(),'data',HydrologydcEnum(),'format','Integer');63 WriteData(fid,'object',obj,'fieldname','spcsediment_head','format','DoubleMat','mattype',1);64 WriteData(fid,'object',obj,'fieldname','sediment_compressibility','format','Double');65 WriteData(fid,'object',obj,'fieldname','sediment_porosity','format','Double');66 WriteData(fid,'object',obj,'fieldname','sediment_thickness','format','Double');67 WriteData(fid,'object',obj,'fieldname','sediment_transmitivity','format','Double');68 WriteData(fid,'object',obj,'fieldname','water_compressibility','format','Double');69 WriteData(fid,'object',obj,'fieldname','water_density','format','Double');70 end % }}}17 function obj = hydrologydc(varargin) % {{{ 18 switch nargin 19 case 0 20 obj=setdefaultparameters(obj); 21 otherwise 22 error('constructor not supported'); 23 end 24 end % }}} 25 function obj = setdefaultparameters(obj) % {{{ 26 27 %Parameters from de Fleurian 2013 28 obj.sediment_compressibility = 1.0e-08; 29 obj.sediment_porosity = .4; 30 obj.sediment_thickness = 20.0; 31 obj.sediment_transmitivity = 8.0e-04; 32 obj.water_compressibility = 5.04e-10; 33 obj.water_density = 1.0e+03; 34 35 end % }}} 36 function md = checkconsistency(obj,md,solution,analyses) % {{{ 37 38 %Early return 39 if ~ismember(HydrologyAnalysisEnum(),analyses), return; end 40 41 md = checkfield(md,'hydrology.spcsediment_head','forcing',1); 42 md = checkfield(md,'hydrology.sediment_compressibility','>',0,'numel',1); 43 md = checkfield(md,'hydrology.sediment_porosity','>',0,'numel',1); 44 md = checkfield(md,'hydrology.sediment_thickness','>',0,'numel',1); 45 md = checkfield(md,'hydrology.sediment_transmitivity','>',0,'numel',1); 46 md = checkfield(md,'hydrology.water_compressibility','>',0,'numel',1); 47 md = checkfield(md,'hydrology.water_density','>',0,'numel',1); 48 end % }}} 49 function disp(obj) % {{{ 50 disp(sprintf(' hydrology Dual Porous Continuum Equivalent parameters:')); 51 52 fielddisplay(obj,'spcsediment_head','sediment water head constraints (NaN means no constraint, meters above MSL)'); 53 fielddisplay(obj,'sediment_compressibility','sediment compressibility (Pa -1)'); 54 fielddisplay(obj,'sediment_porosity','sediment(dimensionless)'); 55 fielddisplay(obj,'sediment_thickness','sediment thickness (meters)'); 56 fielddisplay(obj,'sediment_transmitivity','sediment transmitivity (m2 s-1)'); 57 fielddisplay(obj,'water_compressibility','water compressibility (Pa -1)'); 58 fielddisplay(obj,'water_density','water density (kg m 3)'); 59 60 end % }}} 61 function marshall(obj,fid) % {{{ 62 WriteData(fid,'enum',HydrologyEnum(),'data',HydrologydcEnum(),'format','Integer'); 63 WriteData(fid,'object',obj,'fieldname','spcsediment_head','format','DoubleMat','mattype',1); 64 WriteData(fid,'object',obj,'fieldname','sediment_compressibility','format','Double'); 65 WriteData(fid,'object',obj,'fieldname','sediment_porosity','format','Double'); 66 WriteData(fid,'object',obj,'fieldname','sediment_thickness','format','Double'); 67 WriteData(fid,'object',obj,'fieldname','sediment_transmitivity','format','Double'); 68 WriteData(fid,'object',obj,'fieldname','water_compressibility','format','Double'); 69 WriteData(fid,'object',obj,'fieldname','water_density','format','Double'); 70 end % }}} 71 71 end 72 72 end
Note:
See TracChangeset
for help on using the changeset viewer.
![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)