- Timestamp:
- 10/04/21 06:35:33 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/shared/Elements/PddSurfaceMassBalanceSicopolis.cpp
r23400 r26468 32 32 IssmDouble Pmax = 0.6; 33 33 IssmDouble inv_twelve=1./12.; 34 34 35 35 sconv=(rho_water/rho_ice); //rhow_rain/rhoi 36 36 … … 91 91 IssmDouble coeff5=3.8e-05; // polynomial 92 92 IssmDouble coeff6=6.0e-07; // fit 93 93 94 94 if(tstar>=temp_rain) 95 95 frac_solid = 0.0; … … 100 100 +tstar*(coeff3+tstar*(coeff4+tstar*(coeff5+tstar*coeff6)))); 101 101 } 102 102 103 103 snowfall=snowfall+precip_star*frac_solid*inv_twelve; 104 104 } … … 108 108 if(snowfall<0.0) snowfall=0.0; // correction of 109 109 if(rainfall<0.0) rainfall=0.0; // negative values 110 110 111 111 if(rainfall<=(Pmax*snowfall)){ 112 112 if((rainfall+beta1*pdd)<=(Pmax*snowfall)) { … … 126 126 runoff = smelt+rainfall-Pmax*snowfall; 127 127 } 128 128 129 129 saccu = precip; 130 130
Note:
See TracChangeset
for help on using the changeset viewer.