Changeset 1887 for trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/start2archive.F

Timestamp:

Jan 8, 2018, 11:40:54 PM (7 years ago)

Author:

jvatant

Message:

For more convenience values of upper chemistry
pressure grid "preskim" is stored in the startfi too.
--JVO

File:

• trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/start2archive.F (modified) (2 diffs)

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

@@ -218 +218 @@
        CALL abort
       ENDIF
+
+      ALLOCATE(preskim(nlaykim_up))
       
-      ! Allocate arrays of nlaykim_up size, only if they're present
+      ! Allocate other arrays of nlaykim_up size, only if they're present
       ! The test is on HCN but could be on any as we assume we can't do incomplete chemistry
 
@@ -286 +288 @@
        IF (ierr .NE. NF_NOERR) THEN
           PRINT*, "start2archive: Lecture echoue pour <controle>"
+          CALL abort
+       ENDIF
+
+! load upper chemistry pressure grid from physics start file
+
+      ierr = NF_INQ_VARID (nid1, "preskim", varid)
+      IF (ierr .NE. NF_NOERR) THEN
+       PRINT*, "start2archive: Le champ <preskim> est absent"
+       CALL abort
+      ENDIF
+#ifdef NC_DOUBLE
+       ierr = NF_GET_VAR_DOUBLE(nid1, varid, preskim)
+#else
+      ierr = NF_GET_VAR_REAL(nid1, varid, preskim)
+#endif
+       IF (ierr .NE. NF_NOERR) THEN
+          PRINT*, "start2archive: Lecture echoue pour <preskim>"
           CALL abort
        ENDIF