Ignore:
Timestamp:
Jul 19, 2024, 6:40:44 PM (14 months ago)
Author:
Laurent Fairhead
Message:

Reverting to r4065. Updating fortran standard broke too much stuff. Will do it by smaller chunks
AB, LF

Location:
LMDZ6/trunk/libf/phylmd/Dust
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/Dust/condsurfc.F

    r5075 r5084  
    44     .                     lmt_omnat)
    55      USE dimphy
    6       USE lmdz_netcdf, ONLY: nf_close,nf_noerr,nf_inq_varid,nf_open,nf_nowrite,nf90_get_var
    76      IMPLICIT none
    87!
    … …  
    1110!
    1211      INCLUDE "dimensions.h"
    13 
     12      INCLUDE "netcdf.inc"
     13     
    1414      REAL lmt_bcff(klon), lmt_bcbb(klon),lmt_bc_penner(klon)
    1515      REAL lmt_omff(klon), lmt_ombb(klon)
    … …  
    2424      INTEGER debut(2),epais(2)
    2525!
    26       IF (jour<0 .OR. jour>(360-1)) THEN
    27          IF (jour>(360-1).AND.jour<=367) THEN
     26      IF (jour.LT.0 .OR. jour.GT.(360-1)) THEN
     27         IF (jour.GT.(360-1).AND.jour.LE.367) THEN
    2828           jour=360-1
    2929           print *,'JE: jour changed to jour= ',jour
    … …  
    3535!
    3636      ierr = NF_OPEN ("limitcarbon.nc", NF_NOWRITE, nid1)
    37       if (ierr/=NF_NOERR) then
     37      if (ierr.ne.NF_NOERR) then
    3838        write(6,*)' Pb d''ouverture du fichier limitbc.nc'
    3939        write(6,*)' ierr = ', ierr
    … …  
    4949!
    5050      ierr = NF_INQ_VARID (nid1, "BCFF", nvarid)
    51       ierr = nf90_get_var(nid1, nvarid, lmt_bcff, debut, epais)
     51!nhl #ifdef NC_DOUBLE
     52      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, lmt_bcff)
    5253!      print *,'IERR = ',ierr
    5354!      print *,'NF_NOERR = ',NF_NOERR
    5455!      print *,'debut = ',debut
    5556!      print *,'epais = ',epais
    56       IF (ierr /= NF_NOERR) THEN
     57!nhl #else
     58!nhl       ierr = NF_GET_VARA_REAL (nid1, nvarid, debut, epais, lmt_bcff)
     59!nhl #endif
     60      IF (ierr .NE. NF_NOERR) THEN
    5761         PRINT*, 'Pb de lecture pour les sources BC'
    5862         CALL exit(1)
    … …  
    6165!
    6266      ierr = NF_INQ_VARID (nid1, "BCBB", nvarid)
    63       ierr = nf90_get_var(nid1, nvarid, lmt_bcbb, debut, epais)
    64       IF (ierr /= NF_NOERR) THEN
     67!nhl #ifdef NC_DOUBLE
     68      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, lmt_bcbb)
     69!nhl #else
     70!nhl       ierr = NF_GET_VARA_REAL (nid1, nvarid, debut, epais, lmt_bcbb)
     71!nhl #endif
     72      IF (ierr .NE. NF_NOERR) THEN
    6573         PRINT*, 'Pb de lecture pour les sources BC-biomass'
    6674         CALL exit(1)
    … …  
    6977!
    7078      ierr = NF_INQ_VARID (nid1, "BCBL", nvarid)
    71       ierr = nf90_get_var(nid1, nvarid, lmt_bcbbl, debut, epais)
    72       IF (ierr /= NF_NOERR) THEN
     79!nhl #ifdef NC_DOUBLE
     80      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, lmt_bcbbl)
     81!nhl #else
     82!nhl       ierr = NF_GET_VARA_REAL (nid1, nvarid, debut, epais, lmt_bcbbl)
     83!nhl #endif
     84      IF (ierr .NE. NF_NOERR) THEN
    7385         PRINT*, 'Pb de lecture pour les sources BC low'
    7486         CALL exit(1)
    … …  
    7789!
    7890      ierr = NF_INQ_VARID (nid1, "BCBH", nvarid)
    79       ierr = nf90_get_var (nid1, nvarid, lmt_bcbbh, debut, epais)
    80       IF (ierr /= NF_NOERR) THEN
     91!nhl #ifdef NC_DOUBLE
     92      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, lmt_bcbbh)
     93!nhl #else
     94!nhl       ierr = NF_GET_VARA_REAL (nid1, nvarid, debut, epais, lmt_bcbbh)
     95!nhl #endif
     96      IF (ierr .NE. NF_NOERR) THEN
    8197         PRINT*, 'Pb de lecture pour les sources BC high'
    8298         CALL exit(1)
    … …  
    84100!
    85101      ierr = NF_INQ_VARID (nid1, "TERP", nvarid)
    86       ierr = nf90_get_var (nid1, nvarid, lmt_terp, debut, epais)
    87       IF (ierr /= NF_NOERR) THEN
     102!nhl #ifdef NC_DOUBLE
     103      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, lmt_terp)
     104!nhl #else
     105!nhl       ierr = NF_GET_VARA_REAL (nid1, nvarid, debut, epais, lmt_terp)
     106!nhl #endif
     107      IF (ierr .NE. NF_NOERR) THEN
    88108         PRINT*, 'Pb de lecture pour les sources Terpene'
    89109         CALL exit(1)
    … …  
    92112!
    93113      ierr = NF_INQ_VARID (nid1, "BC_penner", nvarid)
    94       ierr = nf90_get_var (nid1, nvarid, lmt_bc_penner, debut,  epais)
    95       IF (ierr /= NF_NOERR) THEN
     114!nhl #ifdef NC_DOUBLE
     115      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut,
     116     .       epais, lmt_bc_penner)
     117!nhl #else
     118!nhl       ierr = NF_GET_VARA_REAL (nid1, nvarid, debut, epais,
     119!nhl      .       lmt_bc_penner)
     120!nhl #endif
     121      IF (ierr .NE. NF_NOERR) THEN
    96122         PRINT*, 'Pb de lecture pour les sources BC Penner'
    97123         CALL exit(1)
    … …  
    100126!
    101127      ierr = NF_INQ_VARID (nid1, "OMFF", nvarid)
    102       ierr = nf90_get_var (nid1, nvarid, lmt_omff, debut, epais)
    103       IF (ierr /= NF_NOERR) THEN
     128!nhl #ifdef NC_DOUBLE
     129      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, lmt_omff)
     130!nhl #else
     131!nhl       ierr = NF_GET_VARA_REAL (nid1, nvarid, debut, epais, lmt_omff)
     132!nhl #endif
     133      IF (ierr .NE. NF_NOERR) THEN
    104134         PRINT*, 'Pb de lecture pour les sources om-ifossil'
    105135         CALL exit(1)

  • LMDZ6/trunk/libf/phylmd/Dust/condsurfc_new.F

    r5075 r5084  
    66      USE mod_phys_lmdz_para
    77      USE dimphy
    8       USE lmdz_netcdf, ONLY:nf90_get_var,nf_close,nf_noerr,nf_inq_varid,nf_open,nf_nowrite
    98      IMPLICIT none
    109c
    … …  
    1312c
    1413      INCLUDE "dimensions.h"
    15 
     14      INCLUDE "netcdf.inc"
     15     
    1616      REAL lmt_bcff(klon), lmt_bcnff(klon), lmt_bcba(klon)
    1717      REAL lmt_omff(klon), lmt_omnff(klon), lmt_ombb(klon)
    … …  
    3636c
    3737!      IF (jour.LT.0 .OR. jour.GT.(366-1)) THEN
    38       IF (jour<0 .OR. jour>366) THEN
     38      IF (jour.LT.0 .OR. jour.GT.366) THEN
    3939         PRINT*,'Le jour demande n est pas correcte:', jour
    4040         print *,'JE: FORCED TO CONTINUE (emissions have
    … …  
    5858!
    5959      ierr = NF_OPEN ("carbon_emissions.nc", NF_NOWRITE, nid1)
    60       if (ierr/=NF_NOERR) then
     60      if (ierr.ne.NF_NOERR) then
    6161        write(6,*)' Pb d''ouverture du fichier limitbc.nc'
    6262        write(6,*)' ierr = ', ierr
    … …  
    6767!
    6868      ierr = NF_INQ_VARID (nid1, "BCFF", nvarid)
    69       ierr = nf90_get_var (nid1, nvarid, lmt_bcff_glo, debut, epais)
    70       IF (ierr /= NF_NOERR) THEN
     69      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais,
     70     .  lmt_bcff_glo)
     71      IF (ierr .NE. NF_NOERR) THEN
    7172         PRINT*, 'Pb de lecture pour les sources BC'
    7273         CALL exit(1)
    … …  
    7879!
    7980      ierr = NF_INQ_VARID (nid1, "BCNFF", nvarid)
    80       ierr = nf90_get_var (nid1, nvarid, lmt_bcnff_glo, debut, epais)
    81       IF (ierr /= NF_NOERR) THEN
     81      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais,
     82     .    lmt_bcnff_glo)
     83      IF (ierr .NE. NF_NOERR) THEN
    8284         PRINT*, 'Pb de lecture pour les sources BC'
    8385         CALL exit(1)
    … …  
    8789!
    8890      ierr = NF_INQ_VARID (nid1, "BCBBL", nvarid)
    89       ierr = nf90_get_var (nid1, nvarid, lmt_bcbbl_glo, debut, epais)
    90       IF (ierr /= NF_NOERR) THEN
     91      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais,
     92     .  lmt_bcbbl_glo)
     93      IF (ierr .NE. NF_NOERR) THEN
    9194         PRINT*, 'Pb de lecture pour les sources BC low'
    9295         CALL exit(1)
    … …  
    9699!
    97100      ierr = NF_INQ_VARID (nid1, "BCBBH", nvarid)
    98       ierr = nf90_get_var (nid1, nvarid, lmt_bcbbh_glo, debut, epais)
    99       IF (ierr /= NF_NOERR) THEN
     101      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais,
     102     .      lmt_bcbbh_glo)
     103      IF (ierr .NE. NF_NOERR) THEN
    100104         PRINT*, 'Pb de lecture pour les sources BC high'
    101105         CALL exit(1)
    … …  
    105109!
    106110      ierr = NF_INQ_VARID (nid1, "BCBA", nvarid)
    107       ierr = nf90_get_var (nid1, nvarid, lmt_bcba_glo, debut, epais)
    108       IF (ierr /= NF_NOERR) THEN
     111      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais,
     112     .   lmt_bcba_glo)
     113      IF (ierr .NE. NF_NOERR) THEN
    109114         PRINT*, 'Pb de lecture pour les sources BC'
    110115         CALL exit(1)
    … …  
    120125!
    121126      ierr = NF_INQ_VARID (nid1, "OMFF", nvarid)
    122       ierr = nf90_get_var (nid1, nvarid, lmt_omff_glo, debut, epais)
    123       IF (ierr /= NF_NOERR) THEN
     127      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais,
     128     .  lmt_omff_glo)
     129      IF (ierr .NE. NF_NOERR) THEN
    124130         PRINT*, 'Pb de lecture pour les sources OM'
    125131         CALL exit(1)
    … …  
    129135!
    130136      ierr = NF_INQ_VARID (nid1, "OMNFF", nvarid)
    131       ierr = nf90_get_var (nid1, nvarid, lmt_omnff_glo, debut, epais)
    132       IF (ierr /= NF_NOERR) THEN
     137      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais,
     138     .   lmt_omnff_glo)
     139      IF (ierr .NE. NF_NOERR) THEN
    133140         PRINT*, 'Pb de lecture pour les sources OM'
    134141         CALL exit(1)
    … …  
    138145!
    139146      ierr = NF_INQ_VARID (nid1, "OMBBL", nvarid)
    140       ierr = nf90_get_var (nid1, nvarid, lmt_ombbl_glo, debut, epais)
    141       IF (ierr /= NF_NOERR) THEN
     147      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais,
     148     .  lmt_ombbl_glo)
     149      IF (ierr .NE. NF_NOERR) THEN
    142150         PRINT*, 'Pb de lecture pour les sources OM low'
    143151         CALL exit(1)
    … …  
    147155!
    148156      ierr = NF_INQ_VARID (nid1, "OMBBH", nvarid)
    149       ierr = nf90_get_var (nid1, nvarid, lmt_ombbh_glo, debut, epais)
    150       IF (ierr /= NF_NOERR) THEN
     157      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais,
     158     .  lmt_ombbh_glo)
     159      IF (ierr .NE. NF_NOERR) THEN
    151160         PRINT*, 'Pb de lecture pour les sources OM high'
    152161         CALL exit(1)
    … …  
    156165!
    157166      ierr = NF_INQ_VARID (nid1, "OMBA", nvarid)
    158       ierr = nf90_get_var (nid1, nvarid, lmt_omba_glo, debut, epais)
    159       IF (ierr /= NF_NOERR) THEN
     167      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais,
     168     .   lmt_omba_glo)
     169      IF (ierr .NE. NF_NOERR) THEN
    160170         PRINT*, 'Pb de lecture pour les sources OM ship'
    161171         CALL exit(1)
    … …  
    165175!
    166176      ierr = NF_INQ_VARID (nid1, "TERP", nvarid)
    167       ierr = nf90_get_var (nid1, nvarid, lmt_terp_glo, debut, epais)
    168       IF (ierr /= NF_NOERR) THEN
     177      ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais,
     178     .  lmt_terp_glo)
     179      IF (ierr .NE. NF_NOERR) THEN
    169180         PRINT*, 'Pb de lecture pour les sources Terpene'
    170181         CALL exit(1)

  • LMDZ6/trunk/libf/phylmd/Dust/condsurfs.F

    r5075 r5084  
    44     .                     lmt_dmsbio, lmt_h2sbio, lmt_dms, lmt_dmsconc)
    55       USE dimphy
    6        USE lmdz_netcdf, ONLY:nf_close,nf_noerr,nf_inq_varid,nf_open,nf_nowrite,nf90_get_var
    76      IMPLICIT none
    87c
    … …  
    1110c
    1211      INCLUDE "dimensions.h"
     12      INCLUDE "netcdf.inc"
    1313c
    1414      REAL lmt_so2h(klon), lmt_so2b(klon), lmt_so2bb(klon)
    … …  
    2424      INTEGER debut(2),epais(2)
    2525c
    26       IF (jour<0 .OR. jour>(360-1)) THEN
    27          IF ((jour>(360-1)) .AND. (jour<=367)) THEN
     26      IF (jour.LT.0 .OR. jour.GT.(360-1)) THEN
     27         IF ((jour.GT.(360-1)) .AND. (jour.LE.367)) THEN
    2828           jour=360-1
    2929           print *,'JE: jour changed to jour= ',jour
    … …  
    3535c
    3636      ierr = NF_OPEN ("limitsoufre.nc", NF_NOWRITE, nid)
    37       if (ierr/=NF_NOERR) then
     37      if (ierr.ne.NF_NOERR) then
    3838        write(6,*)' Pb d''ouverture du fichier limitsoufre.nc'
    3939        write(6,*)' ierr = ', ierr
    … …  
    4848c
    4949      ierr = NF_INQ_VARID (nid, "VOLC", nvarid)
    50       ierr = nf90_get_var(nid, nvarid, lmt_so2volc, debut, epais)
    51       IF (ierr /= NF_NOERR) THEN
     50cnhl #ifdef NC_DOUBLE
     51      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2volc)
     52cnhl #else
     53cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2volc)
     54cnhl #endif
     55      IF (ierr .NE. NF_NOERR) THEN
    5256         PRINT*, 'Pb de lecture pour les sources so2 volcan'
    5357         CALL exit(1)
    … …  
    5559c
    5660      ierr = NF_INQ_VARID (nid, "ALTI", nvarid)
    57       ierr = nf90_get_var(nid, nvarid, lmt_altvolc, debut, epais)
    58       IF (ierr /= NF_NOERR) THEN
     61cnhl #ifdef NC_DOUBLE
     62      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_altvolc)
     63cnhl #else
     64cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_altvolc)
     65cnhl #endif
     66      IF (ierr .NE. NF_NOERR) THEN
    5967         PRINT*, 'Pb de lecture pour les altitudes volcan'
    6068         CALL exit(1)
    … …  
    6472c
    6573      ierr = NF_INQ_VARID (nid, "SO2ED95L", nvarid)
    66       ierr = nf90_get_var(nid, nvarid, lmt_so2b, debut, epais)
    67       IF (ierr /= NF_NOERR) THEN
     74cnhl #ifdef NC_DOUBLE
     75      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2b)
     76cnhl #else
     77cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2b)
     78cnhl #endif
     79      IF (ierr .NE. NF_NOERR) THEN
    6880         PRINT*, 'Pb de lecture pour les sources so2 edgar low'
    6981         CALL exit(1)
    … …  
    7183c
    7284      ierr = NF_INQ_VARID (nid, "SO2ED95H", nvarid)
    73       ierr = nf90_get_var(nid, nvarid, lmt_so2h, debut, epais)
    74       IF (ierr /= NF_NOERR) THEN
     85cnhl #ifdef NC_DOUBLE
     86      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2h)
     87cnhl #else
     88cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2h)
     89cnhl #endif
     90      IF (ierr .NE. NF_NOERR) THEN
    7591         PRINT*, 'Pb de lecture pour les sources so2 edgar high'
    7692         CALL exit(1)
    … …  
    8096c
    8197      ierr = NF_INQ_VARID (nid, "SO2H", nvarid)
    82       ierr = nf90_get_var(nid, nvarid, lmt_so2h, debut, epais)
    83       IF (ierr /= NF_NOERR) THEN
     98cnhl #ifdef NC_DOUBLE
     99      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2h)
     100cnhl #else
     101cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2h)
     102cnhl #endif
     103      IF (ierr .NE. NF_NOERR) THEN
    84104         PRINT*, 'Pb de lecture pour les sources so2 haut'
    85105         CALL exit(1)
    … …  
    87107c
    88108      ierr = NF_INQ_VARID (nid, "SO2B", nvarid)
    89       ierr = nf90_get_var(nid, nvarid, lmt_so2b, debut, epais)
    90       IF (ierr /= NF_NOERR) THEN
     109cnhl #ifdef NC_DOUBLE
     110      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2b)
     111cnhl #else
     112cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2b)
     113cnhl #endif
     114      IF (ierr .NE. NF_NOERR) THEN
    91115         PRINT*, 'Pb de lecture pour les sources so2 bas'
    92116         CALL exit(1)
    … …  
    96120c
    97121      ierr = NF_INQ_VARID (nid, "SO2BB", nvarid)
    98       ierr = nf90_get_var(nid, nvarid, lmt_so2bb, debut, epais)
    99       IF (ierr /= NF_NOERR) THEN
     122cnhl #ifdef NC_DOUBLE
     123      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2bb)
     124cnhl #else
     125cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2bb)
     126cnhl #endif
     127      IF (ierr .NE. NF_NOERR) THEN
    100128         PRINT*, 'Pb de lecture pour les sources so2 bb'
    101129         CALL exit(1)
    … …  
    103131c
    104132      ierr = NF_INQ_VARID (nid, "SO2BA", nvarid)
    105       ierr = nf90_get_var(nid, nvarid, lmt_so2ba, debut, epais)
    106       IF (ierr /= NF_NOERR) THEN
     133cnhl #ifdef NC_DOUBLE
     134      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2ba)
     135cnhl #else
     136cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2ba)
     137cnhl #endif
     138      IF (ierr .NE. NF_NOERR) THEN
    107139         PRINT*, 'Pb de lecture pour les sources so2 bateau'
    108140         CALL exit(1)
    … …  
    110142c
    111143      ierr = NF_INQ_VARID (nid, "DMSB", nvarid)
    112       ierr = nf90_get_var(nid, nvarid, lmt_dmsbio, debut, epais)
    113       IF (ierr /= NF_NOERR) THEN
     144cnhl #ifdef NC_DOUBLE
     145      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_dmsbio)
     146cnhl #else
     147cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsbio)
     148cnhl #endif
     149      IF (ierr .NE. NF_NOERR) THEN
    114150         PRINT*, 'Pb de lecture pour les sources dms bio'
    115151         CALL exit(1)
    … …  
    117153c
    118154      ierr = NF_INQ_VARID (nid, "H2SB", nvarid)
    119       ierr = nf90_get_var(nid, nvarid, lmt_h2sbio, debut, epais)
    120       IF (ierr /= NF_NOERR) THEN
     155cnhl #ifdef NC_DOUBLE
     156      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_h2sbio)
     157cnhl #else
     158cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_h2sbio)
     159cnhl #endif
     160      IF (ierr .NE. NF_NOERR) THEN
    121161         PRINT*, 'Pb de lecture pour les sources h2s bio'
    122162         CALL exit(1)
    123163      ENDIF
    124164c
    125       IF (flag_dms==1) THEN
     165      IF (flag_dms.EQ.1) THEN
    126166c
    127167      ierr = NF_INQ_VARID (nid, "DMSL", nvarid)
    128       ierr = nf90_get_var(nid, nvarid, lmt_dms, debut, epais)
    129       IF (ierr /= NF_NOERR) THEN
     168cnhl #ifdef NC_DOUBLE
     169      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_dms)
     170cnhl #else
     171cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dms)
     172cnhl #endif
     173      IF (ierr .NE. NF_NOERR) THEN
    130174         PRINT*, 'Pb de lecture pour les sources dms liss'
    131175         CALL exit(1)
    132176      ENDIF
    133177c
    134       ELSEIF (flag_dms==2) THEN
     178      ELSEIF (flag_dms.EQ.2) THEN
    135179c
    136180      ierr = NF_INQ_VARID (nid, "DMSW", nvarid)
    137       ierr = nf90_get_var(nid, nvarid, lmt_dms, debut, epais)
    138       IF (ierr /= NF_NOERR) THEN
     181cnhl #ifdef NC_DOUBLE
     182      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_dms)
     183cnhl #else
     184cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dms)
     185cnhl #endif
     186      IF (ierr .NE. NF_NOERR) THEN
    139187         PRINT*, 'Pb de lecture pour les sources dms wann'
    140188         CALL exit(1)
    141189      ENDIF
    142190c
    143       ELSEIF (flag_dms==3) THEN
     191      ELSEIF (flag_dms.EQ.3) THEN
    144192c
    145193      ierr = NF_INQ_VARID (nid, "DMSC1", nvarid)
    146       ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
    147       IF (ierr /= NF_NOERR) THEN
     194cnhl #ifdef NC_DOUBLE
     195      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
     196cnhl #else
     197cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
     198cnhl #endif
     199      IF (ierr .NE. NF_NOERR) THEN
    148200         PRINT*, 'Pb de lecture pour les sources dmsconc old'
    149201         CALL exit(1)
    150202      ENDIF
    151203c
    152       ELSEIF (flag_dms==4) THEN
     204      ELSEIF (flag_dms.EQ.4) THEN
    153205c
    154206      ierr = NF_INQ_VARID (nid, "DMSC2", nvarid)
    155       ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
    156       IF (ierr /= NF_NOERR) THEN
     207cnhl #ifdef NC_DOUBLE
     208      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
     209cnhl #else
     210cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
     211cnhl #endif
     212      IF (ierr .NE. NF_NOERR) THEN
    157213         PRINT*, 'Pb de lecture pour les sources dms conc 2'
    158214         CALL exit(1)
    159215      ENDIF
    160216c
    161       ELSEIF (flag_dms==5) THEN
     217      ELSEIF (flag_dms.EQ.5) THEN
    162218c
    163219      ierr = NF_INQ_VARID (nid, "DMSC3", nvarid)
    164       ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
    165       IF (ierr /= NF_NOERR) THEN
     220cnhl #ifdef NC_DOUBLE
     221      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
     222cnhl #else
     223cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
     224cnhl #endif
     225      IF (ierr .NE. NF_NOERR) THEN
    166226         PRINT*, 'Pb de lecture pour les sources dms conc 3'
    167227         CALL exit(1)
    168228      ENDIF
    169229c
    170       ELSEIF (flag_dms==6) THEN
     230      ELSEIF (flag_dms.EQ.6) THEN
    171231c
    172232      ierr = NF_INQ_VARID (nid, "DMSC4", nvarid)
    173       ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
    174       IF (ierr /= NF_NOERR) THEN
     233cnhl #ifdef NC_DOUBLE
     234      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
     235cnhl #else
     236cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
     237cnhl #endif
     238      IF (ierr .NE. NF_NOERR) THEN
    175239         PRINT*, 'Pb de lecture pour les sources dms conc 4'
    176240         CALL exit(1)
    177241      ENDIF
    178242c
    179       ELSEIF (flag_dms==7) THEN
     243      ELSEIF (flag_dms.EQ.7) THEN
    180244c
    181245      ierr = NF_INQ_VARID (nid, "DMSC5", nvarid)
    182       ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
    183       IF (ierr /= NF_NOERR) THEN
     246cnhl #ifdef NC_DOUBLE
     247      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
     248cnhl #else
     249cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
     250cnhl #endif
     251      IF (ierr .NE. NF_NOERR) THEN
    184252         PRINT*, 'Pb de lecture pour les sources dms conc 5'
    185253         CALL exit(1)
    186254      ENDIF
    187255c
    188       ELSEIF (flag_dms==8) THEN
     256      ELSEIF (flag_dms.EQ.8) THEN
    189257c
    190258      ierr = NF_INQ_VARID (nid, "DMSC6", nvarid)
    191       ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
    192       IF (ierr /= NF_NOERR) THEN
     259cnhl #ifdef NC_DOUBLE
     260      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
     261cnhl #else
     262cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
     263cnhl #endif
     264      IF (ierr .NE. NF_NOERR) THEN
    193265         PRINT*, 'Pb de lecture pour les sources dms conc 6'
    194266         CALL exit(1)
    195267      ENDIF
    196268c
    197       ELSEIF (flag_dms==9) THEN
     269      ELSEIF (flag_dms.EQ.9) THEN
    198270c
    199271      ierr = NF_INQ_VARID (nid, "DMSC7", nvarid)
    200       ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
    201       IF (ierr /= NF_NOERR) THEN
     272cnhl #ifdef NC_DOUBLE
     273      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
     274cnhl #else
     275cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
     276cnhl #endif
     277      IF (ierr .NE. NF_NOERR) THEN
    202278         PRINT*, 'Pb de lecture pour les sources dms conc 7'
    203279         CALL exit(1)
    204280      ENDIF
    205281c
    206       ELSEIF (flag_dms==10) THEN
     282      ELSEIF (flag_dms.EQ.10) THEN
    207283c
    208284      ierr = NF_INQ_VARID (nid, "DMSC8", nvarid)
    209       ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
    210       IF (ierr /= NF_NOERR) THEN
     285cnhl #ifdef NC_DOUBLE
     286      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
     287cnhl #else
     288cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
     289cnhl #endif
     290      IF (ierr .NE. NF_NOERR) THEN
    211291         PRINT*, 'Pb de lecture pour les sources dms conc 8'
    212292         CALL exit(1)
    … …  
    222302      ierr = NF_CLOSE(nid)
    223303c
    224       IF (flag_dms<=2) THEN
     304      IF (flag_dms.LE.2) THEN
    225305      DO i=1, klon
    226306         lmt_dmsconc(i)=0.0

  • LMDZ6/trunk/libf/phylmd/Dust/condsurfs_new.F

    r5075 r5084  
    99      USE mod_phys_lmdz_para
    1010      USE dimphy
    11       USE lmdz_netcdf, ONLY: nf90_get_var,nf_inq_varid,nf_close,nf_noerr,nf_open,nf_nowrite
    1211      IMPLICIT none
    1312c
    … …  
    1615c
    1716      INCLUDE "dimensions.h"
     17      INCLUDE "netcdf.inc"
    1818c
    1919      REAL lmt_so2b(klon), lmt_so2h(klon), lmt_so2nff(klon)
    … …  
    4040      INTEGER debut(2),epais(2)
    4141c
    42       IF (jour<0 .OR. jour>(366-1)) THEN
     42      IF (jour.LT.0 .OR. jour.GT.(366-1)) THEN
    4343         PRINT*,'Le jour demande n est pas correcte:', jour
    4444         print *,'JE: FORCED TO CONTINUE (emissions have
    … …  
    6262!
    6363      ierr = NF_OPEN ("sulphur_emissions_antro.nc", NF_NOWRITE, nid)
    64       if (ierr/=NF_NOERR) then
     64      if (ierr.ne.NF_NOERR) then
    6565        write(6,*)' Pb d''ouverture du fichier sulphur_emissions_antro'
    6666        write(6,*)' ierr = ', ierr
    … …  
    7272!
    7373      ierr = NF_INQ_VARID (nid, "SO2FF_LOW", nvarid)
    74       ierr = nf90_get_var(nid, nvarid, lmt_so2b_glo, debut, epais)
    75       IF (ierr /= NF_NOERR) THEN
     74      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais,lmt_so2b_glo)
     75      IF (ierr .NE. NF_NOERR) THEN
    7676        PRINT*, 'Pb de lecture pour les sources so2 low'
    7777        print *,'JE klon, jour, debut ,epais ',klon_glo,jour,debut,epais
    … …  
    8484!
    8585      ierr = NF_INQ_VARID (nid, "SO2FF_HIGH", nvarid)
    86       ierr = nf90_get_var(nid, nvarid, lmt_so2h_glo, debut, epais)
    87       IF (ierr /= NF_NOERR) THEN
     86      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais,lmt_so2h_glo)
     87      IF (ierr .NE. NF_NOERR) THEN
    8888        PRINT*, 'Pb de lecture pour les sources so2 high'
    8989        CALL exit(1)
    … …  
    9393!
    9494      ierr = NF_INQ_VARID (nid, "SO2BBH", nvarid)
    95       ierr = nf90_get_var(nid, nvarid,  lmt_so2bb_h_glo, debut, epais)
    96       IF (ierr /= NF_NOERR) THEN
     95      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,
     96     . epais, lmt_so2bb_h_glo)
     97      IF (ierr .NE. NF_NOERR) THEN
    9798        PRINT*, 'Pb de lecture pour les sources so2 BB high'
    9899        CALL exit(1)
    … …  
    102103!
    103104      ierr = NF_INQ_VARID (nid, "SO2BBL", nvarid)
    104       ierr = nf90_get_var(nid, nvarid,  lmt_so2bb_l_glo, debut, epais)
    105       IF (ierr /= NF_NOERR) THEN
     105      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,
     106     . epais, lmt_so2bb_l_glo)
     107      IF (ierr .NE. NF_NOERR) THEN
    106108        PRINT*, 'Pb de lecture pour les sources so2 BB low'
    107109        CALL exit(1)
    … …  
    111113!
    112114      ierr = NF_INQ_VARID (nid, "SO2BA", nvarid)
    113       ierr = nf90_get_var(nid, nvarid, lmt_so2ba_glo, debut, epais)
    114       IF (ierr /= NF_NOERR) THEN
     115      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_so2ba_glo)
     116      IF (ierr .NE. NF_NOERR) THEN
    115117        PRINT*, 'Pb de lecture pour les sources so2 ship'
    116118        CALL exit(1)
    … …  
    120122!
    121123      ierr = NF_INQ_VARID (nid, "SO2NFF", nvarid)
    122       ierr = nf90_get_var(nid, nvarid, lmt_so2nff_glo, debut, epais)
    123       IF (ierr /= NF_NOERR) THEN
     124      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais,
     125     .  lmt_so2nff_glo)
     126      IF (ierr .NE. NF_NOERR) THEN
    124127        PRINT*, 'Pb de lecture pour les sources so2 non FF'
    125128        CALL exit(1)
    … …  
    132135!=======================================================================
    133136      ierr = NF_OPEN ("sulphur_emissions_nat.nc", NF_NOWRITE, nid)
    134       if (ierr/=NF_NOERR) then
     137      if (ierr.ne.NF_NOERR) then
    135138        write(6,*)' Pb d''ouverture du fichier sulphur_emissions_nat'
    136139        write(6,*)' ierr = ', ierr
    … …  
    141144c
    142145      ierr = NF_INQ_VARID (nid, "DMSB", nvarid)
    143       ierr = nf90_get_var(nid, nvarid, lmt_dmsbio_glo, debut, epais)
    144       IF (ierr /= NF_NOERR) THEN
     146      ierr = NF_GET_VARA_DOUBLE (nid, nvarid,debut,epais,lmt_dmsbio_glo)
     147      IF (ierr .NE. NF_NOERR) THEN
    145148         PRINT*, 'Pb de lecture pour les sources dms bio'
    146149         CALL exit(1)
    … …  
    150153c
    151154      ierr = NF_INQ_VARID (nid, "H2SB", nvarid)
    152       ierr = nf90_get_var(nid, nvarid, lmt_h2sbio_glo, debut, epais)
    153       IF (ierr /= NF_NOERR) THEN
     155      ierr = NF_GET_VARA_DOUBLE (nid, nvarid,debut,epais,lmt_h2sbio_glo)
     156      IF (ierr .NE. NF_NOERR) THEN
    154157         PRINT*, 'Pb de lecture pour les sources h2s bio'
    155158         CALL exit(1)
    … …  
    158161c Ocean surface concentration of dms (emissions are computed later)
    159162c
    160       IF (flag_dms==4) THEN
     163      IF (flag_dms.EQ.4) THEN
    161164c
    162165      ierr = NF_INQ_VARID (nid, "DMSC2", nvarid)
    163       ierr = nf90_get_var(nid, nvarid, lmt_dmsconc_glo, debut, epais)
    164       IF (ierr /= NF_NOERR) THEN
     166      ierr = NF_GET_VARA_DOUBLE (nid,nvarid,debut,epais,lmt_dmsconc_glo)
     167      IF (ierr .NE. NF_NOERR) THEN
    165168         PRINT*, 'Pb de lecture pour les sources dms conc 2'
    166169         CALL exit(1)
    … …  
    187190      print *,' Jour = ',jour
    188191      ierr = NF_OPEN ("sulphur_emissions_volc.nc", NF_NOWRITE, nid)
    189       if (ierr/=NF_NOERR) then
     192      if (ierr.ne.NF_NOERR) then
    190193        write(6,*)' Pb d''ouverture du fichier sulphur_emissions_volc'
    191194        write(6,*)' ierr = ', ierr
    … …  
    197200!      ierr = NF_INQ_VARID (nid, "VOLC", nvarid)
    198201      ierr = NF_INQ_VARID (nid, "flx_volc_cont", nvarid)
    199       ierr = nf90_get_var(nid, nvarid, lmt_so2volc_cont_glo, debut, epais)
    200       IF (ierr /= NF_NOERR) THEN
     202      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais,
     203     .                           lmt_so2volc_cont_glo)
     204      IF (ierr .NE. NF_NOERR) THEN
    201205         PRINT*, 'Pb de lecture pour les sources so2 volcan (cont)'
    202206         CALL exit(1)
    … …  
    210214!      ierr = NF_INQ_VARID (nid, "ALTI", nvarid)
    211215      ierr = NF_INQ_VARID (nid, "flx_volc_altcont", nvarid)
    212       ierr = nf90_get_var(nid, nvarid, lmt_altvolc_cont_glo, debut, epais)
    213       IF (ierr /= NF_NOERR) THEN
     216      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais,
     217     .                           lmt_altvolc_cont_glo)
     218      IF (ierr .NE. NF_NOERR) THEN
    214219         PRINT*, 'Pb de lecture pour les altitudes volcan (cont)'
    215220         CALL exit(1)
    … …  
    219224c
    220225      ierr = NF_INQ_VARID (nid, "flx_volc_expl", nvarid)
    221       ierr = nf90_get_var(nid, nvarid, lmt_so2volc_expl_glo, debut, epais)
    222       IF (ierr /= NF_NOERR) THEN
     226      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais,
     227     .                           lmt_so2volc_expl_glo)
     228      IF (ierr .NE. NF_NOERR) THEN
    223229         PRINT*, 'Pb de lecture pour les sources so2 volcan (expl)'
    224230         CALL exit(1)
    … …  
    231237c
    232238      ierr = NF_INQ_VARID (nid, "flx_volc_altexpl", nvarid)
    233       ierr = nf90_get_var(nid, nvarid, lmt_altvolc_expl_glo, debut, epais)
    234       IF (ierr /= NF_NOERR) THEN
     239      ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais,
     240     .                           lmt_altvolc_expl_glo)
     241      IF (ierr .NE. NF_NOERR) THEN
    235242         PRINT*, 'Pb de lecture pour les altitudes volcan'
    236243         CALL exit(1)

  • LMDZ6/trunk/libf/phylmd/Dust/phys_output_write_spl_mod.F90

    r5075 r5084  
    390390    USE ioipsl, ONLY: histend, histsync
    391391    USE iophy, ONLY: set_itau_iophy, histwrite_phy
    392     USE lmdz_netcdf, ONLY: nf90_fill_real
     392    USE netcdf, ONLY: nf90_fill_real
    393393    ! ug Pour les sorties XIOS
    394394    USE lmdz_xios, ONLY: xios_update_calendar, using_xios

  • LMDZ6/trunk/libf/phylmd/Dust/phytracr_spl_mod.F90

    r5075 r5084  
    14411441       call abort_gcm('phytracr_mod', 'pb in ok_chimdust 1 SCDU',1)
    14421442          endif
    1443           if ( (id_codu <= 0) .or. ( id_fine<=0)  ) then
     1443          if ( (id_codu .le. 0) .or. ( id_fine.le.0)  ) then 
    14441444          call abort_gcm('phytracr_mod', 'pb in ok_chimdust 1',1)
    14451445          endif
    … …  
    24372437      ENDDO
    24382438      ENDDO
    2439       IF (iflag_conv==2) THEN
     2439      IF (iflag_conv.EQ.2) THEN
    24402440! Tiedke
    24412441      CALL cltrac_spl(pdtphys,coefh,yu1,yv1,t_seri,tmp_var,  &
    24422442                 aux_var2,paprs,pplay,aux_var3)
    24432443
    2444       ELSE IF (iflag_conv>=3) THEN
     2444      ELSE IF (iflag_conv.GE.3) THEN
    24452445!KE
    24462446      CALL cltrac(pdtphys, coefh,t_seri,tmp_var,aux_var2,paprs,pplay,  &
    … …  
    24942494
    24952495
    2496       IF (iflag_conv>=3) THEN
     2496      IF (iflag_conv.GE.3) THEN
    24972497
    24982498      IF (logitime) THEN
    … …  
    27862786
    27872787
    2788       IF (iflag_conv==2) THEN
     2788      IF (iflag_conv.EQ.2) THEN
    27892789
    27902790      IF (logitime) THEN
    … …  
    28392839      print *,'iflag_conv bef incloud',iflag_conv
    28402840
    2841         IF (iflag_conv==2) THEN
     2841        IF (iflag_conv.EQ.2) THEN
    28422842! Tiedke
    28432843      CALL incloud_scav(.false.,qmin,qmax,masse,henry,kk,prfl,          &
    … …  
    28772877!     .                                  his_dhbclsc,his_dhbccon,tr_seri)
    28782878
    2879         IF (iflag_conv==2) THEN
     2879        IF (iflag_conv.EQ.2) THEN
    28802880! Tiedke
    28812881
    … …  
    29912991!     .                                                 dtrconv,tr_seri)
    29922992! -------------------------------------------------------------     
    2993         IF (iflag_conv==2) THEN
     2993        IF (iflag_conv.EQ.2) THEN
    29942994! Tiedke
    29952995         CALL trconvect(pplay,t_seri,pdtphys,pmfu,pmfd,pen_u,pde_u,  &
    … …  
    30003000         ENDDO
    30013001
    3002         ELSE IF (iflag_conv>=3) THEN
     3002        ELSE IF (iflag_conv.GE.3) THEN
    30033003! KE
    30043004         print *,'JE: KE in phytracr_spl'
    … …  
    31643164
    31653165
    3166        IF (iflag_conv>=3) THEN
     3166       IF (iflag_conv.GE.3) THEN
    31673167       IF (logitime) THEN
    31683168       CALL SYSTEM_CLOCK(COUNT=clock_start)
    … …  
    31953195       ql_incl = ql_incloud_ref
    31963196! choix du lessivage
    3197       IF (iflag_lscav == 3 .OR. iflag_lscav == 4) THEN
     3197      IF (iflag_lscav .EQ. 3 .OR. iflag_lscav .EQ. 4) THEN
    31983198      !IF (.false.) THEN  ! test #DFB (Binta) sans lsc_scav_spl
    31993199        print *,'JE iflag_lscav',iflag_lscav
    … …  
    33623362      CALL satellite_out_spla(jD_cur,jH_cur,pdtphys,rlat,rlon,   &
    33633363                              masque_aqua_cur, masque_terra_cur )
    3364       IF (jH_cur-pdtphys/86400. < 0.) THEN
     3364      IF (jH_cur-pdtphys/86400. .LT. 0.) THEN
    33653365       !new utc day: put in 0 everything
    33663366!JE20150518<<
    … …  
    34703470      ENDDO
    34713471
    3472       IF (jH_cur+pdtphys/86400. >= 1.) THEN
     3472      IF (jH_cur+pdtphys/86400. .GE. 1.) THEN 
    34733473!          print *,'last step of the day'
    34743474          DO i=1,klon
    3475                IF (masque_aqua(i)> 0) THEN
     3475               IF (masque_aqua(i).GT. 0) THEN
    34763476                   aod550_aqua(i)=aod550_aqua(i)/masque_aqua(i)
    34773477                   aod670_aqua(i)=aod670_aqua(i)/masque_aqua(i)
    … …  
    35063506                   aod865_dustsco_aqua(i)= -999.
    35073507               ENDIF
    3508                IF (masque_terra(i)> 0) THEN
     3508               IF (masque_terra(i).GT. 0) THEN
    35093509                   aod550_terra(i)=aod550_terra(i)/masque_terra(i)
    35103510                   aod670_terra(i)=aod670_terra(i)/masque_terra(i)
    … …  
    36353635      fluxss(:)=0.0
    36363636      DO i=1, klon
    3637          IF (iregion_ind(i)>0) THEN           ! LAND
     3637         IF (iregion_ind(i).GT.0) THEN           ! LAND
    36383638           ! SULFUR EMISSIONS
    36393639           fluxh2sff(i)= (lmt_so2ff_l(i)+lmt_so2ff_h(i))*frach2sofso2*  &       
    … …  
    36563656           fluxff(i)=fluxbcff(i)+fluxomff(i)
    36573657         ENDIF
    3658          IF (iregion_bb(i)>0) THEN           ! LAND
     3658         IF (iregion_bb(i).GT.0) THEN           ! LAND
    36593659           ! SULFUR EMISSIONS
    36603660           fluxso2bb(i) =scale_param_bb(iregion_bb(i)) * fracso2emis *  &
    … …  
    45154515      ENDIF
    45164516
    4517       IF (test_sca == 0 ) THEN
     4517      IF (test_sca .EQ. 0 ) THEN
    45184518        ! READ file!!
    45194519        call read_scalenc(filescaleparams,paramname_ind,            &
    … …  
    45564556
    45574557      jH_sca=jH_sca+pdtphys/(24.*3600.)
    4558       IF (jH_sca>(sca_resol)/24.) THEN
     4558      IF (jH_sca.GT.(sca_resol)/24.) THEN
    45594559          test_sca=0
    45604560          jH_sca=jH_ini
    … …  
    45684568      USE mod_grid_phy_lmdz
    45694569      USE mod_phys_lmdz_para
    4570       USE lmdz_netcdf, ONLY:nf_open,nf_close,nf_inq_varid,nf_nowrite,nf_noerr,nf90_get_var
    45714570      IMPLICIT NONE
     4571
     4572      include "netcdf.inc"
    45724573
    45734574      CHARACTER*800 filescaleparams
    … …  
    45884589          !nci=NCOPN(trim(adjustl(filescaleparams)),NCNOWRIT,rcode)
    45894590         ierr = NF_OPEN (trim(adjustl(filescaleparams)),NF_NOWRITE, nid)
    4590           if (ierr == NF_NOERR) THEN
     4591          if (ierr .EQ. NF_NOERR) THEN
    45914592          debutread=step_sca
    45924593          countread=1
    … …  
    45974598            print *,varname
    45984599            ierr = NF_INQ_VARID (nid,trim(adjustl(varname)), nvarid)
    4599             ierr = nf90_get_var (nid, nvarid, auxreal, debutread, countread)
    4600             IF (ierr /= NF_NOERR) THEN
     4600            ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debutread,          &
     4601                         countread, auxreal)
     4602            IF (ierr .NE. NF_NOERR) THEN
    46014603             PRINT*, 'Pb de lecture pour modvalues'
    46024604       print *,'JE  scale_var, step_sca',trim(adjustl(varname)),step_sca

  • LMDZ6/trunk/libf/phylmd/Dust/read_dust.F

    r5075 r5084  
    33      USE mod_grid_phy_lmdz
    44      USE mod_phys_lmdz_para
    5       USE lmdz_netcdf, ONLY:nf90_get_var
    65      IMPLICIT NONE
    76c
    87      INCLUDE "dimensions.h"
    98      INCLUDE "paramet.h"
     9      INCLUDE "netcdf.inc"
    1010c
    1111      INTEGER step, nbjour
    … …  
    4545c
    4646      start(3)=step
    47 
    48       status=nf90_get_var(ncid1,varid1,dust_nc_glo,start,count)
    49 
     47c
     48#ifdef NC_DOUBLE
     49!      status=NF_GET_VARA_DOUBLE(ncid1,varid1,start,count,dust_nc)
     50      status=NF_GET_VARA_DOUBLE(ncid1,varid1,start,count,dust_nc_glo)
     51#else
     52!      status=NF_GET_VARA_REAL(ncid1,varid1,start,count,dust_nc)
     53      status=NF_GET_VARA_REAL(ncid1,varid1,start,count,dust_nc_glo)
     54#endif
     55c
    5056!      call correctbid(iim,jjp1,dust_nc)
    5157      call correctbid(nbp_lon,nbp_lat,dust_nc_glo)

  • LMDZ6/trunk/libf/phylmd/Dust/read_surface.F90

    r5075 r5084  
    1010       USE mod_phys_lmdz_para
    1111       USE iophy
    12        USE lmdz_netcdf, ONLY:nf_inq_varid,nf_noerr,nf90_get_var
     12!       USE netcdf
    1313       IMPLICIT NONE
    1414
     15       INCLUDE "netcdf.inc"
    1516       INCLUDE "dimensions.h"
    1617       INCLUDE "paramet.h"
    … …  
    6465!       print *,'stat,i',status,i,outcycle,aux4s
    6566!       print *,'ifclause',status.NE. NF_NOERR ,outcycle == .false.
    66        IF ((.not.(status/= NF_NOERR) ).and.( .not. outcycle )) THEN
     67       IF ((.not.(status.NE. NF_NOERR) ).and.( .not. outcycle )) THEN
    6768         outcycle=.true.
    6869         latstr=aux4s
    … …  
    7475      varid=NCVID(ncid,latstr,rcode)
    7576
    76           status=nf90_get_var(ncid,varid,lats_glo,startj,endj)
     77#ifdef NC_DOUBLE
     78          status=NF_GET_VARA_DOUBLE(ncid,varid,startj,endj,lats_glo)
     79#else
     80          status=NF_GET_VARA_REAL(ncid,varid,startj,endj,lats_glo)
     81#endif
    7782!      print *,latstr,varid,status,jjp1,rcode
    7883!      IF (status .NE. NF_NOERR) print*,'NOOOOOOO'
    … …  
    108113! Lecture
    109114! -----------------------
    110           status=nf90_get_var(ncid,varid,tmp_dyn_glo,start,count)
     115#ifdef NC_DOUBLE
     116          status=NF_GET_VARA_DOUBLE(ncid,varid,start,count,tmp_dyn_glo)
     117#else
     118          status=NF_GET_VARA_REAL(ncid,varid,start,count,tmp_dyn_glo)
     119#endif
    111120
    112121!      call dump2d(iip1,jjp1,tmp_dyn,'tmp_dyn   ')

  • LMDZ6/trunk/libf/phylmd/Dust/read_vent.F

    r5075 r5084  
    33      USE mod_grid_phy_lmdz
    44      USE mod_phys_lmdz_para
    5       USE lmdz_netcdf, ONLY: nf90_get_var
    65!      USE write_field_phy
    76      IMPLICIT NONE
    … …  
    98c       INCLUDE "dimphy.h"
    109      INCLUDE "paramet.h"
     10      INCLUDE "netcdf.inc"
    1111c
    1212      INTEGER step, nbjour
    … …  
    5151c
    5252      start(3)=step
    53 
    54       status=nf90_get_var(ncidu1,varidu1,u10m_nc_glo,start,count)
    55 
    56       status=nf90_get_var(ncidv1,varidv1,v10m_nc_glo,start,count)
    57 
     53c
     54#ifdef NC_DOUBLE
     55!      status=NF_GET_VARA_DOUBLE(ncidu1,varidu1,start,count,u10m_nc)
     56      status=NF_GET_VARA_DOUBLE(ncidu1,varidu1,start,count,u10m_nc_glo)
     57#else
     58!      status=NF_GET_VARA_REAL(ncidu1,varidu1,start,count,u10m_nc)
     59      status=NF_GET_VARA_REAL(ncidu1,varidu1,start,count,u10m_nc_glo)
     60#endif
     61!       print *,status
     62c
     63#ifdef NC_DOUBLE
     64!      status=NF_GET_VARA_DOUBLE(ncidv1,varidv1,start,count,v10m_nc)
     65      status=NF_GET_VARA_DOUBLE(ncidv1,varidv1,start,count,v10m_nc_glo)
     66#else
     67!      status=NF_GET_VARA_REAL(ncidv1,varidv1,start,count,v10m_nc)
     68      status=NF_GET_VARA_REAL(ncidv1,varidv1,start,count,v10m_nc_glo)
     69#endif
     70c
    5871
    5972!      print *,'beforebidcor u10m_nc', u10m_nc(1,jjp1)
    … …  
    118131      do l=1,nl
    119132         do i=2,iim-1
    120             if(abs(x(i,l))>1.e10) then
     133            if(abs(x(i,l)).gt.1.e10) then
    121134               zz=0.5*(x(i-1,l)+x(i+1,l))
    122135c              print*,'correction ',i,l,x(i,l),zz
Note: See TracChangeset for help on using the changeset viewer.