Changeset 1887 for trunk/LMDZ.TITAN/libf

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

Location:

trunk/LMDZ.TITAN/libf

Files:

• dynphy_lonlat/phytitan/ini_archive.F (modified) (2 diffs)
• dynphy_lonlat/phytitan/start2archive.F (modified) (2 diffs)
• phytitan/comchem_h.F90 (modified) (1 diff)

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

@@ -35 +35 @@
 c=======================================================================
 
-      USE comchem_h, ONLY: nlaykim_up
+      USE comchem_h, ONLY: nlaykim_up, preskim
       USE comsoil_h
 !      use control_mod
@@ -444 +444 @@
 #endif
 
+!---------------------------------------------------------------------
+! Upper chemistry mid-layer pressure preskim() (known from comchem.h)
+!---------------------------------------------------------------------
+      ierr=NF_REDEF (nid) ! Enter NetCDF (re-)define mode
+      ! define variable
+#ifdef NC_DOUBLE
+      ierr=NF_DEF_VAR(nid,"preskim",NF_DOUBLE,1,idim_nlaykim_up,nvarid)
+#else
+      ierr=NF_DEF_VAR(nid,"preskim",NF_FLOAT,1,idim_nlaykim_up,nvarid)
+#endif
+      ierr=NF_PUT_ATT_TEXT (nid,nvarid,"long_name", 34,
+     .                        "Upper chemistry mid-layer pressure")
+      ierr=NF_PUT_ATT_TEXT (nid,nvarid,"units",1,"Pa")
+      ierr=NF_PUT_ATT_TEXT (nid,nvarid,"positive",2,"up")
+      ierr=NF_ENDDEF(nid) ! Leave NetCDF define mode
+      ! write variable
+#ifdef NC_DOUBLE
+      ierr=NF_PUT_VAR_DOUBLE (nid,nvarid,preskim)
+#else
+      ierr=NF_PUT_VAR_REAL (nid,nvarid,preskim)
+#endif
+
 c-----------------------------------------------------------------------
 c  Ecriture aire et coefficients de passage cov. <-> contra. <--> naturel

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

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

@@ -13 +13 @@
 ! NB : For the startfile we use nlaykim_up grid (upper atm) and for outputs we use nlaykim_tot grid (all layers)
 
-   REAL,SAVE,ALLOCATABLE,DIMENSION(:,:) :: zlay_kim  ! Altitude (km) of all chemistry layers
+   REAL,SAVE,ALLOCATABLE,DIMENSION(:) :: preskim  ! Pressure (Pa) of upper chemistry (mid)-layers
+!$OMP_THREADPRIVATE(preskim)
+
+   REAL,SAVE,ALLOCATABLE,DIMENSION(:,:) :: zlay_kim  ! Altitude (km) of all chemistry (mid)-layers
 !$OMP_THREADPRIVATE(zlay_kim)