Changeset 1889

Timestamp:

Jan 9, 2018, 12:26:53 PM (7 years ago)

Author:

jvatant

Message:

Making chemistry handling more flexible - Step 2
+ Added the calculation of the pressure grid in newstart
using Vervack profile in gr_kim_vervack routine
+ Next step : regridding !
--JVO

Location:

trunk/LMDZ.TITAN

Files:

• README (modified) (2 diffs)
• libf/dynphy_lonlat/phytitan/gr_kim_vervack.F90 (added)
• libf/dynphy_lonlat/phytitan/lect_start_archive.F (modified) (8 diffs)
• libf/dynphy_lonlat/phytitan/newstart.F (modified) (7 diffs)
• libf/phytitan/comchem_h.F90 (modified) (1 diff)
• libf/phytitan/iostart.F90 (modified) (2 diffs)
• libf/phytitan/phyredem.F90 (modified) (2 diffs)
• libf/phytitan/tabfi_mod.F90 (modified) (5 diffs)

trunk/LMDZ.TITAN/README

@@ -1343 +1343 @@
 Big modifs of the tracer gestion/init in the physiq with a new query by names (see tracer_h )
 
-== 20/12/2017 - ... : r1871- ... == JVO
+== 20/12/2017 - ... : r1871-1886-1888 ... == JVO
 Management of the chemistry within startfi. Open the way to :
  1) run with chemistry in another resolution than 32x48 !
@@ -1352 +1352 @@
 rather than read in the middle of calchim as before from compo files
 + An "upper_chemistry_layers" dimension is now handled in startfi files
-+ This dimension is also stored in tabfi
 + Added a comchem_h.F90 module for all the stuff related to chemistry
++ In newstart we calculates the pressure grid above GCM top using Vervack profile
 + ...

trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/lect_start_archive.F

@@ -5 +5 @@
 
 !      USE surfdat_h
+      USE comchem_h, ONLY: nlaykimold, preskimold
       USE comsoil_h, ONLY: nsoilmx, layer, mlayer, volcapa, inertiedat
       USE infotrac, ONLY: tname, nqtot
@@ -46 +47 @@
 c------------------------------------
       INTEGER   imold,jmold,lmold,nsoilold,nqold
-
 
 c Variables pour les lectures des fichiers "ini" 
@@ -147 +147 @@
       real, dimension(:,:), allocatable :: emisold
       real, dimension(:,:,:,:), allocatable :: qold
-
+      
       real tab_cntrl(100)
 
@@ -158 +158 @@
       logical :: therminertia_3D=.true. ! flag
 ! therminertia_3D=.true. if thermal inertia is 3D and read from datafile
+
 c Variable intermediaires iutilise pour l'extrapolation verticale 
 c----------------------------------------------------------------
@@ -247 +248 @@
       endif
 
-
       if (nsoilold.ne.nsoilmx) then ! interpolation will be required
         depthinterpol=.true.
       endif
+
+! 1.2.2 find out the # of upper chemistry layers
+
+      ierr= NF_INQ_DIMID(nid,"upper_chemistry_layers",dimid)
+      ierr= NF_INQ_DIMLEN(nid,dimid,nlaykimold)
+
+      ! NB : The regriding, if needed cannot be done here since the new
+      ! pressure grid is only computed at the end of newstart
 
 ! 1.3 Report dimensions
@@ -266 +274 @@
       write(*,*) '  Otherwise, set nsoilmx -in dimphys.h- to: ',nsoilold
       endif
+      write(*,*) "upper_chemistry_layers: ",nlaykimold
       write(*,*) "time lenght: ",timelen
       write(*,*) 
@@ -295 +304 @@
       allocate(mlayerold(nsoilold))
       allocate(qsurfold(imold+1,jmold+1,nqtot))
+      
+      allocate(preskimold(nlaykimold))
 
       allocate(var (imold+1,jmold+1,llm))
