Changeset 2831 for trunk/LMDZ.GENERIC/libf/phystd/calc_cpp_mugaz.F90

Timestamp:

Nov 23, 2022, 4:41:34 PM (2 years ago)

Author:

emillour

Message:

Generic PCM:
Add the possibility to include Venus-like aerosols (triggered by option
aerovenus=.true. in callphys.def); baseline is to use 5 distinct scatterers
but each may be turned on/off (via aerovenus1, aerovenus2, aerovenus2p,
aerovenus3, aerovenusUV flags which may be specified in callphys.def).
GG

File:

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

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

@@ -24 +24 @@
       implicit none
 
+      character(len=20),parameter :: myname="calc_cpp_mugaz"
       real cpp_c   
       real mugaz_c
@@ -64 +65 @@
                ! C2H2 http://encyclopedia.airliquide.com/Encyclopedia.asp?GasID=1
                mugaz_c = mugaz_c + 26.04*gfrac(igas)
+            ! GG MODIF JAN2019
+            elseif(igas.eq.igas_CO)then
+               mugaz_c = mugaz_c + 28.01*gfrac(igas)
+            elseif(igas.eq.igas_OCS)then
+               ! https://pubchem.ncbi.nlm.nih.gov/compound/carbonyl_sulfide
+               mugaz_c = mugaz_c + 60.0751*gfrac(igas)
+            elseif(igas.eq.igas_HCl)then
+               mugaz_c = mugaz_c + 36.46*gfrac(igas)
+            elseif(igas.eq.igas_HF)then
+               mugaz_c = mugaz_c + 20.01*gfrac(igas)
             else
+              if (is_master) then
                print*,'Error in calc_cpp_mugaz: Gas species not recognised!'
-               call abort
+              endif
+              call abort_physic(myname,'Gas species not recognised!',1)
             endif
          endif
@@ -108 +121 @@
                ! http://encyclopedia.airliquide.com/Encyclopedia.asp?GasID=1
                cpp_c   = cpp_c   + 1.575*gfrac(igas)*26.04/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
+            elseif(igas.eq.igas_OCS)then
+               cpp_c   = cpp_c   + 0.6909*gfrac(igas)*60.07/mugaz_c
+            elseif(igas.eq.igas_HCl)then
+               cpp_c   = cpp_c   + 1.7087*gfrac(igas)*36.46/mugaz_c
+            elseif(igas.eq.igas_HF)then
+               cpp_c   = cpp_c   + 1.4553*gfrac(igas)*20.01/mugaz_c
             else
+              if (is_master) then
                print*,'Error in calc_cpp_mugaz: Gas species not recognised!'
-               call abort
+              endif
+              call abort_physic(myname,'Gas species not recognised!',1)
             endif
          endif
@@ -122 +146 @@
       cpp_c = 1000.0*cpp_c
 
-      print*,'Cp in calc_cpp_mugaz is ',cpp_c,'J kg^-1 K^-1'
-      print*,'Mg in calc_cpp_mugaz is ',mugaz_c,'amu'
+      if (is_master) then
+        print*,'Cp in calc_cpp_mugaz is ',cpp_c,'J kg^-1 K^-1'
+        print*,'Mg in calc_cpp_mugaz is ',mugaz_c,'amu'
+      endif
 
       if(((abs(1.-cpp/cpp_c).gt.1.e-6) .or.  &
@@ -146 +172 @@
 
 
-      return
     end subroutine calc_cpp_mugaz