Changeset 2694 for trunk/LMDZ.GENERIC/libf/phystd/genericcommon_h.F90

Timestamp:

Jun 21, 2022, 11:05:32 AM (3 years ago)

Author:

aslmd

Message:

fixed uninitialized variable Lcp in largescale_generic+ typo in index of tendencies after largescale_generic in physiq_mod

File:

• trunk/LMDZ.GENERIC/libf/phystd/genericcommon_h.F90 (modified) (7 diffs)

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

@@ -34 +34 @@
             Tref = 2303.0 !K
             Pref =   1.0E5 !in Pa
+            m = 24.3055
             RLVTT = delta_vapH/(m/1000.)
             metallicity_coeff=1.0*log(10.0)
-            m = 24.3055
+            
         else if (trim(specname) .eq. "Na") then
             print*,"Loading data for Na"
@@ -42 +43 @@
             Tref = 1624.0
             Pref =  1.0E5
+            m = 22.9898
             RLVTT = delta_vapH/(m/1000.)
             metallicity_coeff=0.5*log(10.0)
-            m = 22.9898
+            
         else if (trim(specname) .eq. "Fe") then 
             print*,"Loading data for Fe"
@@ -50 +52 @@
             Tref = 2903.0
             Pref  = 1.0E5
+            m = 55.8450
             RLVTT = delta_vapH/(m/1000.)
             metallicity_coeff=0.0*log(10.0)
-            m = 55.8450
+            
         else if (trim(specname) .eq. "Cr") then
             print*,"Loading data for Cr"
@@ -58 +61 @@
             Tref = 2749.0
             Pref  = 1.0E5
+            m = 51.9961
             RLVTT = delta_vapH/(m/1000.)
             metallicity_coeff=0.0*log(10.0)
-            m = 51.9961
+            
         else if (trim(specname) .eq. "KCl") then 
             print*,"Loading data for KCl"
@@ -68 +72 @@
             metallicity_coeff=0.0*log(10.0)
             m = 74.5498
+            RLVTT = delta_vapH/(m/1000.)
         else if (trim(specname) .eq. "Mn") then 
             print*,"Loading data for Mn"
@@ -75 +80 @@
             metallicity_coeff=1.0*log(10.0)
             m = 54.9380
+            RLVTT = delta_vapH/(m/1000.)
         else if (trim(specname) .eq. "Zn") then 
             print*,"Loading data for Zn"
@@ -82 +88 @@
             metallicity_coeff=1.0*log(10.0)
             m = 65.38
+            RLVTT = delta_vapH/(m/1000.)
         else
             print*,"Unknow species (not in Mg, Fe, Na, KCl, Cr, Mn or Zn)"