Ignore:
Timestamp:
Jul 18, 2024, 4:57:05 PM (5 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/dyn1d
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/dyn1d/mod_1D_amma_read.F90

    r4593 r5073  
    6767      ierr = NF_OPEN(fich_amma,NF_NOWRITE,nid)
    6868      print*,'fich_amma,NF_NOWRITE,nid ',fich_amma,NF_NOWRITE,nid
    69       if (ierr.NE.NF_NOERR) then
     69      if (ierr/=NF_NOERR) then
    7070         write(*,*) 'ERROR: GROS Pb opening forcings nc file '
    7171         write(*,*) NF_STRERROR(ierr)
    … …  
    7474!.......................................................................
    7575      ierr=NF_INQ_DIMID(nid,'lev',rid)
    76       IF (ierr.NE.NF_NOERR) THEN
     76      IF (ierr/=NF_NOERR) THEN
    7777         print*, 'Oh probleme lecture dimension zz'
    7878      ENDIF
    … …  
    8383      print*,'nid,rid',nid,rid
    8484      nt_amma=0
    85       IF (ierr.NE.NF_NOERR) THEN
     85      IF (ierr/=NF_NOERR) THEN
    8686        stop 'probleme lecture dimension sens'
    8787      ENDIF
    … …  
    268268!      call catchaxis(nid,ntime,nlevel,time,z,ierr)
    269269 
    270 #ifdef NC_DOUBLE
    271          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(1),zz)
    272 #else
    273          ierr = NF_GET_VAR_REAL(nid,var3didin(1),zz)
    274 #endif
     270         ierr = nf90_get_var(nid,var3didin(1),zz)
    275271         if(ierr/=NF_NOERR) then
    276272            write(*,*) NF_STRERROR(ierr)
    … …  
    279275!          write(*,*)'lecture z ok',zz
    280276
    281 #ifdef NC_DOUBLE
    282          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(2),temp)
    283 #else
    284          ierr = NF_GET_VAR_REAL(nid,var3didin(2),temp)
    285 #endif
     277         ierr = nf90_get_var(nid,var3didin(2),temp)
    286278         if(ierr/=NF_NOERR) then
    287279            write(*,*) NF_STRERROR(ierr)
    … …  
    290282!          write(*,*)'lecture th ok',temp
    291283
    292 #ifdef NC_DOUBLE
    293          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(3),qv)
    294 #else
    295          ierr = NF_GET_VAR_REAL(nid,var3didin(3),qv)
    296 #endif
     284         ierr = nf90_get_var(nid,var3didin(3),qv)
    297285         if(ierr/=NF_NOERR) then
    298286            write(*,*) NF_STRERROR(ierr)
    … …  
    301289!          write(*,*)'lecture qv ok',qv
    302290 
    303 #ifdef NC_DOUBLE
    304          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(4),u)
    305 #else
    306          ierr = NF_GET_VAR_REAL(nid,var3didin(4),u)
    307 #endif
     291         ierr = nf90_get_var(nid,var3didin(4),u)
    308292         if(ierr/=NF_NOERR) then
    309293            write(*,*) NF_STRERROR(ierr)
    … …  
    312296!          write(*,*)'lecture u ok',u
    313297
    314 #ifdef NC_DOUBLE
    315          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(5),v)
    316 #else
    317          ierr = NF_GET_VAR_REAL(nid,var3didin(5),v)
    318 #endif
     298         ierr = nf90_get_var(nid,var3didin(5),v)
    319299         if(ierr/=NF_NOERR) then
    320300            write(*,*) NF_STRERROR(ierr)
    … …  
    323303!          write(*,*)'lecture v ok',v
    324304
    325 #ifdef NC_DOUBLE
    326          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(6),dw)
    327 #else
    328          ierr = NF_GET_VAR_REAL(nid,var3didin(6),dw)
    329 #endif
     305         ierr = nf90_get_var(nid,var3didin(6),dw)
    330306         if(ierr/=NF_NOERR) then
    331307            write(*,*) NF_STRERROR(ierr)
    … …  
    334310!          write(*,*)'lecture w ok',dw
    335311
    336 #ifdef NC_DOUBLE
    337          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(7),dt)
    338 #else
    339          ierr = NF_GET_VAR_REAL(nid,var3didin(7),dt)
    340 #endif
     312         ierr = nf90_get_var(nid,var3didin(7),dt)
    341313         if(ierr/=NF_NOERR) then
    342314            write(*,*) NF_STRERROR(ierr)
    … …  
    345317!          write(*,*)'lecture dt ok',dt
    346318
    347 #ifdef NC_DOUBLE
    348          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(8),dq)
    349 #else
    350          ierr = NF_GET_VAR_REAL(nid,var3didin(8),dq)
    351 #endif
     319         ierr = nf90_get_var(nid,var3didin(8),dq)
    352320         if(ierr/=NF_NOERR) then
    353321            write(*,*) NF_STRERROR(ierr)
    … …  
    356324!          write(*,*)'lecture dq ok',dq
    357325
    358 #ifdef NC_DOUBLE
    359          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(9),sens)
    360 #else
    361          ierr = NF_GET_VAR_REAL(nid,var3didin(9),sens)
    362 #endif
     326         ierr = nf90_get_var(nid,var3didin(9),sens)
    363327         if(ierr/=NF_NOERR) then
    364328            write(*,*) NF_STRERROR(ierr)
    … …  
    367331!          write(*,*)'lecture sens ok',sens
    368332
    369 #ifdef NC_DOUBLE
    370          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(10),flat)
    371 #else
    372          ierr = NF_GET_VAR_REAL(nid,var3didin(10),flat)
    373 #endif
     333         ierr = nf90_get_var(nid,var3didin(10),flat)
    374334         if(ierr/=NF_NOERR) then
    375335            write(*,*) NF_STRERROR(ierr)
    … …  
    378338!          write(*,*)'lecture flat ok',flat
    379339
    380 #ifdef NC_DOUBLE
    381          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(11),pp)
    382 #else
    383          ierr = NF_GET_VAR_REAL(nid,var3didin(11),pp)
    384 #endif
     340         ierr = nf90_get_var(nid,var3didin(11),pp)
    385341         if(ierr/=NF_NOERR) then
    386342            write(*,*) NF_STRERROR(ierr)
    … …  
    429385
    430386
    431         if (forcing_type.eq.6) then
     387        if (forcing_type==6) then
    432388! Check that initial day of the simulation consistent with AMMA case:
    433        if (annee_ref.ne.2006) then
     389       if (annee_ref/=2006) then
    434390        print*,'Pour AMMA, annee_ref doit etre 2006'
    435391        print*,'Changer annee_ref dans run.def'
    436392        stop
    437393       endif
    438        if (annee_ref.eq.2006 .and. day1.lt.day_ini_amma) then
    439         print*,'AMMA a débuté le 10 juillet 2006',day1,day_ini_amma
     394       if (annee_ref==2006 .and. day1<day_ini_amma) then
     395        print*,'AMMA a d�but� le 10 juillet 2006',day1,day_ini_amma
    440396        print*,'Changer dayref dans run.def'
    441397        stop
    442398       endif
    443        if (annee_ref.eq.2006 .and. day1.gt.day_ini_amma+1) then
     399       if (annee_ref==2006 .and. day1>day_ini_amma+1) then
    444400        print*,'AMMA a fini le 11 juillet'
    445401        print*,'Changer dayref ou nday dans run.def'
    … …  
    464420
    465421       it_amma1=INT(timeit/dt_amma)+1
    466        IF (it_amma1 .EQ. nt_amma) THEN
     422       IF (it_amma1 == nt_amma) THEN
    467423       it_amma2=it_amma1
    468424       ELSE
    … …  
    472428       time_amma2=(it_amma2-1)*dt_amma
    473429
    474        if (it_amma1 .gt. nt_amma) then
     430       if (it_amma1 > nt_amma) then
    475431        write(*,*) 'PB-stop: day, it_amma1, it_amma2, timeit: '            &
    476432     &        ,day,day_ini_amma,it_amma1,it_amma2,timeit/86400.
    … …  
    479435
    480436! time interpolation:
    481        IF (it_amma1 .EQ. it_amma2) THEN
     437       IF (it_amma1 == it_amma2) THEN
    482438          frac=0.
    483439       ELSE

  • LMDZ6/trunk/libf/phylmd/dyn1d/mod_1D_cases_read.F90

    r4593 r5073  
    9999      ierr = NF_OPEN(fich_cas,NF_NOWRITE,nid)
    100100      print*,'fich_cas,NF_NOWRITE,nid ',fich_cas,NF_NOWRITE,nid
    101       if (ierr.NE.NF_NOERR) then
     101      if (ierr/=NF_NOERR) then
    102102         write(*,*) 'ERROR: GROS Pb opening forcings nc file '
    103103         write(*,*) NF_STRERROR(ierr)
    … …  
    106106!.......................................................................
    107107      ierr=NF_INQ_DIMID(nid,'lat',rid)
    108       IF (ierr.NE.NF_NOERR) THEN
     108      IF (ierr/=NF_NOERR) THEN
    109109         print*, 'Oh probleme lecture dimension lat'
    110110      ENDIF
    … …  
    113113!.......................................................................
    114114      ierr=NF_INQ_DIMID(nid,'lon',rid)
    115       IF (ierr.NE.NF_NOERR) THEN
     115      IF (ierr/=NF_NOERR) THEN
    116116         print*, 'Oh probleme lecture dimension lon'
    117117      ENDIF
    … …  
    120120!.......................................................................
    121121      ierr=NF_INQ_DIMID(nid,'lev',rid)
    122       IF (ierr.NE.NF_NOERR) THEN
     122      IF (ierr/=NF_NOERR) THEN
    123123         print*, 'Oh probleme lecture dimension zz'
    124124      ENDIF
    … …  
    129129      print*,'nid,rid',nid,rid
    130130      nt_cas=0
    131       IF (ierr.NE.NF_NOERR) THEN
     131      IF (ierr/=NF_NOERR) THEN
    132132        stop 'probleme lecture dimension sens'
    133133      ENDIF
    … …  
    532532         endif
    533533 
    534 #ifdef NC_DOUBLE
    535          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(1),zz)
    536 #else
    537          ierr = NF_GET_VAR_REAL(nid,var3didin(1),zz)
    538 #endif
     534         ierr = nf90_get_var(nid,var3didin(1),zz)
    539535         if(ierr/=NF_NOERR) then
    540536            write(*,*) NF_STRERROR(ierr)
    … …  
    543539!          write(*,*)'lecture z ok',zz
    544540
    545 #ifdef NC_DOUBLE
    546          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(2),pp)
    547 #else
    548          ierr = NF_GET_VAR_REAL(nid,var3didin(2),pp)
    549 #endif
     541         ierr = nf90_get_var(nid,var3didin(2),pp)
    550542         if(ierr/=NF_NOERR) then
    551543            write(*,*) NF_STRERROR(ierr)
    … …  
    555547
    556548
    557 #ifdef NC_DOUBLE
    558          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(3),temp)
    559 #else
    560          ierr = NF_GET_VAR_REAL(nid,var3didin(3),temp)
    561 #endif
     549         ierr = nf90_get_var(nid,var3didin(3),temp)
    562550         if(ierr/=NF_NOERR) then
    563551            write(*,*) NF_STRERROR(ierr)
    … …  
    566554!          write(*,*)'lecture T ok',temp
    567555
    568 #ifdef NC_DOUBLE
    569          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(4),qv)
    570 #else
    571          ierr = NF_GET_VAR_REAL(nid,var3didin(4),qv)
    572 #endif
     556         ierr = nf90_get_var(nid,var3didin(4),qv)
    573557         if(ierr/=NF_NOERR) then
    574558            write(*,*) NF_STRERROR(ierr)
    … …  
    577561!          write(*,*)'lecture qv ok',qv
    578562 
    579 #ifdef NC_DOUBLE
    580          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(5),rh)
    581 #else
    582          ierr = NF_GET_VAR_REAL(nid,var3didin(5),rh)
    583 #endif
     563         ierr = nf90_get_var(nid,var3didin(5),rh)
    584564         if(ierr/=NF_NOERR) then
    585565            write(*,*) NF_STRERROR(ierr)
    … …  
    588568!          write(*,*)'lecture rh ok',rh
    589569
    590 #ifdef NC_DOUBLE
    591          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(6),theta)
    592 #else
    593          ierr = NF_GET_VAR_REAL(nid,var3didin(6),theta)
    594 #endif
     570         ierr = nf90_get_var(nid,var3didin(6),theta)
    595571         if(ierr/=NF_NOERR) then
    596572            write(*,*) NF_STRERROR(ierr)
    … …  
    599575!          write(*,*)'lecture theta ok',theta
    600576
    601 #ifdef NC_DOUBLE
    602          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(7),rv)
    603 #else
    604          ierr = NF_GET_VAR_REAL(nid,var3didin(7),rv)
    605 #endif
     577         ierr = nf90_get_var(nid,var3didin(7),rv)
    606578         if(ierr/=NF_NOERR) then
    607579            write(*,*) NF_STRERROR(ierr)
    … …  
    610582!          write(*,*)'lecture rv ok',rv
    611583
    612 #ifdef NC_DOUBLE
    613          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(8),u)
    614 #else
    615          ierr = NF_GET_VAR_REAL(nid,var3didin(8),u)
    616 #endif
     584         ierr = nf90_get_var(nid,var3didin(8),u)
    617585         if(ierr/=NF_NOERR) then
    618586            write(*,*) NF_STRERROR(ierr)
    … …  
    621589!          write(*,*)'lecture u ok',u
    622590
    623 #ifdef NC_DOUBLE
    624          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(9),v)
    625 #else
    626          ierr = NF_GET_VAR_REAL(nid,var3didin(9),v)
    627 #endif
     591         ierr = nf90_get_var(nid,var3didin(9),v)
    628592         if(ierr/=NF_NOERR) then
    629593            write(*,*) NF_STRERROR(ierr)
    … …  
    632596!          write(*,*)'lecture v ok',v
    633597
    634 #ifdef NC_DOUBLE
    635          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(10),ug)
    636 #else
    637          ierr = NF_GET_VAR_REAL(nid,var3didin(10),ug)
    638 #endif
     598         ierr = nf90_get_var(nid,var3didin(10),ug)
    639599         if(ierr/=NF_NOERR) then
    640600            write(*,*) NF_STRERROR(ierr)
    … …  
    643603!          write(*,*)'lecture ug ok',ug
    644604
    645 #ifdef NC_DOUBLE
    646          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(11),vg)
    647 #else
    648          ierr = NF_GET_VAR_REAL(nid,var3didin(11),vg)
    649 #endif
     605         ierr = nf90_get_var(nid,var3didin(11),vg)
    650606         if(ierr/=NF_NOERR) then
    651607            write(*,*) NF_STRERROR(ierr)
    … …  
    654610!          write(*,*)'lecture vg ok',vg
    655611
    656 #ifdef NC_DOUBLE
    657          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(12),w)
    658 #else
    659          ierr = NF_GET_VAR_REAL(nid,var3didin(12),w)
    660 #endif
     612         ierr = nf90_get_var(nid,var3didin(12),w)
    661613         if(ierr/=NF_NOERR) then
    662614            write(*,*) NF_STRERROR(ierr)
    … …  
    665617!          write(*,*)'lecture w ok',w
    666618
    667 #ifdef NC_DOUBLE
    668          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(13),du)
    669 #else
    670          ierr = NF_GET_VAR_REAL(nid,var3didin(13),du)
    671 #endif
     619         ierr = nf90_get_var(nid,var3didin(13),du)
    672620         if(ierr/=NF_NOERR) then
    673621            write(*,*) NF_STRERROR(ierr)
    … …  
    676624!          write(*,*)'lecture du ok',du
    677625
    678 #ifdef NC_DOUBLE
    679          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(14),hu)
    680 #else
    681          ierr = NF_GET_VAR_REAL(nid,var3didin(14),hu)
    682 #endif
     626         ierr = nf90_get_var(nid,var3didin(14),hu)
    683627         if(ierr/=NF_NOERR) then
    684628            write(*,*) NF_STRERROR(ierr)
    … …  
    687631!          write(*,*)'lecture hu ok',hu
    688632
    689 #ifdef NC_DOUBLE
    690          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(15),vu)
    691 #else
    692          ierr = NF_GET_VAR_REAL(nid,var3didin(15),vu)
    693 #endif
     633         ierr = nf90_get_var(nid,var3didin(15),vu)
    694634         if(ierr/=NF_NOERR) then
    695635            write(*,*) NF_STRERROR(ierr)
    … …  
    698638!          write(*,*)'lecture vu ok',vu
    699639
    700 #ifdef NC_DOUBLE
    701          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(16),dv)
    702 #else
    703          ierr = NF_GET_VAR_REAL(nid,var3didin(16),dv)
    704 #endif
     640         ierr = nf90_get_var(nid,var3didin(16),dv)
    705641         if(ierr/=NF_NOERR) then
    706642            write(*,*) NF_STRERROR(ierr)
    … …  
    709645!          write(*,*)'lecture dv ok',dv
    710646
    711 #ifdef NC_DOUBLE
    712          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(17),hv)
    713 #else
    714          ierr = NF_GET_VAR_REAL(nid,var3didin(17),hv)
    715 #endif
     647         ierr = nf90_get_var(nid,var3didin(17),hv)
    716648         if(ierr/=NF_NOERR) then
    717649            write(*,*) NF_STRERROR(ierr)
    … …  
    720652!          write(*,*)'lecture hv ok',hv
    721653
    722 #ifdef NC_DOUBLE
    723          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(18),vv)
    724 #else
    725          ierr = NF_GET_VAR_REAL(nid,var3didin(18),vv)
    726 #endif
     654         ierr = nf90_get_var(nid,var3didin(18),vv)
    727655         if(ierr/=NF_NOERR) then
    728656            write(*,*) NF_STRERROR(ierr)
    … …  
    731659!          write(*,*)'lecture vv ok',vv
    732660
    733 #ifdef NC_DOUBLE
    734          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(19),dt)
    735 #else
    736          ierr = NF_GET_VAR_REAL(nid,var3didin(19),dt)
    737 #endif
     661         ierr = nf90_get_var(nid,var3didin(19),dt)
    738662         if(ierr/=NF_NOERR) then
    739663            write(*,*) NF_STRERROR(ierr)
    … …  
    742666!          write(*,*)'lecture dt ok',dt
    743667
    744 #ifdef NC_DOUBLE
    745          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(20),ht)
    746 #else
    747          ierr = NF_GET_VAR_REAL(nid,var3didin(20),ht)
    748 #endif
     668         ierr = nf90_get_var(nid,var3didin(20),ht)
    749669         if(ierr/=NF_NOERR) then
    750670            write(*,*) NF_STRERROR(ierr)
    … …  
    753673!          write(*,*)'lecture ht ok',ht
    754674
    755 #ifdef NC_DOUBLE
    756          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(21),vt)
    757 #else
    758          ierr = NF_GET_VAR_REAL(nid,var3didin(21),vt)
    759 #endif
     675         ierr = nf90_get_var(nid,var3didin(21),vt)
    760676         if(ierr/=NF_NOERR) then
    761677            write(*,*) NF_STRERROR(ierr)
    … …  
    764680!          write(*,*)'lecture vt ok',vt
    765681
    766 #ifdef NC_DOUBLE
    767          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(22),dq)
    768 #else
    769          ierr = NF_GET_VAR_REAL(nid,var3didin(22),dq)
    770 #endif
     682         ierr = nf90_get_var(nid,var3didin(22),dq)
    771683         if(ierr/=NF_NOERR) then
    772684            write(*,*) NF_STRERROR(ierr)
    … …  
    775687!          write(*,*)'lecture dq ok',dq
    776688
    777 #ifdef NC_DOUBLE
    778          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(23),hq)
    779 #else
    780          ierr = NF_GET_VAR_REAL(nid,var3didin(23),hq)
    781 #endif
     689         ierr = nf90_get_var(nid,var3didin(23),hq)
    782690         if(ierr/=NF_NOERR) then
    783691            write(*,*) NF_STRERROR(ierr)
    … …  
    786694!          write(*,*)'lecture hq ok',hq
    787695
    788 #ifdef NC_DOUBLE
    789          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(24),vq)
    790 #else
    791          ierr = NF_GET_VAR_REAL(nid,var3didin(24),vq)
    792 #endif
     696         ierr = nf90_get_var(nid,var3didin(24),vq)
    793697         if(ierr/=NF_NOERR) then
    794698            write(*,*) NF_STRERROR(ierr)
    … …  
    797701!          write(*,*)'lecture vq ok',vq
    798702
    799 #ifdef NC_DOUBLE
    800          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(25),dth)
    801 #else
    802          ierr = NF_GET_VAR_REAL(nid,var3didin(25),dth)
    803 #endif
     703         ierr = nf90_get_var(nid,var3didin(25),dth)
    804704         if(ierr/=NF_NOERR) then
    805705            write(*,*) NF_STRERROR(ierr)
    … …  
    808708!          write(*,*)'lecture dth ok',dth
    809709
    810 #ifdef NC_DOUBLE
    811          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(26),hth)
    812 #else
    813          ierr = NF_GET_VAR_REAL(nid,var3didin(26),hth)
    814 #endif
     710         ierr = nf90_get_var(nid,var3didin(26),hth)
    815711         if(ierr/=NF_NOERR) then
    816712            write(*,*) NF_STRERROR(ierr)
    … …  
    819715!          write(*,*)'lecture hth ok',hth
    820716
    821 #ifdef NC_DOUBLE
    822          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(27),vth)
    823 #else
    824          ierr = NF_GET_VAR_REAL(nid,var3didin(27),vth)
    825 #endif
     717         ierr = nf90_get_var(nid,var3didin(27),vth)
    826718         if(ierr/=NF_NOERR) then
    827719            write(*,*) NF_STRERROR(ierr)
    … …  
    830722!          write(*,*)'lecture vth ok',vth
    831723
    832 #ifdef NC_DOUBLE
    833          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(28),dr)
    834 #else
    835          ierr = NF_GET_VAR_REAL(nid,var3didin(28),dr)
    836 #endif
     724         ierr = nf90_get_var(nid,var3didin(28),dr)
    837725         if(ierr/=NF_NOERR) then
    838726            write(*,*) NF_STRERROR(ierr)
    … …  
    841729!          write(*,*)'lecture dr ok',dr
    842730
    843 #ifdef NC_DOUBLE
    844          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(29),hr)
    845 #else
    846          ierr = NF_GET_VAR_REAL(nid,var3didin(29),hr)
    847 #endif
     731         ierr = nf90_get_var(nid,var3didin(29),hr)
    848732         if(ierr/=NF_NOERR) then
    849733            write(*,*) NF_STRERROR(ierr)
    … …  
    852736!          write(*,*)'lecture hr ok',hr
    853737
    854 #ifdef NC_DOUBLE
    855          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(30),vr)
    856 #else
    857          ierr = NF_GET_VAR_REAL(nid,var3didin(30),vr)
    858 #endif
     738         ierr = nf90_get_var(nid,var3didin(30),vr)
    859739         if(ierr/=NF_NOERR) then
    860740            write(*,*) NF_STRERROR(ierr)
    … …  
    863743!          write(*,*)'lecture vr ok',vr
    864744
    865 #ifdef NC_DOUBLE
    866          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(31),dtrad)
    867 #else
    868          ierr = NF_GET_VAR_REAL(nid,var3didin(31),dtrad)
    869 #endif
     745         ierr = nf90_get_var(nid,var3didin(31),dtrad)
    870746         if(ierr/=NF_NOERR) then
    871747            write(*,*) NF_STRERROR(ierr)
    … …  
    874750!          write(*,*)'lecture dtrad ok',dtrad
    875751
    876 #ifdef NC_DOUBLE
    877          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(32),sens)
    878 #else
    879          ierr = NF_GET_VAR_REAL(nid,var3didin(32),sens)
    880 #endif
     752         ierr = nf90_get_var(nid,var3didin(32),sens)
    881753         if(ierr/=NF_NOERR) then
    882754            write(*,*) NF_STRERROR(ierr)
    … …  
    885757!          write(*,*)'lecture sens ok',sens
    886758
    887 #ifdef NC_DOUBLE
    888          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(33),flat)
    889 #else
    890          ierr = NF_GET_VAR_REAL(nid,var3didin(33),flat)
    891 #endif
     759         ierr = nf90_get_var(nid,var3didin(33),flat)
    892760         if(ierr/=NF_NOERR) then
    893761            write(*,*) NF_STRERROR(ierr)
    … …  
    896764!          write(*,*)'lecture flat ok',flat
    897765
    898 #ifdef NC_DOUBLE
    899          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(34),ts)
    900 #else
    901          ierr = NF_GET_VAR_REAL(nid,var3didin(34),ts)
    902 #endif
     766         ierr = nf90_get_var(nid,var3didin(34),ts)
    903767         if(ierr/=NF_NOERR) then
    904768            write(*,*) NF_STRERROR(ierr)
    … …  
    907771!          write(*,*)'lecture ts ok',ts
    908772
    909 #ifdef NC_DOUBLE
    910          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(35),ustar)
    911 #else
    912          ierr = NF_GET_VAR_REAL(nid,var3didin(35),ustar)
    913 #endif
     773         ierr = nf90_get_var(nid,var3didin(35),ustar)
    914774         if(ierr/=NF_NOERR) then
    915775            write(*,*) NF_STRERROR(ierr)
    … …  
    918778!         write(*,*)'lecture ustar ok',ustar
    919779
    920 #ifdef NC_DOUBLE
    921          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(36),uw)
    922 #else
    923          ierr = NF_GET_VAR_REAL(nid,var3didin(36),uw)
    924 #endif
     780         ierr = nf90_get_var(nid,var3didin(36),uw)
    925781         if(ierr/=NF_NOERR) then
    926782            write(*,*) NF_STRERROR(ierr)
    … …  
    929785!         write(*,*)'lecture uw ok',uw
    930786
    931 #ifdef NC_DOUBLE
    932          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(37),vw)
    933 #else
    934          ierr = NF_GET_VAR_REAL(nid,var3didin(37),vw)
    935 #endif
     787         ierr = nf90_get_var(nid,var3didin(37),vw)
    936788         if(ierr/=NF_NOERR) then
    937789            write(*,*) NF_STRERROR(ierr)
    … …  
    940792!         write(*,*)'lecture vw ok',vw
    941793
    942 #ifdef NC_DOUBLE
    943          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(38),q1)
    944 #else
    945          ierr = NF_GET_VAR_REAL(nid,var3didin(38),q1)
    946 #endif
     794         ierr = nf90_get_var(nid,var3didin(38),q1)
    947795         if(ierr/=NF_NOERR) then
    948796            write(*,*) NF_STRERROR(ierr)
    … …  
    951799!         write(*,*)'lecture q1 ok',q1
    952800
    953 #ifdef NC_DOUBLE
    954          ierr = NF_GET_VAR_DOUBLE(nid,var3didin(39),q2)
    955 #else
    956          ierr = NF_GET_VAR_REAL(nid,var3didin(39),q2)
    957 #endif
     801         ierr = nf90_get_var(nid,var3didin(39),q2)
    958802         if(ierr/=NF_NOERR) then
    959803            write(*,*) NF_STRERROR(ierr)
    … …  
    1081925
    1082926       it_cas1=INT(timeit/pdt_cas)+1
    1083        IF (it_cas1 .EQ. nt_cas) THEN
     927       IF (it_cas1 == nt_cas) THEN
    1084928       it_cas2=it_cas1
    1085929       ELSE
    … …  
    1093937      print *,'time_cas2=',time_cas2
    1094938
    1095        if (it_cas1 .gt. nt_cas) then
     939       if (it_cas1 > nt_cas) then
    1096940        write(*,*) 'PB-stop: day, day_ju_ini_cas,it_cas1, it_cas2, timeit: '            &
    1097941     &        ,day,day_ju_ini_cas,it_cas1,it_cas2,timeit
    … …  
    1100944
    1101945! time interpolation:
    1102        IF (it_cas1 .EQ. it_cas2) THEN
     946       IF (it_cas1 == it_cas2) THEN
    1103947          frac=0.
    1104948       ELSE
Note: See TracChangeset for help on using the changeset viewer.