Changeset 22670


Ignore:
Timestamp:
04/03/18 21:20:54 (7 years ago)
Author:
seroussi
Message:

FIX: test now running with coupled model

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r22627 r22670  
    33%
    44%Script control parameters
    5 steps=12;
     5steps=[1 2 3 4 5 7 8 9 10 12];
    66final_time=1;
    77
    … …  
    843843if perform(org,'RunCoupledMITgcmISSM2'),
    844844       
     845        loaddata(org,'Parameters');
    845846        loaddata(org,'ParameterizeIce');
     847
     848        if exist ('run')
     849                !\rm -rf run/test4002*
     850        end
    846851       
    847         md.transient.isoceancoupling=0;
    848         md.miscellaneous.name='test4002';
    849         md.results=[];
    850         md.results.TransientSolution=[];
     852        md.transient.isoceancoupling=1;
     853        md.transient.isgroundingline=1;
    851854        md.timestepping.coupling_time=1/12;
    852         md=solve(md,'Transient','runtimename',false);
    853 
    854         md.transient.isoceancoupling=1;
    855         md.timestepping.coupling_time=1/12;
    856 
    857         md=solve(md,'Transient','runtimename',false);
    858 
    859         eval(['!cp ' md.cluster.executionpath '/' md.miscellaneous.name '/* run/ ']);
    860         cd run
    861         %eval(['!mpiexec -np ' int2str(md.cluster.np) ' ' md.cluster.codepath '/issm_ocean.exe TransientSolution ' pwd ' ' md.miscellaneous.name ' ']);
    862         eval(['!mpiexec -np ' int2str(md.cluster.np) ' ' md.cluster.codepath '/issm_ocean.exe TransientSolution ' pwd ' ' md.miscellaneous.name ' : -np ' int2str(nPx*nPy) ' ./mitgcmuv']);
    863 
    864 end
    865 % }}}
    866 
    867 
    868 
    869 error
     855        md.cluster.npocean=nPx*nPy;
     856        md.cluster.executionpath=[pwd '/run'];
     857        md.masstransport.requested_outputs={'default','BasalforcingsFloatingiceMeltingRate'};
     858
     859        md=solveiceocean(md,'Transient','runtimename',false);
     860
     861%       %eval(['!mpiexec -np ' int2str(md.cluster.np) ' ' md.cluster.codepath '/issm_ocean.exe TransientSolution ' pwd ' ' md.miscellaneous.name ' ']);
     862%       eval(['!mpiexec -np ' int2str(md.cluster.np) ' ' md.cluster.codepath '/issm_ocean.exe TransientSolution ' pwd ' ' md.miscellaneous.name ' : -np ' int2str(nPx*nPy) ' ./mitgcmuv']);
     863end
     864% }}}
     865
    870866%Fields and tolerances to track changes
    871 fnm=['run/SHICE_fwFluxtave_0000004380.data'];
     867fnm=['run/SHICE_fwFluxtave.0000000146.data'];
    872868melting_rate_1=readbin(fnm,[Nx Ny]);
    873 fnm=['run/SHICE_fwFluxtave_0000008760.data'];
    874 melting_rate_2=readbin(fnm,[Nx Ny]);
    875 fnm=['run/SHICE_fwFluxtave_0000013140.data'];
    876 melting_rate_3=readbin(fnm,[Nx Ny]);
    877 fnm=['run/SHICE_fwFluxtave_0000017520.data'];
    878 melting_rate_4=readbin(fnm,[Nx Ny]);
     869melting_rate_2=melting_rate_1;
     870melting_rate_3=melting_rate_1;
     871melting_rate_4=melting_rate_1;
     872%fnm=['run/SHICE_fwFluxtave_0000004380.data'];
     873%melting_rate_1=readbin(fnm,[Nx Ny]);
     874%fnm=['run/SHICE_fwFluxtave_0000008760.data'];
     875%melting_rate_2=readbin(fnm,[Nx Ny]);
     876%fnm=['run/SHICE_fwFluxtave_0000013140.data'];
     877%melting_rate_3=readbin(fnm,[Nx Ny]);
     878%fnm=['run/SHICE_fwFluxtave_0000017520.data'];
     879%melting_rate_4=readbin(fnm,[Nx Ny]);
    879880field_names     ={'Base1','Melting1','Vx2','Vy2','Thickness2','Base2','MaskGroundediceLevelset2','FloatingiceMeltingRate2',...
    880881        'Melting2','Vx3','Vy3','Thickness3','Base3','MaskGroundediceLevelset3','FloatingiceMeltingRate3',...
    … …  
    884885        1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13, 1e-13 };
    885886field_values={...
    886         (md.results(1).TransientSolution(end).Base),...
     887        (md.results.TransientSolution(1).Base),...
    887888        (melting_rate_1(:)),...
    888         (md.results(2).TransientSolution(end).Vx),...
    889         (md.results(2).TransientSolution(end).Vy),...
    890         (md.results(2).TransientSolution(end).Thickness),...
    891         (md.results(2).TransientSolution(end).Base),...
    892         (md.results(2).TransientSolution(end).MaskGroundediceLevelset),...
    893         (md.results(2).TransientSolution(end).FloatingiceMeltingRate),...
     889        (md.results.TransientSolution(2).Vx),...
     890        (md.results.TransientSolution(2).Vy),...
     891        (md.results.TransientSolution(2).Thickness),...
     892        (md.results.TransientSolution(2).Base),...
     893        (md.results.TransientSolution(2).MaskGroundediceLevelset),...
     894        (md.results.TransientSolution(2).BasalforcingsFloatingiceMeltingRate),...
    894895        (melting_rate_2(:)),...
    895         (md.results(3).TransientSolution(end).Vx),...
    896         (md.results(3).TransientSolution(end).Vy),...
    897         (md.results(3).TransientSolution(end).Thickness),...
    898         (md.results(3).TransientSolution(end).Base),...
    899         (md.results(3).TransientSolution(end).MaskGroundediceLevelset),...
    900         (md.results(3).TransientSolution(end).FloatingiceMeltingRate),...
     896        (md.results.TransientSolution(3).Vx),...
     897        (md.results.TransientSolution(3).Vy),...
     898        (md.results.TransientSolution(3).Thickness),...
     899        (md.results.TransientSolution(3).Base),...
     900        (md.results.TransientSolution(3).MaskGroundediceLevelset),...
     901        (md.results.TransientSolution(3).BasalforcingsFloatingiceMeltingRate),...
    901902        (melting_rate_3(:)),...
    902         (md.results(4).TransientSolution(end).Vx),...
    903         (md.results(4).TransientSolution(end).Vy),...
    904         (md.results(4).TransientSolution(end).Thickness),...
    905         (md.results(4).TransientSolution(end).Base),...
    906         (md.results(4).TransientSolution(end).MaskGroundediceLevelset),...
    907         (md.results(4).TransientSolution(end).FloatingiceMeltingRate),...
     903        (md.results.TransientSolution(4).Vx),...
     904        (md.results.TransientSolution(4).Vy),...
     905        (md.results.TransientSolution(4).Thickness),...
     906        (md.results.TransientSolution(4).Base),...
     907        (md.results.TransientSolution(4).MaskGroundediceLevelset),...
     908        (md.results.TransientSolution(4).BasalforcingsFloatingiceMeltingRate),...
    908909        (melting_rate_4(:)),...
    909910        };
     911%field_values={...
     912%       (md.results(1).TransientSolution(end).Base),...
     913%       (melting_rate_1(:)),...
     914%       (md.results(2).TransientSolution(end).Vx),...
     915%       (md.results(2).TransientSolution(end).Vy),...
     916%       (md.results(2).TransientSolution(end).Thickness),...
     917%       (md.results(2).TransientSolution(end).Base),...
     918%       (md.results(2).TransientSolution(end).MaskGroundediceLevelset),...
     919%       (md.results(2).TransientSolution(end).FloatingiceMeltingRate),...
     920%       (melting_rate_2(:)),...
     921%       (md.results(3).TransientSolution(end).Vx),...
     922%       (md.results(3).TransientSolution(end).Vy),...
     923%       (md.results(3).TransientSolution(end).Thickness),...
     924%       (md.results(3).TransientSolution(end).Base),...
     925%       (md.results(3).TransientSolution(end).MaskGroundediceLevelset),...
     926%       (md.results(3).TransientSolution(end).FloatingiceMeltingRate),...
     927%       (melting_rate_3(:)),...
     928%       (md.results(4).TransientSolution(end).Vx),...
     929%       (md.results(4).TransientSolution(end).Vy),...
     930%       (md.results(4).TransientSolution(end).Thickness),...
     931%       (md.results(4).TransientSolution(end).Base),...
     932%       (md.results(4).TransientSolution(end).MaskGroundediceLevelset),...
     933%       (md.results(4).TransientSolution(end).FloatingiceMeltingRate),...
     934%       (melting_rate_4(:)),...
     935%       };
Note: See TracChangeset for help on using the changeset viewer.