@@ -512 +523 @@
        endif
       endif
-
-c-----------------------------------------------------------------------
-c 3.6 Lecture geopotentiel au sol
+      
+c-----------------------------------------------------------------------
+c 3.6 Read upper chemistry mid-layer pressure
+c-----------------------------------------------------------------------  
+    
+      ierr=NF_INQ_VARID(nid,"preskim",nvarid)
+      IF (ierr .NE. NF_NOERR) THEN
+         PRINT*, "lect_start_archive: Le champ <preskim> est absent"
+         CALL abort
+      ENDIF
+#ifdef NC_DOUBLE
+      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, preskimold)
+#else
+      ierr = NF_GET_VAR_REAL(nid, nvarid, preskimold)
+#endif
+      IF (ierr .NE. NF_NOERR) THEN
+         PRINT*, "lect_start_archive: Lecture echouee pour <preskim>"
+         CALL abort
+      ENDIF
+
+c-----------------------------------------------------------------------
+c 3.7 Lecture geopotentiel au sol
 c-----------------------------------------------------------------------
 c

trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/newstart.F

@@ -17 +17 @@
      &                              is_master
       use infotrac, only: infotrac_init, nqtot, tname
+      USE comchem_h, ONLY: nlaykim_up, nlaykimold
       USE comsoil_h, ONLY: nsoilmx, layer, mlayer, inertiedat
       USE surfdat_h, ONLY: phisfi, albedodat,
@@ -130 +131 @@
       real tab_cntrl(100)
       real tab_cntrl_bis(100)
-
+      
 c variables diverses
 c-------------------
@@ -258 +259 @@
 
       endif
-
 
 c=======================================================================
@@ -536 +536 @@
         CALL lect_start_archive(ngridmx,llm,
      &   date,tsurf,tsoil,emis,q2,
-     &   t,ucov,vcov,ps,teta,phisold_newgrid,q,qsurf,
-     &   surfith,nid)
+     &   t,ucov,vcov,ps,teta,phisold_newgrid,
+     &   q,qsurf,surfith,nid)
         write(*,*) "OK, read start_archive file"
         ! copy soil thermal inertia
@@ -572 +572 @@
       write(*,*) 'q=x : give a specific uniform value to one tracer'
       write(*,*) 'q=profile    : specify a profile for a tracer'
-!      write(*,*) 'ini_q : tracers initialisation for chemistry, water an
-!     $d ice   '
-!      write(*,*) 'ini_q-H2O : tracers initialisation for chemistry and 
-!     $ice '
-!      write(*,*) 'ini_q-iceH2O : tracers initialisation for chemistry on
-!     $ly '
       write(*,*) 'isotherm  : Isothermal Temperatures, wind set to zero'
       write(*,*) 'radequi   : Earth-like radiative equilibrium temperature
@@ -1071 +1065 @@
       CALL inifilr 
       CALL pression(ip1jmp1, ap, bp, ps, p3d)
+      
+
+c=========================================================================
+c  Calcul de la dimension verticale pour la chimie  - JVO 2017
+c  start_archive seulement, la grille verticale pouvant avoir ete modifiee
+c==========================================================================
+
+      IF (choix_1.eq.0) THEN
+
+        ! Calculate the # of upper chemistry layers with the "new" pressure grid
+        ! For this we use Vervack profile for upper atmosphere with dz=10km
+
+        CALL gr_kim_vervack
+
+        WRITE(*,*)
+        WRITE(*,*) " With the compiled vertical grid we found :"
+        WRITE(*,*) " Number of upper chemistry layers =", nlaykim_up
+        
+        ! Regriding is then done, if needed
+        
+        IF (nlaykimold.ne.nlaykim_up) THEN
+
+          WRITE(*,*) " Warning, nlaykimold=", nlaykimold
+          WRITE(*,*) ' which implies that a regriding on upper chemistry
+     & will be performed.'
+          WRITE(*,*)
+          
+!          CALL  regrid_kim
+          
+        ENDIF
+ 
+      endif ! of if (choix_1.eq.0)      
+ 
 
 c-----------------------------------------------------------------------
