| 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) :: q_co2_GCM(ngrid) ! physical point field : Density of co2 in the first layer |
|---|
| 23 | ! REAL, intent(in) :: ps_GCM(ngrid) ! physical point field : Density of co2 in the first layer |
|---|
| 24 | REAL, intent(in) :: global_ave_press_GCM |
|---|
| 25 | REAL, intent(in) :: global_ave_press_new |
|---|
| 26 | REAL, intent(in) :: tendencies_co2_ice_phys_ini(ngrid,nslope) ! physical point field : Evolution of perenial ice over one year |
|---|
| 27 | |
|---|
| 28 | |
|---|
| 29 | ! OUTPUT |
|---|
| 30 | REAL, intent(inout) :: tendencies_co2_ice_phys(ngrid,nslope) ! physical point field : Evolution of perenial ice over one year |
|---|
| 31 | |
|---|
| 32 | |
|---|
| 33 | ! local: |
|---|
| 34 | ! ---- |
|---|
| 35 | |
|---|
| 36 | INTEGER :: i,t,islope |
|---|
| 37 | REAL :: eps, sigma, L, beta, alpha, coef, ave |
|---|
| 38 | |
|---|
| 39 | eps=0.95 |
|---|
| 40 | sigma=5.678E-8 |
|---|
| 41 | L=5.71*10**5 |
|---|
| 42 | beta=3182.48 |
|---|
| 43 | alpha=23.3494 |
|---|
| 44 | |
|---|
| 45 | coef=669*24*3600*eps*sigma/L |
|---|
| 46 | |
|---|
| 47 | print *, "coef", coef |
|---|
| 48 | print *, "global_ave_press_GCM", global_ave_press_GCM |
|---|
| 49 | print *, "global_ave_press_new", global_ave_press_new |
|---|
| 50 | |
|---|
| 51 | ! Evolution of the water ice for each physical point |
|---|
| 52 | do i=1,ngrid |
|---|
| 53 | do islope=1,nslope |
|---|
| 54 | ave=0. |
|---|
| 55 | do t=1,timelen |
|---|
| 56 | |
|---|
| 57 | ! write(*,*)'i,t=',i,t,islope, alpha,beta,ave,vmr_co2_gcm(i,t),vmr_co2_pem(i,t),ps_GCM_2(i,t),global_ave_press_GCM,global_ave_press_new |
|---|
| 58 | ave=ave+(beta/(alpha-log(vmr_co2_gcm(i,t)*ps_GCM_2(i,t)/100)))**4 & |
|---|
| 59 | -(beta/(alpha-log(vmr_co2_pem(i,t)*ps_GCM_2(i,t)*global_ave_press_GCM/global_ave_press_new/100)))**4 |
|---|
| 60 | enddo |
|---|
| 61 | ! print *, "i", i |
|---|
| 62 | ! print *, "tendencies_co2_ice_phys_ini bef", tendencies_co2_ice_phys_ini(i,islope) |
|---|
| 63 | ! tendencies_co2_ice_phys(i,islope)=tendencies_co2_ice_phys_ini(i,islope)+coef*ave/timelen |
|---|
| 64 | tendencies_co2_ice_phys(i,islope)=tendencies_co2_ice_phys_ini(i,islope)-coef*ave/timelen |
|---|
| 65 | |
|---|
| 66 | ! print *, "tendencies after", tendencies_co2_ice_phys(i,islope) |
|---|
| 67 | ! print *, "ave", ave |
|---|
| 68 | ! print *, "timelen", timelen |
|---|
| 69 | ! print *, "vmr_co2_pem(i,t)*ps_GCM_2(i,t)", vmr_co2_pem(i,t)*ps_GCM_2(i,t) |
|---|
| 70 | ! print *, "-------------------" |
|---|
| 71 | enddo |
|---|
| 72 | enddo |
|---|
| 73 | |
|---|
| 74 | |
|---|
| 75 | END SUBROUTINE recomp_tend_co2_slope |
|---|
