Ignore:
Timestamp:
Mar 26, 2009, 1:30:23 PM (15 years ago)
Author:
lguez
Message:

Translated calls using NetCDF 2.4 interface to calls using NetCDF 3.6
Fortran 90 interface.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/read_pstoke0.F

    r940 r1135  
    1717C******************************************************************************
    1818
    19 
     19        use netcdf
    2020       USE dimphy
    2121       IMPLICIT NONE
     
    121121      if (irec .eq. 0) then
    122122
    123             ncidp=NCOPN('phystoke.nc',NCNOWRIT,rcode)
    124 
    125             varidps=NCVID(ncidp,'phis',rcode)
     123            rcode=nf90_open('phystoke.nc',nf90_nowrite,ncidp)
     124
     125            rcode = nf90_inq_varid(ncidp, 'phis', varidps)
    126126            print*,'ncidp,varidps',ncidp,varidps
    127127
    128             varidpl=NCVID(ncidp,'sig_s',rcode)
     128            rcode = nf90_inq_varid(ncidp, 'sig_s', varidpl)
    129129            print*,'ncidp,varidpl',ncidp,varidpl
    130130
    131             varidai=NCVID(ncidp,'aire',rcode)
     131            rcode = nf90_inq_varid(ncidp, 'aire', varidai)
    132132            print*,'ncidp,varidai',ncidp,varidai
    133133
    134                 varidt=NCVID(ncidp,'t',rcode)
     134                rcode = nf90_inq_varid(ncidp, 't', varidt)
    135135                print*,'ncidp,varidt',ncidp,varidt
    136136
    137             varidmfu=NCVID(ncidp,'mfu',rcode)
     137            rcode = nf90_inq_varid(ncidp, 'mfu', varidmfu)
    138138            print*,'ncidp,varidmfu',ncidp,varidmfu
    139139
    140             varidmfd=NCVID(ncidp,'mfd',rcode)
     140            rcode = nf90_inq_varid(ncidp, 'mfd', varidmfd)
    141141            print*,'ncidp,varidmfd',ncidp,varidmfd
    142142
    143             varidenu=NCVID(ncidp,'en_u',rcode)
     143            rcode = nf90_inq_varid(ncidp, 'en_u', varidenu)
    144144            print*,'ncidp,varidenu',ncidp,varidenu
    145145
    146             variddeu=NCVID(ncidp,'de_u',rcode)
     146            rcode = nf90_inq_varid(ncidp, 'de_u', variddeu)
    147147            print*,'ncidp,variddeu',ncidp,variddeu
    148148
    149             varidend=NCVID(ncidp,'en_d',rcode)
     149            rcode = nf90_inq_varid(ncidp, 'en_d', varidend)
    150150            print*,'ncidp,varidend',ncidp,varidend
    151151       
    152             varidded=NCVID(ncidp,'de_d',rcode)
     152            rcode = nf90_inq_varid(ncidp, 'de_d', varidded)
    153153            print*,'ncidp,varidded',ncidp,varidded
    154154       
    155             varidch=NCVID(ncidp,'coefh',rcode)
     155            rcode = nf90_inq_varid(ncidp, 'coefh', varidch)
    156156            print*,'ncidp,varidch',ncidp,varidch
    157157
    158158c Thermiques
    159             varidfmth=NCVID(ncidp,'fm_th',rcode)
     159            rcode = nf90_inq_varid(ncidp, 'fm_th', varidfmth)
    160160            print*,'ncidp,varidfmth',ncidp,varidfmth
    161161
    162             varidenth=NCVID(ncidp,'en_th',rcode)
     162            rcode = nf90_inq_varid(ncidp, 'en_th', varidenth)
    163163            print*,'ncidp,varidenth',ncidp,varidenth
    164164       
    165             varidfi=NCVID(ncidp,'frac_impa',rcode)
     165            rcode = nf90_inq_varid(ncidp, 'frac_impa', varidfi)
    166166            print*,'ncidp,varidfi',ncidp,varidfi
    167167       
    168             varidfn=NCVID(ncidp,'frac_nucl',rcode)
     168            rcode = nf90_inq_varid(ncidp, 'frac_nucl', varidfn)
    169169            print*,'ncidp,varidfn',ncidp,varidfn
    170170       
    171             varidyu1=NCVID(ncidp,'pyu1',rcode)
     171            rcode = nf90_inq_varid(ncidp, 'pyu1', varidyu1)
    172172            print*,'ncidp,varidyu1',ncidp,varidyu1
    173173       
    174             varidyv1=NCVID(ncidp,'pyv1',rcode)
     174            rcode = nf90_inq_varid(ncidp, 'pyv1', varidyv1)
    175175            print*,'ncidp,varidyv1',ncidp,varidyv1
    176176       
    177             varidfts1=NCVID(ncidp,'ftsol1',rcode)
     177            rcode = nf90_inq_varid(ncidp, 'ftsol1', varidfts1)
    178178            print*,'ncidp,varidfts1',ncidp,varidfts1
    179179       
    180             varidfts2=NCVID(ncidp,'ftsol2',rcode)
     180            rcode = nf90_inq_varid(ncidp, 'ftsol2', varidfts2)
    181181            print*,'ncidp,varidfts2',ncidp,varidfts2
    182182         
    183             varidfts3=NCVID(ncidp,'ftsol3',rcode)
     183            rcode = nf90_inq_varid(ncidp, 'ftsol3', varidfts3)
    184184            print*,'ncidp,varidfts3',ncidp,varidfts3
    185185 
    186             varidfts4=NCVID(ncidp,'ftsol4',rcode)
     186            rcode = nf90_inq_varid(ncidp, 'ftsol4', varidfts4)
    187187            print*,'ncidp,varidfts4',ncidp,varidfts4
    188188       
    189             varidpsr1=NCVID(ncidp,'psrf1',rcode)
     189            rcode = nf90_inq_varid(ncidp, 'psrf1', varidpsr1)
    190190            print*,'ncidp,varidpsr1',ncidp,varidpsr1
    191191       
    192             varidpsr2=NCVID(ncidp,'psrf2',rcode)
     192            rcode = nf90_inq_varid(ncidp, 'psrf2', varidpsr2)
    193193            print*,'ncidp,varidpsr2',ncidp,varidpsr2
    194194       
    195             varidpsr3=NCVID(ncidp,'psrf3',rcode)
     195            rcode = nf90_inq_varid(ncidp, 'psrf3', varidpsr3)
    196196            print*,'ncidp,varidpsr3',ncidp,varidpsr3
    197197
    198             varidpsr4=NCVID(ncidp,'psrf4',rcode)
     198            rcode = nf90_inq_varid(ncidp, 'psrf4', varidpsr4)
    199199            print*,'ncidp,varidpsr4',ncidp,varidpsr4
    200200       
Note: See TracChangeset for help on using the changeset viewer.