@@ -1125 +1152 @@
          day_ini=int(date)
       endif
+
 c
       CALL geopot  ( ip1jmp1, teta  , pk , pks,  phis  , phi   )

trunk/LMDZ.TITAN/libf/phytitan/comchem_h.F90

@@ -19 +19 @@
 !$OMP_THREADPRIVATE(zlay_kim)
 
+
+  ! Variable and allocatables for regriding chemistry newstart
+  ! ----------------------------------------------------------
+    
+  INTEGER :: nlaykimold ! Number of upper atm. layers for chemistry in the start_archive file
+  REAL, ALLOCATABLE, DIMENSION(:) :: preskimold ! Pressure grid of upper chemistry in the start_archive file
 
   ! Allocatable arrays for start2archive

trunk/LMDZ.TITAN/libf/phytitan/iostart.F90

@@ -944 +944 @@
                     NF90_PUT_ATT, NF90_NOERR, nf90_strerror, &
                     nf90_inq_dimid, nf90_inquire_dimension, NF90_INQ_VARID
+  USE comchem_h, only: nlaykim_up
   USE comsoil_h, only: nsoilmx
   USE mod_phys_lmdz_para, only: is_master
@@ -997 +998 @@
         ! We know it is an  "mlayer" kind of 1D array
         idim1d=idim3
+      ELSEIF (var_size==nlaykim_up) THEN
+        ! We know it is an  "preskim" kind of 1D array
+        idim1d=idim8
       ELSE 
         PRINT *, "put_var_rgen error : wrong dimension"

trunk/LMDZ.TITAN/libf/phytitan/phyredem.F90

@@ -9 +9 @@
                          alb,ith,pzmea,pzstd,pzsig,pzgam,pzthe)
 ! create physics restart file and write time-independent variables
+  use comchem_h, only: preskim
   use comsoil_h, only: volcapa, mlayer
   use geometry_mod, only: cell_area
@@ -104 +105 @@
   call put_var("soildepth","Soil mid-layer depth",mlayer)
   
+  ! Write the mid-layer upper chemistry pressure
+  call put_var("preskim","Upper chemistry mid-layer pressure",preskim)
+  
   ! Write longitudes
   call put_field("longitude","Longitudes of physics grid",lonfi)

trunk/LMDZ.TITAN/libf/phytitan/tabfi_mod.F90

@@ -55 +55 @@
                            emissiv
       use comsoil_h, only: volcapa
-      use comchem_h, only: nlaykim_up
       use iostart, only: get_var
       use mod_phys_lmdz_para, only: is_parallel
@@ -150 +149 @@
         dtemisice(:)=0 !time scale for snow metamorphism
         volcapa=1000000 ! volumetric heat capacity of subsurface
-! chemistry
-        nlaykim_up=70 ! size of vertical grid for upper chemistry 
         
       ELSE
@@ -207 +204 @@
 ! soil properties
       volcapa = tab_cntrl(tab0+35) ! volumetric heat capacity
-! chemistry
-      nlaykim_up = nint(tab_cntrl(tab0+40)) ! size of vertical grid for upper chemistry
 !-----------------------------------------------------------------------
 !       Save some constants for later use (as routine arguments)
@@ -261 +256 @@
 
       write(*,5) '(35)        volcapa',tab_cntrl(tab0+35),volcapa
-
-      write(*,5) '(40)     nlaykim_up',tab_cntrl(tab0+40),float(nlaykim_up)
 
       write(*,*)
@@ -571 +564 @@
  
       write(*,5) '(35)        volcapa',tab_cntrl(tab0+35),volcapa
-      
-      write(*,5) '(40)     nlaykim_up',tab_cntrl(tab0+40),float(nlaykim_up)
 
       write(*,*)