Ignore:
Timestamp:
Jul 18, 2024, 4:57:05 PM (6 months ago)
Author:
abarral
Message:

Remove all NC_DOUBLE uses outside of lmdz_netcdf.F90 (except in obsolete/, which I hope we'll ditch soon...)
Note: make sure to check convergence at some point, it's possible that we've messed up some when replacing nf_* by nf90_* calls
(lint) replace obsolete logical operators along the way

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

Legend:

Unmodified
Added
Removed

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

    r4593 r5073  
    2424      INTEGER debut(2),epais(2)
    2525!
    26       IF (jour.LT.0 .OR. jour.GT.(360-1)) THEN
    27          IF (jour.GT.(360-1).AND.jour.LE.367) THEN
     26      IF (jour<0 .OR. jour>(360-1)) THEN
     27         IF (jour>(360-1).AND.jour<=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.ne.NF_NOERR) then
     37      if (ierr/=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 !nhl #ifdef NC_DOUBLE
    52       ierr = NF_GET_VARA_DOUBLE (nid1, nvarid, debut, epais, lmt_bcff)
     51      ierr = nf90_get_var(nid1, nvarid, lmt_bcff, debut, epais)
    5352!      print *,'IERR = ',ierr
    5453!      print *,'NF_NOERR = ',NF_NOERR
    5554!      print *,'debut = ',debut
    5655!      print *,'epais = ',epais
    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
     56      IF (ierr /= NF_NOERR) THEN
    6157         PRINT*, 'Pb de lecture pour les sources BC'
    6258         CALL exit(1)
    … …  
    6561!
    6662      ierr = NF_INQ_VARID (nid1, "BCBB", nvarid)
    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
     63      ierr = nf90_get_var(nid1, nvarid, lmt_bcbb, debut, epais)
     64      IF (ierr /= NF_NOERR) THEN
    7365         PRINT*, 'Pb de lecture pour les sources BC-biomass'
    7466         CALL exit(1)
    … …  
    7769!
    7870      ierr = NF_INQ_VARID (nid1, "BCBL", nvarid)
    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
     71      ierr = nf90_get_var(nid1, nvarid, lmt_bcbbl, debut, epais)
     72      IF (ierr /= NF_NOERR) THEN
    8573         PRINT*, 'Pb de lecture pour les sources BC low'
    8674         CALL exit(1)
    … …  
    8977!
    9078      ierr = NF_INQ_VARID (nid1, "BCBH", nvarid)
    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
     79      ierr = nf90_get_var (nid1, nvarid, lmt_bcbbh, debut, epais)
     80      IF (ierr /= NF_NOERR) THEN
    9781         PRINT*, 'Pb de lecture pour les sources BC high'
    9882         CALL exit(1)
    … …  
    10084!
    10185      ierr = NF_INQ_VARID (nid1, "TERP", nvarid)
    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
     86      ierr = nf90_get_var (nid1, nvarid, lmt_terp, debut, epais)
     87      IF (ierr /= NF_NOERR) THEN
    10888         PRINT*, 'Pb de lecture pour les sources Terpene'
    10989         CALL exit(1)
    … …  
    11292!
    11393      ierr = NF_INQ_VARID (nid1, "BC_penner", nvarid)
    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
     94      ierr = nf90_get_var (nid1, nvarid, lmt_bc_penner, debut,  epais)
     95      IF (ierr /= NF_NOERR) THEN
    12296         PRINT*, 'Pb de lecture pour les sources BC Penner'
    12397         CALL exit(1)
    … …  
    126100!
    127101      ierr = NF_INQ_VARID (nid1, "OMFF", nvarid)
    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
     102      ierr = nf90_get_var (nid1, nvarid, lmt_omff, debut, epais)
     103      IF (ierr /= NF_NOERR) THEN
    134104         PRINT*, 'Pb de lecture pour les sources om-ifossil'
    135105         CALL exit(1)

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

    r4593 r5073  
    2424      INTEGER debut(2),epais(2)
    2525c
    26       IF (jour.LT.0 .OR. jour.GT.(360-1)) THEN
    27          IF ((jour.GT.(360-1)) .AND. (jour.LE.367)) THEN
     26      IF (jour<0 .OR. jour>(360-1)) THEN
     27         IF ((jour>(360-1)) .AND. (jour<=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.ne.NF_NOERR) then
     37      if (ierr/=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 cnhl #ifdef NC_DOUBLE
    51       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2volc)
    52 cnhl #else
    53 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2volc)
    54 cnhl #endif
    55       IF (ierr .NE. NF_NOERR) THEN
     50      ierr = nf90_get_var(nid, nvarid, lmt_so2volc, debut, epais)
     51      IF (ierr /= NF_NOERR) THEN
    5652         PRINT*, 'Pb de lecture pour les sources so2 volcan'
    5753         CALL exit(1)
    … …  
    5955c
    6056      ierr = NF_INQ_VARID (nid, "ALTI", nvarid)
    61 cnhl #ifdef NC_DOUBLE
    62       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_altvolc)
    63 cnhl #else
    64 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_altvolc)
    65 cnhl #endif
    66       IF (ierr .NE. NF_NOERR) THEN
     57      ierr = nf90_get_var(nid, nvarid, lmt_altvolc, debut, epais)
     58      IF (ierr /= NF_NOERR) THEN
    6759         PRINT*, 'Pb de lecture pour les altitudes volcan'
    6860         CALL exit(1)
    … …  
    7264c
    7365      ierr = NF_INQ_VARID (nid, "SO2ED95L", nvarid)
    74 cnhl #ifdef NC_DOUBLE
    75       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2b)
    76 cnhl #else
    77 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2b)
    78 cnhl #endif
    79       IF (ierr .NE. NF_NOERR) THEN
     66      ierr = nf90_get_var(nid, nvarid, lmt_so2b, debut, epais)
     67      IF (ierr /= NF_NOERR) THEN
    8068         PRINT*, 'Pb de lecture pour les sources so2 edgar low'
    8169         CALL exit(1)
    … …  
    8371c
    8472      ierr = NF_INQ_VARID (nid, "SO2ED95H", nvarid)
    85 cnhl #ifdef NC_DOUBLE
    86       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2h)
    87 cnhl #else
    88 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2h)
    89 cnhl #endif
    90       IF (ierr .NE. NF_NOERR) THEN
     73      ierr = nf90_get_var(nid, nvarid, lmt_so2h, debut, epais)
     74      IF (ierr /= NF_NOERR) THEN
    9175         PRINT*, 'Pb de lecture pour les sources so2 edgar high'
    9276         CALL exit(1)
    … …  
    9680c
    9781      ierr = NF_INQ_VARID (nid, "SO2H", nvarid)
    98 cnhl #ifdef NC_DOUBLE
    99       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2h)
    100 cnhl #else
    101 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2h)
    102 cnhl #endif
    103       IF (ierr .NE. NF_NOERR) THEN
     82      ierr = nf90_get_var(nid, nvarid, lmt_so2h, debut, epais)
     83      IF (ierr /= NF_NOERR) THEN
    10484         PRINT*, 'Pb de lecture pour les sources so2 haut'
    10585         CALL exit(1)
    … …  
    10787c
    10888      ierr = NF_INQ_VARID (nid, "SO2B", nvarid)
    109 cnhl #ifdef NC_DOUBLE
    110       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2b)
    111 cnhl #else
    112 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2b)
    113 cnhl #endif
    114       IF (ierr .NE. NF_NOERR) THEN
     89      ierr = nf90_get_var(nid, nvarid, lmt_so2b, debut, epais)
     90      IF (ierr /= NF_NOERR) THEN
    11591         PRINT*, 'Pb de lecture pour les sources so2 bas'
    11692         CALL exit(1)
    … …  
    12096c
    12197      ierr = NF_INQ_VARID (nid, "SO2BB", nvarid)
    122 cnhl #ifdef NC_DOUBLE
    123       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2bb)
    124 cnhl #else
    125 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2bb)
    126 cnhl #endif
    127       IF (ierr .NE. NF_NOERR) THEN
     98      ierr = nf90_get_var(nid, nvarid, lmt_so2bb, debut, epais)
     99      IF (ierr /= NF_NOERR) THEN
    128100         PRINT*, 'Pb de lecture pour les sources so2 bb'
    129101         CALL exit(1)
    … …  
    131103c
    132104      ierr = NF_INQ_VARID (nid, "SO2BA", nvarid)
    133 cnhl #ifdef NC_DOUBLE
    134       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_so2ba)
    135 cnhl #else
    136 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_so2ba)
    137 cnhl #endif
    138       IF (ierr .NE. NF_NOERR) THEN
     105      ierr = nf90_get_var(nid, nvarid, lmt_so2ba, debut, epais)
     106      IF (ierr /= NF_NOERR) THEN
    139107         PRINT*, 'Pb de lecture pour les sources so2 bateau'
    140108         CALL exit(1)
    … …  
    142110c
    143111      ierr = NF_INQ_VARID (nid, "DMSB", nvarid)
    144 cnhl #ifdef NC_DOUBLE
    145       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_dmsbio)
    146 cnhl #else
    147 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsbio)
    148 cnhl #endif
    149       IF (ierr .NE. NF_NOERR) THEN
     112      ierr = nf90_get_var(nid, nvarid, lmt_dmsbio, debut, epais)
     113      IF (ierr /= NF_NOERR) THEN
    150114         PRINT*, 'Pb de lecture pour les sources dms bio'
    151115         CALL exit(1)
    … …  
    153117c
    154118      ierr = NF_INQ_VARID (nid, "H2SB", nvarid)
    155 cnhl #ifdef NC_DOUBLE
    156       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_h2sbio)
    157 cnhl #else
    158 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_h2sbio)
    159 cnhl #endif
    160       IF (ierr .NE. NF_NOERR) THEN
     119      ierr = nf90_get_var(nid, nvarid, lmt_h2sbio, debut, epais)
     120      IF (ierr /= NF_NOERR) THEN
    161121         PRINT*, 'Pb de lecture pour les sources h2s bio'
    162122         CALL exit(1)
    163123      ENDIF
    164124c
    165       IF (flag_dms.EQ.1) THEN
     125      IF (flag_dms==1) THEN
    166126c
    167127      ierr = NF_INQ_VARID (nid, "DMSL", nvarid)
    168 cnhl #ifdef NC_DOUBLE
    169       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_dms)
    170 cnhl #else
    171 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dms)
    172 cnhl #endif
    173       IF (ierr .NE. NF_NOERR) THEN
     128      ierr = nf90_get_var(nid, nvarid, lmt_dms, debut, epais)
     129      IF (ierr /= NF_NOERR) THEN
    174130         PRINT*, 'Pb de lecture pour les sources dms liss'
    175131         CALL exit(1)
    176132      ENDIF
    177133c
    178       ELSEIF (flag_dms.EQ.2) THEN
     134      ELSEIF (flag_dms==2) THEN
    179135c
    180136      ierr = NF_INQ_VARID (nid, "DMSW", nvarid)
    181 cnhl #ifdef NC_DOUBLE
    182       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut, epais, lmt_dms)
    183 cnhl #else
    184 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dms)
    185 cnhl #endif
    186       IF (ierr .NE. NF_NOERR) THEN
     137      ierr = nf90_get_var(nid, nvarid, lmt_dms, debut, epais)
     138      IF (ierr /= NF_NOERR) THEN
    187139         PRINT*, 'Pb de lecture pour les sources dms wann'
    188140         CALL exit(1)
    189141      ENDIF
    190142c
    191       ELSEIF (flag_dms.EQ.3) THEN
     143      ELSEIF (flag_dms==3) THEN
    192144c
    193145      ierr = NF_INQ_VARID (nid, "DMSC1", nvarid)
    194 cnhl #ifdef NC_DOUBLE
    195       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
    196 cnhl #else
    197 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
    198 cnhl #endif
    199       IF (ierr .NE. NF_NOERR) THEN
     146      ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
     147      IF (ierr /= NF_NOERR) THEN
    200148         PRINT*, 'Pb de lecture pour les sources dmsconc old'
    201149         CALL exit(1)
    202150      ENDIF
    203151c
    204       ELSEIF (flag_dms.EQ.4) THEN
     152      ELSEIF (flag_dms==4) THEN
    205153c
    206154      ierr = NF_INQ_VARID (nid, "DMSC2", nvarid)
    207 cnhl #ifdef NC_DOUBLE
    208       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
    209 cnhl #else
    210 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
    211 cnhl #endif
    212       IF (ierr .NE. NF_NOERR) THEN
     155      ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
     156      IF (ierr /= NF_NOERR) THEN
    213157         PRINT*, 'Pb de lecture pour les sources dms conc 2'
    214158         CALL exit(1)
    215159      ENDIF
    216160c
    217       ELSEIF (flag_dms.EQ.5) THEN
     161      ELSEIF (flag_dms==5) THEN
    218162c
    219163      ierr = NF_INQ_VARID (nid, "DMSC3", nvarid)
    220 cnhl #ifdef NC_DOUBLE
    221       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
    222 cnhl #else
    223 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
    224 cnhl #endif
    225       IF (ierr .NE. NF_NOERR) THEN
     164      ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
     165      IF (ierr /= NF_NOERR) THEN
    226166         PRINT*, 'Pb de lecture pour les sources dms conc 3'
    227167         CALL exit(1)
    228168      ENDIF
    229169c
    230       ELSEIF (flag_dms.EQ.6) THEN
     170      ELSEIF (flag_dms==6) THEN
    231171c
    232172      ierr = NF_INQ_VARID (nid, "DMSC4", nvarid)
    233 cnhl #ifdef NC_DOUBLE
    234       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
    235 cnhl #else
    236 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
    237 cnhl #endif
    238       IF (ierr .NE. NF_NOERR) THEN
     173      ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
     174      IF (ierr /= NF_NOERR) THEN
    239175         PRINT*, 'Pb de lecture pour les sources dms conc 4'
    240176         CALL exit(1)
    241177      ENDIF
    242178c
    243       ELSEIF (flag_dms.EQ.7) THEN
     179      ELSEIF (flag_dms==7) THEN
    244180c
    245181      ierr = NF_INQ_VARID (nid, "DMSC5", nvarid)
    246 cnhl #ifdef NC_DOUBLE
    247       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
    248 cnhl #else
    249 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
    250 cnhl #endif
    251       IF (ierr .NE. NF_NOERR) THEN
     182      ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
     183      IF (ierr /= NF_NOERR) THEN
    252184         PRINT*, 'Pb de lecture pour les sources dms conc 5'
    253185         CALL exit(1)
    254186      ENDIF
    255187c
    256       ELSEIF (flag_dms.EQ.8) THEN
     188      ELSEIF (flag_dms==8) THEN
    257189c
    258190      ierr = NF_INQ_VARID (nid, "DMSC6", nvarid)
    259 cnhl #ifdef NC_DOUBLE
    260       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
    261 cnhl #else
    262 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
    263 cnhl #endif
    264       IF (ierr .NE. NF_NOERR) THEN
     191      ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
     192      IF (ierr /= NF_NOERR) THEN
    265193         PRINT*, 'Pb de lecture pour les sources dms conc 6'
    266194         CALL exit(1)
    267195      ENDIF
    268196c
    269       ELSEIF (flag_dms.EQ.9) THEN
     197      ELSEIF (flag_dms==9) THEN
    270198c
    271199      ierr = NF_INQ_VARID (nid, "DMSC7", nvarid)
    272 cnhl #ifdef NC_DOUBLE
    273       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
    274 cnhl #else
    275 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
    276 cnhl #endif
    277       IF (ierr .NE. NF_NOERR) THEN
     200      ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
     201      IF (ierr /= NF_NOERR) THEN
    278202         PRINT*, 'Pb de lecture pour les sources dms conc 7'
    279203         CALL exit(1)
    280204      ENDIF
    281205c
    282       ELSEIF (flag_dms.EQ.10) THEN
     206      ELSEIF (flag_dms==10) THEN
    283207c
    284208      ierr = NF_INQ_VARID (nid, "DMSC8", nvarid)
    285 cnhl #ifdef NC_DOUBLE
    286       ierr = NF_GET_VARA_DOUBLE (nid, nvarid, debut,epais,lmt_dmsconc)
    287 cnhl #else
    288 cnhl       ierr = NF_GET_VARA_REAL (nid, nvarid, debut, epais, lmt_dmsconc)
    289 cnhl #endif
    290       IF (ierr .NE. NF_NOERR) THEN
     209      ierr = nf90_get_var(nid, nvarid, lmt_dmsconc, debut, epais)
     210      IF (ierr /= NF_NOERR) THEN
    291211         PRINT*, 'Pb de lecture pour les sources dms conc 8'
    292212         CALL exit(1)
    … …  
    302222      ierr = NF_CLOSE(nid)
    303223c
    304       IF (flag_dms.LE.2) THEN
     224      IF (flag_dms<=2) THEN
    305225      DO i=1, klon
    306226         lmt_dmsconc(i)=0.0

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

    r4593 r5073  
    4646      start(3)=step
    4747c
    48 #ifdef NC_DOUBLE
    4948!      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
     49      status=nf90_get_var(ncid1,varid1,dust_nc_glo,start,count)
    5550c
    5651!      call correctbid(iim,jjp1,dust_nc)

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

    r4593 r5073  
    6565!       print *,'stat,i',status,i,outcycle,aux4s
    6666!       print *,'ifclause',status.NE. NF_NOERR ,outcycle == .false.
    67        IF ((.not.(status.NE. NF_NOERR) ).and.( .not. outcycle )) THEN
     67       IF ((.not.(status/= NF_NOERR) ).and.( .not. outcycle )) THEN
    6868         outcycle=.true.
    6969         latstr=aux4s
    … …  
    7575      varid=NCVID(ncid,latstr,rcode)
    7676
    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
     77          status=nf90_get_var(ncid,varid,lats_glo,startj,endj)
    8278!      print *,latstr,varid,status,jjp1,rcode
    8379!      IF (status .NE. NF_NOERR) print*,'NOOOOOOO'
    … …  
    113109! Lecture
    114110! -----------------------
    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
     111          status=nf90_get_var(ncid,varid,tmp_dyn_glo,start,count)
    120112
    121113!      call dump2d(iip1,jjp1,tmp_dyn,'tmp_dyn   ')

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

    r4593 r5073  
    5151c
    5252      start(3)=step
    53 c
    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
    62 c
    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
    70 c
     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
    7158
    7259!      print *,'beforebidcor u10m_nc', u10m_nc(1,jjp1)
    … …  
    131118      do l=1,nl
    132119         do i=2,iim-1
    133             if(abs(x(i,l)).gt.1.e10) then
     120            if(abs(x(i,l))>1.e10) then
    134121               zz=0.5*(x(i-1,l)+x(i+1,l))
    135122c              print*,'correction ',i,l,x(i,l),zz
Note: See TracChangeset for help on using the changeset viewer.