Changeset 24122
- Timestamp:
- 08/26/19 11:27:09 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test4003.m
r23989 r24122 151 151 !\cp ../MITgcm/build/mitgcmuv run 152 152 !\cp ../MITgcm/install_dngoldberg/verification/shelfice_remeshing/input/* run 153 !\cp ../MITgcm/input_4003/* run 154 !\cp ../MITgcm/input_4003/eedata_uncoupled run/eedata 153 155 cd run 154 156 eval(['!mpirun -np ' int2str(nPx*nPy) ' ./mitgcmuv']); … … 360 362 !\cp ../MITgcm/build/mitgcmuv run 361 363 !\cp ../MITgcm/install_dngoldberg/verification/shelfice_remeshing/input/* run 364 !\cp ../MITgcm/input_4003/* run 362 365 % }}} 363 % {{{ generate MITgcm parameter file data364 fidi=fopen('../MITgcm/input/data','r');365 fido=fopen('run/data','w');366 tline = fgetl(fidi);367 fprintf(fido,'%s\n',tline);368 while 1369 tline = fgetl(fidi);370 if ~ischar(tline), break, end371 %do the change here:372 if strcmpi(tline,' xgOrigin = 0.0,'),373 fprintf(fido,'%s%i%s\n',' xgOrigin = ',xgOrigin,',');374 continue;375 end376 if strcmpi(tline,' ygOrigin = -80.0,'),377 fprintf(fido,'%s%i%s\n',' ygOrigin = ',ygOrigin,',');378 continue;379 end380 if strcmpi(tline,' delX = 20*0.25,'),381 fprintf(fido,'%s%i*%g%s\n',' delX = ',Nx,dLong,',');382 continue;383 end384 if strcmpi(tline,' delY = 20*0.25,'),385 fprintf(fido,'%s%i*%g%s\n',' delY = ',Ny,dLat,',');386 continue;387 end388 if strcmpi(tline,' delZ = 30*30.0,'),389 fprintf(fido,'%s%i*%g%s\n',' delZ = ',Nz,delZ,',');390 continue;391 end392 if strcmpi(tline,' endTime=2592000.,'),393 fprintf(fido,'%s%i%s\n',' endTime= ',endtime,',');394 continue;395 end396 if strcmpi(tline,' deltaT=1200.0,'),397 fprintf(fido,'%s%i%s\n',' deltaT= ',MITgcmDeltaT,',');398 continue;399 end400 if strcmpi(tline,' pChkptFreq=2592000.,'),401 fprintf(fido,'%s%i%s\n',' pChkptFreq= ',endtime,',');402 continue;403 end404 if strcmpi(tline,' taveFreq=2592000.,'),405 fprintf(fido,'%s%i%s\n',' taveFreq= ',outputtime,',');406 continue;407 end408 if strcmpi(tline,' rhoConst=1030.,'),409 fprintf(fido,'%s%i%s\n',' rhoConst= ',rho_water,',');410 continue;411 end412 if strcmpi(tline,' rhoNil=1030.,'),413 fprintf(fido,'%s%i%s\n',' rhoNil= ',rho_water,',');414 continue;415 end416 fprintf(fido,'%s\n',tline);417 end418 %close files419 fclose(fidi);420 fclose(fido);421 % }}}422 423 366 md.timestepping.coupling_time=time_step; 424 367 md.timestepping.time_step=time_step;
Note:
See TracChangeset
for help on using the changeset viewer.