Changeset 1889 for trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/newstart.F

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

File:

• trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/newstart.F (modified) (7 diffs)

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   )