Changeset 2706

Timestamp:

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

Author:

aslmd

Message:

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

Location:

trunk/LMDZ.GENERIC

Files:

• deftank/traceur.def.modern (modified) (1 diff)
• libf/phystd/condensation_generic_mod.F90 (modified) (3 diffs)
• libf/phystd/dyn1d/rcm1d.F (modified) (1 diff)
• libf/phystd/generic_cloud_common_h.F90 (modified) (1 diff)
• libf/phystd/initracer.F90 (modified) (5 diffs)
• libf/phystd/physiq_mod.F90 (modified) (1 diff)
• libf/phystd/tracer_h.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/deftank/traceur.def.modern

@@ -52 +52 @@
 #               is_recomb_qotf ! 1 if tracer recombination is done on-the-fly, else 0
 #                                (if 1, must have is_recomb_qset=0)
-#               is_generic     ! 1 if tracer is generic, else 0
+#               is_condensable  ! 1 if tracer is_condensable, else 0
 #               radius          ! radius value in meter
 #               rho             ! density in kg/m3

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

@@ -8 +8 @@
         use ioipsl_getin_p_mod, only: getin_p !-> to get the metallicity 
         use generic_cloud_common_h, only : RLVTT, cpp, &
-        Psat_generic,Lcpdqsat_generic,specie_parameters
+        Psat_generic,Lcpdqsat_generic,specie_parameters,specie_parameters_table
         USE tracer_h
         IMPLICIT none
@@ -86 +86 @@
         ! Let's loop on tracers 
         do iq=1,nq
-                if((is_generic(iq)==1) .and. (index(noms(iq),"vap") .ne. 0)) then
-                        ! Hyp : vap tracer index comes right before ice tracer index in traceur.def
+                if((is_condensable(iq)==1) .and. (index(noms(iq),"vap") .ne. 0)) then
+                        write(*,*) "the specie ", noms(iq)," is condensable, for generic condensation"
+!                       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
-                        igcm_generic_ice = iq+1
+
+                        ! We look for the corresponding ice traceur (before or after in the list of traceurs, maybe could be generalised to the whole list)
+                        if ((noms(iq)(9:len(trim(noms(iq)))-4) .eq. noms(iq+1)(9:len(trim(noms(iq+1)))-4)) .and. (index(noms(iq+1),"ice") .ne. 0)) then
+                                igcm_generic_ice = iq+1
+                        else if ((iq .gt. 1) .and. (noms(iq)(9:len(trim(noms(iq)))-4) .eq. noms(iq-1)(9:len(trim(noms(iq-1)))-4)) .and. (index(noms(iq-1),"ice") .ne. 0)) then
+                                igcm_generic_ice = iq-1
+                        else
+                                write(*,*) "ERROR : You set a vap traceur but you forgot to set the corresponding ice traceur"
+                        endif
+
                         ! Need to call specie_parameters of the considered specie 
                         call specie_parameters(noms(iq)(9:len(trim(noms(iq)))-4))
+                        !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
 
@@ -130 +146 @@
 
                         Enddo ! k= nlayer, 1, -1
-                endif !(is_generic(iq)==1) .and. (index(noms(iq),"vap") .ne. 0))
+                endif !(is_condensable(iq)==1) .and. (index(noms(iq),"vap") .ne. 0))
         enddo ! iq=1,nq
 

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

@@ -6 +6 @@
       use mod_grid_phy_lmdz, only : regular_lonlat
       use infotrac, only: nqtot, tname
-      use tracer_h, only: noms
+      use tracer_h, only: noms, is_condensable
       use surfdat_h, only: albedodat, phisfi, dryness, watercaptag,
      &                     zmea, zstd, zsig, zgam, zthe,

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)

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

@@ -116 +116 @@
          ALLOCATE(is_recomb_qotf(nqtot))
        ENDIF
-       IF (.NOT. allocated(is_generic)) allocate(is_generic(nq)) !LT
+       IF (.NOT. allocated(is_condensable)) allocate(is_condensable(nq)) !LT
        !! initialization
        alpha_lift(:)     = 0.
@@ -133 +133 @@
        radius(:)=0.
        qext(:)=0.
-       is_generic(:)= 0 !LT
+       is_condensable(:)= 0 !LT
        rho_q(:) = 0. !need to be init here if we want to read it from modern traceur with get_tracdat
 
@@ -400 +400 @@
 
       enddo ! of do iq=1,nq
+
+      ! 3. find condensable traceurs different from h2o and co2
+      do iq=1,nq
+        if ((index(noms(iq),"vap") .ne. 0) .and. (index(noms(iq),"h2o") .eq. 0) .and. (index(noms(iq),"co2") .eq. 0)) then
+          count=count+1
+        endif
+        if ((index(noms(iq),"ice") .ne. 0) .and. (index(noms(iq),"h2o") .eq. 0) .and. (index(noms(iq),"co2") .eq. 0)) then
+          count=count+1
+        endif
+
+      enddo ! of do iq=1,nq
       
       ! check that we identified all tracers:
@@ -430 +441 @@
 
       ! Calculate number of generic tracers
-      ngt = sum(is_generic)
+      ngt = sum(is_condensable)
       write(*,*) 'Number of generic tracer is  ngt = ',ngt
 !     Processing modern traceur options
@@ -608 +619 @@
               is_recomb_qotf(iq)
           end if
-          !option is_generic (LT)
-          if (index(tracline,'is_generic=') /=0) then 
-            read(tracline(index(tracline,'is_generic=') &
-              +len('is_generic='):),*) is_generic(iq)
+          !option is_condensable (LT)
+          if (index(tracline,'is_condensable=') /=0) then 
+            read(tracline(index(tracline,'is_condensable=') &
+              +len('is_condensable='):),*) is_condensable(iq)
             write(*,*) ' Parameter value (traceur.def) :'// &
-              ' is_generic=', is_generic(iq)       
+              ' is_condensable=', is_condensable(iq)       
           else
               write(*,*) ' Parameter value (default)     :'// &
-              ' is_generic=', is_generic(iq)       
+              ' is_condensable=', is_condensable(iq)       
           endif
           !option radius

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

@@ -32 +32 @@
                           alpha_lift, alpha_devil, qextrhor, &
                           igcm_h2o_ice, igcm_h2o_vap, igcm_dustbin, &
-                          igcm_co2_ice, nesp, is_chim, is_generic
+                          igcm_co2_ice, nesp, is_chim, is_condensable
       use time_phylmdz_mod, only: ecritphy, iphysiq, nday
       use phyetat0_mod, only: phyetat0

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

@@ -40 +40 @@
        integer, save, allocatable :: is_recomb_qotf(:) ! 1 if tracer recombination is done on-the-fly, else 0 (if 1, must have is_recomb_qset=0)
 !$OMP THREADPRIVATE(is_recomb,is_recomb_qset,is_recomb_qotf)
-       integer, save, allocatable :: is_generic(:)      ! 1 if tracer is generic, else 0 (added LT)
-!$OMP THREADPRIVATE(is_generic)   !also added by LT
+       integer, save, allocatable :: is_condensable(:)      ! 1 if tracer is generic, else 0 (added LT)
+!$OMP THREADPRIVATE(is_condensable)   !also added by LT
 ! tracer indexes: these are initialized in initracer and should be 0 if the
 !                 corresponding tracer does not exist