Changeset 27599


Ignore:
Timestamp:
02/20/23 18:19:27 (2 years ago)
Author:
dgoldber
Message:

various changes to get coupled stepping working. do not run uncoupled mitgcm

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/test/NightlyRun/test4005.m

    r27591 r27599  
    33%
    44%Script control parameters
    5 steps=1:11;
     5steps=[1:5 7:11];
    66final_time=1;
    77
     
    648648    results=md.results;
    649649
    650     for t=0:time_step:final_time
     650    for t=0:time_step:4*time_step
    651651                 disp(['Year: ' num2str(t)])
    652652
    653653                 %Initialize time steps
    654                  md.timestepping.final_time=(t+1)*time_step;
     654                 md.timestepping.final_time=t+time_step;
    655655                 md.timestepping.time_step=time_step;
    656                  md.timestepping.start_time=t*time_step;;
     656                 md.timestepping.start_time=t;
    657657
    658658                 %Calculate dynamic thinning and write
     
    669669                 % {{{ system call to run MITgcm
    670670                 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
    671677       ds=round(endtime/MITgcmDeltaT);
    672678       ts=round((t+time_step)*y2s/MITgcmDeltaT);
     679                 
    673680                 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                  end
     681                 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
    683690                 cd ..
    684691                 % }}}
    685692
    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'])
    688696
    689697       melt=readbin('./run/melt.data',[Nx,Ny]);
     
    713721
    714722%Fields and tolerances to track changes
    715 fnm=['run/SHICE_fwFluxtave_0000004380.data'];
     723fnm=['run/SHICE_fwFluxtave.0000004380.data'];
    716724melting_rate_1=readbin(fnm,[Nx Ny]);
    717 fnm=['run/SHICE_fwFluxtave_0000008760.data'];
     725fnm=['run/SHICE_fwFluxtave.0000008760.data'];
    718726melting_rate_2=readbin(fnm,[Nx Ny]);
    719 fnm=['run/SHICE_fwFluxtave_0000013140.data'];
     727fnm=['run/SHICE_fwFluxtave.0000013140.data'];
    720728melting_rate_3=readbin(fnm,[Nx Ny]);
    721 fnm=['run/SHICE_fwFluxtave_0000017520.data'];
     729fnm=['run/SHICE_fwFluxtave.0000017520.data'];
    722730melting_rate_4=readbin(fnm,[Nx Ny]);
    723731field_names     ={'Base1','Melting1','Vx2','Vy2','Thickness2','Base2','MaskOceanLevelset2','FloatingiceMeltingRate2',...
Note: See TracChangeset for help on using the changeset viewer.