Ignore:
Timestamp:
08/28/18 09:45:51 (7 years ago)
Author:
Mathieu Morlighem
Message:

merged trunk-jpl and trunk for revision 23187

Location:
issm/trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk

  • issm/trunk/test

  • issm/trunk/test/NightlyRun/test4001.m

    r22758 r23189  
    136136
    137137    %specify computational grid in SIZE.h
    138     fidi=fopen('../MITgcm/code/SIZE.h.bak','r');
    139     fido=fopen('../MITgcm/code/SIZE.h','w');
    140     tline = fgetl(fidi);
    141     fprintf(fido,'%s\n',tline);
    142     while 1
     138    if ~exist('../MITgcm/code/SIZE.h')
     139        fidi=fopen('../MITgcm/code/SIZE.h.bak','r');
     140        fido=fopen('../MITgcm/code/SIZE.h','w');
    143141        tline = fgetl(fidi);
    144         if ~ischar(tline), break, end
    145         %do the change here:
    146         if strcmpi(tline,'     &           sNx =  20,'),
    147             fprintf(fido,'%s%i%s\n','     &           sNx =  ',round(Nx/nPx),',');
    148             continue;
     142        fprintf(fido,'%s\n',tline);
     143        while 1
     144            tline = fgetl(fidi);
     145            if ~ischar(tline), break, end
     146            %do the change here:
     147            if strcmpi(tline,'     &           sNx =  20,'),
     148                fprintf(fido,'%s%i%s\n','     &           sNx =  ',round(Nx/nPx),',');
     149                continue;
     150            end
     151            if strcmpi(tline,'     &           sNy =  20,'),
     152                fprintf(fido,'%s%i%s\n','     &           sNy =  ',round(Ny/nPy),',');
     153                continue;
     154            end
     155            if strcmpi(tline,'     &           nPx =   1,'),
     156                fprintf(fido,'%s%i%s\n','     &           nPx = ',nPx,',');
     157                continue;
     158            end
     159            if strcmpi(tline,'     &           nPy =   2,'),
     160                fprintf(fido,'%s%i%s\n','     &           nPy = ',nPy,',');
     161                continue;
     162            end
     163            fprintf(fido,'%s\n',tline);
    149164        end
    150         if strcmpi(tline,'     &           sNy =  20,'),
    151             fprintf(fido,'%s%i%s\n','     &           sNy =  ',round(Ny/nPy),',');
    152             continue;
    153         end
    154         if strcmpi(tline,'     &           nPx =   1,'),
    155             fprintf(fido,'%s%i%s\n','     &           nPx = ',nPx,',');
    156             continue;
    157         end
    158         if strcmpi(tline,'     &           nPy =   2,'),
    159             fprintf(fido,'%s%i%s\n','     &           nPy = ',nPy,',');
    160             continue;
    161         end
    162         fprintf(fido,'%s\n',tline);
     165        %close  files
     166        fclose(fidi);
     167        fclose(fido);
    163168    end
    164     %close  files
    165     fclose(fidi);
    166     fclose(fido);
    167 
    168          system(['../MITgcm/build.sh generic ' pwd '/../MITgcm']);
     169
     170    system(['../MITgcm/build.sh generic ' pwd '/../MITgcm']);
    169171end
    170172% }}}
     
    542544        md.transient.isgroundingline=1;
    543545        md.transient.isthermal=0;
    544         md.groundingline.migration='SubelementMigration2';
     546        md.groundingline.migration='SubelementMigration';
     547        md.groundingline.melt_interpolation='SubelementMelt2';
     548        md.groundingline.friction_interpolation='SubelementFriction2';
    545549
    546550        md.cluster=generic('name',oshostname(),'np',2);
     
    655659        md.transient.isgroundingline=1;
    656660        md.transient.isthermal=0;
    657         md.groundingline.migration='SubelementMigration2';
     661        md.groundingline.migration='SubelementMigration';
     662        md.groundingline.melt_interpolation='SubelementMelt2';
     663        md.groundingline.friction_interpolation='SubelementFriction2';
    658664
    659665        % }}}
Note: See TracChangeset for help on using the changeset viewer.