Condensable Atmosphere

From krc
(Difference between revisions)
Jump to: navigation, search
(Created page with " define KRC_Cond_Gas(){ if($ARGC == 0){ printf (" Uses the adequate Clausius Clapeyron Coefficients to predict the frost point \n") printf (" Feeds KRC with SatPrA, ...")
 
(Titan)
 
(16 intermediate revisions by one user not shown)
Line 1: Line 1:
 +
The coefficients for the saturation temperature relation for the condensing gas are now input parameters.
 +
KRC uses the Clausius-Clapeyron relation:
  
define KRC_Cond_Gas(){
 
  
if($ARGC == 0){
+
::ln P = a - b/T
    printf (" Uses the adequate Clausius Clapeyron Coefficients to predict the frost point \n")
+
    printf (" Feeds KRC with SatPrA, SatPrB, FANON, AMW \n")
+
    printf (" $1: Body \n")
+
    printf (" Currently Supported: \"Mars\" or \"Pluto\" or \"Titan\" \n")
+
}
+
  
Body  = $1
 
 
if(Body == "Mars"){             
 
  SatPrB = 3182.48
 
  SatPrA = 27.9546
 
  AMW    = 43.54                                                                      #Sylvain's own calculation
 
    FANON  = 0.040                                                                        #Paul R. Mahaffy, Science 2013 (SAM)
 
KPREF  = 1
 
}
 
  
if(Body == "Pluto"){
+
where P is pressure in Pascal and T is temperature in Kelvin.  
  SatPrB = 9.2338
+
  SatPrA = -724.9720
+
  AMW  = 27.
+
  FANON = 0.1
+
KPREF = 0
+
}
+
  
if(Body == "Titan"){
+
This should be useful for Mars, Titan, Pluto etc.
  SatPrB = 9.2338
+
  SatPrA = -724.9720
+
AMW  = 28.8
+
  FANON = 0.016
+
KPREF = 0
+
}
+
  
if(Body != "Mars" && Body != "Titan" && Body != "Pluto") {
+
:SatPrA correspond to the Clausius-Clapeyron coefficient 'a'.
  printf("Only Mars, Pluto, and Titan currently support Condensable Gas defaults\n")
+
 
return(-999.)
+
:SatPrB corresponds to the Clausius-Clapeyron coeffcient 'b'.
} else {
+
 
out        = {}
+
:LVFT = "T"
out.SatPrA = SatPrA
+
 
out.SatPrB = SatPrB
+
Values for CO2 are a = 27.9546 and b = 3182.48 (Provided by HHK, most likely from the Mars book).
out.AMW    = AMW
+
Values for N2 are a = 9.2338 and b = -724.9720 (Derived by SP using this table: [[File:Clausius_Clapeyron.xlsx]]
out.FANON  = FANON
+
 
out.KPREF  = KPREF
+
The proper molecular weight should be input AMW. The Mass-fraction of mean atmosphere that is non-condensing should also be specified with FANON.
return(out)
+
 
}
+
===Mars===
}
+
 
 +
SatPrB = 3182.48
 +
SatPrA = 27.9546
 +
AMW    = 43.54            #SP's own calculation
 +
FANON  = 0.040            #Paul R. Mahaffy, Science 2013 (SAM)
 +
KPREF  = 1
 +
 
 +
[[Image:TFROST_1.png|800px]]
 +
 
 +
out_1 = krc(lat=72.,SatPrA=27.9546,SatPrB=3182.48,AMW=43.54,FANON=0.040,KPREF=1,LVFT="T" )
 +
out_2 = krc(lat=72.,TFROST=151.)
 +
 
 +
===Pluto===
 +
 
 +
  SatPrB = 9.2338
 +
  SatPrA = -724.9720
 +
  AMW    = 27.
 +
  FANON  = 0.1
 +
  KPREF  = 0
 +
 
 +
 
 +
===Titan===
 +
SatPrB = 9.2338
 +
SatPrA = -724.9720
 +
AMW    = 28.8
 +
FANON  = 0.016
 +
KPREF  = 0
 +
 
 +
 
 +
KRC function: KRC_Cond_Gas()
 +
 
 +
 
 +
 
 +
JBARE forces frost free at specified season. Default JBARE=0. (no forcing when = 0.)
 +
 
 +
JBARE corresponds to a CROCUS date, in Ls.
 +
 
 +
Note: the interface facilitates the use of JBARE so that JBARE ~ CROCUS date (Ls), but KRC has a more complex definition.
 +
 
 +
[[Image:TFROST_2.png|800px]]
 +
 
 +
out_1 = krc(lat = 75.)
 +
out_2 = krc(lat = 75, JBARE = 50.)

Latest revision as of 15:37, 14 April 2020

The coefficients for the saturation temperature relation for the condensing gas are now input parameters. KRC uses the Clausius-Clapeyron relation:


ln P = a - b/T


where P is pressure in Pascal and T is temperature in Kelvin.

This should be useful for Mars, Titan, Pluto etc.

SatPrA correspond to the Clausius-Clapeyron coefficient 'a'.
SatPrB corresponds to the Clausius-Clapeyron coeffcient 'b'.
LVFT = "T"

Values for CO2 are a = 27.9546 and b = 3182.48 (Provided by HHK, most likely from the Mars book). Values for N2 are a = 9.2338 and b = -724.9720 (Derived by SP using this table: File:Clausius Clapeyron.xlsx

The proper molecular weight should be input AMW. The Mass-fraction of mean atmosphere that is non-condensing should also be specified with FANON.

[edit] Mars

SatPrB = 3182.48
SatPrA = 27.9546
AMW    = 43.54             #SP's own calculation
FANON  = 0.040             #Paul R. Mahaffy, Science 2013 (SAM)
KPREF  = 1

TFROST 1.png

out_1 = krc(lat=72.,SatPrA=27.9546,SatPrB=3182.48,AMW=43.54,FANON=0.040,KPREF=1,LVFT="T" )
out_2 = krc(lat=72.,TFROST=151.)

[edit] Pluto

 SatPrB = 9.2338
 SatPrA = -724.9720
 AMW    = 27.
 FANON  = 0.1
 KPREF  = 0


[edit] Titan

SatPrB = 9.2338
SatPrA = -724.9720
AMW    = 28.8
FANON  = 0.016
KPREF  = 0


KRC function: KRC_Cond_Gas()


JBARE forces frost free at specified season. Default JBARE=0. (no forcing when = 0.)

JBARE corresponds to a CROCUS date, in Ls.

Note: the interface facilitates the use of JBARE so that JBARE ~ CROCUS date (Ls), but KRC has a more complex definition.

TFROST 2.png

out_1 = krc(lat = 75.)
out_2 = krc(lat = 75, JBARE = 50.)
Personal tools