| 1 | ! |
|---|
| 2 | ! $Id $ |
|---|
| 3 | ! |
|---|
| 4 | SUBROUTINE recomp_tend_co2_slope(tendencies_co2_ice_phys,tendencies_co2_ice_phys_ini,vmr_co2_gcm,vmr_co2_pem,ps_GCM_2,global_ave_press_GCM,global_ave_press_new,timelen,ngrid,nslope) |
|---|
| 5 | |
|---|
| 6 | IMPLICIT NONE |
|---|
| 7 | |
|---|
| 8 | !======================================================================= |
|---|
| 9 | ! |
|---|
| 10 | ! Routine that compute the evolution of the tendencie for co2 ice |
|---|
| 11 | ! |
|---|
| 12 | !======================================================================= |
|---|
| 13 | |
|---|
| 14 | ! arguments: |
|---|
| 15 | ! ---------- |
|---|
| 16 | |
|---|
| 17 | ! INPUT |
|---|
| 18 | INTEGER, intent(in) :: timelen,ngrid,nslope |
|---|
| 19 | REAL, INTENT(in) :: vmr_co2_gcm(ngrid,timelen) ! physical point field : Volume mixing ratio of co2 in the first layer |
|---|
| 20 | REAL, INTENT(in) :: vmr_co2_pem(ngrid,timelen) ! physical point field : Volume mixing ratio of co2 in the first layer |
|---|
| 21 | REAL, intent(in) :: ps_GCM_2(ngrid,timelen) ! physical point field : Surface pressure in the GCM |
|---|
| 22 | REAL, intent(in) :: global_ave_press_GCM |
|---|
| 23 | REAL, intent(in) :: global_ave_press_new |
|---|
| 24 | REAL, intent(in) :: tendencies_co2_ice_phys_ini(ngrid,nslope) ! physical point field : Evolution of perenial ice over one year |
|---|
| 25 | |
|---|
| 26 | ! OUTPUT |
|---|
| 27 | REAL, intent(inout) :: tendencies_co2_ice_phys(ngrid,nslope) ! physical point field : Evolution of perenial ice over one year |
|---|
| 28 | |
|---|
| 29 | ! local: |
|---|
| 30 | ! ---- |
|---|
| 31 | |
|---|
| 32 | INTEGER :: i,t,islope |
|---|
| 33 | REAL :: eps, sigma, L, beta, alpha, coef, ave |
|---|
| 34 | |
|---|
| 35 | eps=0.95 |
|---|
| 36 | sigma=5.678E-8 |
|---|
| 37 | L=5.71*10**5 |
|---|
| 38 | beta=3182.48 |
|---|
| 39 | alpha=23.3494 |
|---|
| 40 | |
|---|
| 41 | coef=669*24*3600*eps*sigma/L |
|---|
| 42 | |
|---|
| 43 | ! Evolution of the water ice for each physical point |
|---|
| 44 | do i=1,ngrid |
|---|
| 45 | do islope=1,nslope |
|---|
| 46 | ave=0. |
|---|
| 47 | if(abs(tendencies_co2_ice_phys(i,islope)).gt.1e-4) then |
|---|
| 48 | do t=1,timelen |
|---|
| 49 | ave=ave+(beta/(alpha-log(vmr_co2_gcm(i,t)*ps_GCM_2(i,t)/100.)))**4 & |
|---|
| 50 | -(beta/(alpha-log(vmr_co2_pem(i,t)*ps_GCM_2(i,t)*global_ave_press_GCM/global_ave_press_new/100.)))**4 |
|---|
| 51 | enddo |
|---|
| 52 | endif |
|---|
| 53 | if(ave.lt.1e-4) ave = 0. |
|---|
| 54 | tendencies_co2_ice_phys(i,islope)=tendencies_co2_ice_phys_ini(i,islope)-coef*ave/timelen |
|---|
| 55 | enddo |
|---|
| 56 | enddo |
|---|
| 57 | |
|---|
| 58 | END SUBROUTINE recomp_tend_co2_slope |
|---|
