Ignore:
Timestamp:
Jul 23, 2024, 7:00:20 AM (4 months ago)
Author:
abarral
Message:

Remove CPP_1D key. It was used once in a single file to wrap a whole internal module -> can't even compile if key is enabled.
(lint) Set NF90_* to lowercase

Location:
LMDZ6/branches/Amaury_dev/libf/phylmd
Files:
13 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/phylmd/create_etat0_unstruct_mod.F90

    r5099 r5100  
    2323  SUBROUTINE init_create_etat0_unstruct
    2424  USE lmdz_xios
    25   USE netcdf, ONLY: NF90_NOWRITE,nf90_close,nf90_noerr,nf90_open
     25  USE netcdf, ONLY: nf90_nowrite,nf90_close,nf90_noerr,nf90_open
    2626  USE mod_phys_lmdz_para
    2727  IMPLICIT NONE
    … …  
    3131    IF (is_omp_master) THEN
    3232
    33       IF (NF90_OPEN("ocean_fraction.nc", NF90_NOWRITE, file_id)==nf90_noerr) THEN
     33      IF (nf90_open("ocean_fraction.nc", nf90_nowrite, file_id)==nf90_noerr) THEN
    3434        CALL xios_set_file_attr("frac_ocean",enabled=.TRUE.)
    3535        CALL xios_set_field_attr("mask",field_ref="frac_ocean_read")
    36         iret=NF90_CLOSE(file_id)
    37       ELSE IF (NF90_OPEN("land_water_0.05.nc", NF90_NOWRITE, file_id)==nf90_noerr) THEN
     36        iret=nf90_close(file_id)
     37      ELSE IF (nf90_open("land_water_0.05.nc", nf90_nowrite, file_id)==nf90_noerr) THEN
    3838        CALL xios_set_file_attr("land_water",name="land_water_0.05",enabled=.TRUE.)
    3939        CALL xios_set_field_attr("mask",field_ref="land_water")
    40         iret=NF90_CLOSE(file_id)
    41       ELSE IF (NF90_OPEN("land_water_0.25.nc", NF90_NOWRITE, file_id)==nf90_noerr) THEN
     40        iret=nf90_close(file_id)
     41      ELSE IF (nf90_open("land_water_0.25.nc", nf90_nowrite, file_id)==nf90_noerr) THEN
    4242        CALL xios_set_file_attr("land_water",name="land_water_0.25",enabled=.TRUE.)
    4343        CALL xios_set_field_attr("mask",field_ref="land_water")
    44         iret=NF90_CLOSE(file_id)
    45       ELSE IF (NF90_OPEN("land_water_0.50.nc", NF90_NOWRITE, file_id)==nf90_noerr) THEN
     44        iret=nf90_close(file_id)
     45      ELSE IF (nf90_open("land_water_0.50.nc", nf90_nowrite, file_id)==nf90_noerr) THEN
    4646        CALL xios_set_file_attr("land_water",name="land_water_0.50",enabled=.TRUE.)
    4747        CALL xios_set_field_attr("mask",field_ref="land_water")
    48         iret=NF90_CLOSE(file_id)
     48        iret=nf90_close(file_id)
    4949      ENDIF
    5050

  • LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/lmdz1d.F90

    r5099 r5100  
    22! $Id$
    33
    4 !#ifdef CPP_1D
    54!#include "../dyn3d/mod_const_mpi.F90"
    65!#include "../dyn3d_common/control_mod.F90"

  • LMDZ6/branches/Amaury_dev/libf/phylmd/grid_noro_m.F90

    r5099 r5100  
    435435! Purpose: Read parameters usually determined with grid_noro from a file.
    436436!===============================================================================
    437   USE netcdf, ONLY: NF90_OPEN,  nf90_inq_dimid, nf90_inquire_dimension,        &
    438         nf90_noerr, NF90_CLOSE, NF90_INQ_VARID, nf90_get_var, NF90_STRERROR,   &
    439         NF90_NOWRITE
     437  USE netcdf, ONLY: nf90_open,  nf90_inq_dimid, nf90_inquire_dimension,        &
     438        nf90_noerr, nf90_close, NF90_INQ_VARID, nf90_get_var, NF90_STRERROR,   &
     439        nf90_nowrite
    440440  IMPLICIT NONE
    441441!-------------------------------------------------------------------------------
    … …  
    471471  masque_lu=ANY(mask/=-99999.); IF(.NOT.masque_lu) mask=0.0
    472472  WRITE(lunout,*)'Masque lu: ',masque_lu
    473   CALL ncerr(NF90_OPEN(fname,NF90_NOWRITE,fid))
     473  CALL ncerr(nf90_open(fname,nf90_nowrite,fid))
    474474  CALL check_dim('x','longitude',x(1:imar))
    475475  CALL check_dim('y','latitude' ,y(1:jmar))
    … …  
    483483  zpic=zmea+2*zstd
    484484  zval=MAX(0.,zmea-2.*zstd)
    485   CALL ncerr(NF90_CLOSE(fid))
     485  CALL ncerr(nf90_close(fid))
    486486  WRITE(lunout,*)'  MEAN ORO:' ,MAXVAL(zmea)
    487487  WRITE(lunout,*)'  ST. DEV.:' ,MAXVAL(zstd)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/iostart.F90

    r5099 r5100  
    3737
    3838    IF (is_mpi_root .AND. is_omp_root) THEN
    39       ierr = NF90_OPEN (filename, NF90_NOWRITE,nid_start)
     39      ierr = nf90_open (filename, nf90_nowrite,nid_start)
    4040      IF (ierr/=nf90_noerr) THEN
    4141        write(6,*)' Pb d''ouverture du fichier '//filename
    … …  
    5454
    5555    IF (is_mpi_root .AND. is_omp_root) THEN
    56         ierr = NF90_CLOSE (nid_start)
     56        ierr = nf90_close (nid_start)
    5757    ENDIF
    5858
    … …  
    312312   
    313313    IF (is_master) THEN
    314       ierr = NF90_CREATE(filename, IOR(NF90_CLOBBER,NF90_64BIT_OFFSET), &
     314      ierr = nf90_create(filename, IOR(nf90_clobber,nf90_64bit_offset), &
    315315                          nid_restart)
    316316      IF (ierr/=nf90_noerr) THEN
    … …  
    320320      ENDIF
    321321
    322       ierr = NF90_PUT_ATT (nid_restart, NF90_GLOBAL, "title","Fichier redemmarage physique")
    323 
    324       ierr = NF90_DEF_DIM (nid_restart, "index", length, idim1)
    325       ierr = NF90_DEF_DIM (nid_restart, "points_physiques", klon_glo, idim2)
    326       ierr = NF90_DEF_DIM (nid_restart, "horizon_vertical", klon_glo*klev, idim3)
    327       ierr = NF90_DEF_DIM (nid_restart, "horizon_klevp1", klon_glo*klevp1, idim4)
    328 
    329 !      ierr = NF90_ENDDEF(nid_restart)
     322      ierr = nf90_put_att (nid_restart, nf90_global, "title","Fichier redemmarage physique")
     323
     324      ierr = nf90_def_dim (nid_restart, "index", length, idim1)
     325      ierr = nf90_def_dim (nid_restart, "points_physiques", klon_glo, idim2)
     326      ierr = nf90_def_dim (nid_restart, "horizon_vertical", klon_glo*klev, idim3)
     327      ierr = nf90_def_dim (nid_restart, "horizon_klevp1", klon_glo*klevp1, idim4)
     328
     329!      ierr = nf90_enddef(nid_restart)
    330330    ENDIF
    331331
    … …  
    338338    INTEGER          :: ierr
    339339
    340     IF (is_master) ierr = NF90_ENDDEF(nid_restart)
     340    IF (is_master) ierr = nf90_enddef(nid_restart)
    341341 
    342342  END SUBROUTINE enddef_restartphy
    … …  
    348348    INTEGER          :: ierr
    349349
    350     IF (is_master) ierr = NF90_CLOSE (nid_restart)
     350    IF (is_master) ierr = nf90_close (nid_restart)
    351351 
    352352  END SUBROUTINE close_restartphy
    … …  
    426426         
    427427!      ierr = NF90_REDEF (nid_restart)
    428       ierr = NF90_DEF_VAR (nid_restart, field_name, nf90_format,(/ idim /),nvarid)
    429       IF (LEN_TRIM(title) > 0) ierr = NF90_PUT_ATT (nid_restart,nvarid,"title", title)
    430 !      ierr = NF90_ENDDEF(nid_restart)
     428      ierr = nf90_def_var (nid_restart, field_name, nf90_format,(/ idim /),nvarid)
     429      IF (LEN_TRIM(title) > 0) ierr = nf90_put_att (nid_restart,nvarid,"title", title)
     430!      ierr = nf90_enddef(nid_restart)
    431431     ENDIF
    432432
    … …  
    454454
    455455       ierr = NF90_INQ_VARID(nid_restart, field_name, nvarid)
    456        ierr = NF90_PUT_VAR(nid_restart,nvarid,RESHAPE(field_glo,(/klon_glo*field_size/)))
     456       ierr = nf90_put_var(nid_restart,nvarid,RESHAPE(field_glo,(/klon_glo*field_size/)))
    457457      ENDIF
    458458   ENDIF
    … …  
    536536!      ierr = NF90_REDEF (nid_restart)
    537537
    538         ierr = NF90_DEF_VAR (nid_restart, var_name, nf90_format,(/ idim1 /),nvarid)
    539         IF (LEN_TRIM(title)>0) ierr = NF90_PUT_ATT (nid_restart,nvarid,"title", title)
    540 !      ierr = NF90_ENDDEF(nid_restart)
     538        ierr = nf90_def_var (nid_restart, var_name, nf90_format,(/ idim1 /),nvarid)
     539        IF (LEN_TRIM(title)>0) ierr = nf90_put_att (nid_restart,nvarid,"title", title)
     540!      ierr = nf90_enddef(nid_restart)
    541541
    542542    ! second pass : write     
    543543      ELSE IF (pass==2) THEN
    544544        ierr = NF90_INQ_VARID(nid_restart, var_name, nvarid)
    545         ierr = NF90_PUT_VAR(nid_restart,nvarid,var)
     545        ierr = nf90_put_var(nid_restart,nvarid,var)
    546546      ENDIF
    547547    ENDIF

  • LMDZ6/branches/Amaury_dev/libf/phylmd/iotd_ecrit.F90

    r5099 r5100  
    109109
    110110
    111            ierr= NF90_PUT_VAR(nid,varid,date,[ntime])
     111           ierr= nf90_put_var(nid,varid,date,[ntime])
    112112
    113113!          print*,'date ',date,ierr,nid
    … …  
    173173
    174174
    175       ierr= NF90_PUT_VAR(nid,varid,zx,corner,edges)
     175      ierr= nf90_put_var(nid,varid,zx,corner,edges)
    176176
    177177      if (ierr/=nf90_noerr) then

  • LMDZ6/branches/Amaury_dev/libf/phylmd/limit_read_mod.F90

    r5099 r5100  
    234234!$OMP MASTER  ! Only master thread
    235235       IF (is_mpi_root) THEN ! Only master processus
    236           ierr = NF90_OPEN ('limit.nc', NF90_NOWRITE, nid)
     236          ierr = nf90_open ('limit.nc', nf90_nowrite, nid)
    237237          IF (ierr /= nf90_noerr) CALL abort_physic(modname,&
    238238               'Pb d''ouverture du fichier de conditions aux limites',1)
    … …  
    240240          !--- WARNING IF CALENDAR IS KNOWN AND DOES NOT MATCH THE ONE OF LMDZ
    241241          ierr=NF90_INQ_VARID(nid, 'TEMPS', nvarid)
    242           ierr=NF90_GET_ATT(nid, nvarid, 'calendar', calendar)
     242          ierr=nf90_get_att(nid, nvarid, 'calendar', calendar)
    243243          IF(ierr==nf90_noerr.AND.calendar/=calend.AND.prt_level>=1) THEN
    244244             WRITE(lunout,*)'BEWARE: gcm and limit.nc calendars differ: '
    … …  
    269269          IF(nn/=klon_glo) CALL abort_physic(modname,abort_message,1)
    270270
    271           ierr = NF90_CLOSE(nid)
     271          ierr = nf90_close(nid)
    272272          IF (ierr /= nf90_noerr) CALL abort_physic(modname,'Pb when closing file', 1)
    273273       END IF ! is_mpi_root
    … …  
    332332       IF (is_mpi_root) THEN ! Only master processus!
    333333
    334           ierr = NF90_OPEN ('limit.nc', NF90_NOWRITE, nid)
     334          ierr = nf90_open ('limit.nc', nf90_nowrite, nid)
    335335          IF (ierr /= nf90_noerr) CALL abort_physic(modname,&
    336336               'Pb d''ouverture du fichier de conditions aux limites',1)
    … …  
    426426
    427427!****************************************************************************************
    428           ierr = NF90_CLOSE(nid)
     428          ierr = nf90_close(nid)
    429429          IF (ierr /= nf90_noerr) CALL abort_physic(modname,'Pb when closing file', 1)
    430430       ENDIF ! is_mpi_root

  • LMDZ6/branches/Amaury_dev/libf/phylmd/limit_slab.F90

    r5099 r5100  
    6161        read_siv=.TRUE.
    6262       
    63         ierr = NF90_OPEN ('limit_slab.nc', NF90_NOWRITE, nid)
     63        ierr = nf90_open ('limit_slab.nc', nf90_nowrite, nid)
    6464        IF (ierr /= nf90_noerr) THEN
    6565            PRINT *,'LIMIT_SLAB file not found'
    … …  
    145145
    146146!****************************************************************************************
    147         ierr = NF90_CLOSE(nid)
     147        ierr = nf90_close(nid)
    148148        IF (ierr /= nf90_noerr) CALL abort_physic(modname,'Pb when closing file', 1)
    149149        END IF ! Read File

  • LMDZ6/branches/Amaury_dev/libf/phylmd/mo_simple_plumes.F90

    r5099 r5100  
    8888    IF (is_mpi_root.AND.is_omp_root) THEN
    8989
    90        iret = nf90_open("MACv2.0-SP_v1.nc", NF90_NOWRITE, ncid)
     90       iret = nf90_open("MACv2.0-SP_v1.nc", nf90_nowrite, ncid)
    9191       IF (iret /= nf90_noerr) THEN
    9292          abort_message='NetCDF File not opened'

  • LMDZ6/branches/Amaury_dev/libf/phylmd/open_climoz_m.F90

    r5099 r5100  
    5353    press_in_cen = press_in_cen * 100.
    5454    nlev = SIZE(press_in_cen)
    55     CALL NF95_INQ_VARID(ncID, "time", varID)
     55    CALL nf95_inq_varid(ncID, "time", varID)
    5656    CALL NF95_GW_VAR(ncid, varid, time_in)
    5757    ntim = SIZE(time_in)

  • LMDZ6/branches/Amaury_dev/libf/phylmd/read_map2D.F90

    r5099 r5100  
    3030! Read variable from file. Done by master process MPI and master thread OpenMP
    3131  IF (is_mpi_root .AND. is_omp_root) THEN
    32      ierr = NF90_OPEN(trim(filename), NF90_NOWRITE, nid)
     32     ierr = nf90_open(trim(filename), nf90_nowrite, nid)
    3333     IF (ierr /= nf90_noerr) CALL write_err_mess('Problem in opening file')
    3434
    … …  
    4141     IF (ierr /= nf90_noerr) CALL write_err_mess('Problem in reading varaiable')
    4242
    43      ierr = NF90_CLOSE(nid)
     43     ierr = nf90_close(nid)
    4444     IF (ierr /= nf90_noerr) CALL write_err_mess('Problem in closing file')
    4545

  • LMDZ6/branches/Amaury_dev/libf/phylmd/readaerosol_mod.F90

    r5099 r5100  
    200200        IF (is_mpi_root) THEN
    201201   
    202           IF (nf90_open(TRIM(file_aerosol), NF90_NOWRITE, ncid) /= nf90_noerr) THEN
    203             CALL check_err( nf90_open(TRIM(file_so4), NF90_NOWRITE, ncid), "pb open "//trim(file_so4) )
     202          IF (nf90_open(TRIM(file_aerosol), nf90_nowrite, ncid) /= nf90_noerr) THEN
     203            CALL check_err( nf90_open(TRIM(file_so4), nf90_nowrite, ncid), "pb open "//trim(file_so4) )
    204204          ENDIF
    205205
    … …  
    351351 
    352352       WRITE(lunout,*) 'reading variable ',TRIM(varname),' in file ', TRIM(fname)
    353        CALL check_err( nf90_open(TRIM(fname), NF90_NOWRITE, ncid), "pb open "//trim(fname) )
     353       CALL check_err( nf90_open(TRIM(fname), nf90_nowrite, ncid), "pb open "//trim(fname) )
    354354
    355355

  • LMDZ6/branches/Amaury_dev/libf/phylmd/regr_horiz_time_climoz_m.F90

    r5099 r5100  
    11MODULE regr_horiz_time_climoz_m
    22
    3   USE interpolation,     ONLY: locate
    4   USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, grid_type, unstructured
    5   USE nrtype,            ONLY: pi
    6   USE netcdf,   ONLY: NF90_CLOBBER, NF90_FLOAT,     NF90_OPEN,  &
    7                       NF90_NOWRITE, nf90_noerr,     NF90_GET_ATT, NF90_GLOBAL
    8   USE netcdf95, ONLY: NF95_DEF_DIM, NF95_INQ_DIMID, NF95_INQUIRE_DIMENSION,    &
    9        NF95_DEF_VAR, NF95_INQ_VARID, NF95_INQUIRE_VARIABLE, &
    10        NF95_OPEN,  NF95_CREATE,  NF95_GET_ATT,   NF95_GW_VAR,  nf95_get_var, &
    11        NF95_CLOSE, NF95_ENDDEF,  NF95_PUT_ATT,  NF95_PUT_VAR, NF95_COPY_ATT
    12   USE print_control_mod, ONLY: lunout
     3  USE interpolation, ONLY : locate
     4  USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat, grid_type, unstructured
     5  USE nrtype, ONLY : pi
     6  USE netcdf, ONLY : nf90_clobber, nf90_float, nf90_open, &
     7          nf90_nowrite, nf90_noerr, nf90_get_att, nf90_global
     8  USE netcdf95, ONLY : nf95_def_dim, nf95_inq_dimid, nf95_inquire_dimension, &
     9          nf95_def_var, nf95_inq_varid, NF95_INQUIRE_VARIABLE, &
     10          NF95_OPEN, NF95_CREATE, NF95_GET_ATT, NF95_GW_VAR, nf95_get_var, &
     11          NF95_CLOSE, NF95_ENDDEF, NF95_PUT_ATT, NF95_PUT_VAR, NF95_COPY_ATT
     12  USE print_control_mod, ONLY : lunout
    1313  USE dimphy
    1414  IMPLICIT NONE
    1515  PRIVATE
    1616  PUBLIC :: regr_horiz_time_climoz
    17   REAL, PARAMETER :: deg2rad=pi/180.
    18   CHARACTER(LEN=13), PARAMETER :: vars_in(2)=['tro3         ','tro3_daylight']
     17  REAL, PARAMETER :: deg2rad = pi / 180.
     18  CHARACTER(LEN = 13), PARAMETER :: vars_in(2) = ['tro3         ', 'tro3_daylight']
    1919
    2020  INTEGER :: nlat_ou, nlon_ou
    2121  REAL, ALLOCATABLE :: latitude_glo(:)
    22 !$OMP THREADPRIVATE(latitude_glo)
     22  !$OMP THREADPRIVATE(latitude_glo)
    2323  INTEGER, ALLOCATABLE :: ind_cell_glo_glo(:)
    24 !$OMP THREADPRIVATE(ind_cell_glo_glo)
     24  !$OMP THREADPRIVATE(ind_cell_glo_glo)
    2525
    2626CONTAINS
    2727
    28 !-------------------------------------------------------------------------------
    29 
    30 SUBROUTINE regr_horiz_time_climoz(read_climoz,interpt)
    31 
    32 !-------------------------------------------------------------------------------
    33 ! Purpose: Regrid horizontally and in time zonal or 3D ozone climatologies.
    34 !   * Read ozone climatology from netcdf file
    35 !   * Regrid it horizontaly to LMDZ grid (quasi-conservative method)
    36 !   * If interpt=T, interpolate linearly in time (one record each day)
    37 !     If interpt=F, keep original time sampling  (14 months).
    38 !   * Save it to a new netcdf file.
    39 !-------------------------------------------------------------------------------
    40 ! Remarks:
    41 !   * Up to 2 variables treated: "tro3" and "tro3_daylight" (if read_climoz=2)
    42 !   * Input fields coordinates: (longitudes, latitudes, pressure_levels, time)
    43 !   * Output grid cells centers coordinates given by [rlonv,] rlatu.
    44 !   * Output grid cells edges   coordinates given by [rlonu,] rlatv.
    45 !   * Input file [longitudes and] latitudes given in degrees.
    46 !   * Input file pressure levels are given in Pa or hPa.
    47 !   * All coordinates variables are stricly monotonic.
    48 !   * Monthly fields are interpolated linearly in time to get daily values.
    49 !   * Fields are known at the middle of the months, so interpolation requires an
    50 !     additional record both for 1st half of january and 2nd half of december:
    51 !     - For a 14-records "climoz.nc": records 1 and 14.
    52 !     - For 12-records files:
    53 !       record 12 of "climoz_m.nc" if available, or record 1  of "climoz.nc".
    54 !       record 1  of "climoz_p.nc" if available, or record 12 of "climoz.nc".
    55 !   * Calendar is taken into account to get one record each day (not 360 always).
    56 !   * Missing values are filled in from sky to ground by copying lowest valid one.
    57 !     Attribute "missing_value" or "_FillValue" must be present in input file.
    58 !-------------------------------------------------------------------------------
    59   USE assert_m,           ONLY: assert
    60   USE cal_tools_m,        ONLY: year_len, mid_month
    61 !!  USE control_mod,        ONLY: anneeref
    62   USE time_phylmdz_mod,   ONLY: annee_ref
    63   USE ioipsl,             ONLY: ioget_year_len, ioget_calendar
    64   USE regr_conserv_m,     ONLY: regr_conserv
    65   USE regr_lint_m,        ONLY: regr_lint
    66   USE regular_lonlat_mod, ONLY: boundslon_reg, boundslat_reg, south, west, east
    67   USE slopes_m,           ONLY: slopes
    68   USE lmdz_xios
    69   USE mod_phys_lmdz_para, ONLY: is_mpi_root, is_master, is_omp_master, gather, gather_mpi, bcast_mpi, klon_mpi
    70   USE geometry_mod, ONLY : latitude_deg, ind_cell_glo
    71   USE mod_grid_phy_lmdz, ONLY: klon_glo
    72 
    73 !-------------------------------------------------------------------------------
    74 ! Arguments:
    75   INTEGER, INTENT(IN) :: read_climoz ! read ozone climatology, 1 or 2
    76 !                         1: read a single ozone climatology used day and night
    77 !                         2: same + read also a daylight climatology
    78   LOGICAL, INTENT(IN) :: interpt     ! TRUE  => daily interpolation
    79                                      ! FALSE => no interpolation (14 months)
    80 !-------------------------------------------------------------------------------
    81 ! Local variables:
    82 
    83 !--- Input files variables
    84   INTEGER :: nlon_in                       ! Number of longitudes
    85   INTEGER :: nlat_in                       ! Number of latitudes
    86   INTEGER :: nlev_in                       ! Number of pressure levels
    87   INTEGER :: nmth_in                       ! Number of months
    88   REAL, ALLOCATABLE:: lon_in(:)           ! Longitudes   (ascending order, rad)
    89   REAL, ALLOCATABLE:: lat_in(:)           ! Latitudes    (ascending order, rad)
    90   REAL, ALLOCATABLE:: lev_in(:)           ! Pressure levels (ascen. order, hPa)
    91   REAL, ALLOCATABLE :: lon_in_edge(:)      ! Longitude intervals edges
    92                                            !              (ascending order,  / )
    93   REAL, ALLOCATABLE :: sinlat_in_edge(:)   ! Sinus of latitude intervals edges
    94                                            !              (ascending order,  / )
    95   LOGICAL :: ldec_lon, ldec_lat, ldec_lev  ! Decreasing order in input file
    96   CHARACTER(LEN=20) :: cal_in              ! Calendar
    97   REAL, ALLOCATABLE :: o3_in3(:,:,:,:,:)   ! Ozone climatologies
    98   REAL, ALLOCATABLE :: o3_in3bis(:,:,:,:,:)   ! Ozone climatologies
    99   REAL, ALLOCATABLE :: o3_in2  (:,:,:,:)   ! Ozone climatologies
    100   REAL, ALLOCATABLE :: o3_in2bis(:,:,:,:,:)   ! Ozone climatologies
    101   ! last index: 1 for the day-night average, 2 for the daylight field.
    102   REAL :: NaN
    103 
    104 !--- Partially or totally regridded variables      (:,:,nlev_in,:,read_climoz)
    105   REAL, ALLOCATABLE :: o3_regr_lon   (:,:,:,:,:) ! (nlon_ou,nlat_in,:,0:13   ,:)
    106   REAL, ALLOCATABLE :: o3_regr_lonlat(:,:,:,:,:) ! (nlon_ou,nlat_ou,:,0:13   ,:)
    107   REAL, ALLOCATABLE :: o3_out3       (:,:,:,:,:) ! (nlon_ou,nlat_ou,:,ntim_ou,:)
    108   REAL, ALLOCATABLE :: o3_out3_glo   (:,:,:,:) !   (nbp_lat,:,ntim_ou,:)
    109   REAL, ALLOCATABLE :: o3_regr_lat     (:,:,:,:) !         (nlat_in,:,0:13   ,:)
    110   REAL, ALLOCATABLE :: o3_out2         (:,:,:,:) !         (nlat_ou,:,ntim_ou,:)
    111   REAL, ALLOCATABLE :: o3_out2_glo     (:,:,:,:) !         (nbp_lat,:,ntim_ou,:)
    112   REAL, ALLOCATABLE :: o3_out          (:,:,:,:) !         (nbp_lat,:,ntim_ou,:)
    113 ! Dimension number  | Interval                | Contains  | For variables:
    114 !   1 (longitude)   | [rlonu(i-1), rlonu(i)]  | rlonv(i)  | all
    115 !   2 (latitude)    | [rlatv(j), rlatv(j-1)]  | rlatu(j)  | all but o3_regr_lon
    116 !   3 (press level) |                         |   lev(k)  | all
    117 ! Note that rlatv(0)=pi/2 and rlatv(nlat_ou)=-pi/2.
    118 ! Dimension 4 is: month number                             (all vars but o3_out)
    119 !                 days elapsed since Jan. 1st 0h at mid-day (o3_out only)
    120   REAL, ALLOCATABLE :: v1(:)
    121 
    122 !--- For NetCDF:
    123   INTEGER :: fID_in_m, fID_in, levID_ou, dimid, vID_in(read_climoz), ntim_ou
    124   INTEGER :: fID_in_p, fID_ou, timID_ou, varid, vID_ou(read_climoz), ndims, ncerr
    125   INTEGER, ALLOCATABLE :: dIDs(:)
    126   CHARACTER(LEN=20) :: cal_ou     !--- Calendar; no time inter => same as input
    127   CHARACTER(LEN=80) :: press_unit !--- Pressure unit
    128   REAL    :: tmidmonth(0:13)      !--- Elapsed days since Jan-1 0h at mid-months
    129                                   ! Additional records 0, 13 for interpolation
    130   REAL, ALLOCATABLE :: tmidday(:) !--- Output times (mid-days since Jan 1st 0h)
    131   LOGICAL :: lprev, lnext         !--- Flags: previous/next files are present
    132   LOGICAL :: l3D, l2D             !--- Flag:  input fields are 3D or zonal
    133   INTEGER :: ii, i, j, k, l, m, dln, ib, ie, iv, dx1, dx2
    134   INTEGER, ALLOCATABLE :: sta(:), cnt(:)
    135   CHARACTER(LEN=80) :: sub, dim_nam, msg
    136   REAL :: null_array(0)
    137   LOGICAL,SAVE :: first=.TRUE.
    138 !$OMP THREADPRIVATE(first) 
    139   REAL, ALLOCATABLE :: test_o3_in(:,:)
    140   REAL, ALLOCATABLE :: test_o3_out(:)
    141 
    142 
    143   IF (grid_type==unstructured) THEN
    144     IF (first) THEN
    145       IF (is_master) THEN
    146         ALLOCATE(latitude_glo(klon_glo))
    147         ALLOCATE(ind_cell_glo_glo(klon_glo))
    148       ELSE
    149         ALLOCATE(latitude_glo(0))
    150         ALLOCATE(ind_cell_glo_glo(0))
     28  !-------------------------------------------------------------------------------
     29
     30  SUBROUTINE regr_horiz_time_climoz(read_climoz, interpt)
     31
     32    !-------------------------------------------------------------------------------
     33    ! Purpose: Regrid horizontally and in time zonal or 3D ozone climatologies.
     34    !   * Read ozone climatology from netcdf file
     35    !   * Regrid it horizontaly to LMDZ grid (quasi-conservative method)
     36    !   * If interpt=T, interpolate linearly in time (one record each day)
     37    !     If interpt=F, keep original time sampling  (14 months).
     38    !   * Save it to a new netcdf file.
     39    !-------------------------------------------------------------------------------
     40    ! Remarks:
     41    !   * Up to 2 variables treated: "tro3" and "tro3_daylight" (if read_climoz=2)
     42    !   * Input fields coordinates: (longitudes, latitudes, pressure_levels, time)
     43    !   * Output grid cells centers coordinates given by [rlonv,] rlatu.
     44    !   * Output grid cells edges   coordinates given by [rlonu,] rlatv.
     45    !   * Input file [longitudes and] latitudes given in degrees.
     46    !   * Input file pressure levels are given in Pa or hPa.
     47    !   * All coordinates variables are stricly monotonic.
     48    !   * Monthly fields are interpolated linearly in time to get daily values.
     49    !   * Fields are known at the middle of the months, so interpolation requires an
     50    !     additional record both for 1st half of january and 2nd half of december:
     51    !     - For a 14-records "climoz.nc": records 1 and 14.
     52    !     - For 12-records files:
     53    !       record 12 of "climoz_m.nc" if available, or record 1  of "climoz.nc".
     54    !       record 1  of "climoz_p.nc" if available, or record 12 of "climoz.nc".
     55    !   * Calendar is taken into account to get one record each day (not 360 always).
     56    !   * Missing values are filled in from sky to ground by copying lowest valid one.
     57    !     Attribute "missing_value" or "_FillValue" must be present in input file.
     58    !-------------------------------------------------------------------------------
     59    USE assert_m, ONLY : assert
     60    USE cal_tools_m, ONLY : year_len, mid_month
     61    !!  USE control_mod,        ONLY: anneeref
     62    USE time_phylmdz_mod, ONLY : annee_ref
     63    USE ioipsl, ONLY : ioget_year_len, ioget_calendar
     64    USE regr_conserv_m, ONLY : regr_conserv
     65    USE regr_lint_m, ONLY : regr_lint
     66    USE regular_lonlat_mod, ONLY : boundslon_reg, boundslat_reg, south, west, east
     67    USE slopes_m, ONLY : slopes
     68    USE lmdz_xios
     69    USE mod_phys_lmdz_para, ONLY : is_mpi_root, is_master, is_omp_master, gather, gather_mpi, bcast_mpi, klon_mpi
     70    USE geometry_mod, ONLY : latitude_deg, ind_cell_glo
     71    USE mod_grid_phy_lmdz, ONLY : klon_glo
     72
     73    !-------------------------------------------------------------------------------
     74    ! Arguments:
     75    INTEGER, INTENT(IN) :: read_climoz ! read ozone climatology, 1 or 2
     76    !                         1: read a single ozone climatology used day and night
     77    !                         2: same + read also a daylight climatology
     78    LOGICAL, INTENT(IN) :: interpt     ! TRUE  => daily interpolation
     79    ! FALSE => no interpolation (14 months)
     80    !-------------------------------------------------------------------------------
     81    ! Local variables:
     82
     83    !--- Input files variables
     84    INTEGER :: nlon_in                       ! Number of longitudes
     85    INTEGER :: nlat_in                       ! Number of latitudes
     86    INTEGER :: nlev_in                       ! Number of pressure levels
     87    INTEGER :: nmth_in                       ! Number of months
     88    REAL, ALLOCATABLE :: lon_in(:)           ! Longitudes   (ascending order, rad)
     89    REAL, ALLOCATABLE :: lat_in(:)           ! Latitudes    (ascending order, rad)
     90    REAL, ALLOCATABLE :: lev_in(:)           ! Pressure levels (ascen. order, hPa)
     91    REAL, ALLOCATABLE :: lon_in_edge(:)      ! Longitude intervals edges
     92    !              (ascending order,  / )
     93    REAL, ALLOCATABLE :: sinlat_in_edge(:)   ! Sinus of latitude intervals edges
     94    !              (ascending order,  / )
     95    LOGICAL :: ldec_lon, ldec_lat, ldec_lev  ! Decreasing order in input file
     96    CHARACTER(LEN = 20) :: cal_in              ! Calendar
     97    REAL, ALLOCATABLE :: o3_in3(:, :, :, :, :)   ! Ozone climatologies
     98    REAL, ALLOCATABLE :: o3_in3bis(:, :, :, :, :)   ! Ozone climatologies
     99    REAL, ALLOCATABLE :: o3_in2  (:, :, :, :)   ! Ozone climatologies
     100    REAL, ALLOCATABLE :: o3_in2bis(:, :, :, :, :)   ! Ozone climatologies
     101    ! last index: 1 for the day-night average, 2 for the daylight field.
     102    REAL :: NaN
     103
     104    !--- Partially or totally regridded variables      (:,:,nlev_in,:,read_climoz)
     105    REAL, ALLOCATABLE :: o3_regr_lon   (:, :, :, :, :) ! (nlon_ou,nlat_in,:,0:13   ,:)
     106    REAL, ALLOCATABLE :: o3_regr_lonlat(:, :, :, :, :) ! (nlon_ou,nlat_ou,:,0:13   ,:)
     107    REAL, ALLOCATABLE :: o3_out3       (:, :, :, :, :) ! (nlon_ou,nlat_ou,:,ntim_ou,:)
     108    REAL, ALLOCATABLE :: o3_out3_glo   (:, :, :, :) !   (nbp_lat,:,ntim_ou,:)
     109    REAL, ALLOCATABLE :: o3_regr_lat     (:, :, :, :) !         (nlat_in,:,0:13   ,:)
     110    REAL, ALLOCATABLE :: o3_out2         (:, :, :, :) !         (nlat_ou,:,ntim_ou,:)
     111    REAL, ALLOCATABLE :: o3_out2_glo     (:, :, :, :) !         (nbp_lat,:,ntim_ou,:)
     112    REAL, ALLOCATABLE :: o3_out          (:, :, :, :) !         (nbp_lat,:,ntim_ou,:)
     113    ! Dimension number  | Interval                | Contains  | For variables:
     114    !   1 (longitude)   | [rlonu(i-1), rlonu(i)]  | rlonv(i)  | all
     115    !   2 (latitude)    | [rlatv(j), rlatv(j-1)]  | rlatu(j)  | all but o3_regr_lon
     116    !   3 (press level) |                         |   lev(k)  | all
     117    ! Note that rlatv(0)=pi/2 and rlatv(nlat_ou)=-pi/2.
     118    ! Dimension 4 is: month number                             (all vars but o3_out)
     119    !                 days elapsed since Jan. 1st 0h at mid-day (o3_out only)
     120    REAL, ALLOCATABLE :: v1(:)
     121
     122    !--- For NetCDF:
     123    INTEGER :: fID_in_m, fID_in, levID_ou, dimid, vID_in(read_climoz), ntim_ou
     124    INTEGER :: fID_in_p, fID_ou, timID_ou, varid, vID_ou(read_climoz), ndims, ncerr
     125    INTEGER, ALLOCATABLE :: dIDs(:)
     126    CHARACTER(LEN = 20) :: cal_ou     !--- Calendar; no time inter => same as input
     127    CHARACTER(LEN = 80) :: press_unit !--- Pressure unit
     128    REAL :: tmidmonth(0:13)      !--- Elapsed days since Jan-1 0h at mid-months
     129    ! Additional records 0, 13 for interpolation
     130    REAL, ALLOCATABLE :: tmidday(:) !--- Output times (mid-days since Jan 1st 0h)
     131    LOGICAL :: lprev, lnext         !--- Flags: previous/next files are present
     132    LOGICAL :: l3D, l2D             !--- Flag:  input fields are 3D or zonal
     133    INTEGER :: ii, i, j, k, l, m, dln, ib, ie, iv, dx1, dx2
     134    INTEGER, ALLOCATABLE :: sta(:), cnt(:)
     135    CHARACTER(LEN = 80) :: sub, dim_nam, msg
     136    REAL :: null_array(0)
     137    LOGICAL, SAVE :: first = .TRUE.
     138    !$OMP THREADPRIVATE(first)
     139    REAL, ALLOCATABLE :: test_o3_in(:, :)
     140    REAL, ALLOCATABLE :: test_o3_out(:)
     141
     142    IF (grid_type==unstructured) THEN
     143      IF (first) THEN
     144        IF (is_master) THEN
     145          ALLOCATE(latitude_glo(klon_glo))
     146          ALLOCATE(ind_cell_glo_glo(klon_glo))
     147        ELSE
     148          ALLOCATE(latitude_glo(0))
     149          ALLOCATE(ind_cell_glo_glo(0))
     150        ENDIF
     151        CALL gather(latitude_deg, latitude_glo)
     152        CALL gather(ind_cell_glo, ind_cell_glo_glo)
    151153      ENDIF
    152       CALL gather(latitude_deg,  latitude_glo)
    153       CALL gather(ind_cell_glo,  ind_cell_glo_glo)
    154     ENDIF
    155   ENDIF
    156    
    157   IF (is_omp_master) THEN
    158     nlat_ou=nbp_lat
    159     nlon_ou=nbp_lon
    160    
    161    !-------------------------------------------------------------------------------
    162     IF (is_mpi_root) THEN
    163       sub="regr_horiz_time_climoz"
    164       WRITE(lunout,*)"Call sequence information: "//TRIM(sub)
    165       CALL assert(read_climoz == 1 .OR. read_climoz == 2, "regr_lat_time_climoz")
    166 
    167       CALL  NF95_OPEN("climoz.nc"  , NF90_NOWRITE, fID_in)
    168       lprev=NF90_OPEN("climoz_m.nc", NF90_NOWRITE, fID_in_m)==nf90_noerr
    169       lnext=NF90_OPEN("climoz_p.nc", NF90_NOWRITE, fID_in_p)==nf90_noerr
    170 
    171       !--- Get coordinates from the input file. Converts lon/lat in radians.
    172       !    Few inversions because "regr_conserv" and gcm need ascending vectors.
    173       CALL NF95_INQ_VARID(fID_in, vars_in(1), varid)
    174       CALL NF95_INQUIRE_VARIABLE(fID_in, varid, dimids=dIDs, ndims=ndims)
    175       l3D=ndims==4; l2D=ndims==3
    176       IF(l3D) WRITE(lunout,*)"Input files contain full 3D ozone fields."
    177       IF(l2D) WRITE(lunout,*)"Input files contain zonal 2D ozone fields."
    178       DO i=1,ndims
    179         CALL NF95_INQUIRE_DIMENSION(fID_in, dIDs(i), name=dim_nam, nclen=dln)
    180         CALL NF95_INQ_VARID(fID_in, dim_nam, varid)
    181         ii=i; IF(l2D) ii=i+1                              !--- ndims==3:NO LONGITUDE
    182         SELECT CASE(ii)
     154    ENDIF
     155
     156    IF (is_omp_master) THEN
     157      nlat_ou = nbp_lat
     158      nlon_ou = nbp_lon
     159
     160      !-------------------------------------------------------------------------------
     161      IF (is_mpi_root) THEN
     162        sub = "regr_horiz_time_climoz"
     163        WRITE(lunout, *)"Call sequence information: " // TRIM(sub)
     164        CALL assert(read_climoz == 1 .OR. read_climoz == 2, "regr_lat_time_climoz")
     165
     166        CALL  NF95_OPEN("climoz.nc", nf90_nowrite, fID_in)
     167        lprev = nf90_open("climoz_m.nc", nf90_nowrite, fID_in_m)==nf90_noerr
     168        lnext = nf90_open("climoz_p.nc", nf90_nowrite, fID_in_p)==nf90_noerr
     169
     170        !--- Get coordinates from the input file. Converts lon/lat in radians.
     171        !    Few inversions because "regr_conserv" and gcm need ascending vectors.
     172        CALL nf95_inq_varid(fID_in, vars_in(1), varid)
     173        CALL NF95_INQUIRE_VARIABLE(fID_in, varid, dimids = dIDs, ndims = ndims)
     174        l3D = ndims==4; l2D = ndims==3
     175        IF(l3D) WRITE(lunout, *)"Input files contain full 3D ozone fields."
     176        IF(l2D) WRITE(lunout, *)"Input files contain zonal 2D ozone fields."
     177        DO i = 1, ndims
     178          CALL nf95_inquire_dimension(fID_in, dIDs(i), name = dim_nam, nclen = dln)
     179          CALL nf95_inq_varid(fID_in, dim_nam, varid)
     180          ii = i; IF(l2D) ii = i + 1                              !--- ndims==3:NO LONGITUDE
     181          SELECT CASE(ii)
    183182          CASE(1)                                         !--- LONGITUDE
    184183            CALL NF95_GW_VAR(fID_in, varid, lon_in)
    185             ldec_lon=lon_in(1)>lon_in(dln); IF(ldec_lon) lon_in=lon_in(dln:1:-1)
    186             nlon_in=dln; lon_in=lon_in*deg2rad
     184            ldec_lon = lon_in(1)>lon_in(dln); IF(ldec_lon) lon_in = lon_in(dln:1:-1)
     185            nlon_in = dln; lon_in = lon_in * deg2rad
    187186          CASE(2)                                         !--- LATITUDE
    188187            CALL NF95_GW_VAR(fID_in, varid, lat_in)
    189             ldec_lat=lat_in(1)>lat_in(dln); IF(ldec_lat) lat_in=lat_in(dln:1:-1)
    190             nlat_in=dln; lat_in=lat_in*deg2rad
     188            ldec_lat = lat_in(1)>lat_in(dln); IF(ldec_lat) lat_in = lat_in(dln:1:-1)
     189            nlat_in = dln; lat_in = lat_in * deg2rad
    191190          CASE(3)                                         !--- PRESSURE LEVELS
    192191            CALL NF95_GW_VAR(fID_in, varid, lev_in)
    193             ldec_lev=lev_in(1)>lev_in(dln); IF(ldec_lev) lev_in=lev_in(dln:1:-1)
    194             nlev_in=dln
     192            ldec_lev = lev_in(1)>lev_in(dln); IF(ldec_lev) lev_in = lev_in(dln:1:-1)
     193            nlev_in = dln
    195194            CALL NF95_GET_ATT(fID_in, varid, "units", press_unit)
    196             k=LEN_TRIM(press_unit)
     195            k = LEN_TRIM(press_unit)
    197196            DO WHILE(ICHAR(press_unit(k:k))==0)
    198               press_unit(k:k)=' '; k=LEN_TRIM(press_unit) !--- REMOVE NULL END CHAR
     197              press_unit(k:k) = ' '; k = LEN_TRIM(press_unit) !--- REMOVE NULL END CHAR
    199198            END DO
    200199            IF(press_unit ==  "Pa") THEN
    201               lev_in = lev_in/100.                        !--- CONVERT TO hPa
     200              lev_in = lev_in / 100.                        !--- CONVERT TO hPa
    202201            ELSE IF(press_unit /= "hPa") THEN
    203               CALL abort_physic(sub, "the only recognized units are Pa and hPa.",1)
     202              CALL abort_physic(sub, "the only recognized units are Pa and hPa.", 1)
    204203            END IF
    205204          CASE(4)                                         !--- TIME
    206             CALL NF95_INQUIRE_DIMENSION(fID_in, dIDs(i), nclen=nmth_in)
    207             cal_in='gregorian'
    208             IF(NF90_GET_ATT(fID_in, varid, 'calendar', cal_in)/=nf90_noerr)        &
    209             WRITE(lunout,*)'WARNING: missing "calendar" attribute for "'//      &
    210             TRIM(dim_nam)//'" in "climoz.nc". Choosing default: "gregorian".'
    211             k=LEN_TRIM(cal_in)
     205            CALL nf95_inquire_dimension(fID_in, dIDs(i), nclen = nmth_in)
     206            cal_in = 'gregorian'
     207            IF(nf90_get_att(fID_in, varid, 'calendar', cal_in)/=nf90_noerr)        &
     208                    WRITE(lunout, *)'WARNING: missing "calendar" attribute for "' // &
     209                            TRIM(dim_nam) // '" in "climoz.nc". Choosing default: "gregorian".'
     210            k = LEN_TRIM(cal_in)
    212211            DO WHILE(ICHAR(cal_in(k:k))==0)
    213               cal_in(k:k)=' '; k=LEN_TRIM(cal_in)         !--- REMOVE NULL END CHAR
    214             END DO
    215         END SELECT
    216       END DO
    217 
    218       !--- Prepare quantities for time interpolation
    219       tmidmonth=mid_month(annee_ref, cal_in)
    220       IF(interpt) THEN
    221         ntim_ou=ioget_year_len(annee_ref)
    222         ALLOCATE(tmidday(ntim_ou))
    223         tmidday=[(REAL(k)-0.5,k=1,ntim_ou)]
    224         CALL ioget_calendar(cal_ou)
     212              cal_in(k:k) = ' '; k = LEN_TRIM(cal_in)         !--- REMOVE NULL END CHAR
     213            END DO
     214          END SELECT
     215        END DO
     216
     217        !--- Prepare quantities for time interpolation
     218        tmidmonth = mid_month(annee_ref, cal_in)
     219        IF(interpt) THEN
     220          ntim_ou = ioget_year_len(annee_ref)
     221          ALLOCATE(tmidday(ntim_ou))
     222          tmidday = [(REAL(k) - 0.5, k = 1, ntim_ou)]
     223          CALL ioget_calendar(cal_ou)
     224        ELSE
     225          ntim_ou = 14
     226          cal_ou = cal_in
     227        END IF
     228      ENDIF
     229
     230      IF (grid_type==unstructured) THEN
     231        CALL bcast_mpi(nlon_in)
     232        CALL bcast_mpi(nlat_in)
     233        CALL bcast_mpi(nlev_in)
     234        CALL bcast_mpi(l3d)
     235        CALL bcast_mpi(tmidmonth)
     236        IF(interpt) CALL bcast_mpi(tmidday)
     237        CALL bcast_mpi(ntim_ou)
     238
     239        IF (is_mpi_root) THEN
     240          CALL xios_set_domain_attr("domain_climoz", nj_glo = nlat_in, nj = nlat_in, jbegin = 0, latvalue_1d = lat_in / deg2rad)
     241          IF (l3D) THEN
     242            CALL xios_set_domain_attr("domain_climoz", ni_glo = nlon_in, ni = nlon_in, ibegin = 0, lonvalue_1d = lon_in / deg2rad)
     243          ELSE
     244            CALL xios_set_domain_attr("domain_climoz", ni_glo = 8, ni = 8, ibegin = 0, lonvalue_1d = (/ 0., 45., 90., 135., 180., 225., 270., 315. /))
     245          ENDIF
     246        ELSE
     247          CALL xios_set_domain_attr("domain_climoz", nj_glo = nlat_in, nj = 0, jbegin = 0, latvalue_1d = null_array)
     248          IF (l3D) THEN
     249            CALL xios_set_domain_attr("domain_climoz", ni_glo = nlon_in, ni = 0, ibegin = 0, lonvalue_1d = null_array)
     250          ELSE
     251            CALL xios_set_domain_attr("domain_climoz", ni_glo = 8, ni = 0, ibegin = 0, lonvalue_1d = null_array)
     252          ENDIF
     253        ENDIF
     254        CALL  xios_set_axis_attr("axis_climoz", n_glo = nlev_in)
     255        CALL  xios_set_axis_attr("time_axis_climoz", n_glo = ntim_ou)
     256        CALL  xios_set_axis_attr("time_axis_climoz", n_glo = ntim_ou)
     257        CALL  xios_set_axis_attr("tr_climoz", n_glo = read_climoz)
     258        CALL  xios_set_field_attr("tro3_out", enabled = .TRUE.)
     259        CALL  xios_set_field_attr("tro3_out", enabled = .TRUE.)
     260
     261        IF (first) THEN
     262          first = .FALSE.
     263          RETURN
     264        ENDIF
     265      ENDIF
     266
     267      IF (is_mpi_root) THEN
     268        !--- Longitudes management:
     269        !    * Need to shift data if the origin of input file longitudes /= -pi
     270        !    * Need to add some margin in longitude to ensure input interval contains
     271        !      all the output intervals => at least one longitudes slice has to be
     272        !      duplicated, possibly more for undersampling.
     273        IF(l3D) THEN
     274          IF (grid_type==unstructured) THEN
     275            dx2 = 0
     276          ELSE
     277            !--- Compute input edges longitudes vector (no end point yet)
     278            ALLOCATE(v1(nlon_in + 1))
     279            v1(1) = (lon_in(nlon_in) + lon_in(1)) / 2. - pi
     280            FORALL(i = 2:nlon_in) v1(i) = (lon_in(i - 1) + lon_in(i)) / 2.
     281            v1(nlon_in + 1) = v1(1) + 2. * pi
     282            DEALLOCATE(lon_in)
     283
     284            !--- Shift input longitudes vector until it contains first output point boundslon_reg(1,west)
     285            v1 = v1 + 2 * pi * REAL(FLOOR((boundslon_reg(1, west) - v1(1)) / (2. * pi)))
     286
     287            !--- Ensure first input longitudes interval contains first output point boundslon_reg(1,west)
     288            dx1 = locate(v1, boundslon_reg(1, west)) - 1
     289            v1 = CSHIFT(v1, SHIFT = dx1, DIM = 1)
     290            v1(nlon_in - dx1 + 2:) = v1(nlon_in - dx1 + 2:) + 2. * pi
     291
     292            !--- Extend input longitudes vector until last interval contains boundslon_reg(nlat_ou,east)
     293            dx2 = 0; DO WHILE(v1(1 + dx2) + 2. * pi<=boundslon_reg(nlon_ou, east)); dx2 = dx2 + 1;
     294            END DO
     295
     296            !--- Final edges longitudes vector (with margin and end point)
     297            ALLOCATE(lon_in_edge(nlon_in + dx2 + 1)); lon_in_edge = [v1, v1(2:1 + dx2) + 2. * pi]
     298            DEALLOCATE(v1)
     299          ENDIF
     300        END IF
     301
     302        !--- Compute sinus of intervals edges latitudes:
     303        ALLOCATE(sinlat_in_edge(nlat_in + 1))
     304        sinlat_in_edge(1) = -1. ; sinlat_in_edge(nlat_in + 1) = 1.
     305        FORALL(j = 2:nlat_in) sinlat_in_edge(j) = SIN((lat_in(j - 1) + lat_in(j)) / 2.)
     306        DEALLOCATE(lat_in)
     307
     308
     309
     310        !--- Check for contiguous years:
     311        ib = 0; ie = 13
     312        IF(nmth_in == 14) THEN; lprev = .FALSE.; lnext = .FALSE.
     313        WRITE(lunout, *)'Using 14 months ozone climatology "climoz.nc"...'
     314        ELSE
     315          IF(lprev) WRITE(lunout, *)'Using "climoz_m.nc" last record (previous year).'
     316          IF(.NOT.lprev) WRITE(lunout, *)"No previous year file ; assuming periodicity."
     317          IF(lnext) WRITE(lunout, *)'Using "climoz_p.nc" first record (next year).'
     318          IF(.NOT.lnext) WRITE(lunout, *)"No next year file ; assuming periodicity."
     319          IF(.NOT.lprev) ib = 1
     320          IF(.NOT.lnext) ie = 12
     321        END IF
     322        ALLOCATE(sta(ndims), cnt(ndims)); sta(:) = 1
     323        IF(l3D) cnt = [nlon_in, nlat_in, nlev_in, 1]
     324        IF(l2D) cnt = [        nlat_in, nlev_in, 1]
     325        IF(l3D) ALLOCATE(o3_in3(nlon_in + dx2, nlat_in, nlev_in, ib:ie, read_climoz))
     326        IF(l2D) ALLOCATE(o3_in2(nlat_in, nlev_in, ib:ie, read_climoz))
     327
     328        !--- Read full current file and one record each available contiguous file
     329        DO iv = 1, read_climoz
     330          CALL nf95_inq_varid(fID_in, vars_in(1), vID_in(iv))
     331          IF(l3D) call NF95_GET_VAR(fID_in, vID_in(iv), o3_in3(1:nlon_in, :, :, 1:12, iv))
     332          IF(l2D) call NF95_GET_VAR(fID_in, vID_in(iv), o3_in2(:, :, 1:12, iv))
     333          IF(lprev) THEN; sta(ndims) = 12
     334          CALL nf95_inq_varid(fID_in_m, vars_in(1), vID_in(iv))
     335          IF(l3D) call NF95_GET_VAR(fID_in_m, vID_in(iv), o3_in3(1:nlon_in, :, :, 0, iv), sta, cnt)
     336          IF(l2d) call NF95_GET_VAR(fID_in_m, vID_in(iv), o3_in2(:, :, 0, iv), sta, cnt)
     337          END IF
     338          IF(lnext) THEN; sta(ndims) = 1
     339          CALL nf95_inq_varid(fID_in_p, vars_in(1), vID_in(iv))
     340          IF(l3D) call NF95_GET_VAR(fID_in_p, vID_in(iv), o3_in3(1:nlon_in, :, :, 13, iv), sta, cnt)
     341          IF(l2D) call NF95_GET_VAR(fID_in_p, vID_in(iv), o3_in2(:, :, 13, iv), sta, cnt)
     342          END IF
     343        END DO
     344        IF(lprev.OR.lnext) DEALLOCATE(sta, cnt)
     345        IF(lprev) CALL NF95_CLOSE(fID_in_m)
     346        IF(lnext) CALL NF95_CLOSE(fID_in_p)
     347
     348        !--- Revert decreasing coordinates vector
     349        IF(l3D) THEN
     350          IF(ldec_lon) o3_in3(1:nlon_in, :, :, :, :) = o3_in3(nlon_in:1:-1, :, :, :, :)
     351          IF(ldec_lat) o3_in3 = o3_in3(:, nlat_in:1:-1, :, :, :)
     352          IF(ldec_lev) o3_in3 = o3_in3(:, :, nlev_in:1:-1, :, :)
     353
     354          IF (grid_type /= unstructured) THEN
     355            !--- Shift values for longitude and duplicate some longitudes slices
     356            o3_in3(1:nlon_in, :, :, :, :) = CSHIFT(o3_in3(1:nlon_in, :, :, :, :), SHIFT = dx1, DIM = 1)
     357            o3_in3(nlon_in + 1:nlon_in + dx2, :, :, :, :) = o3_in3(1:dx2, :, :, :, :)
     358          ENDIF
     359        ELSE
     360          IF(ldec_lat) o3_in2 = o3_in2(nlat_in:1:-1, :, :, :)
     361          IF(ldec_lev) o3_in2 = o3_in2(:, nlev_in:1:-1, :, :)
     362        END IF
     363
     364        !--- Deal with missing values
     365        DO m = 1, read_climoz
     366          WRITE(msg, '(a,i0)')"regr_lat_time_climoz: field Nr.", m
     367          IF(nf90_get_att(fID_in, vID_in(m), "missing_value", NaN)/= nf90_noerr) THEN
     368            IF(nf90_get_att(fID_in, vID_in(m), "_FillValue", NaN)/= nf90_noerr) THEN
     369              WRITE(lunout, *)TRIM(msg) // ": no missing value attribute found."; CYCLE
     370            END IF
     371          END IF
     372          WRITE(lunout, *)TRIM(msg) // ": missing value attribute found."
     373          WRITE(lunout, *)"Trying to fill in NaNs ; a full field would be better."
     374
     375          !--- Check top layer contains no NaNs & search NaNs from top to ground
     376          msg = TRIM(sub) // ": NaNs in top layer !"
     377          IF(l3D) THEN
     378            IF(ANY(o3_in3(:, :, 1, :, m)==NaN)) CALL abort_physic(sub, msg, 1)
     379            DO k = 2, nlev_in
     380              WHERE(o3_in3(:, :, k, :, m)==NaN) o3_in3(:, :, k, :, m) = o3_in3(:, :, k - 1, :, m)
     381            END DO
     382          ELSE
     383            IF(ANY(o3_in2(:, 1, :, m)==NaN)) THEN
     384              WRITE(lunout, *)msg
     385              !--- Fill in latitudes where all values are missing
     386              DO l = 1, nmth_in
     387                !--- Next to south pole
     388                j = 1;       DO WHILE(o3_in2(j, 1, l, m)==NaN); j = j + 1;
     389                END DO
     390                IF(j>1) &
     391                        o3_in2(:j - 1, :, l, m) = SPREAD(o3_in2(j, :, l, m), DIM = 1, ncopies = j - 1)
     392                !--- Next to north pole
     393                j = nlat_in; DO WHILE(o3_in2(j, 1, l, m)==NaN); j = j + 1;
     394                END DO
     395                IF(j<nlat_in) &
     396                        o3_in2(j + 1:, :, l, m) = SPREAD(o3_in2(j, :, l, m), DIM = 1, ncopies = nlat_in - j)
     397              END DO
     398            END IF
     399
     400            !--- Fill in high latitudes missing values
     401            !--- Highest level been filled-in, so has always valid values.
     402            DO k = 2, nlev_in
     403              WHERE(o3_in2(:, k, :, m)==NaN) o3_in2(:, k, :, m) = o3_in2(:, k - 1, :, m)
     404            END DO
     405          END IF
     406        END DO
     407
     408      ENDIF
     409
     410      !=============================================================================
     411      IF(l3D) THEN                                                   !=== 3D FIELDS
     412        !=============================================================================
     413        IF (grid_type==unstructured) THEN
     414          nlat_ou = klon_mpi
     415
     416          IF (is_mpi_root) THEN
     417            ALLOCATE(o3_in3bis(nlon_in, nlat_in, nlev_in, 0:13, read_climoz))
     418            o3_in3bis(:, :, :, ib:ie, :) = o3_in3(1:nlon_in, :, :, ib:ie, :)
     419          ELSE
     420            ALLOCATE(o3_in3bis(0, 0, 0, 0, read_climoz))
     421          ENDIF
     422          ALLOCATE(o3_regr_lonlat(1, nlat_ou, nlev_in, 0:13, read_climoz))
     423
     424          CALL xios_send_field("tro3_in", o3_in3bis(:, :, :, :, :))
     425          CALL xios_recv_field("tro3_out", o3_regr_lonlat(1, :, :, :, :))
     426        ELSE
     427
     428          !--- Regrid in longitude
     429          ALLOCATE(o3_regr_lon(nlon_ou, nlat_in, nlev_in, ie - ib + 1, read_climoz))
     430          CALL regr_conserv(1, o3_in3, xs = lon_in_edge, &
     431                  xt = [boundslon_reg(1, west), boundslon_reg(:, east)], &
     432                  vt = o3_regr_lon, slope = slopes(1, o3_in3, lon_in_edge))
     433          DEALLOCATE(o3_in3)
     434
     435          !--- Regrid in latitude: averaging with respect to SIN(lat) is
     436          !                        equivalent to weighting by COS(lat)
     437          !--- (inverted indices in "o3_regr_lonlat" because "rlatu" is decreasing)
     438          ALLOCATE(o3_regr_lonlat(nlon_ou, nlat_ou, nlev_in, 0:13, read_climoz))
     439          CALL regr_conserv(2, o3_regr_lon, xs = sinlat_in_edge, &
     440                  xt = [- 1., SIN(boundslat_reg(nlat_ou - 1:1:-1, south)), 1.], &
     441                  vt = o3_regr_lonlat(:, nlat_ou:1:- 1, :, ib:ie, :), &
     442                  slope = slopes(2, o3_regr_lon, sinlat_in_edge))
     443          DEALLOCATE(o3_regr_lon)
     444
     445        ENDIF
     446
     447        !--- Duplicate previous/next record(s) if they are not available
     448        IF(.NOT.lprev) o3_regr_lonlat(:, :, :, 0, :) = o3_regr_lonlat(:, :, :, 12, :)
     449        IF(.NOT.lnext) o3_regr_lonlat(:, :, :, 13, :) = o3_regr_lonlat(:, :, :, 1, :)
     450
     451        !--- Regrid in time by linear interpolation:
     452        ALLOCATE(o3_out3(nlon_ou, nlat_ou, nlev_in, ntim_ou, read_climoz))
     453        IF(interpt) CALL regr_lint(4, o3_regr_lonlat, tmidmonth, tmidday, o3_out3)
     454        IF(.NOT.interpt) o3_out3 = o3_regr_lonlat
     455        DEALLOCATE(o3_regr_lonlat)
     456
     457        nlat_ou = nbp_lat
     458        IF (grid_type==unstructured) THEN
     459          CALL xios_send_field('o3_out', o3_out3)
     460          ndims = 3
     461          ALLOCATE(o3_out3_glo(nlat_ou, nlev_in, ntim_ou, read_climoz))
     462          CALL gather_mpi(o3_out3(1, :, :, :, :), o3_out3_glo)
     463        ENDIF
     464
     465        !--- Create the output file and get the variable IDs:
     466        CALL prepare_out(fID_in, nlev_in, ntim_ou, fID_ou, levID_ou, timID_ou, vID_ou, &
     467                ndims, cal_ou)
     468
     469        IF (is_mpi_root) THEN
     470          !--- Write remaining coordinate variables:
     471          CALL NF95_PUT_VAR(fID_ou, levID_ou, lev_in); DEALLOCATE(lev_in)
     472          IF(interpt) CALL NF95_PUT_VAR(fID_ou, timID_ou, tmidday)
     473          IF(.NOT.interpt) CALL NF95_PUT_VAR(fID_ou, timID_ou, tmidmonth)
     474
     475          !--- Write to file (the order of "rlatu" is inverted in the output file):
     476          IF (grid_type==unstructured) THEN
     477
     478            ALLOCATE(o3_out(nlat_ou, nlev_in, ntim_ou, read_climoz))
     479            DO i = 1, klon_glo
     480              o3_out(ind_cell_glo_glo(i), :, :, :) = o3_out3_glo(i, :, :, :)
     481            ENDDO
     482
     483            DO m = 1, read_climoz
     484              CALL NF95_PUT_VAR(fID_ou, vID_ou(m), o3_out(nlat_ou:1:-1, :, :, m))
     485            END DO
     486
     487          ELSE
     488            DO m = 1, read_climoz
     489              CALL NF95_PUT_VAR(fID_ou, vID_ou(m), o3_out3(:, nlat_ou:1:-1, :, :, m))
     490            END DO
     491          ENDIF
     492          CALL NF95_CLOSE(fID_ou)
     493
     494        ENDIF
     495
     496
     497        !=============================================================================
     498      ELSE                                                         !=== ZONAL FIELDS
     499        !=============================================================================
     500
     501        IF (grid_type==unstructured) THEN
     502          nlat_ou = klon_mpi
     503
     504          IF (is_mpi_root) THEN
     505            ALLOCATE(o3_in2bis(8, nlat_in, nlev_in, 0:13, read_climoz))
     506            o3_in2bis(:, :, :, ib:ie, :) = SPREAD(o3_in2, 1, 8)
     507          ELSE
     508            ALLOCATE(o3_in2bis(0, 0, 0, 0, read_climoz))
     509          ENDIF
     510          ALLOCATE(o3_regr_lat(nlat_ou, nlev_in, 0:13, read_climoz))
     511          CALL xios_send_field("tro3_in", o3_in2bis(:, :, :, :, :))
     512          CALL xios_recv_field("tro3_out", o3_regr_lat(:, :, :, :))
     513          IF(.NOT.lprev) o3_regr_lat(:, :, 0, :) = o3_regr_lat(:, :, 12, :)
     514          IF(.NOT.lnext) o3_regr_lat(:, :, 13, :) = o3_regr_lat(:, :, 1, :)
     515
     516        ELSE
     517          !--- Regrid in latitude: averaging with respect to SIN(lat) is
     518          !                        equivalent to weighting by COS(lat)
     519          !--- (inverted indices in "o3_regr_lat" because "rlatu" is decreasing)
     520          ALLOCATE(o3_regr_lat(nlat_ou, nlev_in, 0:13, read_climoz))
     521          CALL regr_conserv(1, o3_in2, xs = sinlat_in_edge, &
     522                  xt = [- 1., SIN(boundslat_reg(nlat_ou - 1:1:-1, south)), 1.], &
     523                  vt = o3_regr_lat(nlat_ou:1:- 1, :, ib:ie, :), &
     524                  slope = slopes(1, o3_in2, sinlat_in_edge))
     525          DEALLOCATE(o3_in2)
     526
     527          !--- Duplicate previous/next record(s) if they are not available
     528          IF(.NOT.lprev) o3_regr_lat(:, :, 0, :) = o3_regr_lat(:, :, 12, :)
     529          IF(.NOT.lnext) o3_regr_lat(:, :, 13, :) = o3_regr_lat(:, :, 1, :)
     530
     531        ENDIF
     532
     533        !--- Regrid in time by linear interpolation:
     534        ALLOCATE(o3_out2(nlat_ou, nlev_in, ntim_ou, read_climoz))
     535        IF(interpt) CALL regr_lint(3, o3_regr_lat, tmidmonth, tmidday, o3_out2)
     536        IF(.NOT.interpt) o3_out2 = o3_regr_lat
     537        DEALLOCATE(o3_regr_lat)
     538
     539        nlat_ou = nbp_lat
     540
     541        IF (grid_type==unstructured) THEN
     542          ndims = 3
     543          ALLOCATE(o3_out2_glo(nlat_ou, nlev_in, ntim_ou, read_climoz))
     544          CALL gather_mpi(o3_out2, o3_out2_glo)
     545        ENDIF
     546
     547        !--- Create the output file and get the variable IDs:
     548        CALL prepare_out(fID_in, nlev_in, ntim_ou, fID_ou, levID_ou, timID_ou, vID_ou, &
     549                ndims, cal_ou)
     550
     551        IF (is_mpi_root) THEN
     552
     553          !--- Write remaining coordinate variables:
     554          CALL NF95_PUT_VAR(fID_ou, levID_ou, lev_in); DEALLOCATE(lev_in)
     555          IF(interpt) CALL NF95_PUT_VAR(fID_ou, timID_ou, tmidday)
     556          IF(.NOT.interpt) CALL NF95_PUT_VAR(fID_ou, timID_ou, tmidmonth)
     557
     558          IF (grid_type==unstructured) THEN
     559
     560            ALLOCATE(o3_out3_glo(nlat_ou, nlev_in, ntim_ou, read_climoz))
     561            DO i = 1, klon_glo
     562              o3_out(ind_cell_glo_glo(i), :, :, :) = o3_out2_glo(i, :, :, :)
     563            ENDDO
     564
     565            DO m = 1, read_climoz
     566              CALL NF95_PUT_VAR(fID_ou, vID_ou(m), o3_out(nlat_ou:1:-1, :, :, m))
     567            END DO
     568          ELSE
     569            !--- Write to file (the order of "rlatu" is inverted in the output file):
     570            DO m = 1, read_climoz
     571              CALL NF95_PUT_VAR(fID_ou, vID_ou(m), o3_out2(nlat_ou:1:-1, :, :, m))
     572            END DO
     573          ENDIF
     574
     575          CALL NF95_CLOSE(fID_ou)
     576
     577        ENDIF
     578
     579        !=============================================================================
     580      END IF
     581      !=============================================================================
     582
     583      IF (is_mpi_root) CALL NF95_CLOSE(fID_in)
     584
     585    ENDIF ! is_omp_master
     586
     587    first = .FALSE.
     588  END SUBROUTINE regr_horiz_time_climoz
     589
     590  !-------------------------------------------------------------------------------
     591
     592
     593  !-------------------------------------------------------------------------------
     594
     595  SUBROUTINE prepare_out(fID_in, nlev_in, ntim_ou, fID_ou, vlevID, vtimID, &
     596          vID_ou, ndims, cal_ou)
     597    !-------------------------------------------------------------------------------
     598    ! Purpose:  This subroutine creates the NetCDF output file, defines
     599    !     dimensions and variables, and writes some of the coordinate variables.
     600    !-------------------------------------------------------------------------------
     601    USE regular_lonlat_mod, ONLY : lon_reg, lat_reg
     602    USE regular_lonlat_mod, ONLY : lon_reg, lat_reg
     603    USE mod_phys_lmdz_para, ONLY : is_mpi_root
     604    USE mod_grid_phy_lmdz, ONLY : klon_glo
     605
     606    !-------------------------------------------------------------------------------
     607    ! Arguments:
     608    INTEGER, INTENT(IN) :: fID_in, nlev_in, ntim_ou
     609    INTEGER, INTENT(OUT) :: fID_ou, vlevID, vtimID
     610    INTEGER, INTENT(OUT) :: vID_ou(:)      ! dim(1/2) 1: O3day&night 2: O3daylight
     611    INTEGER, INTENT(IN) :: ndims          ! fields rank (3 or 4)
     612    CHARACTER(LEN = *), INTENT(IN) :: cal_ou ! calendar
     613    !-------------------------------------------------------------------------------
     614    ! Local variables:
     615    INTEGER :: dlonID, dlatID, dlevID, dtimID, dIDs(4)
     616    INTEGER :: vlonID, vlatID, ncerr, is
     617    REAL, ALLOCATABLE :: latitude_glo_(:)
     618    CHARACTER(LEN = 80) :: sub
     619    INTEGER :: i
     620
     621
     622    !-------------------------------------------------------------------------------
     623
     624    IF (is_mpi_root) THEN
     625      sub = "prepare_out"
     626      WRITE(lunout, *)"CALL sequence information: " // TRIM(sub)
     627      CALL NF95_CREATE("climoz_LMDZ.nc", nf90_clobber, fID_ou)
     628
     629      !--- Dimensions:
     630      IF(ndims==4) &
     631              CALL nf95_def_dim(fID_ou, "rlonv", nlon_ou, dlonID)
     632      CALL nf95_def_dim(fID_ou, "rlatu", nlat_ou, dlatID)
     633      CALL nf95_def_dim(fID_ou, "plev", nlev_in, dlevID)
     634      CALL nf95_def_dim(fID_ou, "time", ntim_ou, dtimID)
     635
     636      !--- Define coordinate variables:
     637      IF(ndims==4) &
     638              CALL nf95_def_var(fID_ou, "rlonv", nf90_float, dlonID, vlonID)
     639      CALL nf95_def_var(fID_ou, "rlatu", nf90_float, dlatID, vlatID)
     640      CALL nf95_def_var(fID_ou, "plev", nf90_float, dlevID, vlevID)
     641      CALL nf95_def_var(fID_ou, "time", nf90_float, dtimID, vtimID)
     642      IF(ndims==4) &
     643              CALL NF95_PUT_ATT(fID_ou, vlonID, "units", "degrees_east")
     644      CALL NF95_PUT_ATT(fID_ou, vlatID, "units", "degrees_north")
     645      CALL NF95_PUT_ATT(fID_ou, vlevID, "units", "millibar")
     646      CALL NF95_PUT_ATT(fID_ou, vtimID, "units", "days since 2000-1-1")
     647      IF(ndims==4) &
     648              CALL NF95_PUT_ATT(fID_ou, vlonID, "standard_name", "longitude")
     649      CALL NF95_PUT_ATT(fID_ou, vlatID, "standard_name", "latitude")
     650      CALL NF95_PUT_ATT(fID_ou, vlevID, "standard_name", "air_pressure")
     651      CALL NF95_PUT_ATT(fID_ou, vtimID, "standard_name", "time")
     652      CALL NF95_PUT_ATT(fID_ou, vlevID, "long_name", "air pressure")
     653      CALL NF95_PUT_ATT(fID_ou, vtimID, "calendar", cal_ou)
     654
     655      !--- Define the main variables:
     656      IF(ndims==3) dIDs(1:3) = [ dlatID, dlevID, dtimID]
     657      IF(ndims==4) dIDs = [dlonID, dlatID, dlevID, dtimID]
     658      CALL nf95_def_var(fID_ou, vars_in(1), nf90_float, dIDs(1:ndims), vID_ou(1))
     659      CALL NF95_PUT_ATT(fID_ou, vID_ou(1), "long_name", "ozone mole fraction")
     660      CALL NF95_PUT_ATT(fID_ou, vID_ou(1), "standard_name", "mole_fraction_of_ozone&
     661              _in_air")
     662      IF(SIZE(vID_ou) == 2) THEN
     663        CALL nf95_def_var(fID_ou, vars_in(2), nf90_float, dIDs(1:ndims), vID_ou(2))
     664        CALL NF95_PUT_ATT(fID_ou, vID_ou(2), "long_name", "ozone mole fraction in da&
     665                ylight")
     666      END IF
     667
     668      !--- Global attributes:
     669      ! The following commands, copying attributes, may fail. That is OK.
     670      ! It should just mean that the attribute is not defined in the input file.
     671      CALL NF95_COPY_ATT(fID_in, nf90_global, "Conventions", fID_ou, nf90_global, ncerr)
     672      CALL handle_err_copy_att("Conventions")
     673      CALL NF95_COPY_ATT(fID_in, nf90_global, "title", fID_ou, nf90_global, ncerr)
     674      CALL handle_err_copy_att("title")
     675      CALL NF95_COPY_ATT(fID_in, nf90_global, "institution", fID_ou, nf90_global, ncerr)
     676      CALL handle_err_copy_att("institution")
     677      CALL NF95_COPY_ATT(fID_in, nf90_global, "source", fID_ou, nf90_global, ncerr)
     678      CALL handle_err_copy_att("source")
     679      CALL NF95_PUT_ATT (fID_ou, nf90_global, "comment", "Regridded for LMDZ")
     680      CALL NF95_ENDDEF(fID_ou)
     681
     682      IF (grid_type==unstructured) THEN
     683        ALLOCATE(latitude_glo_(klon_glo))
     684        DO i = 1, klon_glo
     685          latitude_glo_(ind_cell_glo_glo(i)) = latitude_glo(i)
     686        ENDDO
     687        CALL NF95_PUT_VAR(fID_ou, vlatID, latitude_glo_)
    225688      ELSE
    226         ntim_ou=14
    227         cal_ou=cal_in
    228       END IF
    229     ENDIF
    230 
    231     IF (grid_type==unstructured) THEN
    232       CALL bcast_mpi(nlon_in)
    233       CALL bcast_mpi(nlat_in)
    234       CALL bcast_mpi(nlev_in)
    235       CALL bcast_mpi(l3d)
    236       CALL bcast_mpi(tmidmonth)
    237       IF(interpt) CALL bcast_mpi(tmidday)
    238       CALL bcast_mpi(ntim_ou)
    239 
    240       IF (is_mpi_root) THEN
    241         CALL xios_set_domain_attr("domain_climoz",nj_glo=nlat_in, nj=nlat_in, jbegin=0, latvalue_1d=lat_in/deg2rad)
    242         IF (l3D) THEN
    243           CALL xios_set_domain_attr("domain_climoz",ni_glo=nlon_in, ni=nlon_in, ibegin=0, lonvalue_1d=lon_in/deg2rad)
    244         ELSE
    245           CALL xios_set_domain_attr("domain_climoz",ni_glo=8, ni=8, ibegin=0, lonvalue_1d = (/ 0.,45.,90.,135.,180.,225.,270., 315. /))
    246         ENDIF
    247       ELSE
    248         CALL xios_set_domain_attr("domain_climoz",nj_glo=nlat_in, nj=0, jbegin=0, latvalue_1d=null_array )
    249         IF (l3D) THEN
    250           CALL xios_set_domain_attr("domain_climoz",ni_glo=nlon_in, ni=0, ibegin=0, lonvalue_1d=null_array)
    251         ELSE
    252           CALL xios_set_domain_attr("domain_climoz",ni_glo=8, ni=0, ibegin=0, lonvalue_1d=null_array)
    253         ENDIF
    254       ENDIF
    255       CALL  xios_set_axis_attr("axis_climoz", n_glo=nlev_in)
    256       CALL  xios_set_axis_attr("time_axis_climoz", n_glo=ntim_ou)
    257       CALL  xios_set_axis_attr("time_axis_climoz", n_glo=ntim_ou)
    258       CALL  xios_set_axis_attr("tr_climoz", n_glo=read_climoz)
    259       CALL  xios_set_field_attr("tro3_out", enabled=.TRUE.)
    260       CALL  xios_set_field_attr("tro3_out", enabled=.TRUE.)
    261      
    262       IF (first) THEN
    263         first=.FALSE.
    264         RETURN
     689        !--- Write one of the coordinate variables:
     690        IF(ndims==4) CALL NF95_PUT_VAR(fID_ou, vlonID, lon_reg / deg2rad)
     691        CALL NF95_PUT_VAR(fID_ou, vlatID, lat_reg(nlat_ou:1:-1) / deg2rad)
     692        !    (convert from rad to degrees and sort in ascending order)
    265693      ENDIF
    266694    ENDIF
    267    
    268    
    269     IF (is_mpi_root) THEN     
    270       !--- Longitudes management:
    271       !    * Need to shift data if the origin of input file longitudes /= -pi
    272       !    * Need to add some margin in longitude to ensure input interval contains
    273       !      all the output intervals => at least one longitudes slice has to be
    274       !      duplicated, possibly more for undersampling.
    275       IF(l3D) THEN
    276         IF (grid_type==unstructured) THEN
    277           dx2=0
    278         ELSE
    279           !--- Compute input edges longitudes vector (no end point yet)
    280           ALLOCATE(v1(nlon_in+1))
    281           v1(1)=(lon_in(nlon_in)+lon_in(1))/2.-pi
    282           FORALL(i=2:nlon_in) v1(i)=(lon_in(i-1)+lon_in(i))/2.
    283           v1(nlon_in+1)=v1(1)+2.*pi
    284           DEALLOCATE(lon_in)
    285 
    286           !--- Shift input longitudes vector until it contains first output point boundslon_reg(1,west)
    287           v1=v1+2*pi*REAL(FLOOR((boundslon_reg(1,west)-v1(1))/(2.*pi)))
    288 
    289           !--- Ensure first input longitudes interval contains first output point boundslon_reg(1,west)
    290           dx1=locate(v1,boundslon_reg(1,west))-1
    291           v1=CSHIFT(v1,SHIFT=dx1,DIM=1)
    292           v1(nlon_in-dx1+2:)=v1(nlon_in-dx1+2:)+2.*pi
    293    
    294           !--- Extend input longitudes vector until last interval contains boundslon_reg(nlat_ou,east)
    295           dx2=0; DO WHILE(v1(1+dx2)+2.*pi<=boundslon_reg(nlon_ou,east)); dx2=dx2+1; END DO
    296 
    297           !--- Final edges longitudes vector (with margin and end point)
    298           ALLOCATE(lon_in_edge(nlon_in+dx2+1)); lon_in_edge=[v1,v1(2:1+dx2)+2.*pi]
    299           DEALLOCATE(v1)
    300         ENDIF
    301       END IF
    302 
    303       !--- Compute sinus of intervals edges latitudes:
    304       ALLOCATE(sinlat_in_edge(nlat_in+1))
    305       sinlat_in_edge(1) = -1. ; sinlat_in_edge(nlat_in+1) = 1.
    306       FORALL(j=2:nlat_in) sinlat_in_edge(j)=SIN((lat_in(j-1)+lat_in(j))/2.)
    307       DEALLOCATE(lat_in)
    308 
    309 
    310 
    311       !--- Check for contiguous years:
    312       ib=0; ie=13
    313       IF(nmth_in == 14) THEN; lprev=.FALSE.; lnext=.FALSE.
    314         WRITE(lunout,*)'Using 14 months ozone climatology "climoz.nc"...'
    315       ELSE 
    316         IF(     lprev) WRITE(lunout,*)'Using "climoz_m.nc" last record (previous year).'
    317         IF(.NOT.lprev) WRITE(lunout,*)"No previous year file ; assuming periodicity."
    318         IF(     lnext) WRITE(lunout,*)'Using "climoz_p.nc" first record (next year).'
    319         IF(.NOT.lnext) WRITE(lunout,*)"No next year file ; assuming periodicity."
    320         IF(.NOT.lprev) ib=1
    321         IF(.NOT.lnext) ie=12
    322       END IF
    323       ALLOCATE(sta(ndims),cnt(ndims)); sta(:)=1 
    324       IF(l3D) cnt=[nlon_in,nlat_in,nlev_in,1]
    325       IF(l2D) cnt=[        nlat_in,nlev_in,1] 
    326       IF(l3D) ALLOCATE(o3_in3(nlon_in+dx2,nlat_in,nlev_in,ib:ie,read_climoz))
    327       IF(l2D) ALLOCATE(o3_in2(            nlat_in,nlev_in,ib:ie,read_climoz))
    328 
    329       !--- Read full current file and one record each available contiguous file
    330       DO iv=1,read_climoz
    331         CALL NF95_INQ_VARID(fID_in, vars_in(1), vID_in(iv))
    332         IF(l3D) call NF95_GET_VAR(fID_in, vID_in(iv), o3_in3(1:nlon_in,:,:,1:12,iv))
    333         IF(l2D) call NF95_GET_VAR(fID_in, vID_in(iv), o3_in2(          :,:,1:12,iv))
    334         IF(lprev) THEN; sta(ndims)=12 
    335           CALL NF95_INQ_VARID(fID_in_m, vars_in(1), vID_in(iv))
    336           IF(l3D) call NF95_GET_VAR(fID_in_m,vID_in(iv),o3_in3(1:nlon_in,:,:, 0,iv),sta,cnt)
    337           IF(l2d) call NF95_GET_VAR(fID_in_m,vID_in(iv),o3_in2(          :,:, 0,iv),sta,cnt)
    338         END IF
    339         IF(lnext) THEN; sta(ndims)=1 
    340           CALL NF95_INQ_VARID(fID_in_p, vars_in(1), vID_in(iv))
    341           IF(l3D) call NF95_GET_VAR(fID_in_p,vID_in(iv),o3_in3(1:nlon_in,:,:,13,iv),sta,cnt)
    342           IF(l2D) call NF95_GET_VAR(fID_in_p,vID_in(iv),o3_in2(          :,:,13,iv),sta,cnt)
    343         END IF
    344       END DO
    345       IF(lprev.OR.lnext) DEALLOCATE(sta,cnt)
    346       IF(lprev) CALL NF95_CLOSE(fID_in_m)
    347       IF(lnext) CALL NF95_CLOSE(fID_in_p)
    348 
    349       !--- Revert decreasing coordinates vector
    350       IF(l3D) THEN
    351         IF(ldec_lon) o3_in3(1:nlon_in,:,:,:,:) = o3_in3(nlon_in:1:-1,:,:,:,:)
    352         IF(ldec_lat) o3_in3 = o3_in3(:,nlat_in:1:-1,:,:,:)
    353         IF(ldec_lev) o3_in3 = o3_in3(:,:,nlev_in:1:-1,:,:)
    354        
    355         IF (grid_type /= unstructured) THEN
    356           !--- Shift values for longitude and duplicate some longitudes slices
    357           o3_in3(1:nlon_in,:,:,:,:)=CSHIFT(o3_in3(1:nlon_in,:,:,:,:),SHIFT=dx1,DIM=1)
    358           o3_in3(nlon_in+1:nlon_in+dx2,:,:,:,:)=o3_in3(1:dx2,:,:,:,:)
    359         ENDIF
    360       ELSE
    361         IF(ldec_lat) o3_in2 = o3_in2(  nlat_in:1:-1,:,:,:)
    362         IF(ldec_lev) o3_in2 = o3_in2(  :,nlev_in:1:-1,:,:)
    363       END IF
    364 
    365      !--- Deal with missing values
    366       DO m=1, read_climoz
    367         WRITE(msg,'(a,i0)')"regr_lat_time_climoz: field Nr.",m
    368         IF(NF90_GET_ATT(fID_in,vID_in(m),"missing_value",NaN)/= nf90_noerr) THEN
    369           IF(NF90_GET_ATT(fID_in, vID_in(m),"_FillValue",NaN)/= nf90_noerr) THEN
    370             WRITE(lunout,*)TRIM(msg)//": no missing value attribute found."; CYCLE
    371           END IF
    372         END IF
    373         WRITE(lunout,*)TRIM(msg)//": missing value attribute found."
    374         WRITE(lunout,*)"Trying to fill in NaNs ; a full field would be better."
    375 
    376         !--- Check top layer contains no NaNs & search NaNs from top to ground
    377         msg=TRIM(sub)//": NaNs in top layer !"
    378         IF(l3D) THEN
    379           IF(ANY(o3_in3(:,:,1,:,m)==NaN)) CALL abort_physic(sub,msg,1)
    380           DO k = 2,nlev_in
    381             WHERE(o3_in3(:,:,k,:,m)==NaN) o3_in3(:,:,k,:,m)=o3_in3(:,:,k-1,:,m)
    382           END DO
    383         ELSE
    384           IF(ANY(o3_in2(  :,1,:,m)==NaN)) THEN
    385             WRITE(lunout,*)msg
    386             !--- Fill in latitudes where all values are missing
    387             DO l=1,nmth_in
    388               !--- Next to south pole
    389               j=1;       DO WHILE(o3_in2(j,1,l,m)==NaN); j=j+1; END DO
    390               IF(j>1) &
    391                 o3_in2(:j-1,:,l,m)=SPREAD(o3_in2(j,:,l,m),DIM=1,ncopies=j-1)
    392               !--- Next to north pole
    393               j=nlat_in; DO WHILE(o3_in2(j,1,l,m)==NaN); j=j+1; END DO
    394               IF(j<nlat_in) &
    395                 o3_in2(j+1:,:,l,m)=SPREAD(o3_in2(j,:,l,m),DIM=1,ncopies=nlat_in-j)
    396             END DO
    397           END IF
    398 
    399           !--- Fill in high latitudes missing values
    400           !--- Highest level been filled-in, so has always valid values.
    401           DO k = 2,nlev_in
    402             WHERE(o3_in2(:,k,:,m)==NaN) o3_in2(:,k,:,m)=o3_in2(:,k-1,:,m)
    403           END DO
    404         END IF
    405       END DO
    406 
    407     ENDIF
    408    
    409     !=============================================================================
    410     IF(l3D) THEN                                                   !=== 3D FIELDS
    411     !=============================================================================
    412      IF (grid_type==unstructured) THEN
    413        nlat_ou=klon_mpi
    414        
    415        IF (is_mpi_root) THEN
    416          ALLOCATE(o3_in3bis(nlon_in,nlat_in,nlev_in,0:13,read_climoz))
    417          o3_in3bis(:,:,:,ib:ie,:)=o3_in3(1:nlon_in,:,:,ib:ie,:)
    418        ELSE
    419          ALLOCATE(o3_in3bis(0,0,0,0,read_climoz))
    420        ENDIF
    421        ALLOCATE(o3_regr_lonlat(1, nlat_ou, nlev_in, 0:13, read_climoz))
    422        
    423        CALL xios_send_field("tro3_in",o3_in3bis(:,:,:,:,:))
    424        CALL xios_recv_field("tro3_out",o3_regr_lonlat(1,:,:,:,:))
    425      ELSE
    426          
    427        !--- Regrid in longitude
    428         ALLOCATE(o3_regr_lon(nlon_ou, nlat_in, nlev_in, ie-ib+1, read_climoz))
    429         CALL regr_conserv(1, o3_in3, xs = lon_in_edge,                             &
    430                             xt = [boundslon_reg(1,west),boundslon_reg(:,east)],    &
    431                             vt = o3_regr_lon, slope = slopes(1,o3_in3, lon_in_edge))
    432         DEALLOCATE(o3_in3)
    433 
    434         !--- Regrid in latitude: averaging with respect to SIN(lat) is
    435         !                        equivalent to weighting by COS(lat)
    436         !--- (inverted indices in "o3_regr_lonlat" because "rlatu" is decreasing)
    437         ALLOCATE(o3_regr_lonlat(nlon_ou, nlat_ou, nlev_in, 0:13, read_climoz))
    438         CALL regr_conserv(2, o3_regr_lon, xs = sinlat_in_edge,                     &
    439                         xt = [- 1., SIN(boundslat_reg(nlat_ou-1:1:-1,south)), 1.], &
    440                         vt = o3_regr_lonlat(:,nlat_ou:1:- 1,:,ib:ie,:),            &
    441                    slope = slopes(2,o3_regr_lon, sinlat_in_edge))
    442         DEALLOCATE(o3_regr_lon)
    443 
    444      ENDIF
    445 
    446      !--- Duplicate previous/next record(s) if they are not available
    447      IF(.NOT.lprev) o3_regr_lonlat(:,:,:, 0,:) = o3_regr_lonlat(:,:,:,12,:)
    448      IF(.NOT.lnext) o3_regr_lonlat(:,:,:,13,:) = o3_regr_lonlat(:,:,:, 1,:)
    449      
    450      !--- Regrid in time by linear interpolation:
    451      ALLOCATE(o3_out3(nlon_ou, nlat_ou, nlev_in, ntim_ou, read_climoz))
    452      IF(     interpt) CALL regr_lint(4,o3_regr_lonlat,tmidmonth,tmidday,o3_out3)
    453      IF(.NOT.interpt) o3_out3=o3_regr_lonlat
    454      DEALLOCATE(o3_regr_lonlat)
    455 
    456      nlat_ou=nbp_lat
    457      IF (grid_type==unstructured) THEN
    458        CALL xios_send_field('o3_out',o3_out3)
    459        ndims=3
    460        ALLOCATE(o3_out3_glo(nlat_ou, nlev_in, ntim_ou, read_climoz))
    461        CALL gather_mpi(o3_out3(1,:,:,:,:), o3_out3_glo)
    462      ENDIF
    463 
    464     !--- Create the output file and get the variable IDs:
    465     CALL prepare_out(fID_in,nlev_in,ntim_ou, fID_ou,levID_ou,timID_ou,vID_ou, &
    466                      ndims, cal_ou)
    467 
    468     IF (is_mpi_root) THEN
    469       !--- Write remaining coordinate variables:
    470       CALL NF95_PUT_VAR(fID_ou, levID_ou, lev_in); DEALLOCATE(lev_in)
    471       IF(     interpt) CALL NF95_PUT_VAR(fID_ou, timID_ou, tmidday)
    472       IF(.NOT.interpt) CALL NF95_PUT_VAR(fID_ou, timID_ou, tmidmonth)
    473 
    474       !--- Write to file (the order of "rlatu" is inverted in the output file):
    475         IF (grid_type==unstructured) THEN
    476 
    477           ALLOCATE(o3_out(nlat_ou, nlev_in, ntim_ou, read_climoz))
    478           DO i=1,klon_glo
    479             o3_out(ind_cell_glo_glo(i),:,:,:)=o3_out3_glo(i,:,:,:)
    480           ENDDO
    481 
    482           DO m = 1, read_climoz
    483             CALL NF95_PUT_VAR(fID_ou, vID_ou(m), o3_out(nlat_ou:1:-1,:,:,m))
    484           END DO
    485          
    486         ELSE
    487           DO m = 1, read_climoz
    488             CALL NF95_PUT_VAR(fID_ou, vID_ou(m), o3_out3(:,nlat_ou:1:-1,:,:,m))
    489           END DO
    490       ENDIF
    491       CALL NF95_CLOSE(fID_ou)
    492 
    493 
    494     ENDIF
    495 
    496 
    497     !=============================================================================
    498     ELSE                                                         !=== ZONAL FIELDS
    499     !=============================================================================
    500    
    501      IF (grid_type==unstructured) THEN
    502        nlat_ou=klon_mpi
    503 
    504        IF (is_mpi_root) THEN
    505          ALLOCATE(o3_in2bis(8,nlat_in,nlev_in,0:13,read_climoz))
    506          o3_in2bis(:,:,:,ib:ie,:)=SPREAD(o3_in2,1,8)
    507        ELSE
    508          ALLOCATE(o3_in2bis(0,0,0,0,read_climoz))
    509        ENDIF
    510        ALLOCATE(o3_regr_lat(nlat_ou, nlev_in, 0:13, read_climoz))
    511        CALL xios_send_field("tro3_in",o3_in2bis(:,:,:,:,:))
    512        CALL xios_recv_field("tro3_out",o3_regr_lat(:,:,:,:))
    513        IF(.NOT.lprev) o3_regr_lat(:,:, 0, :) = o3_regr_lat(:,:,12,:)
    514        IF(.NOT.lnext) o3_regr_lat(:,:,13, :) = o3_regr_lat(:,:, 1,:)
    515      
    516      ELSE
    517         !--- Regrid in latitude: averaging with respect to SIN(lat) is
    518         !                        equivalent to weighting by COS(lat)
    519         !--- (inverted indices in "o3_regr_lat" because "rlatu" is decreasing)
    520         ALLOCATE(o3_regr_lat(nlat_ou, nlev_in, 0:13, read_climoz))
    521         CALL regr_conserv(1, o3_in2, xs = sinlat_in_edge,                          &
    522                         xt = [- 1., SIN(boundslat_reg(nlat_ou-1:1:-1,south)), 1.], &
    523                         vt = o3_regr_lat(nlat_ou:1:- 1,:,ib:ie,:),                 &
    524                      slope = slopes(1,o3_in2, sinlat_in_edge))
    525         DEALLOCATE(o3_in2)
    526 
    527         !--- Duplicate previous/next record(s) if they are not available
    528         IF(.NOT.lprev) o3_regr_lat(:,:, 0,:) = o3_regr_lat(:,:,12,:)
    529         IF(.NOT.lnext) o3_regr_lat(:,:,13,:) = o3_regr_lat(:,:, 1,:)
    530 
    531      ENDIF
    532      
    533       !--- Regrid in time by linear interpolation:
    534       ALLOCATE(o3_out2(nlat_ou, nlev_in, ntim_ou, read_climoz))
    535       IF(     interpt) CALL regr_lint(3,o3_regr_lat, tmidmonth, tmidday, o3_out2)
    536       IF(.NOT.interpt) o3_out2=o3_regr_lat
    537       DEALLOCATE(o3_regr_lat)
    538 
    539       nlat_ou=nbp_lat
    540    
    541       IF (grid_type==unstructured) THEN
    542         ndims=3
    543         ALLOCATE(o3_out2_glo(nlat_ou, nlev_in, ntim_ou, read_climoz))
    544         CALL gather_mpi(o3_out2, o3_out2_glo)
    545       ENDIF
    546      
    547       !--- Create the output file and get the variable IDs:
    548       CALL prepare_out(fID_in,nlev_in,ntim_ou, fID_ou,levID_ou,timID_ou,vID_ou, &
    549                          ndims, cal_ou)
    550 
    551       IF (is_mpi_root) THEN
    552      
    553         !--- Write remaining coordinate variables:
    554         CALL NF95_PUT_VAR(fID_ou, levID_ou, lev_in); DEALLOCATE(lev_in)
    555         IF(     interpt) CALL NF95_PUT_VAR(fID_ou, timID_ou, tmidday)
    556         IF(.NOT.interpt) CALL NF95_PUT_VAR(fID_ou, timID_ou, tmidmonth)
    557 
    558         IF (grid_type==unstructured) THEN
    559 
    560           ALLOCATE(o3_out3_glo(nlat_ou, nlev_in, ntim_ou, read_climoz))
    561           DO i=1,klon_glo
    562             o3_out(ind_cell_glo_glo(i),:,:,:)=o3_out2_glo(i,:,:,:)
    563           ENDDO
    564 
    565 
    566           DO m = 1, read_climoz
    567             CALL NF95_PUT_VAR(fID_ou, vID_ou(m), o3_out(nlat_ou:1:-1,:,:,m))
    568           END DO
    569         ELSE
    570           !--- Write to file (the order of "rlatu" is inverted in the output file):
    571           DO m = 1, read_climoz
    572             CALL NF95_PUT_VAR(fID_ou, vID_ou(m), o3_out2(nlat_ou:1:-1,:,:,m))
    573           END DO
    574         ENDIF
    575        
    576         CALL NF95_CLOSE(fID_ou)
    577      
    578       ENDIF
    579 
    580     !=============================================================================
    581     END IF
    582     !=============================================================================
    583 
    584     IF (is_mpi_root) CALL NF95_CLOSE(fID_in)
    585 
    586   ENDIF ! is_omp_master
    587 
    588   first=.FALSE.
    589 END SUBROUTINE regr_horiz_time_climoz
     695
     696  CONTAINS
     697
     698    !-------------------------------------------------------------------------------
     699
     700    SUBROUTINE handle_err_copy_att(att_name)
     701
     702      !-------------------------------------------------------------------------------
     703      USE netcdf, ONLY : nf90_noerr, NF90_strerror
     704      !-------------------------------------------------------------------------------
     705      ! Arguments:
     706      CHARACTER(LEN = *), INTENT(IN) :: att_name
     707      !-------------------------------------------------------------------------------
     708      IF(ncerr /= nf90_noerr) &
     709              WRITE(lunout, *)TRIM(sub) // " prepare_out NF95_COPY_ATT " // TRIM(att_name) // &
     710                      " -- " // TRIM(NF90_strerror(ncerr))
     711
     712    END SUBROUTINE handle_err_copy_att
     713
     714    !-------------------------------------------------------------------------------
     715
     716  END SUBROUTINE prepare_out
     717
     718  !-------------------------------------------------------------------------------
     719
     720END MODULE regr_horiz_time_climoz_m
    590721
    591722!-------------------------------------------------------------------------------
    592 
    593 
    594 !-------------------------------------------------------------------------------
    595 
    596 SUBROUTINE prepare_out(fID_in, nlev_in, ntim_ou, fID_ou, vlevID, vtimID, &
    597                        vID_ou, ndims, cal_ou)
    598 !-------------------------------------------------------------------------------
    599 ! Purpose:  This subroutine creates the NetCDF output file, defines
    600 !     dimensions and variables, and writes some of the coordinate variables.
    601 !-------------------------------------------------------------------------------
    602   USE regular_lonlat_mod, ONLY: lon_reg, lat_reg
    603   USE regular_lonlat_mod, ONLY: lon_reg, lat_reg
    604   USE mod_phys_lmdz_para, ONLY: is_mpi_root
    605   USE mod_grid_phy_lmdz, ONLY: klon_glo
    606 
    607 !-------------------------------------------------------------------------------
    608 ! Arguments:
    609   INTEGER, INTENT(IN)  :: fID_in, nlev_in, ntim_ou
    610   INTEGER, INTENT(OUT) :: fID_ou, vlevID,  vtimID
    611   INTEGER, INTENT(OUT) :: vID_ou(:)      ! dim(1/2) 1: O3day&night 2: O3daylight
    612   INTEGER, INTENT(IN)  :: ndims          ! fields rank (3 or 4)
    613   CHARACTER(LEN=*), INTENT(IN) :: cal_ou ! calendar
    614 !-------------------------------------------------------------------------------
    615 ! Local variables:
    616   INTEGER :: dlonID, dlatID, dlevID, dtimID, dIDs(4)
    617   INTEGER :: vlonID, vlatID, ncerr,  is
    618   REAL,ALLOCATABLE    :: latitude_glo_(:)
    619   CHARACTER(LEN=80) :: sub
    620   INTEGER :: i
    621 
    622 
    623 !-------------------------------------------------------------------------------
    624  
    625   IF (is_mpi_root) THEN 
    626     sub="prepare_out"
    627     WRITE(lunout,*)"CALL sequence information: "//TRIM(sub)
    628     CALL NF95_CREATE("climoz_LMDZ.nc", NF90_clobber, fID_ou)
    629 
    630   !--- Dimensions:
    631     IF(ndims==4) &
    632     CALL NF95_DEF_DIM(fID_ou, "rlonv", nlon_ou, dlonID)
    633     CALL NF95_DEF_DIM(fID_ou, "rlatu", nlat_ou, dlatID)
    634     CALL NF95_DEF_DIM(fID_ou, "plev",  nlev_in, dlevID)
    635     CALL NF95_DEF_DIM(fID_ou, "time",  ntim_ou, dtimID)
    636 
    637     !--- Define coordinate variables:
    638     IF(ndims==4) &
    639     CALL NF95_DEF_VAR(fID_ou, "rlonv", NF90_FLOAT, dlonID, vlonID)
    640     CALL NF95_DEF_VAR(fID_ou, "rlatu", NF90_FLOAT, dlatID, vlatID)
    641     CALL NF95_DEF_VAR(fID_ou, "plev",  NF90_FLOAT, dlevID, vlevID)
    642     CALL NF95_DEF_VAR(fID_ou, "time",  NF90_FLOAT, dtimID, vtimID)
    643     IF(ndims==4) &
    644     CALL NF95_PUT_ATT(fID_ou, vlonID, "units", "degrees_east")
    645     CALL NF95_PUT_ATT(fID_ou, vlatID, "units", "degrees_north")
    646     CALL NF95_PUT_ATT(fID_ou, vlevID, "units", "millibar")
    647     CALL NF95_PUT_ATT(fID_ou, vtimID, "units", "days since 2000-1-1")
    648     IF(ndims==4) &
    649     CALL NF95_PUT_ATT(fID_ou, vlonID, "standard_name", "longitude")
    650     CALL NF95_PUT_ATT(fID_ou, vlatID, "standard_name", "latitude")
    651     CALL NF95_PUT_ATT(fID_ou, vlevID, "standard_name", "air_pressure")
    652     CALL NF95_PUT_ATT(fID_ou, vtimID, "standard_name", "time")
    653     CALL NF95_PUT_ATT(fID_ou, vlevID, "long_name",     "air pressure")
    654     CALL NF95_PUT_ATT(fID_ou, vtimID, "calendar",      cal_ou)
    655 
    656   !--- Define the main variables:
    657     IF(ndims==3) dIDs(1:3) = [ dlatID, dlevID, dtimID]
    658     IF(ndims==4) dIDs=[dlonID, dlatID, dlevID, dtimID]
    659     CALL NF95_DEF_VAR(fID_ou, vars_in(1), NF90_FLOAT, dIDs(1:ndims), vID_ou(1))
    660     CALL NF95_PUT_ATT(fID_ou, vID_ou(1), "long_name", "ozone mole fraction")
    661     CALL NF95_PUT_ATT(fID_ou, vID_ou(1), "standard_name", "mole_fraction_of_ozone&
    662   _in_air")
    663     IF(SIZE(vID_ou) == 2) THEN
    664       CALL NF95_DEF_VAR(fID_ou, vars_in(2), NF90_FLOAT, dIDs(1:ndims), vID_ou(2))
    665       CALL NF95_PUT_ATT(fID_ou, vID_ou(2), "long_name","ozone mole fraction in da&
    666   ylight")
    667     END IF
    668 
    669   !--- Global attributes:
    670   ! The following commands, copying attributes, may fail. That is OK.
    671   ! It should just mean that the attribute is not defined in the input file.
    672     CALL NF95_COPY_ATT(fID_in,NF90_GLOBAL,"Conventions",fID_ou,NF90_GLOBAL, ncerr)
    673     CALL handle_err_copy_att("Conventions")
    674     CALL NF95_COPY_ATT(fID_in,NF90_GLOBAL,"title",      fID_ou,NF90_GLOBAL, ncerr)
    675     CALL handle_err_copy_att("title")
    676     CALL NF95_COPY_ATT(fID_in,NF90_GLOBAL,"institution",fID_ou,NF90_GLOBAL, ncerr)
    677     CALL handle_err_copy_att("institution")
    678     CALL NF95_COPY_ATT(fID_in,NF90_GLOBAL,"source",     fID_ou,NF90_GLOBAL, ncerr)
    679     CALL handle_err_copy_att("source")
    680     CALL NF95_PUT_ATT (fID_ou,NF90_GLOBAL,"comment", "Regridded for LMDZ")
    681     CALL NF95_ENDDEF(fID_ou)
    682 
    683     IF (grid_type==unstructured) THEN
    684       ALLOCATE(latitude_glo_(klon_glo))
    685       DO i=1,klon_glo
    686         latitude_glo_(ind_cell_glo_glo(i))=latitude_glo(i)
    687       ENDDO
    688       CALL NF95_PUT_VAR(fID_ou, vlatID, latitude_glo_)
    689     ELSE
    690       !--- Write one of the coordinate variables:
    691       IF(ndims==4) CALL NF95_PUT_VAR(fID_ou, vlonID, lon_reg/deg2rad)
    692       CALL NF95_PUT_VAR(fID_ou, vlatID, lat_reg(nlat_ou:1:-1)/deg2rad)
    693     !    (convert from rad to degrees and sort in ascending order)
    694     ENDIF
    695   ENDIF
    696  
    697 CONTAINS
    698 
    699 !-------------------------------------------------------------------------------
    700 
    701 SUBROUTINE handle_err_copy_att(att_name)
    702 
    703 !-------------------------------------------------------------------------------
    704   USE netcdf, ONLY: nf90_noerr, NF90_strerror
    705 !-------------------------------------------------------------------------------
    706 ! Arguments:
    707   CHARACTER(LEN=*), INTENT(IN) :: att_name
    708 !-------------------------------------------------------------------------------
    709   IF(ncerr /= nf90_noerr) &
    710     WRITE(lunout,*)TRIM(sub)//" prepare_out NF95_COPY_ATT "//TRIM(att_name)//  &
    711                       " -- "//TRIM(NF90_strerror(ncerr))
    712 
    713 END SUBROUTINE handle_err_copy_att
    714 
    715 !-------------------------------------------------------------------------------
    716 
    717 END SUBROUTINE prepare_out
    718 
    719 !-------------------------------------------------------------------------------
    720 
    721 END MODULE regr_horiz_time_climoz_m
    722 
    723 !-------------------------------------------------------------------------------

  • LMDZ6/branches/Amaury_dev/libf/phylmd/regr_pr_time_av_m.F90

    r5099 r5100  
    113113!-------------------------------------------------------------------------------
    114114  USE dimphy,         ONLY: klon
    115   USE netcdf95,       ONLY: NF95_INQ_VARID, NF95_INQUIRE_VARIABLE, &
    116                             NF95_INQ_DIMID, NF95_INQUIRE_DIMENSION, nf95_get_var
     115  USE netcdf95,       ONLY: nf95_inq_varid, NF95_INQUIRE_VARIABLE, &
     116                            nf95_inq_dimid, nf95_inquire_dimension, nf95_get_var
    117117  USE netcdf,         ONLY: NF90_INQ_VARID, nf90_noerr
    118118  USE assert_m,       ONLY: assert
    … …  
    213213      lPrTfile=lAdjTro.AND.NF90_INQ_VARID(fID,"tropopause_air_pressure",vID)==nf90_noerr
    214214      lO3Tfile=lAdjTro.AND.NF90_INQ_VARID(fID,"tro3_at_tropopause"     ,vID)==nf90_noerr
    215       CALL NF95_INQ_DIMID(fID,"time",vID)
    216       CALL NF95_INQUIRE_DIMENSION(fID,vID,nclen=ntim_in)
     215      CALL nf95_inq_dimid(fID,"time",vID)
     216      CALL nf95_inquire_dimension(fID,vID,nclen=ntim_in)
    217217      linterp=PRESENT(time_in).AND.ntim_in==14
    218218      ALLOCATE(v1(nlon,nlat,nlev_in,n_var))
    … …  
    480480  CHARACTER(LEN=*), INTENT(IN)    :: var
    481481!-------------------------------------------------------------------------------
    482   CALL NF95_INQ_VARID(fID, TRIM(var), vID)
     482  CALL nf95_inq_varid(fID, TRIM(var), vID)
    483483  CALL NF95_INQUIRE_VARIABLE(fID, vID, ndims=n_dim)
    484484  IF(n_dim==2) call NF95_GET_VAR(fID,vID,v(1,:), start=[  1,irec])
    … …  
    511511!-------------------------------------------------------------------------------
    512512  DO i=1,SIZE(nam)
    513     CALL NF95_INQ_VARID(fID, TRIM(nam(i)), vID)
     513    CALL nf95_inq_varid(fID, TRIM(nam(i)), vID)
    514514    CALL NF95_INQUIRE_VARIABLE(fID, vID, ndims=n_dim)
    515515    IF(n_dim==3) call NF95_GET_VAR(fID,vID,v(1,:,:,i), start=[  1,1,irec])
Note: See TracChangeset for help on using the changeset viewer.