Changeset 3003 for LMDZ5/trunk/libf/misc/wxios.F90

Timestamp:

Oct 3, 2017, 4:45:26 PM (7 years ago)

Author:

Laurent Fairhead

Message:

Modifications to the code and xml files to output Ap and B, the coefficients
of the hybrid coordinates as requested by the CMIP6 DataRequest?
LF (with guidance from A. Caubel and S. Senesi)

File:

• LMDZ5/trunk/libf/misc/wxios.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/misc/wxios.F90

@@ -311 +311 @@
     ! Pour déclarer un axe vertical !!!!!!!!!!!!!!!
     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-    SUBROUTINE wxios_add_vaxis(axis_id, axis_size, axis_value)
+    SUBROUTINE wxios_add_vaxis(axis_id, axis_size, axis_value,     &
+                               positif, bnds)
         USE print_control_mod, ONLY : prt_level, lunout
         IMPLICIT NONE
@@ -318 +319 @@
         INTEGER, INTENT(IN) :: axis_size
         REAL, DIMENSION(axis_size), INTENT(IN) :: axis_value
+        CHARACTER (len=*), INTENT(IN), OPTIONAL :: positif
+        REAL, DIMENSION(axis_size, 2), INTENT(IN), OPTIONAL :: bnds
         
 !        TYPE(xios_axisgroup) :: axgroup
@@ -340 +343 @@
         
         ! Ehouarn: New way to declare axis, without axis_group:
-        CALL xios_set_axis_attr(trim(axis_id),n_glo=axis_size,value=axis_value)
+        if (PRESENT(positif) .AND. PRESENT(bnds)) then
+          CALL xios_set_axis_attr(trim(axis_id),n_glo=axis_size,value=axis_value, &
+                                  positive=positif, bounds=bnds)
+        else if (PRESENT(positif)) then
+          CALL xios_set_axis_attr(trim(axis_id),n_glo=axis_size,value=axis_value, &
+                                  positive=positif)
+        else if (PRESENT(bnds)) then
+          CALL xios_set_axis_attr(trim(axis_id),n_glo=axis_size,value=axis_value, &
+                                  bounds=bnds)
+        else
+          CALL xios_set_axis_attr(trim(axis_id),n_glo=axis_size,value=axis_value)
+        endif
 
         !Vérification: