minor edits

- preparation for new release, because I need an updated doi for a new paper
- minor other edits

Author: nschorgh

CITATION.cff

@@ -5,7 +5,7 @@ authors:
     given-names: Norbert
     orcid: https://orcid.org/0000-0002-5821-4066
 title: "Planetary-Code-Collection: Thermal, Ice Evolution, and Exosphere Models for Planetary Surfaces"
-version: 1.2.2
+version: 1.2.3
 doi: 10.5281/zenodo.594268
 url: "https://github.com/nschorgh/Planetary-Code-Collection/"
-date-released: 2024-01-29
+date-released: 2025-04-05

ChangeLog.txt

@@ -1,6 +1,8 @@
 Log of major changes (goes back to times before git version control)
 ---------------------------------------------------------------------
 
+Apr 2025: Release v1.2.3
+
 May 2024: The thermal model is now also described in a paper by Schorghofer & Khatiwala (2024). The model is 20+ years old, but a new Volterra predictor step has been added, which makes it even better.
 
 29 Jan 2024: Release v1.2.2  doi:10.5281/zenodo.10583245

Lunar/CONTENTS.txt

@@ -1,7 +1,7 @@
 Lunar/
 
-Lunar Ice Pump and Diffusion Models
-===================================
+Lunar Vapor Diffusion Models
+============================
 
 oscidea1.f90: (main program)
 	ice pump calculations for one location; uses boundary-value formulation

Lunar/desorptionrates.f90

@@ -1,14 +1,16 @@
 module desorp
+  ! Adsorption parameters
   implicit none
   real(8), parameter :: kBeV=8.617333262e-5 ! [eV/K]
   !real(8), parameter :: nu=1e16
   real(8), parameter :: thetam=1e19
   real(8), parameter :: Eice=0.529, Ec=0.65, W=0.22
-  !real(8), parameter :: Eice=0.65, Ec=0.65, W=0.0
+  !real(8), parameter :: Eice=0.529, Ec=0.65, W=0.0
 end module desorp
 
 
 subroutine sitedistribution(NS,maxE,nArea)
+  ! Site distribution of desorption energies according to Schorghofer (2023)
   use desorp, only: Ec, W
   implicit none
   integer, intent(IN) :: NS
@@ -24,10 +26,10 @@ subroutine sitedistribution(NS,maxE,nArea)
         nArea(k) = 0.
      else
         nArea(k) = 1./W * exp(-(E(k)-Ec)/W)
-        !print *,k,n(k),(ThetaS(k)<=n(k)*thetam)
+        !print *,k,n(k),(ThetaS(k)<=nArea(k)*thetam)
      end if
   enddo
-  !write(*,'(a,1x,f0.3)') 'sum(n)=',sum(n)*dE
+  !write(*,'(a,1x,f0.3)') 'sum(n)=',sum(nArea)*dE
 end subroutine sitedistribution
 
 
@@ -57,12 +59,12 @@ elemental function desorptionrate(T,theta)
      gamma = 1. / ( 1 + (v-1)*(Ec-Eice+W)*b )
      S2 = nu * thetam * exp(-Eice*b) * exp(-gamma*(Ec-Eice)*b) / (1+gamma*W*b)
 
-     ! option 3 (composite)
+     ! option 3
      gamma = 1./v**2
      S1 = nu*thetam * exp(-Eice*b) * &
           & exp(-gamma*(Ec-Eice)*b) / (1 + gamma* W*b)
-     !S = min(S1,S2)
      S = S1
+     !S = min(S1,S2)  ! composite
   end if
 
   desorptionrate = S
@@ -101,6 +103,7 @@ end subroutine desorptionrateXi1
 
 
 elemental function desorptionrateXi1_multi(T,Ek,nArea,ThetaS)
+  ! spectral version of function desorptionrate
   use desorp, only: thetam, kBeV, Eice
   implicit none
   real(8), intent(IN) :: T ! [K]  
@@ -136,7 +139,7 @@ end function desorptionrateXi1_multi
 
 
 subroutine desorptionrate_inverse_Xi(T,S,NS,maxE,ThetaS)
-  ! inverse of function desorptionrateXi_multi (end-member 2)
+  ! inverse of function desorptionrateXi_multi
   use desorp
   implicit none
   integer, intent(IN) :: NS
@@ -162,8 +165,7 @@ end subroutine desorptionrate_inverse_Xi
 
 
 subroutine distribute_pseudo(NS,nArea,maxE,ThetaS)
-  ! redistribute adsorbate mass into highest energies to mimic endmember with
-  ! maximum surface diffusion
+  ! redistribute adsorbate mass into highest energies to mimic endmember 1
   use desorp, only: thetam
   implicit none
   integer, intent(IN) :: NS

Lunar/params.f90

@@ -21,7 +21,7 @@ module params
 
    !real(8), parameter :: g=0.
    real(8), parameter :: g=9.  ! geothermal gradient [K/m]  g = 0.018/0.002
-   real(8) Tm, Ta
+   real(8) Tm, Ta  ! mean temperature and temperature amplitude [K]
 
    ! Static situations
    parameter (Tm=130., Ta=0.)
@@ -33,13 +33,11 @@ module params
    !parameter (Tm=130., Ta=40.)
    !parameter (Tm=250., Ta=100.)
    !parameter (NZ = 200, Deltaz=0.5e-3, maxtime=(100e3+0.1)*secyear, STEPSPERSOL=24) ! 130+/-40K
-   !parameter (NZ = 100, Deltaz=1e-5, maxtime=(1e3+0.1)*secyear, STEPSPERSOL=48)   ! Sin3m 250+/-100K
-   !parameter (NZ = 10, Deltaz=1e-5, maxtime=(100e3+0.1)*secyear, STEPSPERSOL=24)   ! Sin3m 130+/-40K
-   !parameter (NZ = 10, Deltaz=1e-4, maxtime=(10.+0.1)*secyear, STEPSPERSOL=96)   ! endmember2
+   !parameter (NZ = 100, Deltaz=1e-5, maxtime=(1e3+0.1)*secyear, STEPSPERSOL=48) ! Sin3m 250+/-100K
    !real(8), parameter :: dtsec = lunation/STEPSPERSOL ! time step [s]
 
-   ! adjust upper and lower boundary conditions manually below
-   ! adjust weathering rates manually below
+   ! adjust upper and lower boundary conditions manually
+   ! adjust weathering rates manually
 end module params