Changeset 22726
- Timestamp:
- 04/30/18 10:09:55 (7 years ago)
- Location:
- issm/trunk-jpl/test
- Files:
-
- 1 added
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test4002.m
r22723 r22726 846 846 loaddata(org,'Bathymetry'); 847 847 loaddata(org,'IceSheetGeometry'); 848 endtime = round(MITgcmDeltaT * floor(final_time*y2s/MITgcmDeltaT)); 848 endtime = round(MITgcmDeltaT * floor(final_time*y2s/MITgcmDeltaT)); 849 outputtime = round(MITgcmDeltaT * floor(time_step*y2s/MITgcmDeltaT)); 849 850 850 851 % {{{ prepare MITgcm … … 967 968 end 968 969 if strcmpi(tline,' taveFreq=2592000.,'), 969 fprintf(fido,'%s%i%s\n',' taveFreq= ', endtime,',');970 fprintf(fido,'%s%i%s\n',' taveFreq= ',outputtime,','); 970 971 continue; 971 972 end … … 1002 1003 1003 1004 %Fields and tolerances to track changes 1004 fnm=['run/SHICE_fwFluxtave.0000000 146.data'];1005 fnm=['run/SHICE_fwFluxtave.0000000006.data']; 1005 1006 melting_rate_1=readbin(fnm,[Nx Ny]); 1006 melting_rate_2=melting_rate_1; 1007 melting_rate_3=melting_rate_1; 1008 melting_rate_4=melting_rate_1; 1009 %fnm=['run/SHICE_fwFluxtave_0000004380.data']; 1010 %melting_rate_1=readbin(fnm,[Nx Ny]); 1011 %fnm=['run/SHICE_fwFluxtave_0000008760.data']; 1012 %melting_rate_2=readbin(fnm,[Nx Ny]); 1013 %fnm=['run/SHICE_fwFluxtave_0000013140.data']; 1014 %melting_rate_3=readbin(fnm,[Nx Ny]); 1015 %fnm=['run/SHICE_fwFluxtave_0000017520.data']; 1016 %melting_rate_4=readbin(fnm,[Nx Ny]); 1007 fnm=['run/SHICE_fwFluxtave.0000000012.data']; 1008 melting_rate_2=readbin(fnm,[Nx Ny]); 1009 fnm=['run/SHICE_fwFluxtave.0000000018.data']; 1010 melting_rate_3=readbin(fnm,[Nx Ny]); 1011 fnm=['run/SHICE_fwFluxtave.0000000024.data']; 1012 melting_rate_4=readbin(fnm,[Nx Ny]); 1017 1013 field_names ={'Base1','Melting1','Vx2','Vy2','Thickness2','Base2','MaskGroundediceLevelset2','FloatingiceMeltingRate2',... 1018 1014 'Melting2','Vx3','Vy3','Thickness3','Base3','MaskGroundediceLevelset3','FloatingiceMeltingRate3',... … … 1046 1042 (melting_rate_4(:)),... 1047 1043 }; 1048 %field_values={...1049 % (md.results(1).TransientSolution(end).Base),...1050 % (melting_rate_1(:)),...1051 % (md.results(2).TransientSolution(end).Vx),...1052 % (md.results(2).TransientSolution(end).Vy),...1053 % (md.results(2).TransientSolution(end).Thickness),...1054 % (md.results(2).TransientSolution(end).Base),...1055 % (md.results(2).TransientSolution(end).MaskGroundediceLevelset),...1056 % (md.results(2).TransientSolution(end).FloatingiceMeltingRate),...1057 % (melting_rate_2(:)),...1058 % (md.results(3).TransientSolution(end).Vx),...1059 % (md.results(3).TransientSolution(end).Vy),...1060 % (md.results(3).TransientSolution(end).Thickness),...1061 % (md.results(3).TransientSolution(end).Base),...1062 % (md.results(3).TransientSolution(end).MaskGroundediceLevelset),...1063 % (md.results(3).TransientSolution(end).FloatingiceMeltingRate),...1064 % (melting_rate_3(:)),...1065 % (md.results(4).TransientSolution(end).Vx),...1066 % (md.results(4).TransientSolution(end).Vy),...1067 % (md.results(4).TransientSolution(end).Thickness),...1068 % (md.results(4).TransientSolution(end).Base),...1069 % (md.results(4).TransientSolution(end).MaskGroundediceLevelset),...1070 % (md.results(4).TransientSolution(end).FloatingiceMeltingRate),...1071 % (melting_rate_4(:)),...1072 % };
Note:
See TracChangeset
for help on using the changeset viewer.