Changeset 2707 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:

table_tracers_condensable can be used to take back data of condensable tracers

File:

• trunk/LMDZ.GENERIC/libf/phystd/generic_cloud_common_h.F90 (modified) (2 diffs)

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

@@ -102 +102 @@
     !   Load the adequate set of parameters for specname
     !   From a table of traceurs
+    !
+    !   Also set up few others useful parameters, such as epsi=m/mugaz, RLVTT and
+    !   the metallicity_coeff.
+    !   Authors 
+    !   --------
+    !   Noé Clément (2022) 
     !============================================================================
 
         character(*), intent(in) :: specname
-        integer k
+        integer :: ios
         character(len=500):: table_traceurs_line ! table_traceurs_line lines with parameters
 
-        open(117,file='table_traceurs',form='formatted',status='old')
+        open(117,file='table_tracers_condensable',form='formatted',status='old')
 
         read(117,'(A)') table_traceurs_line
 
         do
-            read(117,'(A)') table_traceurs_line
+            read(117,'(A)', iostat=ios) 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
+                if (index(table_traceurs_line,'delta_vapH='   ) /= 0) then
+                    read(table_traceurs_line(index(table_traceurs_line,'delta_vapH=')+len('delta_vapH='):),*) delta_vapH
+                    print*, '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
+                    print*, 'Tref', Tref
                 end if 
                 if (index(table_traceurs_line,'Pref='   ) /= 0) then
                     read(table_traceurs_line(index(table_traceurs_line,'Pref=')+len('Pref='):),*) Pref
+                    print*, 'Pref', Pref
                 end if
                 if (index(table_traceurs_line,'mass='   ) /= 0) then
                     read(table_traceurs_line(index(table_traceurs_line,'mass=')+len('mass='):),*) m
+                    print*, '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
+                   print*, 'metallicity_coeff', metallicity_coeff
+                end if
+
+                ios=1
             end if
-        
+
+            if (ios /= 0) exit
+            
         end do
         ! RLVTT
@@ -142 +154 @@
         !if (is_master)
         close(117)
+
+        RLVTT=delta_vapH/(m/1000.)
+
+        write(*,*) 'RLVTT', RLVTT
+
+        epsi = m/mugaz
 
     end subroutine specie_parameters_table