Changeset 14587
- Timestamp:
- 04/15/13 15:21:20 (12 years ago)
- Location:
- issm/trunk-jpl/src/m/classes
- Files:
-
- 2 edited
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 8 9 10 11 12 13 14 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 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 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 66 67 68 69 70 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.