Changeset 1898

Timestamp:

Jan 26, 2018, 11:29:20 AM (7 years ago)

Author:

jvatant

Message:

Making chemistry handling more flexible - final step for all concerning startfi files
+ Added the vert_regrid_kim.F90 regriding upper fields and chemistry tracers if needed
in the newstart according to the new pressure grid
+ Fix some typos from lasts commit
( e.g. missing init of tracer_surf names in lect_start_archive )
--JVO

Location:

trunk/LMDZ.TITAN

Files:

• README (modified) (2 diffs)
• libf/dynphy_lonlat/phytitan/lect_start_archive.F (modified) (1 diff)
• libf/dynphy_lonlat/phytitan/newstart.F (modified) (4 diffs)
• libf/dynphy_lonlat/phytitan/vert_regrid_kim.F90 (added)
• utilities/prepare_startkim.bash (modified) (1 diff)

trunk/LMDZ.TITAN/README

@@ -1381 +1381 @@
 Big modifs of the tracer gestion/init in the physiq with a new query by names (see tracer_h )
 
-== 20/12/2017 - ... : r1871-86-87-91-94 ... == JVO
+== 12/2017 - 01/2018 : r1871-86-87-91-94-98 == JVO
 Management of the chemistry within startfi. Open the way to :
  1) run with chemistry in another resolution than 32x48 !
@@ -1396 +1396 @@
 with the introduction of gr_kim_vervack.F90 routine
 + For initialisation and loading of chemistry fields in phyetat0 -> chem_settings.F90
-+ ...
++ Upper chemistry fields ( and chemistry tracers if needed ) in newstart are
+ regrided in vert_regrid_kim.F90
 
 == 13/01/18 == JVO

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

@@ -729 +729 @@
 
       DO iq=1,nqtot
+        txt=trim(tname(iq))//"_surf"
         
         write(*,*) "lect_start_archive: loading tracer ",trim(txt)

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

@@ -266 +266 @@
 c=======================================================================
 
-! Initialize global tracer indexes (stored in tracer.h)
-! ... this has to be done before phyetat0
-      call initracer2(nqtot,tname)
-
 ! Take care of arrays in common modules
       ! ALLOCATE ARRAYS in surfdat_h (if not already done, e.g. when using start_archive)
@@ -337 +333 @@
      &                  aire,cu,cv,rad,g,r,cpp,
      &                  1)
+
+! Initialize global tracer indexes (stored in tracer.h)
+! ... this has to be done before phyetat0
+         CALL initracer2(nqtot,tname)
 
         ! Lmodif set to 0 to disable modifications possibility in phyeta0                           
@@ -458 +458 @@
      &                  aire,cu,cv,rad,g,r,cpp,
      &                  1)
+! Initialize global tracer indexes (stored in tracer.h)
+         CALL initracer2(nqtot,tname)
+
       ENDIF
 
@@ -1091 +1094 @@
           ENDIF
           
-          ! Regriding if needed
+          ! Regridind of upper chemistry fields and uppermost GCM
+          ! chemistry tracers if needed
           
           IF (nlaykimold.ne.nlaykim_up) THEN
 
             WRITE(*,*) " Warning, nlaykimold=", nlaykimold
-            WRITE(*,*) ' which implies that a vertical regriding on upper  
-     &chemistry fields will be performed.'
+            WRITE(*,*) ' which implies that a vertical regrid on upper  
+     &chemistry fields will be performed. If the ceiling is lowered I
+     & will also adjust GCM uppermost layers chem. tracers.'
             WRITE(*,*)
           
-!            CALL vert_regrid_kim(ngridmx)
+            CALL vert_regrid_kim(nqtot,q)
 
           ENDIF

trunk/LMDZ.TITAN/utilities/prepare_startkim.bash

@@ -119 +119 @@
         if [ $ilat -eq 1 ] ; then
             ngrid0=0
-            ngrid1=1
+            ngrid1=0
         elif [ $lat -eq 49 ] ; then
             ngrid0=1505