Changeset 2002


Ignore:
Timestamp:
Apr 4, 2014, 2:39:54 PM (11 years ago)
Author:
Ehouarn Millour
Message:

Further cleanup concerning XIOS (mainly about axes being defined as axes and not as groups of axes).
EM

Location:
LMDZ5/trunk
Files:
15 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/DefLists/field_def_lmdz.xml

    r2001 r2002  
    358358    </field_group>
    359359
    360     <field_group id="fields_3D" domain_ref="dom_glo">
     360    <field_group id="fields_3D" domain_ref="dom_glo" axis_ref="presnivs">
    361361        <field id="cldtau"    long_name="Cloud optical thickness"    unit="1" />
    362362        <field id="cldemi"    long_name="Cloud optical emissivity"    unit="1" />

  • LMDZ5/trunk/DefLists/file_def_histLES_lmdz.xml

    r2001 r2002  
    22    <file_group id="defile">
    33        <file id="histLES" name="XhistLES" output_freq="1h" output_level="5" enabled=".FALSE.">
     4           
     5            <!-- VARS 1D -->
     6            <field_group operation="once">
     7              <field field_ref="Ap" level="1" />
     8              <field field_ref="Bp" level="1" />
     9              <field field_ref="Alt" level="1" />
     10            </field_group>
     11       
    412            <!-- VARS 2D -->
    513            <field_group operation="average" freq_op="1ts">
    6                 <field field_ref="Ap" level="1" />
    7                 <field field_ref="Bp" level="1" />
    8                 <field field_ref="Alt" level="1" />
    914                <field field_ref="phis" level="1" />
    1015                <field field_ref="aire" level="1" operation="once" />
    … …  
    2328                <field field_ref="t2m_sic" level="10" />
    2429                <field field_ref="wind10m" level="10" />
    25                 <field field_ref="wind10max" level="10" />
     30                <field field_ref="wind10max" level="10" operation="maximum" />
    2631                <field field_ref="sicf" level="10" />
    2732                <field field_ref="q2m" level="10" />
    … …  
    331336                <field field_ref="lcc" level="10" />
    332337                <field field_ref="wvapp" level="10" />
     338                <field field_ref="ozone_daylight" level="10" />
    333339                <field field_ref="albe_ter" level="10" />
    334340                <field field_ref="albe_lic" level="10" />
    … …  
    361367
    362368            <!-- VARS 3D -->
    363             <field_group operation="average" freq_op="1ts" axis_ref="presnivs_histLES">
     369            <field_group operation="average" freq_op="1ts" axis_ref="presnivs">
    364370                <field field_ref="cldtau" level="10" />
    365371                <field field_ref="cldemi" level="10" />
    … …  
    394400                <field field_ref="rhum" level="10" />
    395401                <field field_ref="ozone" level="10" />
    396                 <field field_ref="ozone_daylight" level="10" />
    397402                <field field_ref="upwd" level="10" />
    398403                <field field_ref="dtphy" level="10" />

  • LMDZ5/trunk/DefLists/file_def_histday_lmdz.xml

    r2001 r2002  
    22    <file_group id="defile">
    33        <file id="histday" name="Xhistday" output_freq="1d" output_level="5" enabled=".FALSE.">
     4           
     5            <!-- VARS 1D -->
     6            <field_group operation="once">
     7              <field field_ref="Ap" level="1" />
     8              <field field_ref="Bp" level="1" />
     9              <field field_ref="Alt" level="1" />
     10            </field_group>
     11       
    412            <!-- VARS 2D -->
    513            <field_group operation="average" freq_op="1ts">
    6                 <field field_ref="Ap" level="1" />
    7                 <field field_ref="Bp" level="1" />
    8                 <field field_ref="Alt" level="1" />
    914                <field field_ref="phis" level="1" />
    10                 <field field_ref="aire" level="1"        operation="once" />
     15                <field field_ref="aire" level="1" operation="once" />
    1116                <field field_ref="contfracATM" level="1" operation="once" />
    1217                <field field_ref="contfracOR" level="1" />
    … …  
    2328                <field field_ref="t2m_sic" level="6" />
    2429                <field field_ref="wind10m" level="1" />
    25                 <field field_ref="wind10max" level="1" />
     30                <field field_ref="wind10max" level="1" operation="maximum" />
    2631                <field field_ref="sicf" level="1" />
    2732                <field field_ref="q2m" level="1" />
    … …  
    331336                <field field_ref="lcc" level="6" />
    332337                <field field_ref="wvapp" level="10" />
     338                <field field_ref="ozone_daylight" level="10" />
    333339                <field field_ref="albe_ter" level="7" />
    334340                <field field_ref="albe_lic" level="7" />
    … …  
    361367
    362368            <!-- VARS 3D -->
    363             <field_group operation="average" freq_op="1ts" axis_ref="presnivs_histday">
     369            <field_group operation="average" freq_op="1ts" axis_ref="presnivs">
    364370                <field field_ref="cldtau" level="5" />
    365371                <field field_ref="cldemi" level="5" />
    … …  
    394400                <field field_ref="rhum" level="5" />
    395401                <field field_ref="ozone" level="10" />
    396                 <field field_ref="ozone_daylight" level="10" />
    397402                <field field_ref="upwd" level="10" />
    398403                <field field_ref="dtphy" level="10" />

  • LMDZ5/trunk/DefLists/file_def_histhf_lmdz.xml

    r2001 r2002  
    22    <file_group id="defile">
    33        <file id="histhf" name="Xhisthf" output_freq="1h" output_level="5" enabled=".FALSE.">
     4           
     5            <!-- VARS 1D -->
     6            <field_group operation="once">
     7              <field field_ref="Ap" level="1" />
     8              <field field_ref="Bp" level="1" />
     9              <field field_ref="Alt" level="1" />
     10            </field_group>
     11       
    412            <!-- VARS 2D -->
    513            <field_group operation="average" freq_op="1ts">
    6                 <field field_ref="Ap" level="1" />
    7                 <field field_ref="Bp" level="1" />
    8                 <field field_ref="Alt" level="1" />
    914                <field field_ref="phis" level="5" />
    1015                <field field_ref="aire" level="10" operation="once" />
    … …  
    2328                <field field_ref="t2m_sic" level="10" />
    2429                <field field_ref="wind10m" level="10" />
    25                 <field field_ref="wind10max" level="10" />
     30                <field field_ref="wind10max" level="10" operation="maximum" />
    2631                <field field_ref="sicf" level="10" />
    2732                <field field_ref="q2m" level="5" />
    … …  
    331336                <field field_ref="lcc" level="10" />
    332337                <field field_ref="wvapp" level="10" />
     338                <field field_ref="ozone_daylight" level="10" />
    333339                <field field_ref="albe_ter" level="7" />
    334340                <field field_ref="albe_lic" level="7" />
    … …  
    361367
    362368            <!-- VARS 3D -->
    363             <field_group operation="average" freq_op="1ts" axis_ref="presnivs_histhf">
     369            <field_group operation="average" freq_op="1ts" axis_ref="presnivs">
    364370                <field field_ref="cldtau" level="10" />
    365371                <field field_ref="cldemi" level="10" />
    … …  
    394400                <field field_ref="rhum" level="10" />
    395401                <field field_ref="ozone" level="10" />
    396                 <field field_ref="ozone_daylight" level="10" />
    397402                <field field_ref="upwd" level="10" />
    398403                <field field_ref="dtphy" level="10" />

  • LMDZ5/trunk/DefLists/file_def_histins_lmdz.xml

    r2001 r2002  
    22    <file_group id="defile">
    33        <file id="histins" name="Xhistins" output_freq="6h" output_level="4" enabled=".FALSE.">
     4           
     5            <!-- VARS 1D -->
     6            <field_group operation="once">
     7              <field field_ref="Ap" level="1" />
     8              <field field_ref="Bp" level="1" />
     9              <field field_ref="Alt" level="1" />
     10            </field_group>
     11       
    412            <!-- VARS 2D -->
    513            <field_group operation="instant" freq_op="1ts">
    6                 <field field_ref="Ap" level="1" />
    7                 <field field_ref="Bp" level="1" />
    8                 <field field_ref="Alt" level="1" />
    914                <field field_ref="phis" level="10" />
    1015                <field field_ref="aire" level="10" operation="once" />
    … …  
    2328                <field field_ref="t2m_sic" level="10" />
    2429                <field field_ref="wind10m" level="1" />
    25                 <field field_ref="wind10max" level="10" />
     30                <field field_ref="wind10max" level="10" operation="maximum" />
    2631                <field field_ref="sicf" level="10" />
    2732                <field field_ref="q2m" level="1" />
    … …  
    331336                <field field_ref="lcc" level="10" />
    332337                <field field_ref="wvapp" level="10" />
     338                <field field_ref="ozone_daylight" level="10" />
    333339                <field field_ref="albe_ter" level="10" />
    334340                <field field_ref="albe_lic" level="10" />
    … …  
    361367
    362368            <!-- VARS 3D -->
    363             <field_group operation="instant" freq_op="1ts" axis_ref="presnivs_histins">
     369            <field_group operation="instant" freq_op="1ts" axis_ref="presnivs">
    364370                <field field_ref="cldtau" level="10" />
    365371                <field field_ref="cldemi" level="10" />
    … …  
    394400                <field field_ref="rhum" level="10" />
    395401                <field field_ref="ozone" level="10" />
    396                 <field field_ref="ozone_daylight" level="10" />
    397402                <field field_ref="upwd" level="10" />
    398403                <field field_ref="dtphy" level="10" />
    … …  
    428433                <field field_ref="wake_deltaq" level="10" />
    429434                <field field_ref="wake_omg" level="10" />
    430                 <field field_ref="wdtrainA" level="10" />
    431                 <field field_ref="wdtrainM" level="10" />
     435                <field field_ref="wdtrainA" level="4" />
     436                <field field_ref="wdtrainM" level="4" />
    432437                <field field_ref="Vprecip" level="10" />
    433438                <field field_ref="ftd" level="10" />

  • LMDZ5/trunk/DefLists/file_def_histmth_lmdz.xml

    r2001 r2002  
    11<file_definition>
    22        <file_group id="defile">
     3        <file id="histmth" name="Xhistmth" output_freq="5d" output_level="10"  enabled=".FALSE.">
    34
    4         <file id="histmth" name="Xhistmth" output_freq="1d" output_level="10"  enabled=".FALSE.">
    5 
     5            <!-- VARS 1D -->
     6            <field_group operation="once">
     7              <field field_ref="Ap" level="1" />
     8              <field field_ref="Bp" level="1" />
     9              <field field_ref="Alt" level="1" />
     10            </field_group>
     11       
    612            <!-- VARS 2D -->
    713            <field_group operation="average" freq_op="1ts">
    8                 <field field_ref="Ap"           level="1" />
    9                 <field field_ref="Bp"           level="1" />
    10                 <field field_ref="Alt"          level="1" />
    1114                <field field_ref="phis"         level="1" />
    1215                <field field_ref="aire"         level="1"   operation="once"/>
    … …  
    2528                <field field_ref="t2m_sic" level="10" />
    2629                <field field_ref="wind10m" level="1" />
    27                 <field field_ref="wind10max" level="10" />
     30                <field field_ref="wind10max" level="10" operation="maximum" />
    2831                <field field_ref="sicf"      level="1" />
    2932                <field field_ref="q2m"       level="1" />
    … …  
    333336                <field field_ref="lcc" level="2" />
    334337                <field field_ref="wvapp" level="2" />
     338                <field field_ref="ozone_daylight" level="2" />
    335339                <field field_ref="albe_ter" level="3" />
    336340                <field field_ref="albe_lic" level="3" />
    … …  
    363367
    364368            <!-- VARS 3D -->
    365             <field_group operation="average" freq_op="1ts" axis_ref="presnivs_histmth">
     369            <field_group operation="average" freq_op="1ts" axis_ref="presnivs">
    366370                <field field_ref="cldtau" level="10" />
    367371                <field field_ref="cldemi" level="10" />
    … …  
    396400                <field field_ref="rhum" level="2" />
    397401                <field field_ref="ozone" level="2" />
    398                 <field field_ref="ozone_daylight" level="2" />
    399402                <field field_ref="upwd" level="2" />
    400403                <field field_ref="dtphy" level="2" />

  • LMDZ5/trunk/DefLists/file_def_histstn_lmdz.xml

    r2001 r2002  
    22    <file_group id="defile">
    33        <file id="histstn" name="Xhiststn" output_freq="1h" output_level="5" enabled=".FALSE.">
     4
     5            <!-- VARS 1D -->
     6            <field_group operation="once">
     7              <field field_ref="Ap" level="1" />
     8              <field field_ref="Bp" level="1" />
     9              <field field_ref="Alt" level="1" />
     10            </field_group>
     11       
    412            <!-- VARS 2D -->
    513            <field_group operation="average" freq_op="1ts">
    6                 <field field_ref="Ap" level="1" />
    7                 <field field_ref="Bp" level="1" />
    8                 <field field_ref="Alt" level="1" />
    914                <field field_ref="phis" level="1" />
    1015                <field field_ref="aire" level="1" operation="once" />
    … …  
    331336                <field field_ref="lcc" level="10" />
    332337                <field field_ref="wvapp" level="10" />
     338                <field field_ref="ozone_daylight" level="10" />
    333339                <field field_ref="albe_ter" level="10" />
    334340                <field field_ref="albe_lic" level="10" />
    … …  
    361367
    362368            <!-- VARS 3D -->
    363             <field_group operation="average" freq_op="1ts" axis_ref="presnivs_histstn">
     369            <field_group operation="average" freq_op="1ts" axis_ref="presnivs">
    364370                <field field_ref="cldtau" level="10" />
    365371                <field field_ref="cldemi" level="10" />
    … …  
    394400                <field field_ref="rhum" level="10" />
    395401                <field field_ref="ozone" level="10" />
    396                 <field field_ref="ozone_daylight" level="10" />
    397402                <field field_ref="upwd" level="10" />
    398403                <field field_ref="dtphy" level="10" />

  • LMDZ5/trunk/DefLists/iodef.xml

    r2001 r2002  
    4040        </domain_definition>
    4141               
    42         <!-- Need to define a grid? Do it here
     42        <!-- Need to define a grid? Do it here -->
    4343        <grid_definition>
    4444          <grid_group id="vertical" axis_ref="presnivs" />
    4545        </grid_definition>
    46         -->
    4746               
    4847        <!-- Define groups of vertical axes -->
    4948        <axis_definition>
    50                 <axis_group id="presnivs" standard_name="Vertical levels" unit="Pa">
    51                 </axis_group>
    52                 <axis_group id="Ahyb" standard_name="Ahyb comp of Hyb Cord" unit="Pa">
    53                 </axis_group>
    54                 <axis_group id="Bhyb" standard_name="Bhyb comp of Hyb Cord" unit="">
    55                 </axis_group>
    56                 <axis_group id="Alt" standard_name="Height approx for scale heigh of 8km at levels" unit="km">
    57                 </axis_group>
    58                 <axis_group id="plev" standard_name="model_level_number" unit="Pa">
    59                 </axis_group>
     49                <axis id="presnivs" standard_name="Vertical levels" unit="Pa">
     50                </axis>
     51                <axis id="Ahyb" standard_name="Ahyb comp of Hyb Cord" unit="Pa">
     52                </axis>
     53                <axis id="Bhyb" standard_name="Bhyb comp of Hyb Cord" unit="">
     54                </axis>
     55                <axis id="Alt" standard_name="Height approx for scale heigh of 8km at levels" unit="km">
     56                </axis>
     57                <axis id="plev" standard_name="model_level_number" unit="Pa">
     58                </axis>
    6059        </axis_definition>
    6160        </context>

  • LMDZ5/trunk/DefLists/iodef_dev.xml

    r1907 r2002  
    22<simulation>
    33    <context id="LMDZ" calendar_type="D360" start_date="1980-01-01 00:00:00">
    4         <!-- Définition des variables -->
    5         <field_definition prec="4" operation="average" freq_op="1ts" enabled=".TRUE.">
     4        <!-- Definition of model variables -->
     5        <field_definition prec="4"
     6                          freq_op="1ts"
     7                          enabled=".TRUE.">
     8           
     9            <!-- 2D variables -->
    610            <field_group id="fields_2D" domain_ref="dom_glo">
    7                 <field id="u"           long_name="Eastward Zonal Wind"                 unit="m/s" />
    8                 <field id="v"           long_name="Northward Meridional Wind"           unit="m/s" />
    9                 <field id="ps"          long_name="Surface Pressure"                    unit="m/s" />
     11                <field id="ps"
     12                       long_name="Surface Pressure"
     13                       unit="m/s" />
    1014            </field_group>
    1115
    12             <field_group id="fields_3D" domain_ref="dom_glo">
    13                 <field id="temperature" long_name="Atmospheric temperature"             unit="K" />
     16            <!-- 3D variables -->
     17            <field_group id="fields_3D"
     18                         domain_ref="dom_glo">
     19                <field id="temperature"
     20                       long_name="Atmospheric temperature"
     21                       unit="K" />
     22                <field id="temp_newton"
     23                       long_name="Relaxation temperature"
     24                       unit="K" />
     25                <field id="u"
     26                       long_name="Eastward Zonal Wind"
     27                       unit="m/s" />
     28                <field id="v"
     29                       long_name="Northward Meridional Wind"
     30                       unit="m/s" />
    1431            </field_group>
    1532        </field_definition>
    1633
    17         <!-- Définition des fichiers de sortie -->
     34        <!-- Définition of output files -->
    1835
    1936        <file_definition>
    20             <file_group id="defile">
    21                 <file id="histins" name="Xhistins" output_freq="1ts" enabled=".TRUE.">
     37            <!-- <file_group id="defile"> -->
     38                <file id="histins"
     39                      name="Xhistins"
     40                      output_freq="12ts"
     41                      enabled=".TRUE."> <!-- NB: output_freq in physics ts-->
     42                   
    2243                    <!-- VARS 2D -->
    23                     <field_group operation="instant" freq_op="1ts">
     44                    <field_group operation="instant"
     45                                 freq_op="1ts">
    2446                        <field field_ref="ps" />
    2547                    </field_group>
     48                    <!-- <field_group field_group_ref="fields_2D"
     49                                 operation="instant" /> doesn't work ?!? -->
    2650
    2751                    <!-- VARS 3D -->
    28                     <field_group operation="instant" freq_op="1ts" axis_ref="presnivs_histins">
     52                    <field_group operation="instant"
     53                                 freq_op="1ts"
     54                                 axis_ref="presnivs">
     55                        <field field_ref="temperature" />
     56                        <field field_ref="temp_newton" operation="once" />
     57                        <field field_ref="u" />
     58                        <field field_ref="v" />
     59                    </field_group>
     60                    <!-- <field_group field_group_ref="fields_3D"
     61                                 operation="instant" /> doesn't work ?!? -->
     62                </file>
     63
     64                <file id="diurnalave"
     65                      name="Xdiurnalave"
     66                      output_freq="1d"
     67                      enabled=".true.">
     68                    <!-- VARS 2D -->
     69                    <field_group operation="average"
     70                                 freq_op="1ts">
     71                        <field field_ref="ps" />
     72                    </field_group>
     73                    <!-- VARS 3D -->
     74                    <field_group operation="average"
     75                                 freq_op="1ts"
     76                                 axis_ref="presnivs">
    2977                        <field field_ref="temperature" />
    3078                        <field field_ref="u" />
    3179                        <field field_ref="v" />
    32                     </field_group>
     80                    </field_group>                     
    3381                </file>
    34             </file_group>
     82            <!-- </file_group> -->
    3583        </file_definition>
    3684
    3785
    38 
    39 
    40         <!-- Définition des domaines et groupes de domaines -->
     86        <!-- Definition of domains and groups of domains -->
    4187        <domain_definition>
    4288            <domain id="dom_glo" data_dim="2" />
    4389        </domain_definition>
    4490       
    45         <!-- Définition des axes verticaux -->
     91        <!-- Definition of vertical axes
     92             NB: these must be defined and set in code via
     93             wxios_add_vaxis -->
    4694        <axis_definition>
    47             <axis_group id="presnivs" standard_name="Vertical levels" unit="Pa">
     95            <!--
     96            <axis_group id="presnivs"
     97                        standard_name="Pseudo-pressure of model vertical levels"
     98                        unit="Pa">
    4899            </axis_group>
    49             <axis_group id="Ahyb" standard_name="Ahyb comp of Hyb Cord" unit="Pa">
    50             </axis_group>
    51             <axis_group id="Bhyb" standard_name="Bhyb comp of Hyb Cord" unit="">
    52             </axis_group>
    53             <axis_group id="Alt" standard_name="Height approx for scale heigh of 8km at levels" unit="Km">
    54             </axis_group>
     100            -->
     101            <axis id="presnivs"
     102                        standard_name="Pseudo-pressure of model vertical levels"
     103                        unit="Pa">
     104            </axis>
    55105        </axis_definition>
    56106    </context>
    … …  
    65115            <variable_group id="parameters" >
    66116                <variable id="using_server" type="boolean">true</variable>
    67                 <variable id="info_level" type="int">100</variable>
     117                <variable id="info_level" type="int">10</variable>
    68118            </variable_group>
    69119        </variable_definition>

  • LMDZ5/trunk/libf/bibio/wxios.F90

    r2001 r2002  
    292292    ! Pour déclarer un axe vertical !!!!!!!!!!!!!!!
    293293    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    294     SUBROUTINE wxios_add_vaxis(axisgroup_id, axis_file, axis_size, axis_value)
    295         IMPLICIT NONE
    296         INCLUDE 'iniprint.h'
    297 
    298         CHARACTER (len=*), INTENT(IN) :: axisgroup_id, axis_file
     294    SUBROUTINE wxios_add_vaxis(axis_id, axis_size, axis_value)
     295        IMPLICIT NONE
     296        INCLUDE 'iniprint.h'
     297
     298        CHARACTER (len=*), INTENT(IN) :: axis_id
    299299        INTEGER, INTENT(IN) :: axis_size
    300300        REAL, DIMENSION(axis_size), INTENT(IN) :: axis_value
    301301       
    302         TYPE(xios_axisgroup) :: axgroup
    303         TYPE(xios_axis) :: ax
    304         CHARACTER(len=50) :: axis_id
    305        
    306         IF (len_trim(axisgroup_id).gt.len(axis_id)) THEN
    307           WRITE(lunout,*) "wxios_add_vaxis: error, size of axis_id too small!!"
    308           WRITE(lunout,*) "     increase it to at least ",len_trim(axisgroup_id)
    309           CALL abort_gcm("wxios_add_vaxis","len(axis_id) too small",1)
    310         ENDIF
    311         axis_id=trim(axisgroup_id)
     302!        TYPE(xios_axisgroup) :: axgroup
     303!        TYPE(xios_axis) :: ax
     304!        CHARACTER(len=50) :: axis_id
     305       
     306!        IF (len_trim(axisgroup_id).gt.len(axis_id)) THEN
     307!          WRITE(lunout,*) "wxios_add_vaxis: error, size of axis_id too small!!"
     308!          WRITE(lunout,*) "     increase it to at least ",len_trim(axisgroup_id)
     309!          CALL abort_gcm("wxios_add_vaxis","len(axis_id) too small",1)
     310!        ENDIF
     311!        axis_id=trim(axisgroup_id)
    312312       
    313313        !On récupère le groupe d'axes qui va bien:
    314         CALL xios_get_axisgroup_handle(axisgroup_id, axgroup)
     314        !CALL xios_get_axisgroup_handle(axisgroup_id, axgroup)
    315315       
    316316        !On ajoute l'axe correspondant à ce fichier:
    317         CALL xios_add_axis(axgroup, ax, TRIM(ADJUSTL(axis_id)))
     317        !CALL xios_add_axis(axgroup, ax, TRIM(ADJUSTL(axis_id)))
    318318       
    319319        !Et on le parametrise:
    320         CALL xios_set_axis_attr_hdl(ax, size=axis_size, value=axis_value)
     320        !CALL xios_set_axis_attr_hdl(ax, size=axis_size, value=axis_value)
     321       
     322        ! Ehouarn: New way to declare axis, without axis_group:
     323        CALL xios_set_axis_attr(trim(axis_id),size=axis_size,value=axis_value)
    321324       
    322325        !Vérification:
    … …  
    515518            !Sinon on se contente de l'activer:
    516519            CALL xios_set_field_attr(fieldname, enabled=.TRUE.)
     520            !NB: This will override an enable=.false. set by a user in the xml file;
     521            !   then the only way to not output the field is by changing its
     522            !   output level
    517523        ENDIF       
    518524       
    519525    END SUBROUTINE wxios_add_field_to_file
    520526   
    521     SUBROUTINE wxios_update_calendar(ito)
    522         INTEGER, INTENT(IN) :: ito
    523         CALL xios_update_calendar(ito)
    524     END SUBROUTINE wxios_update_calendar
    525    
    526     SUBROUTINE wxios_write_2D(fieldname, fdata)
    527         CHARACTER(len=*), INTENT(IN) :: fieldname
    528         REAL, DIMENSION(:,:), INTENT(IN) :: fdata
    529 
    530         CALL xios_send_field(fieldname, fdata)
    531     END SUBROUTINE wxios_write_2D
    532    
    533     SUBROUTINE wxios_write_3D(fieldname, fdata)
    534         CHARACTER(len=*), INTENT(IN) :: fieldname
    535         REAL, DIMENSION(:,:,:), INTENT(IN) :: fdata
    536        
    537         CALL xios_send_field(fieldname, fdata)
    538     END SUBROUTINE wxios_write_3D
     527!    SUBROUTINE wxios_update_calendar(ito)
     528!        INTEGER, INTENT(IN) :: ito
     529!        CALL xios_update_calendar(ito)
     530!    END SUBROUTINE wxios_update_calendar
     531!   
     532!    SUBROUTINE wxios_write_2D(fieldname, fdata)
     533!        CHARACTER(len=*), INTENT(IN) :: fieldname
     534!        REAL, DIMENSION(:,:), INTENT(IN) :: fdata
     535!
     536!        CALL xios_send_field(fieldname, fdata)
     537!    END SUBROUTINE wxios_write_2D
     538   
     539!    SUBROUTINE wxios_write_3D(fieldname, fdata)
     540!        CHARACTER(len=*), INTENT(IN) :: fieldname
     541!        REAL, DIMENSION(:,:,:), INTENT(IN) :: fdata
     542!       
     543!        CALL xios_send_field(fieldname, fdata)
     544!    END SUBROUTINE wxios_write_3D
    539545   
    540546    SUBROUTINE wxios_closedef()

  • LMDZ5/trunk/libf/phydev/iophy.F90

    r1907 r2002  
    340340  USE mod_phys_lmdz_para, only: gather_omp, grid1Dto2D_mpi, &
    341341                                jj_nb, klon_mpi
    342   USE wxios, only: wxios_write_2D
     342  USE xios, only: xios_send_field
    343343
    344344
    … …  
    361361    CALL grid1Dto2D_mpi(buffer_omp,Field2d)
    362362   
    363     CALL wxios_write_2D(field_name, Field2d)
     363    CALL xios_send_field(field_name, Field2d)
    364364!$OMP END MASTER   
    365365
    … …  
    376376  USE mod_phys_lmdz_para, only: gather_omp, grid1Dto2D_mpi, &
    377377                                jj_nb, klon_mpi
    378   USE wxios, only: wxios_write_3D
     378  USE xios, only: xios_send_field
    379379
    380380
    … …  
    401401    CALL grid1Dto2D_mpi(buffer_omp,field3d)
    402402
    403     CALL wxios_write_3D(field_name, Field3d(:,:,1:klev))
     403    CALL xios_send_field(field_name, Field3d(:,:,1:nlev))
    404404!$OMP END MASTER   
    405405

  • LMDZ5/trunk/libf/phydev/physiq.F90

    r1907 r2002  
    2121
    2222#ifdef CPP_XIOS
     23      USE xios, ONLY: xios_update_calendar
    2324      USE wxios, only: wxios_add_vaxis, wxios_set_timestep, wxios_closedef, &
    24                        wxios_update_calendar, histwrite_phy
     25                       histwrite_phy
    2526#endif
    2627
    … …  
    137138#ifdef CPP_XIOS
    138139!XIOS
    139     ! Déclaration de l'axe vertical du fichier:   
    140     CALL wxios_add_vaxis("presnivs", "histins", klev, presnivs)
    141 
    142     !Déclaration du pas de temps:
     140    ! Declare available vertical axes to be used in output files:   
     141    !CALL wxios_add_vaxis("presnivs", "dummy-not-used", klev, presnivs)
     142    CALL wxios_add_vaxis("presnivs", klev, presnivs)
     143
     144    ! Declare time step length (in s):
    143145    CALL wxios_set_timestep(dtime)
    144146
    145     !Finalisation du contexte:
     147    !Finalize the context:
    146148    CALL wxios_closedef()
    147149#endif
    … …  
    187189!$OMP MASTER
    188190    !Increment XIOS time
    189     CALL wxios_update_calendar(itau)
     191    CALL xios_update_calendar(itau)
    190192!$OMP END MASTER
    191193!$OMP BARRIER
    192194
    193     !Send fields to XIOS:
     195    !Send fields to XIOS: (NB these fields must also be defined as
     196    ! <field id="..." /> in iodef.xml to be correctly used
    194197    CALL histwrite_phy("temperature",t)
     198    CALL histwrite_phy("temp_newton",temp_newton)
    195199    CALL histwrite_phy("u",u)
    196200    CALL histwrite_phy("v",v)

  • LMDZ5/trunk/libf/phylmd/iophy.F90

    r2001 r2002  
    883883                                 nid_files
    884884#ifdef CPP_XIOS
    885   USE wxios, only: wxios_write_2D
     885  USE xios, only: xios_send_field
    886886#endif
    887887
    … …  
    952952                     write(lunout,*)'Dans iophy histwrite2D,iff,var%name ',&
    953953                                    iff,trim(var%name)                       
    954                      write(lunout,*)"histwrite2d_phy:.NOT.clef_stations(iff)and iff==iff_beg, call wxios_write_2D"
     954                     write(lunout,*)"histwrite2d_phy:.NOT.clef_stations(iff)and iff==iff_beg, call xios_send_field"
    955955                  endif
    956                   CALL wxios_write_2D(var%name, Field2d)
     956                  CALL xios_send_field(var%name, Field2d)
    957957                  firstx=.false.
    958958               ENDIF
    … …  
    968968!                        IF (iff == iff_beg) THEN
    969969!                          if (prt_level >= 10) then
    970 !                            write(lunout,*)"histwrite2d_phy: .NOT.clef_stations(iff) and iff==iff_beg, call wxios_write_2D"
     970!                            write(lunout,*)"histwrite2d_phy: .NOT.clef_stations(iff) and iff==iff_beg, call xios_send_field"
    971971!                          endif
    972 !                          CALL wxios_write_2D(var%name, Field2d)
     972!                          CALL xios_send_field(var%name, Field2d)
    973973!                        ENDIF
    974974!#endif
    … …  
    10191019                                 nid_files
    10201020#ifdef CPP_XIOS
    1021   USE wxios, only: wxios_write_3D
     1021  USE xios, only: xios_send_field
    10221022#endif
    10231023
    … …  
    10851085                  write(lunout,*)'Dans iophy, histwrite3D iff nlev klev firstx', &
    10861086                                  iff,nlev,klev, firstx                       
    1087                   write(lunout,*)'histwrite3d_phy: call wxios_write_3D for ', &
     1087                  write(lunout,*)'histwrite3d_phy: call xios_send_field for ', &
    10881088                                  trim(var%name), ' with iim jjm nlevx = ', &
    10891089                                  iim,jj_nb,nlevx
    10901090                endif
    1091                 CALL wxios_write_3D(var%name, Field3d(:,:,1:nlevx))
     1091                CALL xios_send_field(var%name, Field3d(:,:,1:nlevx))
    10921092                            firstx=.false.
    10931093              ENDIF
    … …  
    11031103!#ifdef CPP_XIOS
    11041104!                        IF (iff == 1) THEN
    1105 !                              CALL wxios_write_3D(var%name, Field3d(:,:,1:klev))
     1105!                              CALL xios_send_field(var%name, Field3d(:,:,1:klev))
    11061106!                        ENDIF
    11071107!#endif
    … …  
    11491149                                is_sequential, klon_mpi_begin, klon_mpi_end, &
    11501150                                jj_nb, klon_mpi
    1151   USE wxios, only: wxios_write_2D
     1151  USE xios, only: xios_send_field
    11521152
    11531153
    … …  
    11841184
    11851185
    1186         CALL wxios_write_2D(field_name, Field2d)
     1186        CALL xios_send_field(field_name, Field2d)
    11871187
    11881188    ELSE
    … …  
    12201220                                is_sequential, klon_mpi_begin, klon_mpi_end, &
    12211221                                jj_nb, klon_mpi
    1222   USE wxios, only: wxios_write_3D
     1222  USE xios, only: xios_send_field
    12231223
    12241224
    … …  
    12541254        ALLOCATE(index3d(iim*jj_nb*nlev))
    12551255        ALLOCATE(fieldok(iim*jj_nb,nlev))
    1256         CALL wxios_write_3D(field_name, Field3d(:,:,1:nlev))
     1256        CALL xios_send_field(field_name, Field3d(:,:,1:nlev))
    12571257                       
    12581258    ELSE

  • LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

    r2001 r2002  
    316316    endif
    317317   if (iff.le.6) then
    318     CALL wxios_add_vaxis("presnivs", phys_out_filenames(iff), &
     318    CALL wxios_add_vaxis("presnivs", &
    319319            levmax(iff) - levmin(iff) + 1, presnivs(levmin(iff):levmax(iff)))
    320     CALL wxios_add_vaxis("Ahyb", phys_out_filenames(iff), &
     320    CALL wxios_add_vaxis("Ahyb", &
    321321            levmax(iff) - levmin(iff) + 1, Ahyb)
    322     CALL wxios_add_vaxis("Bhyb", phys_out_filenames(iff), &
     322    CALL wxios_add_vaxis("Bhyb", &
    323323            levmax(iff) - levmin(iff) + 1, Bhyb)
    324     CALL wxios_add_vaxis("Alt", phys_out_filenames(iff), &
     324    CALL wxios_add_vaxis("Alt", &
    325325            levmax(iff) - levmin(iff) + 1, Alt)
    326326   else
    327327    ! NMC files
    328     CALL wxios_add_vaxis("plev", phys_out_filenames(iff), &
     328    CALL wxios_add_vaxis("plev", &
    329329            levmax(iff) - levmin(iff) + 1, rlevSTD(levmin(iff):levmax(iff)))
    330330   endif

  • LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

    r1971 r2002  
    230230#ifdef CPP_XIOS
    231231    ! ug Pour les sorties XIOS
    232     USE wxios, only: wxios_update_calendar, wxios_closedef
     232    USE xios, ONLY: xios_update_calendar
     233    USE wxios, only: wxios_closedef
    233234#endif
    234235    USE phys_cal_mod, only : mth_len
    … …  
    292293       IF (vars_defined) THEN
    293294          if (prt_level >= 10) then
    294              write(lunout,*)"phys_output_write: call wxios_update_calendar, itau_w=",itau_w
     295             write(lunout,*)"phys_output_write: call xios_update_calendar, itau_w=",itau_w
    295296          endif
    296           CALL wxios_update_calendar(itau_w)
     297          CALL xios_update_calendar(itau_w)
    297298       END IF
    298299       !$OMP END MASTER
Note: See TracChangeset for help on using the changeset viewer.