source: issm/oecreview/Archive/24307-24683/ISSM-24567-24568.diff@ 24684

Last change on this file since 24684 was 24684, checked in by Mathieu Morlighem, 5 years ago

CHG: added new review
File size: 7.2 KB

  • ../trunk-jpl/test/NightlyRun/test247.m

     
     1%Test Name: SquareShelfTranIspddIsdeltaO18pdNoInterpSSA2d(
     2md=triangle(model(),'../Exp/Square.exp',60000.);
     3md=setmask(md,'all','');
     4md=parameterize(md,'../Par/SquareShelf.par');
     5
     6%md.verbose=verbose('all');
     7
     8% Use of ispdd and isdelta18o methods
     9md.smb = SMBd18opdd();
     10md.smb.isd18opd=1;
     11
     12% Add temperature, precipitation and delta18o needed to measure the surface mass balance
     13%  creating delta18o
     14load '../Data/delta18o.data'
     15md.smb.delta18o=delta18o;
     16
     17% creating Present day  temperatures
     18% Same temperature over the all region:
     19tmonth(1:12)=238.15+20.;
     20for imonth=0:11
     21        md.smb.temperatures_presentday(1:md.mesh.numberofvertices,imonth+1)=tmonth(imonth+1);
     22        % Time for the last line:
     23        md.smb.temperatures_presentday(md.mesh.numberofvertices+1,imonth+1)=((imonth+1)/12);
     24end
     25md.smb = initialize(md.smb,md);
     26
     27% creating initialization and spc temperatures initialization and
     28% spc
     29md.thermal.spctemperature=mean(md.smb.temperatures_presentday(1:md.mesh.numberofvertices,1:12),2)-10; %-10*ones(md.mesh.numberofvertices,1);
     30md.initialization.temperature=md.thermal.spctemperature; %md.smb.temperatures_presentday(1:md.mesh.numberofvertices,1);
     31
     32% creating precipitation
     33for imonth=0:11
     34        md.smb.precipitations_presentday(1:md.mesh.numberofvertices,imonth+1)=-0.4*10^(-6)*md.mesh.y+0.5;
     35        % Time for the last line:
     36        md.smb.precipitations_presentday(md.mesh.numberofvertices+1,imonth+1)=((imonth+1)/12);
     37end
     38
     39md.smb.s0p = (max(md.geometry.surface,0));
     40md.smb.s0t = (max(md.geometry.surface,0));
     41md.smb.issetpddfac = 1;
     42md.smb.pddfac_snow =4.3;
     43md.smb.pddfac_ice = 8.3;
     44
     45md=extrude(md,5,1.2);
     46md=setflowequation(md,'HO','all');
     47md.settings.results_on_nodes={'Temperature','Waterfraction','Enthalpy'};
     48
     49md.thermal.isenthalpy=1;
     50md.thermal.isdynamicbasalspc=1;
     51md.thermal.fe = 'P1xP2';
     52%md.thermal.spctemperature(find(md.mesh.vertexonbase)) = NaN;
     53
     54md.initialization.waterfraction= zeros(md.mesh.numberofvertices,1);
     55md.initialization.watercolumn  = zeros(md.mesh.numberofvertices,1);
     56md.initialization.enthalpy = zeros(md.mesh.numberofvertices,1);
     57md.thermal.isdrainicecolumn=0;
     58
     59md=solve(md,'thermal');
     60
     61md.initialization.temperature=md.results.ThermalSolution.Temperature;
     62md.initialization.enthalpy=md.results.ThermalSolution.Enthalpy;
     63%x = find(md.initialization.temperature<210);
     64%md.initialization.temperature(x) = 268;
     65md.materials.rheology_B = cuffey(md.initialization.temperature);
     66
     67% Friction
     68TEMP = zeros(md.mesh.numberofvertices,1);
     69TEMP(md.mesh.elements) = md.initialization.temperature(:,1:6);
     70
     71temperature =TEMP;
     72pressure = md.initialization.pressure;%md.materials.rho_ice*md.constants.g*md.geometry.thickness;
     73Tm = md.materials.meltingpoint-md.materials.beta*pressure;
     74
     75md.friction.coefficient = md.friction.coefficient;%averaging(md,md.friction.coefficient,2);
     76md.friction=frictionjosh(md.friction);
     77md.friction.pressure_adjusted_temperature = temperature - Tm;
     78md.friction.gamma= 5;
     79
     80
     81% time steps and resolution
     82md.timestepping.time_step=0.5;
     83md.settings.output_frequency=1;
     84md.timestepping.final_time=2;
     85md.timestepping.interp_forcings=0;
     86
     87%md.transient.requested_outputs={'default','IceVolumeAboveFloatation','IceVolume','TemperaturePDD','SmbMonthlytemperatures','SmbPrecipitation'};
     88md.transient.requested_outputs={'default','IceVolumeAboveFloatation','IceVolume','TemperaturePDD'};
     89md=setflowequation(md,'SSA','all');
     90md.cluster=generic('name',oshostname(),'np',1); % 3 for the cluster
     91md=solve(md,'Transient');
     92
     93field_names = {...
     94        'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','Enthalpy1','SmbMassBalance1', ...
     95        'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','Enthalpy2','SmbMassBalance2', ...
     96        'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','Enthalpy3','SmbMassBalance3', ...
     97        'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','Temperature4','Enthalpy4','SmbMassBalance4'};
     98field_tolerances={...
     99        1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13,...
     100        1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13,...
     101        1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13,...
     102        1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-13};
     103field_values={...
     104        (md.results.TransientSolution(1).Vx),...
     105        (md.results.TransientSolution(1).Vy),...
     106        (md.results.TransientSolution(1).Vz),...
     107        (md.results.TransientSolution(1).Vel),...
     108        (md.results.TransientSolution(1).Pressure),...
     109        (md.results.TransientSolution(1).Base),...
     110        (md.results.TransientSolution(1).Surface),...
     111        (md.results.TransientSolution(1).Thickness),...
     112        (md.results.TransientSolution(1).Temperature),...
     113        (md.results.TransientSolution(1).Enthalpy),...
     114        (md.results.TransientSolution(1).SmbMassBalance),...
     115        (md.results.TransientSolution(2).Vx),...
     116        (md.results.TransientSolution(2).Vy),...
     117        (md.results.TransientSolution(2).Vz),...
     118        (md.results.TransientSolution(2).Vel),...
     119        (md.results.TransientSolution(2).Pressure),...
     120        (md.results.TransientSolution(2).Base),...
     121        (md.results.TransientSolution(2).Surface),...
     122        (md.results.TransientSolution(2).Thickness),...
     123        (md.results.TransientSolution(2).Temperature),...
     124        (md.results.TransientSolution(2).Enthalpy),...
     125        (md.results.TransientSolution(2).SmbMassBalance),...
     126        (md.results.TransientSolution(3).Vx),...
     127        (md.results.TransientSolution(3).Vy),...
     128        (md.results.TransientSolution(3).Vz),...
     129        (md.results.TransientSolution(3).Vel),...
     130        (md.results.TransientSolution(3).Pressure),...
     131        (md.results.TransientSolution(3).Base),...
     132        (md.results.TransientSolution(3).Surface),...
     133        (md.results.TransientSolution(3).Thickness),...
     134        (md.results.TransientSolution(3).Temperature),...
     135        (md.results.TransientSolution(3).Enthalpy),...
     136        (md.results.TransientSolution(3).SmbMassBalance),...
     137        (md.results.TransientSolution(4).Vx),...
     138        (md.results.TransientSolution(4).Vy),...
     139        (md.results.TransientSolution(4).Vz),...
     140        (md.results.TransientSolution(4).Vel),...
     141        (md.results.TransientSolution(4).Pressure),...
     142        (md.results.TransientSolution(4).Base),...
     143        (md.results.TransientSolution(4).Surface),...
     144        (md.results.TransientSolution(4).Thickness),...
     145        (md.results.TransientSolution(4).Temperature),...
     146        (md.results.TransientSolution(4).Enthalpy),...
     147        (md.results.TransientSolution(4).SmbMassBalance),...
     148        };

  • ../trunk-jpl/test/Archives/Archive247.arch 

    Cannot display: file marked as a binary type.
svn:mime-type = application/octet-stream
Note: See TracBrowser for help on using the repository browser.