make_profile_rc

src/CMakeLists.txt

@@ -24,6 +24,7 @@ set(RADIATE_SOURCES
 
 set(ADIABAT_SOURCES
   adiabat/clima_adiabat_general.f90
+  adiabat/clima_adiabat_rc.f90
   adiabat/clima_adiabat.f90
 )
 

src/adiabat/clima_adiabat.f90

@@ -14,7 +14,7 @@ module clima_adiabat
 
     ! settings and free parameters
     integer :: nz
-    real(dp) :: P_top = 1.0e-2_dp !! (dynes/cm2)
+    real(dp) :: P_top = 1.0_dp !! (dynes/cm2)
     real(dp) :: T_trop = 180.0_dp !! (T)
     real(dp), allocatable :: RH(:) !! relative humidity (ng)
     !> If .true., then any function that calls `make_column` will
@@ -105,6 +105,8 @@ module clima_adiabat
     procedure :: surface_temperature => AdiabatClimate_surface_temperature
     procedure :: surface_temperature_column => AdiabatClimate_surface_temperature_column
     procedure :: surface_temperature_bg_gas => AdiabatClimate_surface_temperature_bg_gas
+    ! Test
+    procedure :: make_profile_rc => AdiabatClimate_make_profile_rc
     ! Utilities
     procedure :: set_ocean_solubility_fcn => AdiabatClimate_set_ocean_solubility_fcn
     procedure :: to_regular_grid => AdiabatClimate_to_regular_grid
@@ -765,6 +767,64 @@ subroutine fcn(n_, x_, fvec_, iflag_)
     end subroutine
   end function
 
+  subroutine AdiabatClimate_make_profile_rc(self, P_i_surf, T_surf, T, err)
+    use clima_const, only: k_boltz, N_avo
+    use clima_adiabat_rc, only: make_profile_rc
+    class(AdiabatClimate), intent(inout) :: self
+    real(dp), intent(in) :: P_i_surf(:) !! dynes/cm^2
+    real(dp), intent(in) :: T_surf, T(:)
+    character(:), allocatable, intent(out) :: err
+
+    real(dp), allocatable :: P_e(:), z_e(:), f_i_e(:,:), lapse_rate_e(:)
+    real(dp), allocatable :: density(:)
+    integer :: i, j
+
+    allocate(P_e(2*self%nz+1),z_e(2*self%nz+1),f_i_e(2*self%nz+1,self%sp%ng),lapse_rate_e(2*self%nz+1))
+    allocate(density(self%nz))
+
+    if (size(P_i_surf) /= self%sp%ng) then
+      err = "P_i_surf has the wrong dimension"
+      return
+    endif
+    if (size(T) /= self%nz) then
+      err = "T has the wrong dimension"
+    endif
+
+    self%T_surf = T_surf
+    self%T = T
+    call make_profile_rc(self%T_surf, self%T, P_i_surf, &
+                         self%sp, self%nz, self%planet_mass, &
+                         self%planet_radius, self%P_top, self%RH, &
+                         self%rtol, self%atol, &
+                         self%ocean_fcns, self%ocean_args_p, &
+                         P_e, z_e, f_i_e, lapse_rate_e, &
+                         self%N_surface, self%N_ocean, &
+                         err)
+    if (allocated(err)) return
+
+    self%T_surf = self%T_surf
+    self%P_surf = P_e(1)
+    do i = 1,self%nz
+      self%P(i) = P_e(2*i)
+      self%z(i) = z_e(2*i)
+      self%dz(i) = z_e(2*i+1) - z_e(2*i-1)
+      do j =1,self%sp%ng
+        self%f_i(i,j) = f_i_e(2*i,j)
+      enddo
+    enddo
+
+    density = self%P/(k_boltz*self%T)
+    do j =1,self%sp%ng
+      self%densities(:,j) = self%f_i(:,j)*density(:)
+    enddo
+
+    do i = 1,self%sp%ng
+      ! mol/cm^2 in atmosphere
+      self%N_atmos(i) = sum(density*self%f_i(:,i)*self%dz)/N_avo
+    enddo
+
+  end subroutine
+
   !> Sets a function for describing how gases dissolve in a liquid ocean.
   subroutine AdiabatClimate_set_ocean_solubility_fcn(self, species, fcn, err)
     class(AdiabatClimate), intent(inout) :: self

src/adiabat/clima_adiabat_general.f90

@@ -9,6 +9,7 @@ module clima_adiabat_general
 
   public :: make_profile, make_column
   public :: OceanFunction
+  public :: DrySpeciesType, CondensingSpeciesType
 
   type :: OceanFunction
     procedure(ocean_solubility_fcn), nopass, pointer :: fcn => null()