source: LMDZ6/trunk/libf/misc/wxios.F90 @ 3165

Last change on this file since 3165 was 3165, checked in by Laurent Fairhead, 6 years ago

Update to DR 1.00.21 and modifications following IPSCL6.0.15 workflow tests

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 23.4 KB
RevLine 
[2509]1! $Id$
[1825]2#ifdef CPP_XIOS
3MODULE wxios
4    USE xios
5    USE iaxis
6    USE iaxis_attr
7    USE icontext_attr
8    USE idate
9    USE idomain_attr
10    USE ifield_attr
11    USE ifile_attr
12    USE ixml_tree
13
14    !Variables disponibles pendant toute l'execution du programme:
15   
16    INTEGER, SAVE :: g_comm
17    CHARACTER(len=100), SAVE :: g_ctx_name
18    TYPE(xios_context), SAVE :: g_ctx
19!$OMP THREADPRIVATE(g_comm,g_cts_name,g_ctx)
[1852]20    LOGICAL, SAVE :: g_flag_xml = .FALSE.
21    CHARACTER(len=100) :: g_field_name = "nofield"
22!$OMP THREADPRIVATE(g_flag_xml,g_field_name)
[2271]23    REAL :: missing_val_omp
24    REAL :: missing_val
25!$OMP THREADPRIVATE(missing_val)
[1825]26
[2741]27#ifdef XIOS1
28#error "XIOS v1 no longer supported, use XIOS v2."
29#endif
30
[1825]31    CONTAINS
32   
33    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
34    !   36day => 36d etc     !!!!!!!!!!!!!!!!!!!!
35    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
36   
37    SUBROUTINE reformadate(odate, ndate)
38        CHARACTER(len=*), INTENT(IN) :: odate
[2509]39        TYPE(xios_duration) :: ndate
[1825]40       
41        INTEGER :: i = 0
[1852]42         !!!!!!!!!!!!!!!!!!
43         ! Pour XIOS:
44         !  year : y
45         !  month : mo
46         !  day : d
47         !  hour : h
48         !  minute : mi
49         !  second : s
50         !!!!!!!!!!!!!!!!!!
51
[1825]52        i = INDEX(odate, "day")
53        IF (i > 0) THEN
[2509]54            read(odate(1:i-1),*) ndate%day
[1825]55        END IF
[1852]56
57        i = INDEX(odate, "hr")
58        IF (i > 0) THEN
[2509]59            read(odate(1:i-1),*) ndate%hour
[1852]60        END IF
61
62        i = INDEX(odate, "mth")
63        IF (i > 0) THEN
[2509]64            read(odate(1:i-1),*) ndate%month
[1852]65        END IF
[1825]66       
[1852]67        !IF (prt_level >= 10) WRITE(lunout,*) "Xios. ", odate, " => ", ndate
[1825]68    END SUBROUTINE reformadate
69   
70    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
71    !   ave(X) => average etc     !!!!!!!!!!!!!!!
72    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
73   
74    CHARACTER(len=7) FUNCTION reformaop(op)
75        CHARACTER(len=*), INTENT(IN) :: op
76       
77        INTEGER :: i = 0
78        reformaop = "average"
79       
80        IF (op.EQ."inst(X)") THEN
81            reformaop = "instant"
82        END IF
83       
84        IF (op.EQ."once") THEN
85            reformaop = "once"
86        END IF
87       
88        IF (op.EQ."t_max(X)") THEN
89            reformaop = "maximum"
90        END IF
91       
92        IF (op.EQ."t_min(X)") THEN
93            reformaop = "minimum"
94        END IF
95       
[1852]96        !IF (prt_level >= 10) WRITE(lunout,*) "Xios. ", op, " => ", reformaop
[1825]97    END FUNCTION reformaop
98
99    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
100    ! Routine d'initialisation      !!!!!!!!!!!!!
101    !     A lancer juste après mpi_init !!!!!!!!!!!!!
102    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
103
[2055]104    SUBROUTINE wxios_init(xios_ctx_name, locom, outcom, type_ocean)
[2509]105        USE print_control_mod, ONLY : prt_level, lunout
[1852]106        IMPLICIT NONE
107
[1825]108      CHARACTER(len=*), INTENT(IN) :: xios_ctx_name
[1852]109      INTEGER, INTENT(IN), OPTIONAL :: locom
110      INTEGER, INTENT(OUT), OPTIONAL :: outcom
[2055]111      CHARACTER(len=6), INTENT(IN), OPTIONAL :: type_ocean
[1852]112
113   
[1825]114        TYPE(xios_context) :: xios_ctx
[1852]115        INTEGER :: xios_comm
[1825]116
[1897]117        IF (prt_level >= 10) WRITE(lunout,*) "wxios_init: Initialization"
[1852]118
119
120
121        IF (PRESENT(locom)) THEN
[1897]122          CALL xios_initialize(xios_ctx_name, local_comm = locom, return_comm = xios_comm )
123          IF (prt_level >= 10) WRITE(lunout,*) "wxios_init: ctx=",xios_ctx_name," local_comm=",locom,", return_comm=",xios_comm
[1852]124        ELSE
[1897]125          CALL xios_initialize(xios_ctx_name, return_comm = xios_comm )
126          IF (prt_level >= 10) WRITE(lunout,*) "wxios_init: ctx=",xios_ctx_name," return_comm=",xios_comm
[1852]127        END IF
[1825]128       
[1852]129        IF (PRESENT(outcom)) THEN
[1897]130          outcom = xios_comm
131          IF (prt_level >= 10) WRITE(lunout,*) "wxios_init: ctx=",xios_ctx_name," outcom=",outcom
[1852]132        END IF
[1825]133       
134        !Enregistrement des variables globales:
135        g_comm = xios_comm
136        g_ctx_name = xios_ctx_name
137       
[2055]138        ! Si couple alors init fait dans cpl_init
139        IF (.not. PRESENT(type_ocean)) THEN
[2054]140            CALL wxios_context_init()
[2055]141        ENDIF
[2054]142
[1825]143    END SUBROUTINE wxios_init
144
[1852]145    SUBROUTINE wxios_context_init()
[2509]146        USE print_control_mod, ONLY : prt_level, lunout
147!        USE mod_phys_lmdz_mpi_data, ONLY : COMM_LMDZ_PHY
[1852]148        IMPLICIT NONE
149
150        TYPE(xios_context) :: xios_ctx
151
[2509]152!$OMP MASTER
[1852]153        !Initialisation du contexte:
154        CALL xios_context_initialize(g_ctx_name, g_comm)
155        CALL xios_get_handle(g_ctx_name, xios_ctx)    !Récupération
156        CALL xios_set_current_context(xios_ctx)            !Activation
157        g_ctx = xios_ctx
158
[2001]159        IF (prt_level >= 10) THEN
160          WRITE(lunout,*) "wxios_context_init: Current context is ",trim(g_ctx_name)
161          WRITE(lunout,*) "     now call xios_solve_inheritance()"
162        ENDIF
[1852]163        !Une première analyse des héritages:
164        CALL xios_solve_inheritance()
[2509]165!$OMP END MASTER
[1852]166    END SUBROUTINE wxios_context_init
167
[1825]168    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
169    ! Routine de paramétrisation !!!!!!!!!!!!!!!!!!
170    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
171
[2095]172    SUBROUTINE wxios_set_cal(pasdetemps, calendrier, annee, mois, jour, heure, ini_an, ini_mois, ini_jour, ini_heure)
[2509]173        USE print_control_mod, ONLY : prt_level, lunout
[1852]174        IMPLICIT NONE
175
[1825]176     !Paramètres:
177     CHARACTER(len=*), INTENT(IN) :: calendrier
[2095]178     INTEGER, INTENT(IN) :: annee, mois, jour, ini_an, ini_mois, ini_jour
179     REAL, INTENT(IN) :: pasdetemps, heure, ini_heure
[1825]180     
181     !Variables:
182     CHARACTER(len=80) :: abort_message
183     CHARACTER(len=19) :: date
184     INTEGER :: njour = 1
185     
186     !Variables pour xios:
[2509]187     TYPE(xios_duration) :: mdtime
[1825]188     !REAL(kind = 8) :: year=0, month=0, day=0, hour=0, minute=0, second=0
189     
[2509]190        mdtime%second=pasdetemps
[1825]191
192        !Réglage du calendrier:
193        SELECT CASE (calendrier)
194            CASE('earth_360d')
[2509]195                CALL xios_define_calendar("D360")
196                IF (prt_level >= 10) WRITE(lunout,*) 'wxios_set_cal: Calendrier terrestre a 360 jours/an'
197            CASE('earth_365d')
198                CALL xios_define_calendar("NoLeap")
199                IF (prt_level >= 10) WRITE(lunout,*) 'wxios_set_cal: Calendrier terrestre a 365 jours/an'
200            CASE('gregorian')
201                CALL xios_define_calendar("Gregorian")
202                IF (prt_level >= 10) WRITE(lunout,*) 'wxios_set_cal: Calendrier gregorien'
203            CASE DEFAULT
204                abort_message = 'wxios_set_cal: Mauvais choix de calendrier'
205                CALL abort_gcm('Gcm:Xios',abort_message,1)
206        END SELECT
[1825]207       
[2095]208        !Formatage de la date d'origine:
[2509]209        WRITE(date, "(i4.4,'-',i2.2,'-',i2.2,' ',i2.2,':00:00')") annee, mois, jour, int(heure)
[1825]210       
[2095]211        IF (prt_level >= 10) WRITE(lunout,*) "wxios_set_cal: Time origin: ", date
[2509]212        CALL xios_set_time_origin(xios_date(annee,mois,jour,int(heure),0,0))
[2095]213
214        !Formatage de la date de debut:
215
216        WRITE(date, "(i4.4,'-',i2.2,'-',i2.2,' ',i2.2,':00:00')") ini_an, ini_mois, ini_jour, int(ini_heure)
[1825]217       
[2095]218        IF (prt_level >= 10) WRITE(lunout,*) "wxios_set_cal: Start date: ", date
219       
[2509]220        CALL xios_set_start_date(xios_date(ini_an,ini_mois,ini_jour,int(ini_heure),0,0))
[2095]221       
[1825]222        !Et enfin,le pas de temps:
223        CALL xios_set_timestep(mdtime)
[1897]224        IF (prt_level >= 10) WRITE(lunout,*) "wxios_set_cal: ts=",mdtime
[1825]225    END SUBROUTINE wxios_set_cal
226
[1852]227    SUBROUTINE wxios_set_timestep(ts)
228        REAL, INTENT(IN) :: ts
[2509]229        TYPE(xios_duration) :: mdtime     
[1852]230
[2509]231        mdtime%timestep = ts
232
[1852]233        CALL xios_set_timestep(mdtime)
234    END SUBROUTINE wxios_set_timestep
235
[1825]236    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
237    ! Pour initialiser un domaine !!!!!!!!!!!!!!!!!!!!
238    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1852]239    SUBROUTINE wxios_domain_param(dom_id, is_sequential, ni, nj, ni_glo, nj_glo,        &
240                                    ibegin, iend, ii_begin, ii_end, jbegin, jend,       &
241                                    data_ni, data_ibegin, data_iend,                    &
[1897]242                                    io_lat, io_lon,is_south_pole,mpi_rank)
[1825]243         
[1897]244
[2509]245        USE print_control_mod, ONLY : prt_level, lunout
[1852]246        IMPLICIT NONE
247
[1897]248        CHARACTER(len=*),INTENT(IN) :: dom_id ! domain identifier
249        LOGICAL,INTENT(IN) :: is_sequential ! flag
250        INTEGER,INTENT(IN) :: ni ! local MPI domain number of longitudes
251        INTEGER,INTENT(IN) :: nj ! local MPI domain number of latitudes
252        INTEGER,INTENT(IN) :: ni_glo ! global grid number of longitudes
253        INTEGER,INTENT(IN) :: nj_glo ! global grid number of latitudes
254        INTEGER,INTENT(IN) :: ibegin ! start index, on global grid, of local MPI domain
255        INTEGER,INTENT(IN) :: iend ! end index, on global grid, of local MPI domain
256        INTEGER,INTENT(IN) :: ii_begin ! i index at which local data starts (first row)
257        INTEGER,INTENT(IN) :: ii_end ! i index at which local data ends (last row)
258        INTEGER,INTENT(IN) :: jbegin ! start index, on global grid, of local MPI domain
259        INTEGER,INTENT(IN) :: jend ! end index, on global grid, of local MPI domain
260        INTEGER,INTENT(IN) :: data_ni
261        INTEGER,INTENT(IN) :: data_ibegin
262        INTEGER,INTENT(IN) :: data_iend
263        REAL,INTENT(IN) :: io_lat(:) ! latitudes (of global grid)
264        REAL,INTENT(IN) :: io_lon(:) ! longitudes (of global grid)
265        logical,intent(in) :: is_south_pole ! does this process include the south pole?
266        integer,intent(in) :: mpi_rank ! rank of process
[1825]267       
268        TYPE(xios_domain) :: dom
[1852]269        LOGICAL :: boool
[1825]270       
[1852]271        !Masque pour les problèmes de recouvrement MPI:
272        LOGICAL :: mask(ni,nj)
[1825]273       
274        !On récupère le handle:
275        CALL xios_get_domain_handle(dom_id, dom)
276       
[1897]277        IF (prt_level >= 10) THEN
278          WRITE(lunout,*) "wxios_domain_param: mpirank=",mpi_rank," ni:",ni," ni_glo:", ni_glo, " nj:", nj, " nj_glo:", nj_glo
279          WRITE(lunout,*) "wxios_domain_param: mpirank=",mpi_rank," ibegin:",ibegin," iend:", iend, " jbegin:", jbegin, " jend:", jend
280          WRITE(lunout,*) "wxios_domain_param: mpirank=",mpi_rank," ii_begin:",ii_begin," ii_end:", ii_end
281          WRITE(lunout,*) "wxios_domain_param: mpirank=",mpi_rank," Size io_lon:", SIZE(io_lon(ibegin:iend)), " io_lat:", SIZE(io_lat(jbegin:jend))
282        ENDIF
[1825]283       
284        !On parametrise le domaine:
[2509]285        CALL xios_set_domain_attr_hdl(dom, ni_glo=ni_glo, ibegin=ibegin-1, ni=ni, type="rectilinear")
286        CALL xios_set_domain_attr_hdl(dom, nj_glo=nj_glo, jbegin=jbegin-1, nj=nj, data_dim=2)
287        CALL xios_set_domain_attr_hdl(dom, lonvalue_1d=io_lon(ibegin:iend), latvalue_1d=io_lat(jbegin:jend))
[3165]288        !On definit un axe de latitudes pour les moyennes zonales
289        IF (xios_is_valid_axis("axis_lat")) THEN
290           CALL xios_set_axis_attr( "axis_lat", n_glo=nj_glo, n=nj, begin=jbegin-1, value=io_lat(jbegin:jend))
291        ENDIF
292
[1852]293        IF (.NOT.is_sequential) THEN
294            mask(:,:)=.TRUE.
[1897]295            if (ii_begin>1) mask(1:ii_begin-1,1) = .FALSE.
296            if (ii_end<ni) mask(ii_end+1:ni,nj) = .FALSE.
297            ! special case for south pole
298            if ((ii_end.eq.1).and.(is_south_pole)) mask(1:ni,nj)=.true.
299            IF (prt_level >= 10) THEN
300              WRITE(lunout,*) "wxios_domain_param: mpirank=",mpi_rank," mask(:,1)=",mask(:,1)
301              WRITE(lunout,*) "wxios_domain_param: mpirank=",mpi_rank," mask(:,nj)=",mask(:,nj)
302            ENDIF
[2509]303            CALL xios_set_domain_attr_hdl(dom, mask_2d=mask)
[1852]304        END IF
305
[1825]306         CALL xios_is_defined_domain_attr_hdl(dom,ni_glo=boool)
307        !Vérification:
308        IF (xios_is_valid_domain(dom_id)) THEN
[1897]309            IF (prt_level >= 10) WRITE(lunout,*) "wxios_domain_param: Domain initialized: ", trim(dom_id), boool
[1825]310        ELSE
[1897]311            IF (prt_level >= 10) WRITE(lunout,*) "wxios_domain_param: Invalid domain: ", trim(dom_id)
[1825]312        END IF
313    END SUBROUTINE wxios_domain_param
314   
315    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
316    ! Pour déclarer un axe vertical !!!!!!!!!!!!!!!
317    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[3003]318    SUBROUTINE wxios_add_vaxis(axis_id, axis_size, axis_value,     &
319                               positif, bnds)
[2509]320        USE print_control_mod, ONLY : prt_level, lunout
[1852]321        IMPLICIT NONE
322
[2002]323        CHARACTER (len=*), INTENT(IN) :: axis_id
[1852]324        INTEGER, INTENT(IN) :: axis_size
[1825]325        REAL, DIMENSION(axis_size), INTENT(IN) :: axis_value
[3003]326        CHARACTER (len=*), INTENT(IN), OPTIONAL :: positif
327        REAL, DIMENSION(axis_size, 2), INTENT(IN), OPTIONAL :: bnds
[1825]328       
[2002]329!        TYPE(xios_axisgroup) :: axgroup
330!        TYPE(xios_axis) :: ax
331!        CHARACTER(len=50) :: axis_id
[1825]332       
[2002]333!        IF (len_trim(axisgroup_id).gt.len(axis_id)) THEN
334!          WRITE(lunout,*) "wxios_add_vaxis: error, size of axis_id too small!!"
335!          WRITE(lunout,*) "     increase it to at least ",len_trim(axisgroup_id)
336!          CALL abort_gcm("wxios_add_vaxis","len(axis_id) too small",1)
337!        ENDIF
338!        axis_id=trim(axisgroup_id)
[1825]339       
340        !On récupère le groupe d'axes qui va bien:
[2002]341        !CALL xios_get_axisgroup_handle(axisgroup_id, axgroup)
[1825]342       
343        !On ajoute l'axe correspondant à ce fichier:
[2002]344        !CALL xios_add_axis(axgroup, ax, TRIM(ADJUSTL(axis_id)))
[1825]345       
346        !Et on le parametrise:
[2002]347        !CALL xios_set_axis_attr_hdl(ax, size=axis_size, value=axis_value)
[1825]348       
[2002]349        ! Ehouarn: New way to declare axis, without axis_group:
[3003]350        if (PRESENT(positif) .AND. PRESENT(bnds)) then
351          CALL xios_set_axis_attr(trim(axis_id),n_glo=axis_size,value=axis_value, &
352                                  positive=positif, bounds=bnds)
353        else if (PRESENT(positif)) then
354          CALL xios_set_axis_attr(trim(axis_id),n_glo=axis_size,value=axis_value, &
355                                  positive=positif)
356        else if (PRESENT(bnds)) then
357          CALL xios_set_axis_attr(trim(axis_id),n_glo=axis_size,value=axis_value, &
358                                  bounds=bnds)
359        else
360          CALL xios_set_axis_attr(trim(axis_id),n_glo=axis_size,value=axis_value)
361        endif
[2741]362
[1825]363        !Vérification:
[1852]364        IF (xios_is_valid_axis(TRIM(ADJUSTL(axis_id)))) THEN
[1897]365            IF (prt_level >= 10) WRITE(lunout,*) "wxios_add_vaxis: Axis created: ", TRIM(ADJUSTL(axis_id))
[1825]366        ELSE
[2001]367            WRITE(lunout,*) "wxios_add_vaxis: Invalid axis: ", TRIM(ADJUSTL(axis_id))
[1825]368        END IF
369
370    END SUBROUTINE wxios_add_vaxis
371   
372   
373    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
374    ! Pour déclarer un fichier  !!!!!!!!!!!!!!!!!!!
375    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
376    SUBROUTINE wxios_add_file(fname, ffreq, flvl)
[2509]377        USE print_control_mod, ONLY : prt_level, lunout
[1852]378        IMPLICIT NONE
379
[1825]380        CHARACTER(len=*), INTENT(IN) :: fname
381        CHARACTER(len=*), INTENT(IN) :: ffreq
382        INTEGER, INTENT(IN) :: flvl
383       
384        TYPE(xios_file) :: x_file
385        TYPE(xios_filegroup) :: x_fg
[2509]386        TYPE(xios_duration) :: nffreq
[1825]387       
[1852]388        !On regarde si le fichier n'est pas défini par XML:
389        IF (.NOT.xios_is_valid_file(fname)) THEN
390            !On créé le noeud:
391            CALL xios_get_filegroup_handle("defile", x_fg)
392            CALL xios_add_file(x_fg, x_file, fname)
[1825]393       
[1852]394            !On reformate la fréquence:
395            CALL reformadate(ffreq, nffreq)
[1825]396       
[1852]397            !On configure:
398            CALL xios_set_file_attr_hdl(x_file, name="X"//fname,&
[2509]399                output_freq=nffreq, output_level=flvl, enabled=.TRUE.)
400       
[1852]401            IF (xios_is_valid_file("X"//fname)) THEN
[2001]402                IF (prt_level >= 10) THEN
403                  WRITE(lunout,*) "wxios_add_file: New file: ", "X"//fname
[2509]404                  WRITE(lunout,*) "wxios_add_file: output_freq=",nffreq,"; output_lvl=",flvl
[2001]405                ENDIF
[1852]406            ELSE
[2001]407                WRITE(lunout,*) "wxios_add_file: Error, invalid file: ", "X"//trim(fname)
[2509]408                WRITE(lunout,*) "wxios_add_file: output_freq=",nffreq,"; output_lvl=",flvl
[1852]409            END IF
[1825]410        ELSE
[2001]411            IF (prt_level >= 10) THEN
412              WRITE(lunout,*) "wxios_add_file: File ",trim(fname), " défined using XML."
413            ENDIF
414            ! Ehouarn: add an enable=.true. on top of xml definitions... why???
415            CALL xios_set_file_attr(fname, enabled=.TRUE.)
[1825]416        END IF
417    END SUBROUTINE wxios_add_file
418   
419    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1852]420    ! Pour créer un champ      !!!!!!!!!!!!!!!!!!!!
[1825]421    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
422    SUBROUTINE wxios_add_field(fieldname, fieldgroup, fieldlongname, fieldunit)
[1897]423        USE netcdf, only: nf90_fill_real
[1852]424
425        IMPLICIT NONE
426        INCLUDE 'iniprint.h'
[1825]427       
428        CHARACTER(len=*), INTENT(IN) :: fieldname
429        TYPE(xios_fieldgroup), INTENT(IN) :: fieldgroup
430        CHARACTER(len=*), INTENT(IN) :: fieldlongname
431        CHARACTER(len=*), INTENT(IN) :: fieldunit
432       
433        TYPE(xios_field) :: field
434        CHARACTER(len=10) :: newunit
435        REAL(KIND=8) :: def
436       
437        !La valeur par défaut des champs non définis:
438        def = nf90_fill_real
439       
440        IF (fieldunit .EQ. " ") THEN
441            newunit = "-"
442        ELSE
443            newunit = fieldunit
444        ENDIF
445       
446        !On ajoute le champ:
447        CALL xios_add_field(fieldgroup, field, fieldname)
[1897]448        !IF (prt_level >= 10) WRITE(lunout,*) "wxios_add_field: ",fieldname,fieldgroup, fieldlongname, fieldunit
[1825]449       
450        !On rentre ses paramètres:
451        CALL xios_set_field_attr_hdl(field, standard_name=fieldlongname, unit=newunit, default_value=def)
[1897]452        IF (prt_level >= 10) WRITE(lunout,*) "wxios_add_field: Field ",trim(fieldname), "cree:"
453        IF (prt_level >= 10) WRITE(lunout,*) "wxios_add_field: long_name=",trim(fieldlongname),"; unit=",trim(newunit),";  default_value=",nf90_fill_real
[1825]454
455    END SUBROUTINE wxios_add_field
456   
457    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[1852]458    ! Pour déclarer un champ      !!!!!!!!!!!!!!!!!
[1825]459    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[2137]460    SUBROUTINE wxios_add_field_to_file(fieldname, fdim, fid, fname, fieldlongname, fieldunit, field_level, op, nam_axvert)
[2509]461        USE print_control_mod, ONLY : prt_level, lunout
[1852]462        IMPLICIT NONE
463
[1825]464        CHARACTER(len=*), INTENT(IN) :: fieldname
465        INTEGER, INTENT(IN)          :: fdim, fid
466        CHARACTER(len=*), INTENT(IN) :: fname
467        CHARACTER(len=*), INTENT(IN) :: fieldlongname
468        CHARACTER(len=*), INTENT(IN) :: fieldunit
469        INTEGER, INTENT(IN)          :: field_level
470        CHARACTER(len=*), INTENT(IN) :: op
471       
[2001]472        CHARACTER(len=20) :: axis_id ! Ehouarn: dangerous...
[2137]473        CHARACTER(len=20), INTENT(IN), OPTIONAL :: nam_axvert
[1825]474        CHARACTER(len=100) :: operation
475        TYPE(xios_file) :: f
476        TYPE(xios_field) :: field
477        TYPE(xios_fieldgroup) :: fieldgroup
[2509]478        TYPE(xios_duration) :: freq_op
[2741]479
[1852]480        LOGICAL :: bool=.FALSE.
481        INTEGER :: lvl =0
[1825]482       
483       
[2001]484        ! Ajout Abd pour NMC:
485        IF (fid.LE.6) THEN
486          axis_id="presnivs"
487        ELSE
488          axis_id="plev"
489        ENDIF
[2137]490 
491        IF (PRESENT(nam_axvert)) THEN
492           axis_id=nam_axvert
493           print*,'nam_axvert=',axis_id
494        ENDIF
[1825]495       
496        !on prépare le nom de l'opération:
497        operation = reformaop(op)
498       
499       
500        !On selectionne le bon groupe de champs:
501        IF (fdim.EQ.2) THEN
[2001]502          CALL xios_get_fieldgroup_handle("fields_2D", fieldgroup)
[1825]503        ELSE
504          CALL xios_get_fieldgroup_handle("fields_3D", fieldgroup)
505        ENDIF
506       
507        !On regarde si le champ à déjà été créé ou non:
[1852]508        IF (xios_is_valid_field(fieldname) .AND. .NOT. g_field_name == fieldname) THEN
509            !Si ce champ existe via XML (ie, dès le premier passage, ie g_field_name != fieldname) alors rien d'autre à faire
[1897]510            IF (prt_level >= 10) WRITE(lunout,*) "wxios_add_field_to_file: Field ", trim(fieldname), "exists via XML"
[1852]511            g_flag_xml = .TRUE.
512            g_field_name = fieldname
513
514        ELSE IF (.NOT. g_field_name == fieldname) THEN
515            !Si premier pssage et champ indéfini, alors on le créé
516
[1897]517            IF (prt_level >= 10) WRITE(lunout,*) "wxios_add_field_to_file: Field ", trim(fieldname), "does not exist"
[1825]518           
519            !On le créé:
520            CALL wxios_add_field(fieldname,  fieldgroup, fieldlongname, fieldunit)
521            IF (xios_is_valid_field(fieldname)) THEN
[1897]522                IF (prt_level >= 10) WRITE(lunout,*) "wxios_add_field_to_file: Field ", trim(fieldname), "created"
[1825]523            ENDIF
524
[1852]525            g_flag_xml = .FALSE.
526            g_field_name = fieldname
527
528        END IF
529
530        IF (.NOT. g_flag_xml) THEN
531            !Champ existe déjà, mais pas XML, alors on l'ajoute
532            !On ajoute le champ:
533            CALL xios_get_file_handle(fname, f)
534            CALL xios_add_fieldtofile(f, field)
535           
536           
537            !L'operation, sa frequence:
[2509]538            freq_op%timestep=1
539            CALL xios_set_field_attr_hdl(field, field_ref=fieldname, operation=TRIM(ADJUSTL(operation)), freq_op=freq_op, prec=4)
[1852]540
541           
542            !On rentre ses paramètres:
543            CALL xios_set_field_attr_hdl(field, level=field_level, enabled=.TRUE.)
544           
545            IF (fdim.EQ.2) THEN
546                !Si c'est un champ 2D:
[1897]547                IF (prt_level >= 10) THEN
548                  WRITE(lunout,*) "wxios_add_field_to_file: 2D Field ", trim(fieldname), " in ", "X"//trim(fname) ," configured with:"
549                  WRITE(lunout,*) "wxios_add_field_to_file: op=", TRIM(ADJUSTL(operation))
550                  WRITE(lunout,*) "wxios_add_field_to_file: freq_op=1ts","; lvl=",field_level
551                ENDIF
[1852]552            ELSE
553                !Si 3D :
554                !On ajoute l'axe vertical qui va bien:
555                CALL xios_set_field_attr_hdl(field, axis_ref=TRIM(ADJUSTL(axis_id)))
556               
[1897]557                IF (prt_level >= 10) THEN
558                  WRITE(lunout,*) "wxios_add_field_to_file: 3D Field",trim(fieldname), " in ", "X"//trim(fname), "configured with:"
559                  WRITE(lunout,*) "wxios_add_field_to_file: freq_op=1ts","; lvl=",field_level
560                  WRITE(lunout,*) "wxios_add_field_to_file: axis=",TRIM(ADJUSTL(axis_id))
561                ENDIF
[1852]562            END IF
[1825]563       
564        ELSE
[1852]565            !Sinon on se contente de l'activer:
566            CALL xios_set_field_attr(fieldname, enabled=.TRUE.)
[2002]567            !NB: This will override an enable=.false. set by a user in the xml file;
568            !   then the only way to not output the field is by changing its
569            !   output level
[1852]570        ENDIF       
[1825]571       
572    END SUBROUTINE wxios_add_field_to_file
573   
[2002]574!    SUBROUTINE wxios_update_calendar(ito)
575!        INTEGER, INTENT(IN) :: ito
576!        CALL xios_update_calendar(ito)
577!    END SUBROUTINE wxios_update_calendar
578!   
579!    SUBROUTINE wxios_write_2D(fieldname, fdata)
580!        CHARACTER(len=*), INTENT(IN) :: fieldname
581!        REAL, DIMENSION(:,:), INTENT(IN) :: fdata
582!
583!        CALL xios_send_field(fieldname, fdata)
584!    END SUBROUTINE wxios_write_2D
[1825]585   
[2002]586!    SUBROUTINE wxios_write_3D(fieldname, fdata)
587!        CHARACTER(len=*), INTENT(IN) :: fieldname
588!        REAL, DIMENSION(:,:,:), INTENT(IN) :: fdata
589!       
590!        CALL xios_send_field(fieldname, fdata)
591!    END SUBROUTINE wxios_write_3D
[1825]592   
593    SUBROUTINE wxios_closedef()
594        CALL xios_close_context_definition()
[2095]595!        CALL xios_update_calendar(0)
[1825]596    END SUBROUTINE wxios_closedef
597   
598    SUBROUTINE wxios_close()
599        CALL xios_context_finalize()
600         CALL xios_finalize()
601     END SUBROUTINE wxios_close
602END MODULE wxios
603#endif
Note: See TracBrowser for help on using the repository browser.