Changeset 26933


Ignore:
Timestamp:
03/11/22 10:41:47 (3 years ago)
Author:
aleahsommers
Message:

BUG: fixed PMPheat term in meltrate

File:
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp

    r26696 r26933  
    381381        dpressure_water[0] = rho_water*g*(dh[0] - dbed[0]);
    382382                dpressure_water[1] = rho_water*g*(dh[1] - dbed[1]);
    383                 PMPheat=-CT*CW*conductivity*(dh[0]*dpressure_water[0]+dh[1]*dpressure_water[1]);
     383                PMPheat=CT*CW*conductivity*rho_water*(dh[0]*dpressure_water[0]+dh[1]*dpressure_water[1]);
    384384
    385385        meltrate = 1/latentheat*(G+frictionheat+rho_water*g*conductivity*(dh[0]*dh[0]+dh[1]*dh[1])-PMPheat);
    … …  
    629629           dpressure_water[0] = rho_water*g*(dh[0] - dbed[0]);
    630630                dpressure_water[1] = rho_water*g*(dh[1] - dbed[1]);
    631                 PMPheat=-CT*CW*conductivity*(dh[0]*dpressure_water[0]+dh[1]*dpressure_water[1]);
     631                PMPheat=CT*CW*conductivity*rho_water*(dh[0]*dpressure_water[0]+dh[1]*dpressure_water[1]);
    632632
    633633                meltrate = 1/latentheat*(G+frictionheat+rho_water*g*conductivity*(dh[0]*dh[0]+dh[1]*dh[1])-PMPheat);
Note: See TracChangeset for help on using the changeset viewer.