Changeset 27601


Ignore:
Timestamp:
02/21/23 07:04:32 (2 years ago)
Author:
seroussi
Message:

CHG: tests for coupled run not working yet so removed

File:
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/test/NightlyRun/test4005.m

    r27599 r27601  
    666666                 system(['cp run/shelfice_dmdt.bin run/shelfice_dmdt_' int2str(t) '.bin'])
    667667
    668 
    669668                 % {{{ system call to run MITgcm
    670669                 cd run
    671 
    672                  newline = [' niter0 = ' num2str(t*y2s/MITgcmDeltaT)];
     670       newline = [' niter0 = ' num2str(t*y2s/MITgcmDeltaT)];
    673671                 command=['!sed "s/.*niter0.*/' newline '/" data > data.temp; mv data.temp data'];
    674672                 eval(command)
    675 
    676673
    677674       ds=round(endtime/MITgcmDeltaT);
    … …  
    696693
    697694       melt=readbin('./run/melt.data',[Nx,Ny]);
    698 
    699695       md.basalforcings.floatingice_melting_rate=-melt(:)*y2s/rho_ice;
    700696                 md=solve(md,'Transient');
    … …  
    720716% }}}
    721717
    722 %Fields and tolerances to track changes
    723 fnm=['run/SHICE_fwFluxtave.0000004380.data'];
    724 melting_rate_1=readbin(fnm,[Nx Ny]);
    725 fnm=['run/SHICE_fwFluxtave.0000008760.data'];
    726 melting_rate_2=readbin(fnm,[Nx Ny]);
    727 fnm=['run/SHICE_fwFluxtave.0000013140.data'];
    728 melting_rate_3=readbin(fnm,[Nx Ny]);
    729 fnm=['run/SHICE_fwFluxtave.0000017520.data'];
    730 melting_rate_4=readbin(fnm,[Nx Ny]);
    731 field_names     ={'Base1','Melting1','Vx2','Vy2','Thickness2','Base2','MaskOceanLevelset2','FloatingiceMeltingRate2',...
    732     'Melting2','Vx3','Vy3','Thickness3','Base3','MaskOceanLevelset3','FloatingiceMeltingRate3',...
    733     'Melting3','Vx4','Vy4','Thickness4','Base4','MaskOceanLevelset4','FloatingiceMeltingRate4','Melting4'};
    734 field_tolerances={2e-13,1e-13,...
    735     8e-06,7e-06,2e-07,2e-08,3e-08,1e-13,1e-13,...
    736     8e-06,7e-06,4e-07,3e-08,5e-08,1e-13,6e-08,...
    737     8e-06,7e-06,5e-07,4e-08,8e-08,6e-08,3e-07};
    738 field_values={...
    739     (md.results(1).TransientSolution(end).Base),...
    740     (melting_rate_1(:)),...
    741     (md.results(2).TransientSolution(end).Vx),...
    742     (md.results(2).TransientSolution(end).Vy),...
    743     (md.results(2).TransientSolution(end).Thickness),...
    744     (md.results(2).TransientSolution(end).Base),...
    745     (md.results(2).TransientSolution(end).MaskOceanLevelset),...
    746     (md.results(2).TransientSolution(end).FloatingiceMeltingRate),...
    747     (melting_rate_2(:)),...
    748     (md.results(3).TransientSolution(end).Vx),...
    749     (md.results(3).TransientSolution(end).Vy),...
    750     (md.results(3).TransientSolution(end).Thickness),...
    751     (md.results(3).TransientSolution(end).Base),...
    752     (md.results(3).TransientSolution(end).MaskOceanLevelset),...
    753     (md.results(3).TransientSolution(end).FloatingiceMeltingRate),...
    754     (melting_rate_3(:)),...
    755     (md.results(4).TransientSolution(end).Vx),...
    756     (md.results(4).TransientSolution(end).Vy),...
    757     (md.results(4).TransientSolution(end).Thickness),...
    758     (md.results(4).TransientSolution(end).Base),...
    759     (md.results(4).TransientSolution(end).MaskOceanLevelset),...
    760     (md.results(4).TransientSolution(end).FloatingiceMeltingRate),...
    761     (melting_rate_4(:)),...
    762     };
     718%%Fields and tolerances to track changes
     719%fnm=['run/SHICE_fwFluxtave.0000004380.data'];
     720%melting_rate_1=readbin(fnm,[Nx Ny]);
     721%fnm=['run/SHICE_fwFluxtave.0000008760.data'];
     722%melting_rate_2=readbin(fnm,[Nx Ny]);
     723%fnm=['run/SHICE_fwFluxtave.0000013140.data'];
     724%melting_rate_3=readbin(fnm,[Nx Ny]);
     725%fnm=['run/SHICE_fwFluxtave.0000017520.data'];
     726%melting_rate_4=readbin(fnm,[Nx Ny]);
     727%field_names     ={'Base1','Melting1','Vx2','Vy2','Thickness2','Base2','MaskOceanLevelset2','FloatingiceMeltingRate2',...
     728%    'Melting2','Vx3','Vy3','Thickness3','Base3','MaskOceanLevelset3','FloatingiceMeltingRate3',...
     729%    'Melting3','Vx4','Vy4','Thickness4','Base4','MaskOceanLevelset4','FloatingiceMeltingRate4','Melting4'};
     730%field_tolerances={2e-13,1e-13,...
     731%    8e-06,7e-06,2e-07,2e-08,3e-08,1e-13,1e-13,...
     732%    8e-06,7e-06,4e-07,3e-08,5e-08,1e-13,6e-08,...
     733%    8e-06,7e-06,5e-07,4e-08,8e-08,6e-08,3e-07};
     734%field_values={...
     735%    (md.results(1).TransientSolution(end).Base),...
     736%    (melting_rate_1(:)),...
     737%    (md.results(2).TransientSolution(end).Vx),...
     738%    (md.results(2).TransientSolution(end).Vy),...
     739%    (md.results(2).TransientSolution(end).Thickness),...
     740%    (md.results(2).TransientSolution(end).Base),...
     741%    (md.results(2).TransientSolution(end).MaskOceanLevelset),...
     742%    (md.results(2).TransientSolution(end).FloatingiceMeltingRate),...
     743%    (melting_rate_2(:)),...
     744%    (md.results(3).TransientSolution(end).Vx),...
     745%    (md.results(3).TransientSolution(end).Vy),...
     746%    (md.results(3).TransientSolution(end).Thickness),...
     747%    (md.results(3).TransientSolution(end).Base),...
     748%    (md.results(3).TransientSolution(end).MaskOceanLevelset),...
     749%    (md.results(3).TransientSolution(end).FloatingiceMeltingRate),...
     750%    (melting_rate_3(:)),...
     751%    (md.results(4).TransientSolution(end).Vx),...
     752%    (md.results(4).TransientSolution(end).Vy),...
     753%    (md.results(4).TransientSolution(end).Thickness),...
     754%    (md.results(4).TransientSolution(end).Base),...
     755%    (md.results(4).TransientSolution(end).MaskOceanLevelset),...
     756%    (md.results(4).TransientSolution(end).FloatingiceMeltingRate),...
     757%    (melting_rate_4(:)),...
     758%    };
Note: See TracChangeset for help on using the changeset viewer.