Changeset 1895 for trunk/LMDZ.TITAN/libf

Timestamp:

Jan 16, 2018, 10:53:17 AM (7 years ago)

Author:

jvatant

Message:

Follow-up of r1894 - fix a missing init in phyetat0
--JVO

Location:

trunk/LMDZ.TITAN/libf

Files:

• dynphy_lonlat/phytitan/start2archive.F (modified) (1 diff)
• phytitan/chem_settings.F90 (modified) (3 diffs)

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

@@ -204 +204 @@
 ! Get value of the "upper_chemistry_layers" dimension from physics start file
 
+      ! NB : nlaykim_up is recalculated in phyetat0 ... this is redundant for newstart
+      ! If you want to do it properly put ykim_upS in a specific separate header
+      ! Anyway it should be safe. -- JVO 18
+
       ierr = NF_INQ_DIMID(nid1,"upper_chemistry_layers",varid)
       IF (ierr .NE. NF_NOERR) THEN

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

@@ -19 +19 @@
 USE comchem_h
 
-USE iostart, only: get_field, get_var
+USE iostart, only: get_field, get_var, inquire_dimension_length
+USE callkeys_mod, only : callchim
 
 IMPLICIT NONE
@@ -42 +43 @@
 !======================================================================
 
-  ! 0. Allocates arrays in comchem_h
+  ! 0. Start by reading how many layers of upper chemistry there are
+  
+  nlaykim_up=inquire_dimension_length("upper_chemistry_layers")
+  
+  ! 1. Allocates arrays in comchem_h
+  
   CALL ini_comchem_h(ngrid)
 
-  ! 1. Load upper chemistry pressure grid
+  ! 2. Load upper chemistry pressure grid
   
   CALL get_var("preskim",preskim,found)
@@ -54 +60 @@
                maxval(preskim), minval(preskim) 
  
-  ! 2. Load upper chemistry composition fields
+  ! 3. Inquire ( and load ) upper chemistry composition fields
  
   CALL get_field("H_up",ykim_up(1,:,:),found,indextime)
   IF (.NOT.found) THEN
+  
     ! We assume we can't do uncomplete chemistry
     WRITE(*,*) "chem_settings: No upper chemistry fields."
+    
+    IF ( callchim ) THEN ! if chemistry we must have the upper fields !
+      CALL abort_physic(modname,"Failed loading uppper chemistry fields, whereas callchim set to true !",1)
+    ENDIF
+    
+    DEALLOCATE(ykim_up) ! it will be useless
+    
   ELSE