Changeset 27599
- Timestamp:
- 02/20/23 18:19:27 (2 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test4005.m
r27591 r27599 3 3 % 4 4 %Script control parameters 5 steps= 1:11;5 steps=[1:5 7:11]; 6 6 final_time=1; 7 7 … … 648 648 results=md.results; 649 649 650 for t=0:time_step: final_time650 for t=0:time_step:4*time_step 651 651 disp(['Year: ' num2str(t)]) 652 652 653 653 %Initialize time steps 654 md.timestepping.final_time= (t+1)*time_step;654 md.timestepping.final_time=t+time_step; 655 655 md.timestepping.time_step=time_step; 656 md.timestepping.start_time=t *time_step;;656 md.timestepping.start_time=t; 657 657 658 658 %Calculate dynamic thinning and write … … 669 669 % {{{ system call to run MITgcm 670 670 cd run 671 672 newline = [' niter0 = ' num2str(t*y2s/MITgcmDeltaT)]; 673 command=['!sed "s/.*niter0.*/' newline '/" data > data.temp; mv data.temp data']; 674 eval(command) 675 676 671 677 ds=round(endtime/MITgcmDeltaT); 672 678 ts=round((t+time_step)*y2s/MITgcmDeltaT); 679 673 680 eval(['!mpirun -np ' int2str(nPx*nPy) ' ./mitgcmuv']); 674 eval(['!\mv STDERR.0000 STDERR_' myint2str(t ) '.data'])675 eval(['!\mv STDOUT.0000 STDOUT_' myint2str(t ) '.data'])676 eval(['!\cp hFacC.data hFacC_' myint2str(t ) '.data'])677 eval(['!\cp icetopo.bin icetopo_' myint2str(t ) '.data'])678 for fld={'S','T','U','V','Eta', ...679 'SHICE_heatFluxtave','SHICE_fwFluxtave'}680 eval(['!\mv ' fld{1} '.' myint2str(ds,10) '.data ' ...681 fld{1} '_' myint2str(t,10) '.data'])682 end681 eval(['!\mv STDERR.0000 STDERR_' myint2str(t/time_step) '.data']) 682 eval(['!\mv STDOUT.0000 STDOUT_' myint2str(t/time_step) '.data']) 683 eval(['!\cp hFacC.data hFacC_' myint2str(t/time_step) '.data']) 684 eval(['!\cp icetopo.bin icetopo_' myint2str(t/time_step) '.data']) 685 % for fld={'S','T','U','V','Eta', ... 686 % 'SHICE_heatFluxtave','SHICE_fwFluxtave'} 687 % eval(['!\mv ' fld{1} '.' myint2str(ds,10) '.data ' ... 688 % fld{1} '_' myint2str(t,10) '.data']) 689 % end 683 690 cd .. 684 691 % }}} 685 692 686 system(['cp run/SHICE_fwFluxtave.' myint2str((t+1)*2160,10) '.data run/melt.data']) 687 system(['mv run/SHICE_fwFluxtave.' myint2str((t+1)*2160,10) '.data run/melt' int2str(t) '.data']) 693 system(['cp run/SHICE_fwFluxtave.' myint2str((t+time_step)*y2s/MITgcmDeltaT,10) '.data run/melt.data']) 694 % system(['cp run/SHICE_fwFluxtave.' myint2str((t+1)*2160,10) '.data run/melt.data']) 695 % system(['mv run/SHICE_fwFluxtave.' myint2str((t+1)*2160,10) '.data run/melt' int2str(t) '.data']) 688 696 689 697 melt=readbin('./run/melt.data',[Nx,Ny]); … … 713 721 714 722 %Fields and tolerances to track changes 715 fnm=['run/SHICE_fwFluxtave _0000004380.data'];723 fnm=['run/SHICE_fwFluxtave.0000004380.data']; 716 724 melting_rate_1=readbin(fnm,[Nx Ny]); 717 fnm=['run/SHICE_fwFluxtave _0000008760.data'];725 fnm=['run/SHICE_fwFluxtave.0000008760.data']; 718 726 melting_rate_2=readbin(fnm,[Nx Ny]); 719 fnm=['run/SHICE_fwFluxtave _0000013140.data'];727 fnm=['run/SHICE_fwFluxtave.0000013140.data']; 720 728 melting_rate_3=readbin(fnm,[Nx Ny]); 721 fnm=['run/SHICE_fwFluxtave _0000017520.data'];729 fnm=['run/SHICE_fwFluxtave.0000017520.data']; 722 730 melting_rate_4=readbin(fnm,[Nx Ny]); 723 731 field_names ={'Base1','Melting1','Vx2','Vy2','Thickness2','Base2','MaskOceanLevelset2','FloatingiceMeltingRate2',...
Note:
See TracChangeset
for help on using the changeset viewer.