Changeset 20900


Ignore:
Timestamp:
07/13/16 14:25:48 (9 years ago)
Author:
Mathieu Morlighem
Message:

CHG: yts is now consistent with md.constants.yts

Location:
issm/trunk-jpl/src/m/solve
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/m/solve/loadresultsfrombuffer.js

    r20299 r20900  
    1414
    1515                //load results:
    16                 structure=parseresultsfrombuffer(buffer,buffersize);
     16                structure=parseresultsfrombuffer(md,buffer,buffersize);
    1717
    1818                //load structure onto results:

  • issm/trunk-jpl/src/m/solve/loadresultsfromdisk.m

    r20821 r20900  
    3131
    3232        %load results onto model
    33         structure=parseresultsfromdisk(filename,~md.settings.io_gather);
     33        structure=parseresultsfromdisk(md,filename,~md.settings.io_gather);
    3434        if isempty(fieldnames(structure)),
    3535                error(['No result found in binary file ' filename '. Check for solution crash.']);

  • issm/trunk-jpl/src/m/solve/loadresultsfromdisk.py

    r17480 r20900  
    2828
    2929                #load results onto model
    30                 structure=parseresultsfromdisk(filename,not md.settings.io_gather)
     30                structure=parseresultsfromdisk(md,filename,not md.settings.io_gather)
    3131                if not len(structure):
    3232                        raise RuntimeError("No result found in binary file '%s'. Check for solution crash." % filename)

  • issm/trunk-jpl/src/m/solve/parseresultsfrombuffer.js

    r19831 r20900  
    1 function parseresultsfrombuffer(buffer,buffersize){ //{{{
    2 //PARSERESULTSFROMDISK - ...
    3 //
    4 //   Usage:
    5 //      results=parseresultsfrombuffer(buffer,buffersize);
     1function parseresultsfrombuffer(md,buffer,buffersize){ //{{{
    62
    73        //Open file
    … …  
    117
    128        //Read fields until the end of the file.
    13         result  = ReadData(fid);
     9        result  = ReadData(fid,md);
    1410
    1511        if (MapIsEmpty(result))throw Error('no results found in binary memory buffer');
    … …  
    5753
    5854                //read next result
    59                 result  = ReadData(fid);
     55                result  = ReadData(fid,md);
    6056        }
    6157        return results;
    6258} // }}}
    63 function ReadData(fid){ //{{{
     59function ReadData(fid,md){ //{{{
    6460//READDATA - ...
    6561//
    6662//   Usage:
    67 //      field=ReadData(fid)
     63//      field=ReadData(fid,md)
    6864
    6965        //read field
    … …  
    8884
    8985                //Process units here FIXME: this should not be done here!
    90                 var yts=365.0*24.0*3600.0;
     86                var yts=md.constants.yts;
    9187                if (fieldname == 'BalancethicknessThickeningRate') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;
    9288                else if (fieldname == 'HydrologyWaterVx') for (var i=0;i<field.length;i++)field[i]= field[i]*yts;

  • issm/trunk-jpl/src/m/solve/parseresultsfromdisk.m

    r20212 r20900  
    1 function results=parseresultsfromdisk(filename,iosplit)
    2 %PARSERESULTSFROMDISK - ...
    3 %
    4 %   Usage:
    5 %      results=parseresultsfromdisk(filename,iosplit)
     1function results=parseresultsfromdisk(md,filename,iosplit)
    62
    73if iosplit,
    8         results=parseresultsfromdiskiosplit(filename);
     4        results=parseresultsfromdiskiosplit(md,filename);
    95else
    10         results=parseresultsfromdiskioserial(filename);
    11 end
    12 
    13 
    14 function results=parseresultsfromdiskioserial(filename) % {{{
    15 %PARSERESULTSFROMDISK - ...
    16 %
    17 %   Usage:
    18 %      results=parseresultsfromdiskioserial(filename)
     6        results=parseresultsfromdiskioserial(md,filename);
     7end
     8
     9
     10function results=parseresultsfromdiskioserial(md,filename) % {{{
    1911
    2012%Open file
    … …  
    2618
    2719%Read fields until the end of the file.
    28 result  = ReadData(fid);
     20result  = ReadData(fid,md);
    2921if isempty(result), error(['no results found in binary file ' filename]); end
    3022check_nomoresteps=0;
    … …  
    6658        %read next result
    6759        try,
    68                 result  = ReadData(fid);
     60                result  = ReadData(fid,md);
    6961        catch me,
    7062                disp('WARNING: file corrupted, trying partial recovery');
    … …  
    7668fclose(fid);
    7769% }}}
    78 function results=parseresultsfromdiskiosplit(filename) % {{{
    79 %PARSERESULTSFROMDISKIOSPLIT - ...
    80 %
    81 %   Usage:
    82 %      results=parseresultsfromdiskiosplit(filename)
     70function results=parseresultsfromdiskiosplit(md,filename) % {{{
    8371
    8472%Open file
    … …  
    117105%third pass, this time to read the real information
    118106fseek(fid,0,-1); %rewind
    119 result=ReadData(fid);
     107result=ReadData(fid,md);
    120108while ~isempty(result),
    121109
    … …  
    131119        %read next result
    132120        try,
    133                 result=ReadData(fid);
     121                result=ReadData(fid,md);
    134122        catch me,
    135123                disp('WARNING: file corrupted, results partial recovery');
    … …  
    142130fclose(fid);
    143131        % }}}
    144 function result=ReadData(fid) % {{{
    145 %READDATA - ...
    146 %
    147 %   Usage:
    148 %      field=ReadData(fid)
     132function result=ReadData(fid,md) % {{{
    149133
    150134%read field
    … …  
    175159
    176160        %Process units here FIXME: this should not be done here!
    177         yts=365.0*24.0*3600.0;
     161        yts=md.constants.yts;
    178162        if strcmp(fieldname,'BalancethicknessThickeningRate'),
    179163                field = field*yts;

  • issm/trunk-jpl/src/m/solve/parseresultsfromdisk.py

    r20222 r20900  
    55import MatlabFuncs as m
    66
    7 def parseresultsfromdisk(filename,iosplit):
    8         """
    9         PARSERESULTSFROMDISK - ...
    10 
    11            Usage:
    12               results=parseresultsfromdisk(filename,iosplit)
    13         """
    14 
     7def parseresultsfromdisk(md,filename,iosplit):
    158        if iosplit:
    16                 results=parseresultsfromdiskiosplit(filename)
     9                results=parseresultsfromdiskiosplit(md,filename)
    1710        else:
    18                 results=parseresultsfromdiskioserial(filename)
     11                results=parseresultsfromdiskioserial(md,filename)
    1912
    2013        return results
    2114
    22 def parseresultsfromdiskioserial(filename):    # {{{
    23         """
    24         PARSERESULTSFROMDISK - ...
    25          
    26             Usage:
    27                results=parseresultsfromdiskioserial(filename)
    28         """
     15def parseresultsfromdiskioserial(md,filename):    # {{{
    2916
    3017        #Open file
    … …  
    3926
    4027        #Read fields until the end of the file.
    41         result=ReadData(fid)
     28        result=ReadData(fid,md)
    4229
    4330        counter=0
    … …  
    9077
    9178                #read next result
    92                 result=ReadData(fid)
     79                result=ReadData(fid,md)
    9380
    9481        fid.close()
    … …  
    9683        return results
    9784        # }}}
    98 def parseresultsfromdiskiosplit(filename):    # {{{
    99         """
    100         PARSERESULTSFROMDISKIOSPLIT - ...
    101          
    102             Usage:
    103                results=parseresultsfromdiskiosplit(filename)
    104         """
     85def parseresultsfromdiskiosplit(md,filename):    # {{{
    10586
    10687        #Open file
    … …  
    141122        #third pass, this time to read the real information
    142123        fid.seek(0)    #rewind
    143         result=ReadData(fid)
     124        result=ReadData(fid,md)
    144125        while result:
    145126
    … …  
    156137
    157138                #read next result
    158                 result=ReadData(fid)
     139                result=ReadData(fid,md)
    159140
    160141        #close file
    … …  
    163144        return results
    164145        # }}}
    165 def ReadData(fid):    # {{{
     146def ReadData(fid,md):    # {{{
    166147        """
    167148        READDATA - ...
    168149         
    169150            Usage:
    170                field=ReadData(fid)
     151               field=ReadData(fid,md)
    171152        """
    172153
    … …  
    195176
    196177                #Process units here FIXME: this should not be done here!
    197                 yts=365.0*24.0*3600.0
     178                yts=md.constants.yts
    198179                if m.strcmp(fieldname,'BalancethicknessThickeningRate'):
    199180                        field = field*yts
    … …  
    224205                elif m.strcmp(fieldname,'CalvingCalvingrate'):
    225206                        field = field*yts
    226 
    227207
    228208                result=OrderedDict()
Note: See TracChangeset for help on using the changeset viewer.