Changeset 3653 for trunk/LMDZ.GENERIC/libf/phystd/initracer.F90

Timestamp:

Feb 27, 2025, 10:10:41 AM (4 months ago)

Author:

gmilcareck

Message:

Generic PCM:
Adding molecular weight for each gas used by the model.
Updating molecular weight for each tracers (to be consistent with gases).
GM

File:

• trunk/LMDZ.GENERIC/libf/phystd/initracer.F90 (modified) (5 diffs)

trunk/LMDZ.GENERIC/libf/phystd/initracer.F90

@@ -227 +227 @@
         if (noms(iq).eq."co2") then
           igcm_co2=iq
-          mmol(igcm_co2)=44.
+          mmol(igcm_co2)=44.010
           count=count+1
 !          write(*,*) 'co2: count=',count
@@ -233 +233 @@
         if (noms(iq).eq."co2_ice") then
           igcm_co2_ice=iq
-          mmol(igcm_co2_ice)=44.
+          mmol(igcm_co2_ice)=44.010
           count=count+1
 !          write(*,*) 'co2_ice: count=',count
@@ -239 +239 @@
         if (noms(iq).eq."h2o_vap") then
           igcm_h2o_vap=iq
-          mmol(igcm_h2o_vap)=18.
+          mmol(igcm_h2o_vap)=18.01528
           count=count+1
 !          write(*,*) 'h2o_vap: count=',count
@@ -245 +245 @@
         if (noms(iq).eq."h2o_ice") then
           igcm_h2o_ice=iq
-          mmol(igcm_h2o_ice)=18.
+          mmol(igcm_h2o_ice)=18.01528
           count=count+1
 !          write(*,*) 'h2o_ice: count=',count
@@ -251 +251 @@
         if (noms(iq).eq."co") then
           igcm_co=iq
-          mmol(igcm_co)=28.
+          mmol(igcm_co)=28.010
           count=count+1
         endif
         if (noms(iq).eq."o") then
           igcm_o=iq
-          mmol(igcm_o)=16.
+          mmol(igcm_o)=15.9994
           count=count+1
         endif
         if (noms(iq).eq."o1d") then
           igcm_o1d=iq
-          mmol(igcm_o1d)=16.
+          mmol(igcm_o1d)=15.9994
           count=count+1
         endif
         if (noms(iq).eq."o2") then
           igcm_o2=iq
-          mmol(igcm_o2)=32.
+          mmol(igcm_o2)=31.998
           count=count+1
         endif
         if (noms(iq).eq."o3") then
           igcm_o3=iq
-          mmol(igcm_o3)=48.
+          mmol(igcm_o3)=47.9982
           count=count+1
         endif
         if (noms(iq).eq."h") then
           igcm_h=iq
-          mmol(igcm_h)=1.
+          mmol(igcm_h)=1.00784
           count=count+1
         endif
         if (noms(iq).eq."h2") then
           igcm_h2=iq
-          mmol(igcm_h2)=2.
+          mmol(igcm_h2)=2.016
           count=count+1
         endif
         if (noms(iq).eq."oh") then
           igcm_oh=iq
-          mmol(igcm_oh)=17.
+          mmol(igcm_oh)=17.008
           count=count+1
         endif
         if (noms(iq).eq."ho2") then
           igcm_ho2=iq
-          mmol(igcm_ho2)=33.
+          mmol(igcm_ho2)=33.0067
           count=count+1
         endif
         if (noms(iq).eq."h2o2") then
           igcm_h2o2=iq
-          mmol(igcm_h2o2)=34.
+          mmol(igcm_h2o2)=34.0147
           count=count+1
         endif
         if (noms(iq).eq."n2") then
           igcm_n2=iq
-          mmol(igcm_n2)=28.
+          mmol(igcm_n2)=28.0134
           count=count+1
         endif
         if (noms(iq).eq."ch4") then
           igcm_ch4=iq
-          mmol(igcm_ch4)=16.
+          mmol(igcm_ch4)=16.0425
           count=count+1
         endif
         if (noms(iq).eq."ar") then
           igcm_ar=iq
