|
Last change
on this file since 2794 was
2794,
checked in by llange, 2 years ago
|
|
MARS PEM:
- Add a PEMETAT0 that read "startfi_pem.nc"
- Add the soil in the model: soil temperature, thermal properties, ice table
- Add a routine that compute CO2 + H2O adsorption
- Minor corrections in PEM.F90
LL
|
|
File size:
968 bytes
|
| Line | |
|---|
| 1 | ! Update surface ice distribution to iterate to steady state if requested |
|---|
| 2 | if(ice_update)then |
|---|
| 3 | |
|---|
| 4 | do ig=1,ngrid |
|---|
| 5 | |
|---|
| 6 | delta_ice = (qsurf(ig,igcm_h2o_ice)-ice_initial(ig)) |
|---|
| 7 | |
|---|
| 8 | ! add multiple years of evolution |
|---|
| 9 | qsurf_hist(ig,igcm_h2o_ice) = qsurf_hist(ig,igcm_h2o_ice) + delta_ice*icetstep |
|---|
| 10 | |
|---|
| 11 | ! if ice has gone -ve, set to zero |
|---|
| 12 | if(qsurf_hist(ig,igcm_h2o_ice).lt.0.0)then |
|---|
| 13 | qsurf_hist(ig,igcm_h2o_ice) = 0.0 |
|---|
| 14 | endif |
|---|
| 15 | |
|---|
| 16 | ! if ice is seasonal, set to zero (NEW) |
|---|
| 17 | if(ice_min(ig).lt.0.01)then |
|---|
| 18 | qsurf_hist(ig,igcm_h2o_ice) = 0.0 |
|---|
| 19 | endif |
|---|
| 20 | |
|---|
| 21 | enddo |
|---|
| 22 | |
|---|
| 23 | ! enforce ice conservation |
|---|
| 24 | ice_tot= SUM(qsurf_hist(:,igcm_h2o_ice)*cell_area(:) )/SUM(cell_area(:)) |
|---|
| 25 | qsurf_hist(:,igcm_h2o_ice) = qsurf_hist(:,igcm_h2o_ice)*(icesrf/ice_tot) |
|---|
| 26 | |
|---|
| 27 | endif |
|---|
Note: See
TracBrowser
for help on using the repository browser.
