source: trunk/LMDZ.COMMON/libf/evolution/adsorption_mod.F90 @ 3241

Last change on this file since 3241 was 3206, checked in by jbclement, 17 months ago

PEM:
Cleanings of unused variables/arguments and bad type conversions.
JBC
File size: 18.0 KB
Line 
1module adsorption_mod
2
3  implicit none
4
5  LOGICAL adsorption_pem ! True by default, to compute adsorption/desorption. Read in  pem.def
6  real, save, allocatable :: co2_adsorbded_phys(:,:,:)  ! co2 that is in the regolith [kg/m^2]
7  real, save, allocatable :: h2o_adsorbded_phys(:,:,:)  ! h2o that is in the regolith [kg/m^2]
8
9  contains
10
11!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
12!!!
13!!! Purpose: Compute CO2 and H2O adsorption, following the methods from Zent & Quinn 1995, Jackosky et al., 1997
14!!!
15!!! Author: LL, 01/2023
16!!!
17!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
18  subroutine ini_adsorption_h_PEM(ngrid,nslope,nsoilmx_PEM)
19
20  implicit none
21  integer,intent(in) :: ngrid ! number of atmospheric columns
22  integer,intent(in) :: nslope ! number of slope within a mesh
23  integer,intent(in) :: nsoilmx_PEM ! number of soil layer in the PEM
24    allocate(co2_adsorbded_phys(ngrid,nsoilmx_PEM,nslope))
25    allocate(h2o_adsorbded_phys(ngrid,nsoilmx_PEM,nslope))
26  end subroutine ini_adsorption_h_PEM
27
28!!! -----------------------------------------------
29
30  subroutine end_adsorption_h_PEM
31
32  implicit none
33    if (allocated(co2_adsorbded_phys)) deallocate(co2_adsorbded_phys)
34    if (allocated(h2o_adsorbded_phys)) deallocate(h2o_adsorbded_phys)
35  end subroutine end_adsorption_h_PEM
36
37!!! -----------------------------------------------
38
39  subroutine regolith_adsorption(ngrid,nslope,nsoil_PEM,timelen,tend_h2oglaciers,tend_co2glaciers,waterice,co2ice,tsoil_PEM,TI_PEM,ps,q_co2,q_h2o, &
40                                m_h2o_completesoil,delta_mh2oreg, m_co2_completesoil,delta_mco2reg)
41
42! inputs
43 INTEGER,INTENT(IN) :: ngrid, nslope, nsoil_PEM,timelen ! size dimension: physics x subslope x soil x timeseries
44 REAL,INTENT(IN) :: tend_h2oglaciers(ngrid,nslope),tend_co2glaciers(ngrid,nslope) !tendancies on the glaciers [1]
45 REAL,INTENT(IN) :: waterice(ngrid,nslope)              ! water ice at the surface [kg/m^2]
46 REAL,INTENT(IN) :: co2ice(ngrid,nslope)                ! co2 ice at the surface [kg/m^2]
47 REAL,INTENT(IN) :: TI_PEM(ngrid,nsoil_PEM,nslope)      ! Soil Thermal inertia (J/K/^2/s^1/2)
48 REAL,INTENT(IN) :: tsoil_PEM(ngrid,nsoil_PEM,nslope)   ! Soil temperature (K)
49 REAL,INTENT(IN) :: ps(ngrid,timelen)                   ! Average surface pressure [Pa]
50 REAL,INTENT(IN) :: q_co2(ngrid,timelen)                ! Mass mixing ratio of co2 in the first layer (kg/kg)
51 REAL,INTENT(IN) :: q_h2o(ngrid,timelen)                ! Mass mixing ratio of H2o in the first layer (kg/kg)
52
53! outputs
54 REAL,INTENT(INOUT) :: m_h2o_completesoil(ngrid,nsoil_PEM,nslope) ! Density of h2o adsorbed (kg/m^3)(ngrid,nsoil_PEM,nslope)     
55 REAL,INTENT(OUT) :: delta_mh2oreg(ngrid)                         ! Difference density of h2o adsorbed (kg/m^3)
56
57 REAL,INTENT(INOUT) :: m_co2_completesoil(ngrid,nsoil_PEM,nslope)   ! Density of co2 adsorbed (kg/m^3)
58 REAL,INTENT(OUT) :: delta_mco2reg(ngrid)                            ! Difference density of co2 adsorbed (kg/m^3)
59 
60! local variables
61 REAL :: theta_h2o_adsorbded(ngrid,nsoil_PEM,nslope) ! Fraction of the pores occupied by H2O molecules
62! -------------
63
64
65! Compute H2O adsorption, then CO2 adsorption
66
67  call regolith_h2oadsorption(ngrid,nslope,nsoil_PEM,timelen,tend_h2oglaciers,tend_co2glaciers,waterice,co2ice,ps,q_co2,q_h2o,tsoil_PEM,TI_PEM, &
68                                   theta_h2o_adsorbded,m_h2o_completesoil,delta_mh2oreg)
69
70
71  call regolith_co2adsorption(ngrid,nslope,nsoil_PEM,timelen,tend_h2oglaciers,tend_co2glaciers,waterice,co2ice,ps,q_co2,q_h2o, &
72                                                          tsoil_PEM,TI_PEM,m_co2_completesoil,delta_mco2reg)
73
74   RETURN
75  end subroutine
76
77!------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
78
79  subroutine regolith_h2oadsorption(ngrid,nslope,nsoil_PEM,timelen,tend_h2oglaciers,tend_co2glaciers,waterice,co2ice,ps,q_co2,q_h2o,tsoil_PEM,TI_PEM, &
80                                   theta_h2o_adsorbded,m_h2o_completesoil,delta_mreg)
81
82!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
83!!!
84!!! Purpose: Compute H2O adsorption, following the methods from Jackosky et al., 1997
85!!!
86!!! Author: LL, 01/2023
87!!!
88!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
89
90use comsoil_h_PEM,         only: layer_PEM, index_breccia
91use comslope_mod,          only: subslope_dist, def_slope_mean
92use vertical_layers_mod,   only: ap,bp
93use constants_marspem_mod, only: alpha_clap_h2o, beta_clap_h2o, m_h2o, m_co2,m_noco2, rho_regolith
94
95#ifndef CPP_STD
96    use comcstfi_h,   only: pi
97#else
98    use comcstfi_mod, only: pi
99#endif
100
101 implicit none
102! inputs
103 INTEGER,INTENT(IN) :: ngrid, nslope, nsoil_PEM,timelen ! Size dimension
104 REAL,INTENT(IN) :: tend_h2oglaciers(ngrid,nslope),tend_co2glaciers(ngrid,nslope) ! Tendencies on the glaciers ()
105 REAL,INTENT(IN) :: waterice(ngrid,nslope)              ! Water ice at the surface [kg/m^2]
106 REAL,INTENT(IN) :: co2ice(ngrid,nslope)                ! CO2 ice at the surface [kg/m^2]
107 REAL,INTENT(IN) :: ps(ngrid,timelen)                   ! Surface pressure (Pa)         
108 REAL,INTENT(IN) :: q_co2(ngrid,timelen)                ! Mass mixing ratio of co2 in the first layer (kg/kg)
109 REAL,INTENT(IN) :: q_h2o(ngrid,timelen)                ! Mass mixing ratio of H2o in the first layer (kg/kg)
110 REAL,INTENT(IN) :: TI_PEM(ngrid,nsoil_PEM,nslope)      ! Soil Thermal inertia (J/K/^2/s^1/2)
111 REAL,INTENT(IN) :: tsoil_PEM(ngrid,nsoil_PEM,nslope)   ! Soil temperature (K)
112
113! outputs
114 REAL,INTENT(INOUT) :: m_h2o_completesoil(ngrid,nsoil_PEM,nslope) ! Density of h2o adsorbed (kg/m^3)(ngrid,nsoil_PEM,nslope)     
115 REAL,INTENT(OUT) :: theta_h2o_adsorbded(ngrid,nsoil_PEM,nslope)  ! Fraction of the pores occupied by H2O molecules
116 REAL,INTENT(OUT) :: delta_mreg(ngrid)                            ! Difference density of h2o adsorbed (kg/m^3)
117
118! constants
119 REAL :: Ko =  1.57e-8            ! Jackosky et al. 1997
120 REAL :: e = 2573.9               ! Jackosky et al. 1997
121 REAL :: mu = 0.48                ! Jackosky et al. 1997
122 real :: m_theta = 2.84e-7        ! Mass of h2o per m^2 absorbed Jackosky et al. 1997
123! real :: as = 18.9e3             ! Specific area, Buhler & Piqueux 2021
124 real :: as = 9.48e4              ! Specific area, Zent
125 real ::  inertie_thresold = 800. ! TI > 800 means cementation
126
127 ! local variables
128 REAL :: deltam_reg_complete(ngrid,index_breccia,nslope) ! Difference in the mass per slope and soil layer (kg/m^3)
129 real :: K                        ! Used to compute theta
130 integer ig, iloop, islope, it    ! For loops
131 INTEGER :: ispermanent_co2glaciers(ngrid,nslope)      ! Check if the co2 glacier is permanent
132 INTEGER :: ispermanent_h2oglaciers(ngrid,nslope)      ! Check if the h2o glacier is permanent
133 REAL :: deltam_reg_slope(ngrid,nslope)                ! Difference density of h2o adsorbed per slope (kg/m^3)
134 REAL :: dm_h2o_regolith_slope(ngrid,nsoil_PEM,nslope) ! Elementary h2o mass adsorded per mesh per slope
135 real :: A,B                                           ! Used to compute the mean mass above the surface
136 real :: p_sat                                         ! Saturated vapor pressure of ice
137 real,allocatable :: mass_mean(:,:)                    ! Mean mass above the surface
138 real,allocatable :: zplev_mean(:,:)                   ! Pressure above the surface
139 real,allocatable :: pvapor(:,:)                       ! Partial pressure above the surface
140 real, allocatable :: pvapor_avg(:)                    ! Yearly averaged
141
142! 0. Some initializations
143
144
145     allocate(mass_mean(ngrid,timelen))
146     allocate(zplev_mean(ngrid,timelen))
147     allocate(pvapor(ngrid,timelen))
148     allocate(pvapor_avg(ngrid))
149     A =(1/m_co2 - 1/m_noco2)
150     B=1/m_noco2
151     theta_h2o_adsorbded(:,:,:) = 0.
152     dm_h2o_regolith_slope(:,:,:) = 0.
153
154#ifndef CPP_STD
155
156!0.1 Look at perennial ice
157  do ig = 1,ngrid
158    do islope = 1,nslope
159     if((abs(tend_h2oglaciers(ig,islope)).gt.1e-5).and.(abs(waterice(ig,islope)).gt.0)) then
160        ispermanent_h2oglaciers(ig,islope) = 1
161     else
162        ispermanent_h2oglaciers(ig,islope) = 0
163     endif
164
165     if((abs(tend_co2glaciers(ig,islope)).gt.1e-5).and.(abs(co2ice(ig,islope)).gt.0)) then
166        ispermanent_co2glaciers(ig,islope) = 1
167     else
168        ispermanent_co2glaciers(ig,islope) = 0
169     endif
170    enddo
171   enddo
172
173!   0.2 Compute the partial pressure of vapor
174!a. the molecular mass into the column
175     do ig = 1,ngrid
176       mass_mean(ig,:) = 1/(A*q_co2(ig,:) +B)
177     enddo
178
179
180! b. pressure level
181     do it = 1,timelen
182       do ig = 1,ngrid
183         zplev_mean(ig,it) = ap(1) + bp(1)*ps(ig,it)
184       enddo
185     enddo
186! c. Vapor pressure
187     pvapor(:,:) = mass_mean(:,:)/m_h2o*q_h2o(:,:)*zplev_mean(:,:)
188     pvapor_avg(:) = sum(pvapor(:,:),2)/timelen
189#endif
190     deallocate(pvapor)
191     deallocate(zplev_mean)
192     deallocate(mass_mean)
193#ifndef CPP_STD
194
195! 1. we compute the mass of H2O adsorded in each layer of the meshes 
196
197 do ig = 1,ngrid
198  do islope = 1,nslope
199    do iloop = 1,index_breccia
200      K = Ko*exp(e/tsoil_PEM(ig,iloop,islope))
201      if(TI_PEM(ig,iloop,islope).lt.inertie_thresold)  then
202        theta_h2o_adsorbded(ig,iloop,islope) = (K*pvapor_avg(ig)/(1+K*pvapor_avg(ig)))**mu
203      else
204        p_sat =exp(beta_clap_h2o/tsoil_PEM(ig,iloop,islope) +alpha_clap_h2o) ! we assume fixed temperature in the ice ... not really good but ...
205        theta_h2o_adsorbded(ig,iloop,islope) = (K*p_sat/(1+K*p_sat))**mu
206      endif
207        dm_h2o_regolith_slope(ig,iloop,islope) = as*theta_h2o_adsorbded(ig,iloop,islope)*m_theta*rho_regolith
208   enddo
209  enddo
210 enddo
211
212! 2. Check the exchange between the atmosphere and the regolith
213
214  do ig = 1,ngrid
215   delta_mreg(ig) = 0.
216   do islope = 1,nslope
217    deltam_reg_slope(ig,islope) = 0.
218    do iloop = 1,index_breccia
219       if((TI_PEM(ig,iloop,islope).lt.inertie_thresold).and.(ispermanent_h2oglaciers(ig,islope).eq.0).and.(ispermanent_co2glaciers(ig,islope).eq.0)) then
220             deltam_reg_complete(ig,iloop,islope) = (dm_h2o_regolith_slope(ig,iloop,islope) - m_h2o_completesoil(ig,iloop,islope)) &
221                                                   *(layer_PEM(iloop+1) - layer_PEM(iloop))
222       else ! NO EXCHANGE AS ICE BLOCK THE DYNAMIC!
223             deltam_reg_complete(ig,iloop,islope) = 0.
224       endif
225       deltam_reg_slope(ig,islope) = deltam_reg_slope(ig,islope) + deltam_reg_complete(ig,iloop,islope)
226    enddo
227   delta_mreg(ig) = delta_mreg(ig) + deltam_reg_slope(ig,islope)*subslope_dist(ig,islope)/cos(pi*def_slope_mean(islope)/180.)
228   enddo
229  enddo
230   m_h2o_completesoil(:,:,:) = dm_h2o_regolith_slope(:,:,:)
231
232 RETURN
233#endif
234  end subroutine
235
236!------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
237!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
238!!!
239!!! Purpose: Compute CO2  following the methods from Zent & Quinn 1995
240!!!
241!!! Author: LL, 01/2023
242!!!
243!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
244
245  SUBROUTINE regolith_co2adsorption(ngrid,nslope,nsoil_PEM,timelen,tend_h2oglaciers,tend_co2glaciers,waterice,co2ice,ps,q_co2,q_h2o,&
246                                   tsoil_PEM,TI_PEM,m_co2_completesoil,delta_mreg)
247
248use comsoil_h_PEM,         only: layer_PEM, index_breccia, index_breccia
249use comslope_mod,          only: subslope_dist, def_slope_mean
250use vertical_layers_mod,   only: ap, bp
251use constants_marspem_mod, only: m_co2, m_noco2, rho_regolith
252
253#ifndef CPP_STD
254    use comcstfi_h,   only: pi
255#else
256    use comcstfi_mod, only: pi
257#endif
258
259      IMPLICIT NONE
260! Inputs: 
261 INTEGER,INTENT(IN) :: ngrid, nslope, nsoil_PEM,timelen             ! Size dimension
262 REAL,INTENT(IN) :: ps(ngrid,timelen)                               ! Average surface pressure [Pa]
263 REAL,INTENT(IN) :: tsoil_PEM(ngrid,nsoil_PEM,nslope)               ! Mean Soil Temperature [K]
264 REAL,INTENT(IN) :: TI_PEM(ngrid,nsoil_PEM,nslope)                  ! Mean Thermal Inertia [USI]
265 REAL,INTENT(IN) :: tend_h2oglaciers(ngrid,nslope),tend_co2glaciers(ngrid,nslope) ! Tendencies on the glaciers ()
266 REAL,INTENT(IN) :: q_co2(ngrid,timelen),q_h2o(ngrid,timelen)       ! Mass mixing ratio of co2 and h2o in the first layer (kg/kg)
267 REAL,INTENT(IN) :: waterice(ngrid,nslope)                          ! Water ice at the surface [kg/m^2]
268 REAL,INTENT(IN) :: co2ice(ngrid,nslope)                            ! CO2 ice at the surface [kg/m^2]
269
270! Outputs:
271 REAL,INTENT(INOUT) :: m_co2_completesoil(ngrid,nsoil_PEM,nslope)   ! Density of co2 adsorbed (kg/m^3)
272 REAL,INTENT(OUT) :: delta_mreg(ngrid)                              ! Difference density of co2 adsorbed (kg/m^3)
273 
274! Constants:
275
276 REAL :: alpha = 7.512e-6         ! Zent & Quinn 1995
277 REAL :: beta =  -1541.5          ! Zent & Quinn 1995
278 REAL ::  inertie_thresold = 800. ! TI > 800 means cementation
279 real :: m_theta = 4.27e-7        ! Mass of co2 per m^2 absorbed
280! real :: as = 18.9e3             ! Specific area, Buhler & Piqueux 2021
281 real :: as = 9.48e4              ! Same as previous but from zent
282! Local         
283 real :: A,B                                             ! Used to compute the mean mass above the surface
284 INTEGER :: ig,islope,iloop,it                           ! For loops
285 REAL :: dm_co2_regolith_slope(ngrid,nsoil_PEM,nslope)   ! Elementary mass adsorded per mesh per slope
286 INTEGER :: ispermanent_co2glaciers(ngrid,nslope)        ! Check if the co2 glacier is permanent
287 INTEGER :: ispermanent_h2oglaciers(ngrid,nslope)        ! Check if the h2o glacier is permanent
288 REAL :: deltam_reg_complete(ngrid,index_breccia,nslope) ! Difference in the mass per slope and soil layer (kg/m^3)
289 REAL :: deltam_reg_slope(ngrid,nslope)                  ! Difference in the mass per slope  (kg/m^3)
290 REAL :: m_h2o_adsorbed(ngrid,nsoil_PEM,nslope)          ! Density of CO2 adsorbed (kg/m^3)
291 REAL :: theta_h2o_adsorbed(ngrid,nsoil_PEM,nslope)      ! Fraction of the pores occupied by H2O molecules
292 REAL :: delta_mh2o(ngrid)                               ! Difference density of h2o adsorbed (kg/m^3)
293!timelen array are allocated because heavy ...
294 real,allocatable :: mass_mean(:,:)                      ! Mean mass above the surface
295 real,allocatable :: zplev_mean(:,:)                     ! Pressure above the surface
296 real,allocatable :: pco2(:,:)                           ! Partial pressure above the surface
297 real, allocatable :: pco2_avg(:)                        ! Yearly averaged
298
299! 0. Some initializations
300
301     allocate(mass_mean(ngrid,timelen))
302     allocate(zplev_mean(ngrid,timelen))
303     allocate(pco2(ngrid,timelen))
304     allocate(pco2_avg(ngrid))
305
306
307 
308      m_h2o_adsorbed(:,:,:) = 0.
309      A =(1/m_co2 - 1/m_noco2)
310      B=1/m_noco2
311
312     dm_co2_regolith_slope(:,:,:) = 0
313     delta_mreg(:) = 0.
314
315#ifndef CPP_STD
316
317!0.1 Look at perennial ice
318  do ig = 1,ngrid
319    do islope = 1,nslope
320     if((abs(tend_h2oglaciers(ig,islope)).gt.1e-5).and.(abs(waterice(ig,islope)).gt.0)) then
321        ispermanent_h2oglaciers(ig,islope) = 1
322     else
323        ispermanent_h2oglaciers(ig,islope) = 0
324     endif
325
326     if((abs(tend_co2glaciers(ig,islope)).gt.1e-5).and.(abs(co2ice(ig,islope)).gt.0)) then
327        ispermanent_co2glaciers(ig,islope) = 1
328     else
329        ispermanent_co2glaciers(ig,islope) = 0
330     endif
331    enddo
332   enddo
333
334!   0.2  Compute the partial pressure of CO2
335!a. the molecular mass into the column
336     do ig = 1,ngrid
337       mass_mean(ig,:) = 1./(A*q_co2(ig,:) +B)
338     enddo
339
340! b. pressure level
341     do it = 1,timelen
342       do ig = 1,ngrid
343         zplev_mean(ig,it) = ap(1) + bp(1)*ps(ig,it)
344       enddo
345     enddo
346
347! c. Vapor pressure
348     pco2(:,:) = mass_mean(:,:)/m_co2*q_co2(:,:)*zplev_mean(:,:)
349     pco2_avg(:) = sum(pco2(:,:),2)/timelen
350
351     deallocate(zplev_mean)
352     deallocate(mass_mean)
353     deallocate(pco2)
354
355
356! 1. Compute the fraction of the pores occupied by H2O
357
358 call regolith_h2oadsorption(ngrid,nslope,nsoil_PEM,timelen,tend_h2oglaciers,tend_co2glaciers,waterice,co2ice,ps,q_co2,q_h2o,tsoil_PEM,TI_PEM, &
359                                   theta_h2o_adsorbed, m_h2o_adsorbed,delta_mh2o)
360
361! 2.  we compute the mass of co2 adsorded in each layer of the meshes 
362
363 do ig = 1,ngrid
364  do islope = 1,nslope
365    do iloop = 1,index_breccia
366     if((TI_PEM(ig,iloop,islope).lt.inertie_thresold).and.(ispermanent_h2oglaciers(ig,islope).eq.0).and.(ispermanent_co2glaciers(ig,islope).eq.0)) then
367     dm_co2_regolith_slope(ig,iloop,islope) = as*rho_regolith*m_theta*(1-theta_h2o_adsorbed(ig,iloop,islope))*alpha*pco2_avg(ig)/ &
368                                             (alpha*pco2_avg(ig)+sqrt(tsoil_PEM(ig,iloop,islope))*exp(beta/tsoil_PEM(ig,iloop,islope)))
369     else
370        if(abs(m_co2_completesoil(ig,iloop,islope)).lt.(1e-10)) then !!! we are at first call
371          dm_co2_regolith_slope(ig,iloop,islope) = as*rho_regolith*m_theta*(1-theta_h2o_adsorbed(ig,iloop,islope))*alpha*pco2_avg(ig) &
372                                                 /(alpha*pco2_avg(ig)+sqrt(tsoil_PEM(ig,iloop,islope))*exp(beta/tsoil_PEM(ig,iloop,islope)))
373        else ! no change: permanent ice stick the atoms of CO2
374          dm_co2_regolith_slope(ig,iloop,islope) = m_co2_completesoil(ig,iloop,islope)
375        endif
376     endif
377    enddo
378 enddo
379enddo
380
381! 3. Check the exchange between the atmosphere and the regolith
382
383  do ig = 1,ngrid
384   delta_mreg(ig) = 0.
385   do islope = 1,nslope
386    deltam_reg_slope(ig,islope) = 0.
387    do iloop = 1,index_breccia
388       if((TI_PEM(ig,iloop,islope).lt.inertie_thresold).and.(ispermanent_h2oglaciers(ig,islope).eq.0).and.(ispermanent_co2glaciers(ig,islope).eq.0)) then
389             deltam_reg_complete(ig,iloop,islope) = (dm_co2_regolith_slope(ig,iloop,islope) - m_co2_completesoil(ig,iloop,islope)) &
390                                                   *(layer_PEM(iloop+1) - layer_PEM(iloop))
391       else ! NO EXCHANGE AS ICE BLOCK THE DYNAMIC!
392             deltam_reg_complete(ig,iloop,islope) = 0.
393       endif
394       deltam_reg_slope(ig,islope) = deltam_reg_slope(ig,islope) + deltam_reg_complete(ig,iloop,islope)
395    enddo
396   delta_mreg(ig) = delta_mreg(ig) + deltam_reg_slope(ig,islope)*subslope_dist(ig,islope)/cos(pi*def_slope_mean(islope)/180.)
397   enddo
398  enddo
399  m_co2_completesoil(:,:,:) = dm_co2_regolith_slope(:,:,:)
400
401!=======================================================================
402      RETURN
403#endif
404      END
405
406
407   end module
Note: See TracBrowser for help on using the repository browser.