Changeset 1898 for trunk/LMDZ.TITAN/libf

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/libf/dynphy_lonlat/phytitan

Files:

• lect_start_archive.F (modified) (1 diff)
• newstart.F (modified) (4 diffs)
• vert_regrid_kim.F90 (added)

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