Changeset 3629 for LMDZ5/branches/IPSLCM5A2.1/libf/misc/wxios.F90

Timestamp:

Feb 10, 2020, 9:54:26 AM (4 years ago)

Author:

acozic

Message:

Add new grid, new axis and new variables for cmip protocole and dr2xml

• field_group id="coord_hyb"
• grid_ref="klevp1_bnds"
• grid_ref="klev_bnds"
• domain id="greordered"
• axis id="axis_lat"
• axis id="bnds"
• axis id="klevp1"
• axis id="klev"

File:

• LMDZ5/branches/IPSLCM5A2.1/libf/misc/wxios.F90 (modified) (4 diffs)

LMDZ5/branches/IPSLCM5A2.1/libf/misc/wxios.F90

@@ -342 +342 @@
         CALL xios_set_domain_attr_hdl(dom, nj_glo=nj_glo, jbegin=jbegin-1, nj=nj, data_dim=2)
         CALL xios_set_domain_attr_hdl(dom, lonvalue_1d=io_lon(ibegin:iend), latvalue_1d=io_lat(jbegin:jend))
+        IF (xios_is_valid_axis("axis_lat")) THEN
+           CALL xios_set_axis_attr( "axis_lat", n_glo=nj_glo, n=nj, begin=jbegin-1, value=io_lat(jbegin:jend))
+        ENDIF
 #endif
         IF (.NOT.is_sequential) THEN
@@ -372 +375 @@
     ! 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
@@ -379 +383 @@
         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
@@ -404 +410 @@
         CALL xios_set_axis_attr(trim(axis_id),size=axis_size,value=axis_value)
 #else
-        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
 #endif        
         !Vérification: