| 1 | ! |
|---|
| 2 | ! $Id $ |
|---|
| 3 | ! |
|---|
| 4 | SUBROUTINE evol_co2_ice_s_slope(qsurf,tendencies_co2_ice_phys,& |
|---|
| 5 | iim_input,jjm_input,ngrid,cell_area,STOPPING,nslope) |
|---|
| 6 | |
|---|
| 7 | USE temps_mod_evol, ONLY: dt_pem |
|---|
| 8 | |
|---|
| 9 | IMPLICIT NONE |
|---|
| 10 | |
|---|
| 11 | !======================================================================= |
|---|
| 12 | ! |
|---|
| 13 | ! Routine that compute the evolution of the water ice |
|---|
| 14 | ! |
|---|
| 15 | !======================================================================= |
|---|
| 16 | |
|---|
| 17 | ! arguments: |
|---|
| 18 | ! ---------- |
|---|
| 19 | |
|---|
| 20 | ! INPUT |
|---|
| 21 | |
|---|
| 22 | INTEGER, intent(in) :: iim_input,jjm_input, ngrid,nslope ! # of grid points along longitude/latitude/ total |
|---|
| 23 | ! REAL, intent(in) :: tendencies_h2o_ice_phys(ngrid) ! physical point field : Evolution of perenial ice over one year |
|---|
| 24 | REAL, intent(in) :: cell_area(ngrid) |
|---|
| 25 | |
|---|
| 26 | ! OUTPUT |
|---|
| 27 | REAL, INTENT(INOUT) :: qsurf(ngrid,nslope) ! physical point field : Previous and actual density of water ice |
|---|
| 28 | LOGICAL :: STOPPING |
|---|
| 29 | REAL, intent(inout) :: tendencies_co2_ice_phys(ngrid,nslope) ! physical point field : Evolution of perenial ice over one year |
|---|
| 30 | |
|---|
| 31 | |
|---|
| 32 | ! local: |
|---|
| 33 | ! ---- |
|---|
| 34 | |
|---|
| 35 | INTEGER :: i,j,ig0,islope ! loop variable |
|---|
| 36 | ! REAL :: not_budget, budget |
|---|
| 37 | REAL :: pos_tend, neg_tend, real_coefficient,negative_part |
|---|
| 38 | REAL :: new_tendencies(ngrid) |
|---|
| 39 | |
|---|
| 40 | STOPPING=.false. |
|---|
| 41 | |
|---|
| 42 | |
|---|
| 43 | ! Evolution of the water ice for each physical point |
|---|
| 44 | do i=1,ngrid |
|---|
| 45 | do islope=1,nslope |
|---|
| 46 | qsurf(i,islope)=qsurf(i,islope)+tendencies_co2_ice_phys(i,islope)*dt_pem |
|---|
| 47 | if (qsurf(i,islope).lt.0) then |
|---|
| 48 | qsurf(i,islope)=0. |
|---|
| 49 | tendencies_co2_ice_phys(i,islope)=0. |
|---|
| 50 | endif |
|---|
| 51 | enddo |
|---|
| 52 | enddo |
|---|
| 53 | |
|---|
| 54 | |
|---|
| 55 | END SUBROUTINE evol_co2_ice_s_slope |
|---|
