Ignore:
Timestamp:
Jul 20, 2024, 4:30:31 PM (3 months ago)
Author:
abarral
Message:

Remove all managed netcdf77 usage (excluding external: rrtm, ecrad)
Remove call to netcdf.inc
Replace USE lmdz_netcdf by USE netcdf
Replace lmdz_netcdf.F90 by lmdz_netcdf_format.F90

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/old_1DUTILS_read_interp.h

    r5086 r5088  
    146146!program reading forcings of the TWP-ICE experiment
    147147
    148         use lmdz_netcdf, ONLY: nf_open,nf_nowrite,nf_noerr,nf_strerror,nf_inq_varid,nf90_get_var,&
    149             nf_inq_dimid,nf_inq_dimlen
     148        use netcdf, ONLY: nf90_open,nf90_nowrite,nf90_noerr,nf90_strerror,nf90_inq_varid,nf90_get_var,&
     149            nf90_inq_dimid,nf90_inquire_dimension
    150150
    151151
    … …  
    177177      integer var3didin(nbvar3d)
    178178
    179       ierr = NF_OPEN(fich_twpice,NF_NOWRITE,nid)
    180       if (ierr.NE.NF_NOERR) then
     179      ierr = nf90_open(fich_twpice,nf90_nowrite,nid)
     180      if (ierr.NE.nf90_noerr) then
    181181         write(*,*) 'ERROR: Pb opening forcings cdf file '
    182          write(*,*) NF_STRERROR(ierr)
     182         write(*,*) nf90_strerror(ierr)
    183183         stop ""
    184184      endif
    185185
    186       ierr=NF_INQ_VARID(nid,"lat",var3didin(1))
    187          if(ierr/=NF_NOERR) then
    188            write(*,*) NF_STRERROR(ierr)
     186      ierr=nf90_inq_varid(nid,"lat",var3didin(1))
     187         if(ierr/=nf90_noerr) then
     188           write(*,*) nf90_strerror(ierr)
    189189           stop 'lat'
    190190         endif
    191191     
    192        ierr=NF_INQ_VARID(nid,"lon",var3didin(2))
    193          if(ierr/=NF_NOERR) then
    194            write(*,*) NF_STRERROR(ierr)
     192       ierr=nf90_inq_varid(nid,"lon",var3didin(2))
     193         if(ierr/=nf90_noerr) then
     194           write(*,*) nf90_strerror(ierr)
    195195           stop 'lon'
    196196         endif
    197197
    198        ierr=NF_INQ_VARID(nid,"alt",var3didin(3))
    199          if(ierr/=NF_NOERR) then
    200            write(*,*) NF_STRERROR(ierr)
     198       ierr=nf90_inq_varid(nid,"alt",var3didin(3))
     199         if(ierr/=nf90_noerr) then
     200           write(*,*) nf90_strerror(ierr)
    201201           stop 'alt'
    202202         endif
    203203
    204       ierr=NF_INQ_VARID(nid,"phis",var3didin(4))
    205          if(ierr/=NF_NOERR) then
    206            write(*,*) NF_STRERROR(ierr)
     204      ierr=nf90_inq_varid(nid,"phis",var3didin(4))
     205         if(ierr/=nf90_noerr) then
     206           write(*,*) nf90_strerror(ierr)
    207207           stop 'phis'
    208208         endif
    209209
    210       ierr=NF_INQ_VARID(nid,"T",var3didin(5))
    211          if(ierr/=NF_NOERR) then
    212            write(*,*) NF_STRERROR(ierr)
     210      ierr=nf90_inq_varid(nid,"T",var3didin(5))
     211         if(ierr/=nf90_noerr) then
     212           write(*,*) nf90_strerror(ierr)
    213213           stop 'T'
    214214         endif
    215215
    216       ierr=NF_INQ_VARID(nid,"q",var3didin(6))
    217          if(ierr/=NF_NOERR) then
    218            write(*,*) NF_STRERROR(ierr)
     216      ierr=nf90_inq_varid(nid,"q",var3didin(6))
     217         if(ierr/=nf90_noerr) then
     218           write(*,*) nf90_strerror(ierr)
    219219           stop 'q'
    220220         endif
    221221
    222       ierr=NF_INQ_VARID(nid,"u",var3didin(7))
    223          if(ierr/=NF_NOERR) then
    224            write(*,*) NF_STRERROR(ierr)
     222      ierr=nf90_inq_varid(nid,"u",var3didin(7))
     223         if(ierr/=nf90_noerr) then
     224           write(*,*) nf90_strerror(ierr)
    225225           stop 'u'
    226226         endif
    227227
    228       ierr=NF_INQ_VARID(nid,"v",var3didin(8))
    229          if(ierr/=NF_NOERR) then
    230            write(*,*) NF_STRERROR(ierr)
     228      ierr=nf90_inq_varid(nid,"v",var3didin(8))
     229         if(ierr/=nf90_noerr) then
     230           write(*,*) nf90_strerror(ierr)
    231231           stop 'v'
    232232         endif
    233233
    234       ierr=NF_INQ_VARID(nid,"omega",var3didin(9))
    235          if(ierr/=NF_NOERR) then
    236            write(*,*) NF_STRERROR(ierr)
     234      ierr=nf90_inq_varid(nid,"omega",var3didin(9))
     235         if(ierr/=nf90_noerr) then
     236           write(*,*) nf90_strerror(ierr)
    237237           stop 'omega'
    238238         endif
    239239
    240       ierr=NF_INQ_VARID(nid,"div",var3didin(10))
    241          if(ierr/=NF_NOERR) then
    242            write(*,*) NF_STRERROR(ierr)
     240      ierr=nf90_inq_varid(nid,"div",var3didin(10))
     241         if(ierr/=nf90_noerr) then
     242           write(*,*) nf90_strerror(ierr)
    243243           stop 'div'
    244244         endif
    245245
    246       ierr=NF_INQ_VARID(nid,"T_adv_h",var3didin(11))
    247          if(ierr/=NF_NOERR) then
    248            write(*,*) NF_STRERROR(ierr)
     246      ierr=nf90_inq_varid(nid,"T_adv_h",var3didin(11))
     247         if(ierr/=nf90_noerr) then
     248           write(*,*) nf90_strerror(ierr)
    249249           stop 'T_adv_h'
    250250         endif
    251251
    252       ierr=NF_INQ_VARID(nid,"T_adv_v",var3didin(12))
    253          if(ierr/=NF_NOERR) then
    254            write(*,*) NF_STRERROR(ierr)
     252      ierr=nf90_inq_varid(nid,"T_adv_v",var3didin(12))
     253         if(ierr/=nf90_noerr) then
     254           write(*,*) nf90_strerror(ierr)
    255255           stop 'T_adv_v'
    256256         endif
    257257
    258       ierr=NF_INQ_VARID(nid,"q_adv_h",var3didin(13))
    259          if(ierr/=NF_NOERR) then
    260            write(*,*) NF_STRERROR(ierr)
     258      ierr=nf90_inq_varid(nid,"q_adv_h",var3didin(13))
     259         if(ierr/=nf90_noerr) then
     260           write(*,*) nf90_strerror(ierr)
    261261           stop 'q_adv_h'
    262262         endif
    263263
    264       ierr=NF_INQ_VARID(nid,"q_adv_v",var3didin(14))
    265          if(ierr/=NF_NOERR) then
    266            write(*,*) NF_STRERROR(ierr)
     264      ierr=nf90_inq_varid(nid,"q_adv_v",var3didin(14))
     265         if(ierr/=nf90_noerr) then
     266           write(*,*) nf90_strerror(ierr)
    267267           stop 'q_adv_v'
    268268         endif
    269269
    270       ierr=NF_INQ_VARID(nid,"s",var3didin(15))
    271          if(ierr/=NF_NOERR) then
    272            write(*,*) NF_STRERROR(ierr)
     270      ierr=nf90_inq_varid(nid,"s",var3didin(15))
     271         if(ierr/=nf90_noerr) then
     272           write(*,*) nf90_strerror(ierr)
    273273           stop 's'
    274274         endif
    275275
    276       ierr=NF_INQ_VARID(nid,"s_adv_h",var3didin(16))
    277          if(ierr/=NF_NOERR) then
    278            write(*,*) NF_STRERROR(ierr)
     276      ierr=nf90_inq_varid(nid,"s_adv_h",var3didin(16))
     277         if(ierr/=nf90_noerr) then
     278           write(*,*) nf90_strerror(ierr)
    279279           stop 's_adv_h'
    280280         endif
    281281   
    282       ierr=NF_INQ_VARID(nid,"s_adv_v",var3didin(17))
    283          if(ierr/=NF_NOERR) then
    284            write(*,*) NF_STRERROR(ierr)
     282      ierr=nf90_inq_varid(nid,"s_adv_v",var3didin(17))
     283         if(ierr/=nf90_noerr) then
     284           write(*,*) nf90_strerror(ierr)
    285285           stop 's_adv_v'
    286286         endif
    287287
    288       ierr=NF_INQ_VARID(nid,"p_srf_aver",var3didin(18))
    289          if(ierr/=NF_NOERR) then
    290            write(*,*) NF_STRERROR(ierr)
     288      ierr=nf90_inq_varid(nid,"p_srf_aver",var3didin(18))
     289         if(ierr/=nf90_noerr) then
     290           write(*,*) nf90_strerror(ierr)
    291291           stop 'p_srf_aver'
    292292         endif
    293293
    294       ierr=NF_INQ_VARID(nid,"p_srf_center",var3didin(19))
    295          if(ierr/=NF_NOERR) then
    296            write(*,*) NF_STRERROR(ierr)
     294      ierr=nf90_inq_varid(nid,"p_srf_center",var3didin(19))
     295         if(ierr/=nf90_noerr) then
     296           write(*,*) nf90_strerror(ierr)
    297297           stop 'p_srf_center'
    298298         endif
    299299
    300       ierr=NF_INQ_VARID(nid,"T_srf",var3didin(20))
    301          if(ierr/=NF_NOERR) then
    302            write(*,*) NF_STRERROR(ierr)
     300      ierr=nf90_inq_varid(nid,"T_srf",var3didin(20))
     301         if(ierr/=nf90_noerr) then
     302           write(*,*) nf90_strerror(ierr)
    303303           stop 'T_srf'
    304304         endif
    … …  
    315315         
    316316         ierr = NF90_GET_VAR(nid,var3didin(1),lat)
    317          if(ierr/=NF_NOERR) then
    318             write(*,*) NF_STRERROR(ierr)
     317         if(ierr/=nf90_noerr) then
     318            write(*,*) nf90_strerror(ierr)
    319319            stop "getvarup"
    320320         endif
    … …  
    322322
    323323         ierr = NF90_GET_VAR(nid,var3didin(2),lon)
    324          if(ierr/=NF_NOERR) then
    325             write(*,*) NF_STRERROR(ierr)
     324         if(ierr/=nf90_noerr) then
     325            write(*,*) nf90_strerror(ierr)
    326326            stop "getvarup"
    327327         endif
    … …  
    329329 
    330330         ierr = NF90_GET_VAR(nid,var3didin(3),alt)
    331          if(ierr/=NF_NOERR) then
    332             write(*,*) NF_STRERROR(ierr)
     331         if(ierr/=nf90_noerr) then
     332            write(*,*) nf90_strerror(ierr)
    333333            stop "getvarup"
    334334         endif
    … …  
    336336 
    337337         ierr = NF90_GET_VAR(nid,var3didin(4),phis)
    338          if(ierr/=NF_NOERR) then
    339             write(*,*) NF_STRERROR(ierr)
     338         if(ierr/=nf90_noerr) then
     339            write(*,*) nf90_strerror(ierr)
    340340            stop "getvarup"
    341341         endif
    … …  
    343343         
    344344         ierr = NF90_GET_VAR(nid,var3didin(5),T)
    345          if(ierr/=NF_NOERR) then
    346             write(*,*) NF_STRERROR(ierr)
     345         if(ierr/=nf90_noerr) then
     346            write(*,*) nf90_strerror(ierr)
    347347            stop "getvarup"
    348348         endif
    … …  
    350350
    351351         ierr = NF90_GET_VAR(nid,var3didin(6),q)
    352          if(ierr/=NF_NOERR) then
    353             write(*,*) NF_STRERROR(ierr)
     352         if(ierr/=nf90_noerr) then
     353            write(*,*) nf90_strerror(ierr)
    354354            stop "getvarup"
    355355         endif
    … …  
    362362       enddo
    363363         ierr = NF90_GET_VAR(nid,var3didin(7),u)
    364          if(ierr/=NF_NOERR) then
    365             write(*,*) NF_STRERROR(ierr)
     364         if(ierr/=nf90_noerr) then
     365            write(*,*) nf90_strerror(ierr)
    366366            stop "getvarup"
    367367         endif
    … …  
    369369
    370370         ierr = NF90_GET_VAR(nid,var3didin(8),v)
    371          if(ierr/=NF_NOERR) then
    372             write(*,*) NF_STRERROR(ierr)
     371         if(ierr/=nf90_noerr) then
     372            write(*,*) nf90_strerror(ierr)
    373373            stop "getvarup"
    374374         endif
    … …  
    376376
    377377         ierr = NF90_GET_VAR(nid,var3didin(9),omega)
    378          if(ierr/=NF_NOERR) then
    379             write(*,*) NF_STRERROR(ierr)
     378         if(ierr/=nf90_noerr) then
     379            write(*,*) nf90_strerror(ierr)
    380380            stop "getvarup"
    381381         endif
    … …  
    389389
    390390         ierr = NF90_GET_VAR(nid,var3didin(10),div)
    391          if(ierr/=NF_NOERR) then
    392             write(*,*) NF_STRERROR(ierr)
     391         if(ierr/=nf90_noerr) then
     392            write(*,*) nf90_strerror(ierr)
    393393            stop "getvarup"
    394394         endif
    … …  
    396396
    397397         ierr = NF90_GET_VAR(nid,var3didin(11),T_adv_h)
    398          if(ierr/=NF_NOERR) then
    399             write(*,*) NF_STRERROR(ierr)
     398         if(ierr/=nf90_noerr) then
     399            write(*,*) nf90_strerror(ierr)
    400400            stop "getvarup"
    401401         endif
    … …  
    410410
    411411         ierr = NF90_GET_VAR(nid,var3didin(12),T_adv_v)
    412          if(ierr/=NF_NOERR) then
    413             write(*,*) NF_STRERROR(ierr)
     412         if(ierr/=nf90_noerr) then
     413            write(*,*) nf90_strerror(ierr)
    414414            stop "getvarup"
    415415         endif
    … …  
    423423
    424424         ierr = NF90_GET_VAR(nid,var3didin(13),q_adv_h)
    425          if(ierr/=NF_NOERR) then
    426             write(*,*) NF_STRERROR(ierr)
     425         if(ierr/=nf90_noerr) then
     426            write(*,*) nf90_strerror(ierr)
    427427            stop "getvarup"
    428428         endif
    … …  
    437437
    438438         ierr = NF90_GET_VAR(nid,var3didin(14),q_adv_v)
    439          if(ierr/=NF_NOERR) then
    440             write(*,*) NF_STRERROR(ierr)
     439         if(ierr/=nf90_noerr) then
     440            write(*,*) nf90_strerror(ierr)
    441441            stop "getvarup"
    442442         endif
    … …  
    451451
    452452         ierr = NF90_GET_VAR(nid,var3didin(15),s)
    453          if(ierr/=NF_NOERR) then
    454             write(*,*) NF_STRERROR(ierr)
     453         if(ierr/=nf90_noerr) then
     454            write(*,*) nf90_strerror(ierr)
    455455            stop "getvarup"
    456456         endif
    457457
    458458         ierr = NF90_GET_VAR(nid,var3didin(16),s_adv_h)
    459          if(ierr/=NF_NOERR) then
    460             write(*,*) NF_STRERROR(ierr)
     459         if(ierr/=nf90_noerr) then
     460            write(*,*) nf90_strerror(ierr)
    461461            stop "getvarup"
    462462         endif
    463463
    464464         ierr = NF90_GET_VAR(nid,var3didin(17),s_adv_v)
    465          if(ierr/=NF_NOERR) then
    466             write(*,*) NF_STRERROR(ierr)
     465         if(ierr/=nf90_noerr) then
     466            write(*,*) nf90_strerror(ierr)
    467467            stop "getvarup"
    468468         endif
    469469
    470470         ierr = NF90_GET_VAR(nid,var3didin(18),p_srf_aver)
    471          if(ierr/=NF_NOERR) then
    472             write(*,*) NF_STRERROR(ierr)
     471         if(ierr/=nf90_noerr) then
     472            write(*,*) nf90_strerror(ierr)
    473473            stop "getvarup"
    474474         endif
    475475
    476476         ierr = NF90_GET_VAR(nid,var3didin(19),p_srf_center)
    477          if(ierr/=NF_NOERR) then
    478             write(*,*) NF_STRERROR(ierr)
     477         if(ierr/=nf90_noerr) then
     478            write(*,*) nf90_strerror(ierr)
    479479            stop "getvarup"
    480480         endif
    481481
    482482         ierr = NF90_GET_VAR(nid,var3didin(20),T_srf)
    483          if(ierr/=NF_NOERR) then
    484             write(*,*) NF_STRERROR(ierr)
     483         if(ierr/=nf90_noerr) then
     484            write(*,*) nf90_strerror(ierr)
    485485            stop "getvarup"
    486486         endif
    … …  
    492492         subroutine catchaxis(nid,ttm,llm,time,lev,ierr)
    493493
    494          use lmdz_netcdf, ONLY: nf_open,nf_nowrite,nf_noerr,nf_strerror,nf_inq_varid,nf90_get_var,&
    495             nf_inq_dimid,nf_inq_dimlen
     494         use netcdf, ONLY: nf90_open,nf90_nowrite,nf90_noerr,nf90_strerror,nf90_inq_varid,nf90_get_var,&
     495            nf90_inq_dimid,nf90_inquire_dimension
    496496
    497497         implicit none
    … …  
    507507! Control & lecture on dimensions
    508508! ===============================
    509          ierr=NF_INQ_DIMID(nid,"time",timedimin)
    510          ierr=NF_INQ_VARID(nid,"time",timevar)
    511          if (ierr.NE.NF_NOERR) then
     509         ierr=nf90_inq_dimid(nid,"time",timedimin)
     510         ierr=nf90_inq_varid(nid,"time",timevar)
     511         if (ierr.NE.nf90_noerr) then
    512512            write(*,*) 'ERROR: Field <time> is missing'
    513513            stop "" 
    514514         endif
    515          ierr=NF_INQ_DIMLEN(nid,timedimin,timelen)
    516 
    517          ierr=NF_INQ_DIMID(nid,"lev",levdimin)
    518          ierr=NF_INQ_VARID(nid,"lev",levvar)
    519          if (ierr.NE.NF_NOERR) then
     515         ierr=nf90_inquire_dimension(nid,timedimin,len=timelen)
     516
     517         ierr=nf90_inq_dimid(nid,"lev",levdimin)
     518         ierr=nf90_inq_varid(nid,"lev",levvar)
     519         if (ierr.NE.nf90_noerr) then
    520520             write(*,*) 'ERROR: Field <lev> is lacking'
    521521             stop ""
    522522         endif
    523          ierr=NF_INQ_DIMLEN(nid,levdimin,levlen)
     523         ierr=nf90_inquire_dimension(nid,levdimin,len=levlen)
    524524
    525525         if((timelen/=ttm).or.(levlen/=llm)) then
    … …  
    21702170
    21712171
    2172       use lmdz_netcdf, ONLY: nf_open,nf_nowrite,nf_noerr,nf_strerror,nf_inq_varid,nf90_get_var,&
    2173             nf_inq_dimid,nf_inq_dimlen
     2172      use netcdf, ONLY: nf90_open,nf90_nowrite,nf90_noerr,nf90_strerror,nf90_inq_varid,nf90_get_var,&
     2173            nf90_inq_dimid,nf90_inquire_dimension
    21742174      implicit none
    21752175
    … …  
    21902190      integer var3didin(nbvar3d)
    21912191
    2192       ierr = NF_OPEN(fich_fire,NF_NOWRITE,nid)
    2193       if (ierr.NE.NF_NOERR) then
     2192      ierr = nf90_open(fich_fire,nf90_nowrite,nid)
     2193      if (ierr.NE.nf90_noerr) then
    21942194         write(*,*) 'ERROR: Pb opening forcings nc file '
    2195          write(*,*) NF_STRERROR(ierr)
     2195         write(*,*) nf90_strerror(ierr)
    21962196         stop ""
    21972197      endif
    21982198
    21992199
    2200        ierr=NF_INQ_VARID(nid,"zz",var3didin(1))
    2201          if(ierr/=NF_NOERR) then
    2202            write(*,*) NF_STRERROR(ierr)
     2200       ierr=nf90_inq_varid(nid,"zz",var3didin(1))
     2201         if(ierr/=nf90_noerr) then
     2202           write(*,*) nf90_strerror(ierr)
    22032203           stop 'lev'
    22042204         endif
    22052205
    22062206
    2207       ierr=NF_INQ_VARID(nid,"thetal",var3didin(2))
    2208          if(ierr/=NF_NOERR) then
    2209            write(*,*) NF_STRERROR(ierr)
     2207      ierr=nf90_inq_varid(nid,"thetal",var3didin(2))
     2208         if(ierr/=nf90_noerr) then
     2209           write(*,*) nf90_strerror(ierr)
    22102210           stop 'temp'
    22112211         endif
    22122212
    2213       ierr=NF_INQ_VARID(nid,"qt",var3didin(3))
    2214          if(ierr/=NF_NOERR) then
    2215            write(*,*) NF_STRERROR(ierr)
     2213      ierr=nf90_inq_varid(nid,"qt",var3didin(3))
     2214         if(ierr/=nf90_noerr) then
     2215           write(*,*) nf90_strerror(ierr)
    22162216           stop 'qv'
    22172217         endif
    22182218
    2219       ierr=NF_INQ_VARID(nid,"u",var3didin(4))
    2220          if(ierr/=NF_NOERR) then
    2221            write(*,*) NF_STRERROR(ierr)
     2219      ierr=nf90_inq_varid(nid,"u",var3didin(4))
     2220         if(ierr/=nf90_noerr) then
     2221           write(*,*) nf90_strerror(ierr)
    22222222           stop 'u'
    22232223         endif
    22242224
    2225       ierr=NF_INQ_VARID(nid,"v",var3didin(5))
    2226          if(ierr/=NF_NOERR) then
    2227            write(*,*) NF_STRERROR(ierr)
     2225      ierr=nf90_inq_varid(nid,"v",var3didin(5))
     2226         if(ierr/=nf90_noerr) then
     2227           write(*,*) nf90_strerror(ierr)
    22282228           stop 'v'
    22292229         endif
    22302230
    2231       ierr=NF_INQ_VARID(nid,"tke",var3didin(6))
    2232          if(ierr/=NF_NOERR) then
    2233            write(*,*) NF_STRERROR(ierr)
     2231      ierr=nf90_inq_varid(nid,"tke",var3didin(6))
     2232         if(ierr/=nf90_noerr) then
     2233           write(*,*) nf90_strerror(ierr)
    22342234           stop 'tke'
    22352235         endif
    22362236
    2237       ierr=NF_INQ_VARID(nid,"ugeo",var3didin(7))
    2238          if(ierr/=NF_NOERR) then
    2239            write(*,*) NF_STRERROR(ierr)
     2237      ierr=nf90_inq_varid(nid,"ugeo",var3didin(7))
     2238         if(ierr/=nf90_noerr) then
     2239           write(*,*) nf90_strerror(ierr)
    22402240           stop 'ug'
    22412241         endif
    22422242
    2243       ierr=NF_INQ_VARID(nid,"vgeo",var3didin(8))
    2244          if(ierr/=NF_NOERR) then
    2245            write(*,*) NF_STRERROR(ierr)
     2243      ierr=nf90_inq_varid(nid,"vgeo",var3didin(8))
     2244         if(ierr/=nf90_noerr) then
     2245           write(*,*) nf90_strerror(ierr)
    22462246           stop 'vg'
    22472247         endif
    22482248     
    2249       ierr=NF_INQ_VARID(nid,"wls",var3didin(9))
    2250          if(ierr/=NF_NOERR) then
    2251            write(*,*) NF_STRERROR(ierr)
     2249      ierr=nf90_inq_varid(nid,"wls",var3didin(9))
     2250         if(ierr/=nf90_noerr) then
     2251           write(*,*) nf90_strerror(ierr)
    22522252           stop 'wls'
    22532253         endif
    22542254
    2255       ierr=NF_INQ_VARID(nid,"dqtdx",var3didin(10))
    2256          if(ierr/=NF_NOERR) then
    2257            write(*,*) NF_STRERROR(ierr)
     2255      ierr=nf90_inq_varid(nid,"dqtdx",var3didin(10))
     2256         if(ierr/=nf90_noerr) then
     2257           write(*,*) nf90_strerror(ierr)
    22582258           stop 'dqtdx'
    22592259         endif
    22602260
    2261       ierr=NF_INQ_VARID(nid,"dqtdy",var3didin(11))
    2262          if(ierr/=NF_NOERR) then
    2263            write(*,*) NF_STRERROR(ierr)
     2261      ierr=nf90_inq_varid(nid,"dqtdy",var3didin(11))
     2262         if(ierr/=nf90_noerr) then
     2263           write(*,*) nf90_strerror(ierr)
    22642264           stop 'dqtdy'
    22652265      endif
    22662266
    2267       ierr=NF_INQ_VARID(nid,"dqtdt",var3didin(12))
    2268          if(ierr/=NF_NOERR) then
    2269            write(*,*) NF_STRERROR(ierr)
     2267      ierr=nf90_inq_varid(nid,"dqtdt",var3didin(12))
     2268         if(ierr/=nf90_noerr) then
     2269           write(*,*) nf90_strerror(ierr)
    22702270           stop 'dqtdt'
    22712271      endif
    22722272
    2273       ierr=NF_INQ_VARID(nid,"thl_rad",var3didin(13))
    2274          if(ierr/=NF_NOERR) then
    2275            write(*,*) NF_STRERROR(ierr)
     2273      ierr=nf90_inq_varid(nid,"thl_rad",var3didin(13))
     2274         if(ierr/=nf90_noerr) then
     2275           write(*,*) nf90_strerror(ierr)
    22762276           stop 'thl_rad'
    22772277      endif
    … …  
    22802280 
    22812281         ierr = NF90_GET_VAR(nid,var3didin(1),zz)
    2282          if(ierr/=NF_NOERR) then
    2283             write(*,*) NF_STRERROR(ierr)
     2282         if(ierr/=nf90_noerr) then
     2283            write(*,*) nf90_strerror(ierr)
    22842284            stop "getvarup"
    22852285         endif
    … …  
    22872287
    22882288         ierr = NF90_GET_VAR(nid,var3didin(2),thl)
    2289          if(ierr/=NF_NOERR) then
    2290             write(*,*) NF_STRERROR(ierr)
     2289         if(ierr/=nf90_noerr) then
     2290            write(*,*) nf90_strerror(ierr)
    22912291            stop "getvarup"
    22922292         endif
    … …  
    22942294
    22952295         ierr = NF90_GET_VAR(nid,var3didin(3),qt)
    2296          if(ierr/=NF_NOERR) then
    2297             write(*,*) NF_STRERROR(ierr)
     2296         if(ierr/=nf90_noerr) then
     2297            write(*,*) nf90_strerror(ierr)
    22982298            stop "getvarup"
    22992299         endif
    … …  
    23012301 
    23022302         ierr = NF90_GET_VAR(nid,var3didin(4),u)
    2303          if(ierr/=NF_NOERR) then
    2304             write(*,*) NF_STRERROR(ierr)
     2303         if(ierr/=nf90_noerr) then
     2304            write(*,*) nf90_strerror(ierr)
    23052305            stop "getvarup"
    23062306         endif
    … …  
    23082308
    23092309         ierr = NF90_GET_VAR(nid,var3didin(5),v)
    2310          if(ierr/=NF_NOERR) then
    2311             write(*,*) NF_STRERROR(ierr)
     2310         if(ierr/=nf90_noerr) then
     2311            write(*,*) nf90_strerror(ierr)
    23122312            stop "getvarup"
    23132313         endif
    … …  
    23152315
    23162316         ierr = NF90_GET_VAR(nid,var3didin(6),tke)
    2317          if(ierr/=NF_NOERR) then
    2318             write(*,*) NF_STRERROR(ierr)
     2317         if(ierr/=nf90_noerr) then
     2318            write(*,*) nf90_strerror(ierr)
    23192319            stop "getvarup"
    23202320         endif
    … …  
    23222322
    23232323         ierr = NF90_GET_VAR(nid,var3didin(7),ug)
    2324          if(ierr/=NF_NOERR) then
    2325             write(*,*) NF_STRERROR(ierr)
     2324         if(ierr/=nf90_noerr) then
     2325            write(*,*) nf90_strerror(ierr)
    23262326            stop "getvarup"
    23272327         endif
    … …  
    23292329
    23302330         ierr = NF90_GET_VAR(nid,var3didin(8),vg)
    2331          if(ierr/=NF_NOERR) then
    2332             write(*,*) NF_STRERROR(ierr)
     2331         if(ierr/=nf90_noerr) then
     2332            write(*,*) nf90_strerror(ierr)
    23332333            stop "getvarup"
    23342334         endif
    … …  
    23362336
    23372337         ierr = NF90_GET_VAR(nid,var3didin(9),wls)
    2338          if(ierr/=NF_NOERR) then
    2339             write(*,*) NF_STRERROR(ierr)
     2338         if(ierr/=nf90_noerr) then
     2339            write(*,*) nf90_strerror(ierr)
    23402340            stop "getvarup"
    23412341         endif
    … …  
    23432343
    23442344         ierr = NF90_GET_VAR(nid,var3didin(10),dqtdx)
    2345          if(ierr/=NF_NOERR) then
    2346             write(*,*) NF_STRERROR(ierr)
     2345         if(ierr/=nf90_noerr) then
     2346            write(*,*) nf90_strerror(ierr)
    23472347            stop "getvarup"
    23482348         endif
    … …  
    23502350
    23512351         ierr = NF90_GET_VAR(nid,var3didin(11),dqtdy)
    2352          if(ierr/=NF_NOERR) then
    2353             write(*,*) NF_STRERROR(ierr)
     2352         if(ierr/=nf90_noerr) then
     2353            write(*,*) nf90_strerror(ierr)
    23542354            stop "getvarup"
    23552355         endif
    … …  
    23572357
    23582358         ierr = NF90_GET_VAR(nid,var3didin(12),dqtdt)
    2359          if(ierr/=NF_NOERR) then
    2360             write(*,*) NF_STRERROR(ierr)
     2359         if(ierr/=nf90_noerr) then
     2360            write(*,*) nf90_strerror(ierr)
    23612361            stop "getvarup"
    23622362         endif
    … …  
    23642364
    23652365         ierr = NF90_GET_VAR(nid,var3didin(13),thl_rad)
    2366          if(ierr/=NF_NOERR) then
    2367             write(*,*) NF_STRERROR(ierr)
     2366         if(ierr/=nf90_noerr) then
     2367            write(*,*) nf90_strerror(ierr)
    23682368            stop "getvarup"
    23692369         endif
    … …  
    23802380!program reading initial profils and forcings of the Dice case study
    23812381
    2382       use lmdz_netcdf, ONLY: nf_open,nf_nowrite,nf_noerr,nf_strerror,nf_inq_varid,nf90_get_var,&
    2383             nf_inq_dimid,nf_inq_dimlen
     2382      use netcdf, ONLY: nf90_open,nf90_nowrite,nf90_noerr,nf90_strerror,nf90_inq_varid,nf90_get_var,&
     2383            nf90_inq_dimid,nf90_inquire_dimension
    23842384
    23852385      implicit none
    … …  
    24072407
    24082408      pzero=100000.
    2409       ierr = NF_OPEN(fich_dice,NF_NOWRITE,nid)
    2410       if (ierr.NE.NF_NOERR) then
     2409      ierr = nf90_open(fich_dice,nf90_nowrite,nid)
     2410      if (ierr.NE.nf90_noerr) then
    24112411         write(*,*) 'ERROR: Pb opening forcings nc file '
    2412          write(*,*) NF_STRERROR(ierr)
     2412         write(*,*) nf90_strerror(ierr)
    24132413         stop ""
    24142414      endif
    24152415
    24162416
    2417        ierr=NF_INQ_VARID(nid,"height",var3didin(1))
    2418          if(ierr/=NF_NOERR) then
    2419            write(*,*) NF_STRERROR(ierr)
     2417       ierr=nf90_inq_varid(nid,"height",var3didin(1))
     2418         if(ierr/=nf90_noerr) then
     2419           write(*,*) nf90_strerror(ierr)
    24202420           stop 'height'
    24212421         endif
    24222422
    2423        ierr=NF_INQ_VARID(nid,"pf",var3didin(11))
    2424          if(ierr/=NF_NOERR) then
    2425            write(*,*) NF_STRERROR(ierr)
     2423       ierr=nf90_inq_varid(nid,"pf",var3didin(11))
     2424         if(ierr/=nf90_noerr) then
     2425           write(*,*) nf90_strerror(ierr)
    24262426           stop 'pf'
    24272427         endif
    24282428
    2429       ierr=NF_INQ_VARID(nid,"theta",var3didin(12))
    2430          if(ierr/=NF_NOERR) then
    2431            write(*,*) NF_STRERROR(ierr)
     2429      ierr=nf90_inq_varid(nid,"theta",var3didin(12))
     2430         if(ierr/=nf90_noerr) then
     2431           write(*,*) nf90_strerror(ierr)
    24322432           stop 'theta'
    24332433         endif
    24342434
    2435       ierr=NF_INQ_VARID(nid,"qv",var3didin(13))
    2436          if(ierr/=NF_NOERR) then
    2437            write(*,*) NF_STRERROR(ierr)
     2435      ierr=nf90_inq_varid(nid,"qv",var3didin(13))
     2436         if(ierr/=nf90_noerr) then
     2437           write(*,*) nf90_strerror(ierr)
    24382438           stop 'qv'
    24392439         endif
    24402440
    2441       ierr=NF_INQ_VARID(nid,"u",var3didin(14))
    2442          if(ierr/=NF_NOERR) then
    2443            write(*,*) NF_STRERROR(ierr)
     2441      ierr=nf90_inq_varid(nid,"u",var3didin(14))
     2442         if(ierr/=nf90_noerr) then
     2443           write(*,*) nf90_strerror(ierr)
    24442444           stop 'u'
    24452445         endif
    24462446
    2447       ierr=NF_INQ_VARID(nid,"v",var3didin(15))
    2448          if(ierr/=NF_NOERR) then
    2449            write(*,*) NF_STRERROR(ierr)
     2447      ierr=nf90_inq_varid(nid,"v",var3didin(15))
     2448         if(ierr/=nf90_noerr) then
     2449           write(*,*) nf90_strerror(ierr)
    24502450           stop 'v'
    24512451         endif
    24522452
    2453       ierr=NF_INQ_VARID(nid,"o3mmr",var3didin(16))
    2454          if(ierr/=NF_NOERR) then
    2455            write(*,*) NF_STRERROR(ierr)
     2453      ierr=nf90_inq_varid(nid,"o3mmr",var3didin(16))
     2454         if(ierr/=nf90_noerr) then
     2455           write(*,*) nf90_strerror(ierr)
    24562456           stop 'o3'
    24572457         endif
    24582458
    2459       ierr=NF_INQ_VARID(nid,"shf",var3didin(2))
    2460          if(ierr/=NF_NOERR) then
    2461            write(*,*) NF_STRERROR(ierr)
     2459      ierr=nf90_inq_varid(nid,"shf",var3didin(2))
     2460         if(ierr/=nf90_noerr) then
     2461           write(*,*) nf90_strerror(ierr)
    24622462           stop 'shf'
    24632463         endif
    24642464
    2465       ierr=NF_INQ_VARID(nid,"lhf",var3didin(3))
    2466          if(ierr/=NF_NOERR) then
    2467            write(*,*) NF_STRERROR(ierr)
     2465      ierr=nf90_inq_varid(nid,"lhf",var3didin(3))
     2466         if(ierr/=nf90_noerr) then
     2467           write(*,*) nf90_strerror(ierr)
    24682468           stop 'lhf'
    24692469         endif
    24702470     
    2471       ierr=NF_INQ_VARID(nid,"lwup",var3didin(4))
    2472          if(ierr/=NF_NOERR) then
    2473            write(*,*) NF_STRERROR(ierr)
     2471      ierr=nf90_inq_varid(nid,"lwup",var3didin(4))
     2472         if(ierr/=nf90_noerr) then
     2473           write(*,*) nf90_strerror(ierr)
    24742474           stop 'lwup'
    24752475         endif
    24762476
    2477       ierr=NF_INQ_VARID(nid,"swup",var3didin(5))
    2478          if(ierr/=NF_NOERR) then
    2479            write(*,*) NF_STRERROR(ierr)
     2477      ierr=nf90_inq_varid(nid,"swup",var3didin(5))
     2478         if(ierr/=nf90_noerr) then
     2479           write(*,*) nf90_strerror(ierr)
    24802480           stop 'dqtdx'
    24812481         endif
    24822482
    2483       ierr=NF_INQ_VARID(nid,"Tg",var3didin(6))
    2484          if(ierr/=NF_NOERR) then
    2485            write(*,*) NF_STRERROR(ierr)
     2483      ierr=nf90_inq_varid(nid,"Tg",var3didin(6))
     2484         if(ierr/=nf90_noerr) then
     2485           write(*,*) nf90_strerror(ierr)
    24862486           stop 'Tg'
    24872487      endif
    24882488
    2489       ierr=NF_INQ_VARID(nid,"ustar",var3didin(7))
    2490          if(ierr/=NF_NOERR) then
    2491            write(*,*) NF_STRERROR(ierr)
     2489      ierr=nf90_inq_varid(nid,"ustar",var3didin(7))
     2490         if(ierr/=nf90_noerr) then
     2491           write(*,*) nf90_strerror(ierr)
    24922492           stop 'ustar'
    24932493      endif
    24942494
    2495       ierr=NF_INQ_VARID(nid,"psurf",var3didin(8))
    2496          if(ierr/=NF_NOERR) then
    2497            write(*,*) NF_STRERROR(ierr)
     2495      ierr=nf90_inq_varid(nid,"psurf",var3didin(8))
     2496         if(ierr/=nf90_noerr) then
     2497           write(*,*) nf90_strerror(ierr)
    24982498           stop 'psurf'
    24992499      endif
    25002500
    2501       ierr=NF_INQ_VARID(nid,"Ug",var3didin(9))
    2502          if(ierr/=NF_NOERR) then
    2503            write(*,*) NF_STRERROR(ierr)
     2501      ierr=nf90_inq_varid(nid,"Ug",var3didin(9))
     2502         if(ierr/=nf90_noerr) then
     2503           write(*,*) nf90_strerror(ierr)
    25042504           stop 'Ug'
    25052505      endif
    25062506
    2507       ierr=NF_INQ_VARID(nid,"Vg",var3didin(10))
    2508          if(ierr/=NF_NOERR) then
    2509            write(*,*) NF_STRERROR(ierr)
     2507      ierr=nf90_inq_varid(nid,"Vg",var3didin(10))
     2508         if(ierr/=nf90_noerr) then
     2509           write(*,*) nf90_strerror(ierr)
    25102510           stop 'Vg'
    25112511      endif
    25122512
    2513       ierr=NF_INQ_VARID(nid,"hadvT",var3didin(17))
    2514          if(ierr/=NF_NOERR) then
    2515            write(*,*) NF_STRERROR(ierr)
     2513      ierr=nf90_inq_varid(nid,"hadvT",var3didin(17))
     2514         if(ierr/=nf90_noerr) then
     2515           write(*,*) nf90_strerror(ierr)
    25162516           stop 'hadvT'
    25172517      endif
    25182518
    2519       ierr=NF_INQ_VARID(nid,"hadvq",var3didin(18))
    2520          if(ierr/=NF_NOERR) then
    2521            write(*,*) NF_STRERROR(ierr)
     2519      ierr=nf90_inq_varid(nid,"hadvq",var3didin(18))
     2520         if(ierr/=nf90_noerr) then
     2521           write(*,*) nf90_strerror(ierr)
    25222522           stop 'hadvq'
    25232523      endif
    25242524
    2525       ierr=NF_INQ_VARID(nid,"hadvu",var3didin(19))
    2526          if(ierr/=NF_NOERR) then
    2527            write(*,*) NF_STRERROR(ierr)
     2525      ierr=nf90_inq_varid(nid,"hadvu",var3didin(19))
     2526         if(ierr/=nf90_noerr) then
     2527           write(*,*) nf90_strerror(ierr)
    25282528           stop 'hadvu'
    25292529      endif
    25302530
    2531       ierr=NF_INQ_VARID(nid,"hadvv",var3didin(20))
    2532          if(ierr/=NF_NOERR) then
    2533            write(*,*) NF_STRERROR(ierr)
     2531      ierr=nf90_inq_varid(nid,"hadvv",var3didin(20))
     2532         if(ierr/=nf90_noerr) then
     2533           write(*,*) nf90_strerror(ierr)
    25342534           stop 'hadvv'
    25352535      endif
    25362536
    2537       ierr=NF_INQ_VARID(nid,"w",var3didin(21))
    2538          if(ierr/=NF_NOERR) then
    2539            write(*,*) NF_STRERROR(ierr)
     2537      ierr=nf90_inq_varid(nid,"w",var3didin(21))
     2538         if(ierr/=nf90_noerr) then
     2539           write(*,*) nf90_strerror(ierr)
    25402540           stop 'w'
    25412541      endif
    25422542
    2543       ierr=NF_INQ_VARID(nid,"omega",var3didin(22))
    2544          if(ierr/=NF_NOERR) then
    2545            write(*,*) NF_STRERROR(ierr)
     2543      ierr=nf90_inq_varid(nid,"omega",var3didin(22))
     2544         if(ierr/=nf90_noerr) then
     2545           write(*,*) nf90_strerror(ierr)
    25462546           stop 'omega'
    25472547      endif
    … …  
    25502550 
    25512551         ierr = NF90_GET_VAR(nid,var3didin(1),zz)
    2552          if(ierr/=NF_NOERR) then
    2553             write(*,*) NF_STRERROR(ierr)
     2552         if(ierr/=nf90_noerr) then
     2553            write(*,*) nf90_strerror(ierr)
    25542554            stop "getvarup"
    25552555         endif
    … …  
    25572557 
    25582558         ierr = NF90_GET_VAR(nid,var3didin(11),pres)
    2559          if(ierr/=NF_NOERR) then
    2560             write(*,*) NF_STRERROR(ierr)
     2559         if(ierr/=nf90_noerr) then
     2560            write(*,*) nf90_strerror(ierr)
    25612561            stop "getvarup"
    25622562         endif
    … …  
    25642564
    25652565         ierr = NF90_GET_VAR(nid,var3didin(12),th)
    2566          if(ierr/=NF_NOERR) then
    2567             write(*,*) NF_STRERROR(ierr)
     2566         if(ierr/=nf90_noerr) then
     2567            write(*,*) nf90_strerror(ierr)
    25682568            stop "getvarup"
    25692569         endif
    … …  
    25742574
    25752575         ierr = NF90_GET_VAR(nid,var3didin(13),qv)
    2576          if(ierr/=NF_NOERR) then
    2577             write(*,*) NF_STRERROR(ierr)
     2576         if(ierr/=nf90_noerr) then
     2577            write(*,*) nf90_strerror(ierr)
    25782578            stop "getvarup"
    25792579         endif
    … …  
    25812581 
    25822582         ierr = NF90_GET_VAR(nid,var3didin(14),u)
    2583          if(ierr/=NF_NOERR) then
    2584             write(*,*) NF_STRERROR(ierr)
     2583         if(ierr/=nf90_noerr) then
     2584            write(*,*) nf90_strerror(ierr)
    25852585            stop "getvarup"
    25862586         endif
    … …  
    25882588
    25892589         ierr = NF90_GET_VAR(nid,var3didin(15),v)
    2590          if(ierr/=NF_NOERR) then
    2591             write(*,*) NF_STRERROR(ierr)
     2590         if(ierr/=nf90_noerr) then
     2591            write(*,*) nf90_strerror(ierr)
    25922592            stop "getvarup"
    25932593         endif
    … …  
    25952595
    25962596         ierr = NF90_GET_VAR(nid,var3didin(16),o3)
    2597          if(ierr/=NF_NOERR) then
    2598             write(*,*) NF_STRERROR(ierr)
     2597         if(ierr/=nf90_noerr) then
     2598            write(*,*) nf90_strerror(ierr)
    25992599            stop "getvarup"
    26002600         endif
    … …  
    26022602
    26032603         ierr = NF90_GET_VAR(nid,var3didin(2),shf)
    2604          if(ierr/=NF_NOERR) then
    2605             write(*,*) NF_STRERROR(ierr)
     2604         if(ierr/=nf90_noerr) then
     2605            write(*,*) nf90_strerror(ierr)
    26062606            stop "getvarup"
    26072607         endif
    … …  
    26092609
    26102610         ierr = NF90_GET_VAR(nid,var3didin(3),lhf)
    2611          if(ierr/=NF_NOERR) then
    2612             write(*,*) NF_STRERROR(ierr)
     2611         if(ierr/=nf90_noerr) then
     2612            write(*,*) nf90_strerror(ierr)
    26132613            stop "getvarup"
    26142614         endif
    … …  
    26162616
    26172617         ierr = NF90_GET_VAR(nid,var3didin(4),lwup)
    2618          if(ierr/=NF_NOERR) then
    2619             write(*,*) NF_STRERROR(ierr)
     2618         if(ierr/=nf90_noerr) then
     2619            write(*,*) nf90_strerror(ierr)
    26202620            stop "getvarup"
    26212621         endif
    … …  
    26232623
    26242624         ierr = NF90_GET_VAR(nid,var3didin(5),swup)
    2625          if(ierr/=NF_NOERR) then
    2626             write(*,*) NF_STRERROR(ierr)
     2625         if(ierr/=nf90_noerr) then
     2626            write(*,*) nf90_strerror(ierr)
    26272627            stop "getvarup"
    26282628         endif
    … …  
    26302630
    26312631         ierr = NF90_GET_VAR(nid,var3didin(6),tg)
    2632          if(ierr/=NF_NOERR) then
    2633             write(*,*) NF_STRERROR(ierr)
     2632         if(ierr/=nf90_noerr) then
     2633            write(*,*) nf90_strerror(ierr)
    26342634            stop "getvarup"
    26352635         endif
    … …  
    26372637
    26382638         ierr = NF90_GET_VAR(nid,var3didin(7),ustar)
    2639          if(ierr/=NF_NOERR) then
    2640             write(*,*) NF_STRERROR(ierr)
     2639         if(ierr/=nf90_noerr) then
     2640            write(*,*) nf90_strerror(ierr)
    26412641            stop "getvarup"
    26422642         endif
    … …  
    26442644
    26452645         ierr = NF90_GET_VAR(nid,var3didin(8),psurf)
    2646          if(ierr/=NF_NOERR) then
    2647             write(*,*) NF_STRERROR(ierr)
     2646         if(ierr/=nf90_noerr) then
     2647            write(*,*) nf90_strerror(ierr)
    26482648            stop "getvarup"
    26492649         endif
    … …  
    26512651
    26522652         ierr = NF90_GET_VAR(nid,var3didin(9),ug)
    2653          if(ierr/=NF_NOERR) then
    2654             write(*,*) NF_STRERROR(ierr)
     2653         if(ierr/=nf90_noerr) then
     2654            write(*,*) nf90_strerror(ierr)
    26552655            stop "getvarup"
    26562656         endif
    … …  
    26582658
    26592659         ierr = NF90_GET_VAR(nid,var3didin(10),vg)
    2660          if(ierr/=NF_NOERR) then
    2661             write(*,*) NF_STRERROR(ierr)
     2660         if(ierr/=nf90_noerr) then
     2661            write(*,*) nf90_strerror(ierr)
    26622662            stop "getvarup"
    26632663         endif
    … …  
    26652665
    26662666         ierr = NF90_GET_VAR(nid,var3didin(17),hadvt)
    2667          if(ierr/=NF_NOERR) then
    2668             write(*,*) NF_STRERROR(ierr)
     2667         if(ierr/=nf90_noerr) then
     2668            write(*,*) nf90_strerror(ierr)
    26692669            stop "getvarup"
    26702670         endif
    … …  
    26722672
    26732673         ierr = NF90_GET_VAR(nid,var3didin(18),hadvq)
    2674          if(ierr/=NF_NOERR) then
    2675             write(*,*) NF_STRERROR(ierr)
     2674         if(ierr/=nf90_noerr) then
     2675            write(*,*) nf90_strerror(ierr)
    26762676            stop "getvarup"
    26772677         endif
    … …  
    26792679
    26802680         ierr = NF90_GET_VAR(nid,var3didin(19),hadvu)
    2681          if(ierr/=NF_NOERR) then
    2682             write(*,*) NF_STRERROR(ierr)
     2681         if(ierr/=nf90_noerr) then
     2682            write(*,*) nf90_strerror(ierr)
    26832683            stop "getvarup"
    26842684         endif
    … …  
    26862686
    26872687         ierr = NF90_GET_VAR(nid,var3didin(20),hadvv)
    2688          if(ierr/=NF_NOERR) then
    2689             write(*,*) NF_STRERROR(ierr)
     2688         if(ierr/=nf90_noerr) then
     2689            write(*,*) nf90_strerror(ierr)
    26902690            stop "getvarup"
    26912691         endif
    … …  
    26932693
    26942694         ierr = NF90_GET_VAR(nid,var3didin(21),w)
    2695          if(ierr/=NF_NOERR) then
    2696             write(*,*) NF_STRERROR(ierr)
     2695         if(ierr/=nf90_noerr) then
     2696            write(*,*) nf90_strerror(ierr)
    26972697            stop "getvarup"
    26982698         endif
    … …  
    27002700
    27012701         ierr = NF90_GET_VAR(nid,var3didin(22),omega)
    2702          if(ierr/=NF_NOERR) then
    2703             write(*,*) NF_STRERROR(ierr)
     2702         if(ierr/=nf90_noerr) then
     2703            write(*,*) nf90_strerror(ierr)
    27042704            stop "getvarup"
    27052705         endif
    … …  
    27142714!program reading initial profils and forcings of the Gabls4 case study
    27152715
    2716       use lmdz_netcdf, ONLY: nf_open,nf_nowrite,nf_noerr,nf_strerror,nf_inq_varid,nf90_get_var,&
    2717             nf_inq_dimid,nf_inq_dimlen
     2716      use netcdf, ONLY: nf90_open,nf90_nowrite,nf90_noerr,nf90_strerror,nf90_inq_varid,nf90_get_var,&
     2717            nf90_inq_dimid,nf90_inquire_dimension
    27182718
    27192719      implicit none
    … …  
    27432743      integer var3didin(nbvar3d)
    27442744
    2745       ierr = NF_OPEN(fich_gabls4,NF_NOWRITE,nid)
    2746       if (ierr.NE.NF_NOERR) then
     2745      ierr = nf90_open(fich_gabls4,nf90_nowrite,nid)
     2746      if (ierr.NE.nf90_noerr) then
    27472747         write(*,*) 'ERROR: Pb opening forcings nc file '
    2748          write(*,*) NF_STRERROR(ierr)
     2748         write(*,*) nf90_strerror(ierr)
    27492749         stop ""
    27502750      endif
    27512751
    27522752
    2753        ierr=NF_INQ_VARID(nid,"height",var3didin(1))
    2754          if(ierr/=NF_NOERR) then
    2755            write(*,*) NF_STRERROR(ierr)
     2753       ierr=nf90_inq_varid(nid,"height",var3didin(1))
     2754         if(ierr/=nf90_noerr) then
     2755           write(*,*) nf90_strerror(ierr)
    27562756           stop 'height'
    27572757         endif
    27582758
    2759       ierr=NF_INQ_VARID(nid,"depth_sn",var3didin(2))
    2760          if(ierr/=NF_NOERR) then
    2761            write(*,*) NF_STRERROR(ierr)
     2759      ierr=nf90_inq_varid(nid,"depth_sn",var3didin(2))
     2760         if(ierr/=nf90_noerr) then
     2761           write(*,*) nf90_strerror(ierr)
    27622762           stop 'depth_sn'
    27632763      endif
    27642764
    2765       ierr=NF_INQ_VARID(nid,"Ug",var3didin(3))
    2766          if(ierr/=NF_NOERR) then
    2767            write(*,*) NF_STRERROR(ierr)
     2765      ierr=nf90_inq_varid(nid,"Ug",var3didin(3))
     2766         if(ierr/=nf90_noerr) then
     2767           write(*,*) nf90_strerror(ierr)
    27682768           stop 'Ug'
    27692769      endif
    27702770
    2771       ierr=NF_INQ_VARID(nid,"Vg",var3didin(4))
    2772          if(ierr/=NF_NOERR) then
    2773            write(*,*) NF_STRERROR(ierr)
     2771      ierr=nf90_inq_varid(nid,"Vg",var3didin(4))
     2772         if(ierr/=nf90_noerr) then
     2773           write(*,*) nf90_strerror(ierr)
    27742774           stop 'Vg'
    27752775      endif
    2776        ierr=NF_INQ_VARID(nid,"pf",var3didin(5))
    2777          if(ierr/=NF_NOERR) then
    2778            write(*,*) NF_STRERROR(ierr)
     2776       ierr=nf90_inq_varid(nid,"pf",var3didin(5))
     2777         if(ierr/=nf90_noerr) then
     2778           write(*,*) nf90_strerror(ierr)
    27792779           stop 'pf'
    27802780         endif
    27812781
    2782       ierr=NF_INQ_VARID(nid,"theta",var3didin(6))
    2783          if(ierr/=NF_NOERR) then
    2784            write(*,*) NF_STRERROR(ierr)
     2782      ierr=nf90_inq_varid(nid,"theta",var3didin(6))
     2783         if(ierr/=nf90_noerr) then
     2784           write(*,*) nf90_strerror(ierr)
    27852785           stop 'theta'
    27862786         endif
    27872787
    2788       ierr=NF_INQ_VARID(nid,"tempe",var3didin(7))
    2789          if(ierr/=NF_NOERR) then
    2790            write(*,*) NF_STRERROR(ierr)
     2788      ierr=nf90_inq_varid(nid,"tempe",var3didin(7))
     2789         if(ierr/=nf90_noerr) then
     2790           write(*,*) nf90_strerror(ierr)
    27912791           stop 'tempe'
    27922792         endif
    27932793
    2794       ierr=NF_INQ_VARID(nid,"qv",var3didin(8))
    2795          if(ierr/=NF_NOERR) then
    2796            write(*,*) NF_STRERROR(ierr)
     2794      ierr=nf90_inq_varid(nid,"qv",var3didin(8))
     2795         if(ierr/=nf90_noerr) then
     2796           write(*,*) nf90_strerror(ierr)
    27972797           stop 'qv'
    27982798         endif
    27992799
    2800       ierr=NF_INQ_VARID(nid,"u",var3didin(9))
    2801          if(ierr/=NF_NOERR) then
    2802            write(*,*) NF_STRERROR(ierr)
     2800      ierr=nf90_inq_varid(nid,"u",var3didin(9))
     2801         if(ierr/=nf90_noerr) then
     2802           write(*,*) nf90_strerror(ierr)
    28032803           stop 'u'
    28042804         endif
    28052805
    2806       ierr=NF_INQ_VARID(nid,"v",var3didin(10))
    2807          if(ierr/=NF_NOERR) then
    2808            write(*,*) NF_STRERROR(ierr)
     2806      ierr=nf90_inq_varid(nid,"v",var3didin(10))
     2807         if(ierr/=nf90_noerr) then
     2808           write(*,*) nf90_strerror(ierr)
    28092809           stop 'v'
    28102810         endif
    28112811
    2812       ierr=NF_INQ_VARID(nid,"hadvT",var3didin(11))
    2813          if(ierr/=NF_NOERR) then
    2814            write(*,*) NF_STRERROR(ierr)
     2812      ierr=nf90_inq_varid(nid,"hadvT",var3didin(11))
     2813         if(ierr/=nf90_noerr) then
     2814           write(*,*) nf90_strerror(ierr)
    28152815           stop 'hadvt'
    28162816         endif
    28172817
    2818       ierr=NF_INQ_VARID(nid,"hadvQ",var3didin(12))
    2819          if(ierr/=NF_NOERR) then
    2820            write(*,*) NF_STRERROR(ierr)
     2818      ierr=nf90_inq_varid(nid,"hadvQ",var3didin(12))
     2819         if(ierr/=nf90_noerr) then
     2820           write(*,*) nf90_strerror(ierr)
    28212821           stop 'hadvq'
    28222822      endif
    28232823
    2824       ierr=NF_INQ_VARID(nid,"Tsnow",var3didin(14))
    2825          if(ierr/=NF_NOERR) then
    2826            write(*,*) NF_STRERROR(ierr)
     2824      ierr=nf90_inq_varid(nid,"Tsnow",var3didin(14))
     2825         if(ierr/=nf90_noerr) then
     2826           write(*,*) nf90_strerror(ierr)
    28272827           stop 'tsnow'
    28282828      endif
    28292829
    2830       ierr=NF_INQ_VARID(nid,"snow_density",var3didin(15))
    2831          if(ierr/=NF_NOERR) then
    2832            write(*,*) NF_STRERROR(ierr)
     2830      ierr=nf90_inq_varid(nid,"snow_density",var3didin(15))
     2831         if(ierr/=nf90_noerr) then
     2832           write(*,*) nf90_strerror(ierr)
    28332833           stop 'snow_density'
    28342834      endif
    28352835
    2836       ierr=NF_INQ_VARID(nid,"Tg",var3didin(16))
    2837          if(ierr/=NF_NOERR) then
    2838            write(*,*) NF_STRERROR(ierr)
     2836      ierr=nf90_inq_varid(nid,"Tg",var3didin(16))
     2837         if(ierr/=nf90_noerr) then
     2838           write(*,*) nf90_strerror(ierr)
    28392839           stop 'Tg'
    28402840      endif
    … …  
    28452845 
    28462846         ierr = NF90_GET_VAR(nid,var3didin(1),zz_i)
    2847          if(ierr/=NF_NOERR) then
    2848             write(*,*) NF_STRERROR(ierr)
     2847         if(ierr/=nf90_noerr) then
     2848            write(*,*) nf90_strerror(ierr)
    28492849            stop "getvarup"
    28502850         endif
    28512851 
    28522852         ierr = NF90_GET_VAR(nid,var3didin(2),depth_sn)
    2853          if(ierr/=NF_NOERR) then
    2854             write(*,*) NF_STRERROR(ierr)
     2853         if(ierr/=nf90_noerr) then
     2854            write(*,*) nf90_strerror(ierr)
    28552855            stop "getvarup"
    28562856         endif
    28572857 
    28582858         ierr = NF90_GET_VAR(nid,var3didin(3),ug_i)
    2859          if(ierr/=NF_NOERR) then
    2860             write(*,*) NF_STRERROR(ierr)
     2859         if(ierr/=nf90_noerr) then
     2860            write(*,*) nf90_strerror(ierr)
    28612861            stop "getvarup"
    28622862         endif
    28632863 
    28642864         ierr = NF90_GET_VAR(nid,var3didin(4),vg_i)
    2865          if(ierr/=NF_NOERR) then
    2866             write(*,*) NF_STRERROR(ierr)
     2865         if(ierr/=nf90_noerr) then
     2866            write(*,*) nf90_strerror(ierr)
    28672867            stop "getvarup"
    28682868         endif
    28692869 
    28702870         ierr = NF90_GET_VAR(nid,var3didin(5),pf_i)
    2871          if(ierr/=NF_NOERR) then
    2872             write(*,*) NF_STRERROR(ierr)
     2871         if(ierr/=nf90_noerr) then
     2872            write(*,*) nf90_strerror(ierr)
    28732873            stop "getvarup"
    28742874         endif
    28752875
    28762876         ierr = NF90_GET_VAR(nid,var3didin(6),th_i)
    2877          if(ierr/=NF_NOERR) then
    2878             write(*,*) NF_STRERROR(ierr)
     2877         if(ierr/=nf90_noerr) then
     2878            write(*,*) nf90_strerror(ierr)
    28792879            stop "getvarup"
    28802880         endif
    28812881
    28822882         ierr = NF90_GET_VAR(nid,var3didin(7),t_i)
    2883          if(ierr/=NF_NOERR) then
    2884             write(*,*) NF_STRERROR(ierr)
     2883         if(ierr/=nf90_noerr) then
     2884            write(*,*) nf90_strerror(ierr)
    28852885            stop "getvarup"
    28862886         endif
    28872887
    28882888         ierr = NF90_GET_VAR(nid,var3didin(8),qv_i)
    2889          if(ierr/=NF_NOERR) then
    2890             write(*,*) NF_STRERROR(ierr)
     2889         if(ierr/=nf90_noerr) then
     2890            write(*,*) nf90_strerror(ierr)
    28912891            stop "getvarup"
    28922892         endif
    28932893 
    28942894         ierr = NF90_GET_VAR(nid,var3didin(9),u_i)
    2895          if(ierr/=NF_NOERR) then
    2896             write(*,*) NF_STRERROR(ierr)
     2895         if(ierr/=nf90_noerr) then
     2896            write(*,*) nf90_strerror(ierr)
    28972897            stop "getvarup"
    28982898         endif
    28992899 
    29002900         ierr = NF90_GET_VAR(nid,var3didin(10),v_i)
    2901          if(ierr/=NF_NOERR) then
    2902             write(*,*) NF_STRERROR(ierr)
     2901         if(ierr/=nf90_noerr) then
     2902            write(*,*) nf90_strerror(ierr)
    29032903            stop "getvarup"
    29042904         endif
    29052905 
    29062906         ierr = NF90_GET_VAR(nid,var3didin(11),hadvt_i)
    2907          if(ierr/=NF_NOERR) then
    2908             write(*,*) NF_STRERROR(ierr)
     2907         if(ierr/=nf90_noerr) then
     2908            write(*,*) nf90_strerror(ierr)
    29092909            stop "getvarup"
    29102910         endif
    29112911 
    29122912         ierr = NF90_GET_VAR(nid,var3didin(12),hadvq_i)
    2913          if(ierr/=NF_NOERR) then
    2914             write(*,*) NF_STRERROR(ierr)
     2913         if(ierr/=nf90_noerr) then
     2914            write(*,*) nf90_strerror(ierr)
    29152915            stop "getvarup"
    29162916         endif
    29172917 
    29182918         ierr = NF90_GET_VAR(nid,var3didin(14),tsnow)
    2919          if(ierr/=NF_NOERR) then
    2920             write(*,*) NF_STRERROR(ierr)
     2919         if(ierr/=nf90_noerr) then
     2920            write(*,*) nf90_strerror(ierr)
    29212921            stop "getvarup"
    29222922         endif
    29232923 
    29242924         ierr = NF90_GET_VAR(nid,var3didin(15),snow_dens)
    2925          if(ierr/=NF_NOERR) then
    2926             write(*,*) NF_STRERROR(ierr)
     2925         if(ierr/=nf90_noerr) then
     2926            write(*,*) nf90_strerror(ierr)
    29272927            stop "getvarup"
    29282928         endif
    29292929
    29302930         ierr = NF90_GET_VAR(nid,var3didin(16),tg)
    2931          if(ierr/=NF_NOERR) then
    2932             write(*,*) NF_STRERROR(ierr)
     2931         if(ierr/=nf90_noerr) then
     2932            write(*,*) nf90_strerror(ierr)
    29332933            stop "getvarup"
    29342934         endif
Note: See TracChangeset for help on using the changeset viewer.