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

Timestamp:

Dec 20, 2017, 12:23:42 PM (8 years ago)

Author:

jvatant

Message:

Making chemistry handling more flexible - Step 1
+ Enable upper_chemistry_dimension in the startfi files
+ Added a comchem_h for all the chemistry related stuff
+ start2archive adapted to these modifs : next step newstart !
--JVO

Location:

trunk/LMDZ.TITAN/libf/phytitan

Files:

• comchem_h.F90 (added)
• iostart.F90 (modified) (5 diffs)
• tabfi_mod.F90 (modified) (5 diffs)

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

@@ -15 +15 @@
     INTEGER,SAVE :: idim6 ! "nlayer" dimension
     INTEGER,SAVE :: idim7 ! "Time" dimension
+    INTEGER,SAVE :: idim8 ! "upper_chemistry_layers" dimension
     INTEGER,SAVE :: timeindex ! current time index (for time-dependent fields)
-!$OMP THREADPRIVATE(idim1,idim2,idim3,idim4,idim5,idim6,idim7,timeindex)
+!$OMP THREADPRIVATE(idim1,idim2,idim3,idim4,idim5,idim6,idim7,idim8,idim9,timeindex)
     INTEGER,PARAMETER :: length=100 ! size of tab_cntrl array
     
@@ -468 +469 @@
   USE tracer_h, only: nqtot_p
   USE comsoil_h, only: nsoilmx
+  USE comchem_h, only: nlaykim_up
 
   IMPLICIT NONE
@@ -556 +558 @@
       ENDIF
 
+      ierr=NF90_DEF_DIM(nid_restart,"upper_chemistry_layers",nlaykim_up,idim8)
+      IF (ierr/=NF90_NOERR) THEN
+        write(*,*)'open_restartphy: problem defining upper_chemistry_layers dimension '
+        write(*,*)trim(nf90_strerror(ierr))
+        CALL ABORT
+      ENDIF
+      
       ierr=NF90_ENDDEF(nid_restart)
       IF (ierr/=NF90_NOERR) THEN
@@ -634 +643 @@
   USE dimphy
   USE comsoil_h, only: nsoilmx
+  USE comchem_h, only: nlaykim_up
   USE mod_grid_phy_lmdz
   USE mod_phys_lmdz_para
@@ -807 +817 @@
             ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
                               (/idim2,idim3,idim7/),nvarid)
+#endif
+           if (ierr.ne.NF90_NOERR) then
+              write(*,*)"put_field_rgen error: failed to define "//trim(field_name)
+              write(*,*)trim(nf90_strerror(ierr))
+            endif
+            IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
+            ierr=NF90_ENDDEF(nid_restart)
+          endif
+          ! Write the variable
+          ierr=NF90_PUT_VAR(nid_restart,nvarid,field_glo,&
+                            start=(/1,1,timeindex/))
+
+        endif ! of if (.not.present(time))
+
+      ELSE IF (field_size==nlaykim_up) THEN
+        ! input is a 2D "upper chemistry" array
+        if (.not.present(time)) then ! for a time-independent field
+          ierr = NF90_REDEF(nid_restart)
+#ifdef NC_DOUBLE
+          ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
+                            (/idim2,idim8/),nvarid)
+#else
+          ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
+                            (/idim2,idim8/),nvarid)
+#endif
+          if (ierr.ne.NF90_NOERR) then
+            write(*,*)"put_field_rgen error: failed to define "//trim(field_name)
+            write(*,*)trim(nf90_strerror(ierr))
+          endif
+          IF (LEN_TRIM(title) > 0) ierr=NF90_PUT_ATT(nid_restart,nvarid,"title",title)
+          ierr = NF90_ENDDEF(nid_restart)
+          ierr = NF90_PUT_VAR(nid_restart,nvarid,field_glo)
+        else
+          ! check if the variable has already been defined:
+          ierr=NF90_INQ_VARID(nid_restart,field_name,nvarid)
+          if (ierr/=NF90_NOERR) then ! variable not found, define it
+            ierr=NF90_REDEF(nid_restart)
+#ifdef NC_DOUBLE
+            ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_DOUBLE,&
+                              (/idim2,idim8,idim7/),nvarid)
+#else
+            ierr=NF90_DEF_VAR(nid_restart,field_name,NF90_FLOAT,&
+                              (/idim2,idim8,idim7/),nvarid)
 #endif
            if (ierr.ne.NF90_NOERR) then

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

@@ -55 +55 @@
                            emissiv
       use comsoil_h, only: volcapa
+      use comchem_h, only: nlaykim_up
       use iostart, only: get_var
       use mod_phys_lmdz_para, only: is_parallel
@@ -149 +150 @@
         dtemisice(:)=0 !time scale for snow metamorphism
         volcapa=1000000 ! volumetric heat capacity of subsurface
-
+! chemistry
+        nlaykim_up=70 ! size of vertical grid for upper chemistry 
+        
       ELSE
 !-----------------------------------------------------------------------
@@ -204 +207 @@
 ! soil properties
       volcapa = tab_cntrl(tab0+35) ! volumetric heat capacity
+! chemistry
+      nlaykim_up = nint(tab_cntrl(tab0+40)) ! size of vertical grid for upper chemistry
 !-----------------------------------------------------------------------
 !       Save some constants for later use (as routine arguments)
@@ -256 +261 @@
 
       write(*,5) '(35)        volcapa',tab_cntrl(tab0+35),volcapa
+
+      write(*,5) '(40)     nlaykim_up',tab_cntrl(tab0+40),float(nlaykim_up)
 
       write(*,*)
@@ -564 +571 @@
  
       write(*,5) '(35)        volcapa',tab_cntrl(tab0+35),volcapa
+      
+      write(*,5) '(40)     nlaykim_up',tab_cntrl(tab0+40),float(nlaykim_up)
 
       write(*,*)