Changeset 3669 for trunk/LMDZ.PLUTO/libf/phypluto/calc_cpp_mugaz.F90

Timestamp:

Mar 5, 2025, 2:13:17 PM (4 months ago)

Author:

afalco

Message:

Pluto: add C2H4 as gas.
AF

File:

• trunk/LMDZ.PLUTO/libf/phypluto/calc_cpp_mugaz.F90 (modified) (3 diffs)

trunk/LMDZ.PLUTO/libf/phypluto/calc_cpp_mugaz.F90

@@ -63 +63 @@
                mugaz_c = mugaz_c + 30.07*gfrac(igas)
             elseif(igas.eq.igas_C2H2)then
-               ! C2H2 http://encyclopedia.airliquide.com/Encyclopedia.asp?GasID=1
+               ! C2H2 https://encyclopedia.airliquide.com/fr/acetylene
                mugaz_c = mugaz_c + 26.04*gfrac(igas)
             ! GG MODIF JAN2019
+            elseif(igas.eq.igas_C2H4)then
+               mugaz_c = mugaz_c + 28.054*gfrac(igas)
             elseif(igas.eq.igas_CO)then
                mugaz_c = mugaz_c + 28.01*gfrac(igas)
@@ -91 +93 @@
             ! ignore variable gas in cpp calculation
          else
-            ! all values at 300 K from Engineering Toolbox
+            ! values at minimum temperature from Engineering Toolbox
+            ! (better for Pluto)
             if(igas.eq.igas_CO2)then
-               !cpp_c   = cpp_c   + 0.744*gfrac(igas) ! @ ~210 K (better for
-               !Mars conditions)
-               cpp_c   = cpp_c   + 0.846*gfrac(igas)*44.01/mugaz_c
+            ! https://www.engineeringtoolbox.com/carbon-dioxide-d_974.html
+               cpp_c   = cpp_c   + 0.709*gfrac(igas)*44.01/mugaz_c
             elseif(igas.eq.igas_N2)then
-               cpp_c   = cpp_c   + 1.040*gfrac(igas)*28.01/mugaz_c
+               cpp_c   = cpp_c   + 1.039*gfrac(igas)*28.01/mugaz_c
             elseif(igas.eq.igas_H2)then
                cpp_c   = cpp_c   + 14.31*gfrac(igas)*2.01/mugaz_c
@@ -109 +111 @@
                cpp_c   = cpp_c   + 1.003*gfrac(igas)*34.08/mugaz_c ! from wikipedia...
             elseif(igas.eq.igas_CH4)then
-               cpp_c   = cpp_c   + 2.226*gfrac(igas)*16.04/mugaz_c
+               cpp_c   = cpp_c   + 2.087*gfrac(igas)*16.04/mugaz_c
             elseif(igas.eq.igas_NH3)then
                cpp_c   = cpp_c   + 2.175*gfrac(igas)*17.03/mugaz_c
                print*,'WARNING, cpp for NH3 may be for liquid'
             elseif(igas.eq.igas_C2H6)then
-               ! C2H6
                ! http://encyclopedia.airliquide.com/Encyclopedia.asp?GasID=28
-               cpp_c   = cpp_c   + 1.763*gfrac(igas)*30.07/mugaz_c
+               cpp_c   = cpp_c   + 1.535*gfrac(igas)*30.07/mugaz_c
             elseif(igas.eq.igas_C2H2)then
-               ! C2H2
-               ! http://encyclopedia.airliquide.com/Encyclopedia.asp?GasID=1
+               ! https://encyclopedia.airliquide.com/fr/acetylene
                cpp_c   = cpp_c   + 1.575*gfrac(igas)*26.04/mugaz_c
             !!!!! MODIF GG JAN 2019  (check source values !!)
+            elseif(igas.eq.igas_C2H4)then
+               ! https://www.engineeringtoolbox.com/ethylene-ethene-C2H4-properties-d_2104.html
+               cpp_c   = cpp_c   + 1.53*gfrac(igas)*28.054/mugaz_c
+            !!!!! MODIF GG JAN 2019  (check source values !!)
             elseif(igas.eq.igas_CO)then
-               cpp_c   = cpp_c   + 1.0425*gfrac(igas)*28.01/mugaz_c
+               cpp_c   = cpp_c   + 1.039*gfrac(igas)*28.01/mugaz_c
             elseif(igas.eq.igas_OCS)then
                cpp_c   = cpp_c   + 0.6909*gfrac(igas)*60.07/mugaz_c