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