Changeset 3233 for trunk/LMDZ.GENERIC/libf/phystd/dyn1d

Timestamp:

Feb 23, 2024, 10:22:20 AM (9 months ago)

Author:

gmilcareck

Message:

Add the possibility to use a fixed vertical molar fraction profile for the
collision-induced absorption like the correlated-k.

File:

• trunk/LMDZ.GENERIC/libf/phystd/dyn1d/kcm1d.F90 (modified) (6 diffs)

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/kcm1d.F90

@@ -8 +8 @@
   use comsaison_h, only: mu0, fract, dist_star
   use planete_mod
-  use callkeys_mod, only: pceil, tstrat, tracer, global1d
+  use callkeys_mod, only: pceil, tstrat, tracer, global1d, &
+                          varspec, varspec_data, nvarlayer
   use inifis_mod, only: inifis
   use aerosol_mod, only: iniaerosol
@@ -20 +21 @@
   use geometry_mod, only: init_geometry
   use dimphy, only : init_dimphy
+  use gases_h, only: ngasmx
 
   implicit none
@@ -87 +89 @@
   real int_dtauv(1,llm,L_NSPECTV)
   real Eatmtot
+  
+  ! For fixed variable molar mass
+  real, dimension(:),allocatable,save :: p_var
+  real, dimension(:),allocatable,save :: mu_var
+  real, dimension(:,:),allocatable,save :: frac_var
 
   integer ierr
@@ -116 +123 @@
   ALLOCATE(fract(1))
   ALLOCATE(glat(1))
+
+! Initialize fixed variable mu
+!        ~~~~~~~~~~~~~~~~~~~~~~~~~~~~   
+
+  if(varspec) then
+    IF (.NOT.ALLOCATED(p_var))  ALLOCATE(p_var(nvarlayer))
+    IF (.NOT.ALLOCATED(mu_var))  ALLOCATE(mu_var(nvarlayer))
+    IF (.NOT.ALLOCATED(frac_var))  ALLOCATE(frac_var(nvarlayer,ngasmx))
+    p_var = 0.0
+    mu_var = 0.0
+    frac_var = 0.0
+    dt_file=TRIM(varspec_data)
+    open(34,file=dt_file,form='formatted',status='old',iostat=ios)
+    if (ios.ne.0) then        ! file not found
+      write(*,*) 'Error from varspec'
+      write(*,*) 'Data file ',trim(varspec_data),' not found.'
+      write(*,*) 'Check that the data is in your run repository.'
+      call abort_physic !a verifier
+    else
+      do k=1,nvarlayer
+        read(34,*) p_var(k), mu_var(k),frac_var(k,1:ngasmx)
+        !The order of columns in frac_var must correspond to order of molecules gases.def 
+        !The format of your file must be:
+        ! pressure(k) molar_mass(k), molar_fraction_of_gas_1(k), molar_fraction_of_gas_2(k),..., molar_fraction_of_gas_ngasmx(k)
+      enddo
+    endif
+    close(34)
+  endif
 
   !  Load parameters from "run.def" 
@@ -386 +421 @@
           int_dtaui,int_dtauv,                            &
           tau_col,cloudfrac,totcloudfrac,                 &
-          .false.,firstcall,lastcall)
+          .false.,p_var,frac_var,firstcall,lastcall)
 
      !write(*,*) 'BASE 3'
@@ -438 +473 @@
        int_dtaui,int_dtauv,                                    &
        tau_col,cloudfrac,totcloudfrac,                         &
-       .false.,firstcall,lastcall)
+       .false.,p_var,frac_var,firstcall,lastcall)