Changeset 1886 for trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan

Timestamp:

Jan 7, 2018, 5:35:22 PM (7 years ago)

Author:

jvatant

Message:

Some adjustments to r1871 on start_archive managment of chemistry dimension.
--JVO

Location:

trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan

Files:

• ini_archive.F (modified) (3 diffs)
• start2archive.F (modified) (5 diffs)
• write_archive.F (modified) (1 diff)

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

@@ -35 +35 @@
 c=======================================================================
 
+      USE comchem_h, ONLY: nlaykim_up
       USE comsoil_h
 !      use control_mod
@@ -80 +81 @@
       INTEGER idim_tim
       INTEGER idim_nsoilmx ! "subsurface_layers" dimension ID #
+      INTEGER idim_nlaykim_up ! "upper_chemistry_layers" dimension ID #
       INTEGER nid,nvarid
       real sig_s(llm),s(llm)
@@ -162 +164 @@
       ierr = NF_DEF_DIM (nid, "altitude", llm, idim_llm)
       ierr = NF_DEF_DIM (nid,"subsurface_layers",nsoilmx,idim_nsoilmx)
+      ierr = NF_DEF_DIM (nid,"upper_chemistry_layers",
+     .                        nlaykim_up,idim_nlaykim_up)
 
       ierr = NF_DEF_DIM (nid,"index", length, idim_index)

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

@@ -111 +111 @@
       INTEGER Lmodif
 
-      LOGICAL nokim
+      LOGICAL kim
 
       REAL ptotal
@@ -226 +226 @@
         PRINT*, "start2archive: Missing field(s) for upper chemistry ...
      . I presume they're all absent !"
-        nokim=.TRUE.
+        kim=.FALSE.
       ELSE
         PRINT*,"start2archive: I found a field for upper chemistry ...
@@ -232 +232 @@
         ! Allocates upper chemistry fields in comchem_h on physical and scalar grid
         CALL allokim_start2archive(ngridmx,ip1jmp1)
+        kim=.TRUE.
       ENDIF
 
@@ -366 +367 @@
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,tankCH4,tankCH4S)
 
-      IF (nokim .eqv. .FALSE.) THEN ! NB : fields are in comchem_h
+      IF (kim) THEN ! NB : fields are in comchem_h
          call gr_fi_dyn(nlaykim_up,ngridmx,iip1,jjp1,H,H_S)
          call gr_fi_dyn(nlaykim_up,ngridmx,iip1,jjp1,H2,H2_S)
@@ -560 +561 @@
 c-----------------------------------------------------------------
 
-      IF (nokim .eqv. .FALSE.) THEN
+      IF (kim) THEN
          call write_archive(nid,ntime,'H_up',
      .              'H in upper atmosphere','kg/kg',4,H_S)

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

@@ -175 +175 @@
           ierr=NF_INQ_DIMID(nid,"upper_chemistry_layers",id(3))
           if (ierr.ne.NF_NOERR) then
-           write(*,*)"write_archive: dimension 
-     &                <upper_chemistry_layers> is missing !!!"
+           write(*,*)"write_archive: dimension", 
+     &               " <upper_chemistry_layers> is missing !!!"
            call abort
           endif