Changeset 2706 for trunk/LMDZ.GENERIC/libf/phystd/generic_cloud_common_h.F90

Timestamp:

Jun 21, 2022, 11:05:38 AM (2 years ago)

Author:

aslmd

Message:

is_condensable has replaced is_generic, some changes in getting back tracers names

File:

• trunk/LMDZ.GENERIC/libf/phystd/generic_cloud_common_h.F90 (modified) (1 diff)

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

@@ -94 +94 @@
         epsi = m/mugaz
     end subroutine specie_parameters
+
+
+    subroutine specie_parameters_table(specname)
+        
+        implicit none
+    !============================================================================
+    !   Load the adequate set of parameters for specname
+    !   From a table of traceurs
+    !============================================================================
+
+        character(*), intent(in) :: specname
+        integer k
+        character(len=500):: table_traceurs_line ! table_traceurs_line lines with parameters
+
+        open(117,file='table_traceurs',form='formatted',status='old')
+
+        read(117,'(A)') table_traceurs_line
+
+        do
+            read(117,'(A)') table_traceurs_line
+            
+            if (index(table_traceurs_line,specname) /= 0) then
+
+                write(*,*) table_traceurs_line
+                
+                if (index(table_traceurs_line,'deltavapH='   ) /= 0) then
+                    read(table_traceurs_line(index(table_traceurs_line,'deltavapH=')+len('deltavapH='):),*) delta_vapH
+                    write(*,*) 'delta_vapH ', delta_vapH
+                end if
+                if (index(table_traceurs_line,'Tref='   ) /= 0) then
+                    read(table_traceurs_line(index(table_traceurs_line,'Tref=')+len('Tref='):),*) Tref
+                end if 
+                if (index(table_traceurs_line,'Pref='   ) /= 0) then
+                    read(table_traceurs_line(index(table_traceurs_line,'Pref=')+len('Pref='):),*) Pref
+                end if
+                if (index(table_traceurs_line,'mass='   ) /= 0) then
+                    read(table_traceurs_line(index(table_traceurs_line,'mass=')+len('mass='):),*) m
+                end if
+                if (index(table_traceurs_line,'metallicity_coeff='   ) /= 0) then
+                   read(table_traceurs_line(index(table_traceurs_line,'metallicity_coeff=')+len('metallicity_coeff='):),*) metallicity_coeff
+                end if
+            end if
+        
+        end do
+        ! RLVTT
+
+        !if (is_master)
+        close(117)
+
+    end subroutine specie_parameters_table
 
     subroutine Psat_generic(T,p,metallicity,psat,qsat)