Changeset 1895 for trunk/LMDZ.TITAN/libf/phytitan

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

File:

• trunk/LMDZ.TITAN/libf/phytitan/chem_settings.F90 (modified) (3 diffs)

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