[2779] | 1 | ! |
---|
| 2 | ! $Id $ |
---|
| 3 | ! |
---|
| 4 | SUBROUTINE compute_tendencies_slope(tendencies_h2o_ice,min_h2o_ice_Y1,& |
---|
| 5 | min_h2o_ice_Y2,iim_input,jjm_input,ngrid,tendencies_h2o_ice_phys,nslope) |
---|
| 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 ,nslope ! # of grid points along longitude/latitude/ total |
---|
| 22 | REAL, intent(in) , dimension(iim_input+1,jjm_input+1,nslope):: 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,nslope):: 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,nslope) :: tendencies_h2o_ice ! LON x LAT field : difference between the minima = evolution of perenial ice |
---|
| 27 | REAL, intent(out) , dimension(ngrid,nslope) :: 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,islope ! 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 | DO j=1,jjm_input+1 |
---|
| 41 | DO i = 1, iim_input |
---|
| 42 | DO islope = 1, nslope |
---|
| 43 | tendencies_h2o_ice(i,j,islope)=min_h2o_ice_Y2(i,j,islope)-min_h2o_ice_Y1(i,j,islope) |
---|
| 44 | enddo |
---|
| 45 | ENDDO |
---|
| 46 | ENDDO |
---|
| 47 | |
---|
[2794] | 48 | print *, "jjm_input+1", jjm_input+1 |
---|
| 49 | print *, "iim_input+1", iim_input+1 |
---|
| 50 | print *, "nslope+1", nslope+1 |
---|
| 51 | |
---|
[2779] | 52 | ! If the difference is too small; there is no evolution |
---|
| 53 | DO j=1,jjm_input+1 |
---|
| 54 | DO i = 1, iim_input |
---|
| 55 | DO islope = 1, nslope |
---|
[2794] | 56 | ! print *, "tendencies_h2o_ice(i,j,islope)LAAA", tendencies_h2o_ice(i,j,islope) |
---|
[2779] | 57 | if(abs(tendencies_h2o_ice(i,j,islope)).LT.1.0E-10) then |
---|
| 58 | tendencies_h2o_ice(i,j,islope)=0. |
---|
| 59 | endif |
---|
[2794] | 60 | ! print *, "tendencies_h2o_ice(i,j,islope)HERE", tendencies_h2o_ice(i,j,islope) |
---|
[2779] | 61 | enddo |
---|
| 62 | ENDDO |
---|
| 63 | ENDDO |
---|
| 64 | |
---|
| 65 | |
---|
| 66 | ! We reorganise the difference on the physical grid |
---|
| 67 | tendencies_h2o_ice_phys(1,:)=tendencies_h2o_ice(1,1,:) |
---|
| 68 | |
---|
| 69 | ig0 = 2 |
---|
| 70 | DO j=2,jjm_input |
---|
| 71 | DO i = 1, iim_input |
---|
| 72 | tendencies_h2o_ice_phys(ig0,:) =tendencies_h2o_ice(i,j,:) |
---|
| 73 | ig0= ig0 + 1 |
---|
| 74 | ENDDO |
---|
| 75 | ENDDO |
---|
| 76 | |
---|
| 77 | tendencies_h2o_ice_phys(ig0,:) = tendencies_h2o_ice(1,jjm_input+1,:) |
---|
| 78 | |
---|
[2794] | 79 | ! print *, "tendencies_h2o_ice_physze", tendencies_h2o_ice_phys(:,:) |
---|
[2779] | 80 | |
---|
[2794] | 81 | |
---|
[2779] | 82 | END SUBROUTINE compute_tendencies_slope |
---|
| 83 | |
---|
| 84 | |
---|
| 85 | |
---|
| 86 | |
---|
| 87 | |
---|