-          mmol(igcm_ar)=40.
+          mmol(igcm_ar)=39.948
           count=count+1
         endif
         if (noms(iq).eq."n") then
           igcm_n=iq
-          mmol(igcm_n)=14.
+          mmol(igcm_n)=14.0067
           count=count+1
         endif
         if (noms(iq).eq."no") then
           igcm_no=iq
-          mmol(igcm_no)=30.
+          mmol(igcm_no)=30.0061
           count=count+1
         endif
         if (noms(iq).eq."no2") then
           igcm_no2=iq
-          mmol(igcm_no2)=46.
+          mmol(igcm_no2)=46.0055 
           count=count+1
         endif
         if (noms(iq).eq."n2d") then
           igcm_n2d=iq
-          mmol(igcm_n2d)=28.
+          mmol(igcm_n2d)=28.0134
           count=count+1
         endif
         if (noms(iq).eq."ch3") then
           igcm_ch3=iq
-          mmol(igcm_ch3)=15.
+          mmol(igcm_ch3)=15.034
           count=count+1
         endif
         if (noms(iq).eq."ch") then
           igcm_ch=iq
-          mmol(igcm_ch)=13.
+          mmol(igcm_ch)=13.0186
           count=count+1
         endif
         if (noms(iq).eq."3ch2") then
           igcm_3ch2=iq
-          mmol(igcm_3ch2)=14.
+          mmol(igcm_3ch2)=14.0266
           count=count+1
         endif
         if (noms(iq).eq."1ch2") then
           igcm_1ch2=iq
-          mmol(igcm_1ch2)=14.
+          mmol(igcm_1ch2)=14.0266
           count=count+1
         endif
         if (noms(iq).eq."cho") then
           igcm_cho=iq
-          mmol(igcm_cho)=29.
+          mmol(igcm_cho)=29.018
           count=count+1
         endif
         if (noms(iq).eq."ch2o") then
           igcm_ch2o=iq
-          mmol(igcm_ch2o)=30.
+          mmol(igcm_ch2o)=30.031
           count=count+1
         endif
         if (noms(iq).eq."ch3o") then
           igcm_ch3o=iq
-          mmol(igcm_ch3o)=31.
+          mmol(igcm_ch3o)=31.0339
           count=count+1
         endif
         if (noms(iq).eq."c") then
           igcm_c=iq
-          mmol(igcm_c)=12.
+          mmol(igcm_c)=12.0107
           count=count+1
         endif
         if (noms(iq).eq."c2") then
           igcm_c2=iq
-          mmol(igcm_c2)=24.
+          mmol(igcm_c2)=24.0214 
           count=count+1
         endif
         if (noms(iq).eq."c2h") then
           igcm_c2h=iq
-          mmol(igcm_c2h)=25.
+          mmol(igcm_c2h)=25.0293
           count=count+1
         endif
         if (noms(iq).eq."c2h2") then
           igcm_c2h2=iq
-          mmol(igcm_c2h2)=26.
+          mmol(igcm_c2h2)=26.038
           count=count+1
         endif
         if (noms(iq).eq."c2h3") then
           igcm_c2h3=iq
-          mmol(igcm_c2h3)=27.
+          mmol(igcm_c2h3)=27.0452
           count=count+1
         endif
         if (noms(iq).eq."c2h4") then
           igcm_c2h4=iq
-          mmol(igcm_c2h4)=28.
+          mmol(igcm_c2h4)=28.053
           count=count+1
         endif
         if (noms(iq).eq."c2h6") then
           igcm_c2h6=iq
-          mmol(igcm_c2h6)=30.
+          mmol(igcm_c2h6)=30.069
           count=count+1
         endif
         if (noms(iq).eq."ch2co") then
           igcm_ch2co=iq
-          mmol(igcm_ch2co)=42.
+          mmol(igcm_ch2co)=42.0367
           count=count+1
         endif
         if (noms(iq).eq."ch3co") then
           igcm_ch3co=iq
-          mmol(igcm_ch3co)=43.
+          mmol(igcm_ch3co)=43.0446
           count=count+1
         endif