Changeset 2707

Timestamp:

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

Author:

aslmd

Message:

table_tracers_condensable can be used to take back data of condensable tracers

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• condensation_generic_mod.F90 (modified) (2 diffs)
• generic_cloud_common_h.F90 (modified) (2 diffs)

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

@@ -87 +87 @@
         do iq=1,nq
                 if((is_condensable(iq)==1) .and. (index(noms(iq),"vap") .ne. 0)) then
-                        write(*,*) "the specie ", noms(iq)," is condensable, for generic condensation"
+                        write(*,*) "There is a specie which is condensable, for generic condensation : ", noms(iq)
+
 !                       Let's get the index of our tracers (we look for igcm _generic_vap and igcm_generic_ice)
                         ! tname_ice = trim(noms(iq)(1:len(tname_ice)-3))//"ice" 
                         ! print*,trim(adjustl(trim(noms(iq))(9:len(trim(noms(iq)))-4))) !testing here, should go away
-                        print*,noms(iq)(9:len(trim(noms(iq)))-4)
                         ! stop
+
                         igcm_generic_vap=iq
 
@@ -108 +109 @@
                         !call specie_parameters_table(noms(iq)(9:len(trim(noms(iq)))-4))
                         
-                        Lcp=RLVTT/cpp !need to be init here, otherwise we don't know RLVTT yet
+                        Lcp=RLVTT/cpp ! need to be init here
 
                         !  Vertical loop (from top to bottom)

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