Changeset 19548


Ignore:
Timestamp:
09/10/15 11:11:26 (10 years ago)
Author:
Mathieu Morlighem
Message:

BUG: do isthermal even if issmb=0

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/cores/transient_core.cpp

    r19545 r19548  
    8686                femmodel->parameters->SetParam(save_results,SaveResultsEnum);
    8787
    88                 if(issmb && isthermal && domaintype==Domain3DEnum){
    89                         femmodel->parameters->FindParam(&smb_model,SmbEnum);
    90                         if(VerboseSolution()) _printf0_("   computing thermal regime\n");
    91                         if(smb_model==SMBpddEnum) ResetBoundaryConditions(femmodel,ThermalAnalysisEnum);
    92                         if(smb_model==SMBd18opddEnum) ResetBoundaryConditions(femmodel,ThermalAnalysisEnum);
     88                if(isthermal && domaintype==Domain3DEnum){
     89                        if(issmb){
     90                                femmodel->parameters->FindParam(&smb_model,SmbEnum);
     91                                if(VerboseSolution()) _printf0_("   computing thermal regime\n");
     92                                if(smb_model==SMBpddEnum) ResetBoundaryConditions(femmodel,ThermalAnalysisEnum);
     93                                if(smb_model==SMBd18opddEnum) ResetBoundaryConditions(femmodel,ThermalAnalysisEnum);
     94                        }
    9395                        thermal_core(femmodel);
    9496                }
Note: See TracChangeset for help on using the changeset viewer.