| 1 | ! |
|---|
| 2 | ! $Id $ |
|---|
| 3 | ! |
|---|
| 4 | SUBROUTINE compute_tendencies(tendencies_h2o_ice,min_h2o_ice_Y1,& |
|---|
| 5 | min_h2o_ice_Y2,iim_input,jjm_input,ngrid,tendencies_h2o_ice_phys) |
|---|
| 6 | |
|---|
| 7 | IMPLICIT NONE |
|---|
| 8 | |
|---|
| 9 | |
|---|
| 10 | !======================================================================= |
|---|
| 11 | ! |
|---|
| 12 | ! Compute the tendencies of the evolution of water ice over the years |
|---|
| 13 | ! |
|---|
| 14 | !======================================================================= |
|---|
| 15 | |
|---|
| 16 | ! arguments: |
|---|
| 17 | ! ---------- |
|---|
| 18 | |
|---|
| 19 | ! INPUT |
|---|
| 20 | |
|---|
| 21 | INTEGER, intent(in) :: iim_input,jjm_input,ngrid ! # of grid points along longitude/latitude/ total |
|---|
| 22 | REAL, intent(in) , dimension(iim_input+1,jjm_input+1):: min_h2o_ice_Y1 ! LON x LAT field : minimum of water ice at each point for the first year |
|---|
| 23 | REAL, intent(in) , dimension(iim_input+1,jjm_input+1):: min_h2o_ice_Y2 ! LON x LAT field : minimum of water ice at each point for the second year |
|---|
| 24 | |
|---|
| 25 | ! OUTPUT |
|---|
| 26 | REAL, intent(out) , dimension(iim_input+1,jjm_input+1) :: tendencies_h2o_ice ! LON x LAT field : difference between the minima = evolution of perenial ice |
|---|
| 27 | REAL, intent(out) , dimension(ngrid) :: tendencies_h2o_ice_phys ! physical point field : difference between the minima = evolution of perenial ice |
|---|
| 28 | |
|---|
| 29 | ! local: |
|---|
| 30 | ! ------ |
|---|
| 31 | |
|---|
| 32 | INTEGER :: i,j,ig0 ! loop variable |
|---|
| 33 | |
|---|
| 34 | !======================================================================= |
|---|
| 35 | |
|---|
| 36 | |
|---|
| 37 | ! We compute the difference |
|---|
| 38 | tendencies_h2o_ice(:,:)=min_h2o_ice_Y2(:,:)-min_h2o_ice_Y1(:,:) |
|---|
| 39 | |
|---|
| 40 | ! If the difference is too small; there is no evolution |
|---|
| 41 | DO j=1,jjm_input+1 |
|---|
| 42 | DO i = 1, iim_input |
|---|
| 43 | if(abs(tendencies_h2o_ice(i,j)).LT.1.0E-10) then |
|---|
| 44 | tendencies_h2o_ice(i,j)=0. |
|---|
| 45 | endif |
|---|
| 46 | ENDDO |
|---|
| 47 | ENDDO |
|---|
| 48 | |
|---|
| 49 | |
|---|
| 50 | ! We reorganise the difference on the physical grid |
|---|
| 51 | tendencies_h2o_ice_phys(1)=tendencies_h2o_ice(1,1) |
|---|
| 52 | |
|---|
| 53 | ig0 = 2 |
|---|
| 54 | DO j=2,jjm_input |
|---|
| 55 | DO i = 1, iim_input |
|---|
| 56 | tendencies_h2o_ice_phys(ig0) =tendencies_h2o_ice(i,j) |
|---|
| 57 | ig0= ig0 + 1 |
|---|
| 58 | ENDDO |
|---|
| 59 | ENDDO |
|---|
| 60 | |
|---|
| 61 | tendencies_h2o_ice_phys(ig0) = tendencies_h2o_ice(1,jjm_input+1) |
|---|
| 62 | |
|---|
| 63 | |
|---|
| 64 | END SUBROUTINE compute_tendencies |
|---|
| 65 | |
|---|
| 66 | |
|---|
| 67 | |
|---|
| 68 | |
|---|
| 69 | |
|---|
