Changeset 1135 for LMDZ4


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

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

Location:
LMDZ4/branches/LMDZ4-dev/libf
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ4/branches/LMDZ4-dev/libf/dyn3d/guide.F

    r1046 r1135  
    33!
    44      subroutine guide(itau,ucov,vcov,teta,q,masse,ps)
     5
     6      use netcdf
    57
    68      IMPLICIT NONE
    … …  
    225227c lecture d'un fichier netcdf pour determiner le nombre de niveaux
    226228         if (guide_modele) then
    227            if (ncidpl.eq.-99) ncidpl=NCOPN('apbp.nc',NCNOWRIT,rcod)
     229           if (ncidpl.eq.-99) rcod=nf90_open('apbp.nc',Nf90_NOWRITe,
     230     $           ncidpl)
    228231         else
    229232         if (guide_u) then
    230            if (ncidpl.eq.-99) ncidpl=NCOPN('u.nc',NCNOWRIT,rcod)
     233           if (ncidpl.eq.-99) rcod=nf90_open('u.nc',Nf90_NOWRITe,ncidpl)
    231234         endif
    232235c
    233236         if (guide_v) then
    234            if (ncidpl.eq.-99) ncidpl=NCOPN('v.nc',NCNOWRIT,rcod)
     237           if (ncidpl.eq.-99) rcod=nf90_open('v.nc',nf90_nowrite,ncidpl)
    235238         endif
    236239c
    237240         if (guide_T) then
    238            if (ncidpl.eq.-99) ncidpl=NCOPN('T.nc',NCNOWRIT,rcod)
     241           if (ncidpl.eq.-99) rcod=nf90_open('T.nc',nf90_nowrite,ncidpl)
    239242         endif
    240243c
    241244         if (guide_Q) then
    242            if (ncidpl.eq.-99) ncidpl=NCOPN('hur.nc',NCNOWRIT,rcod)
     245           if (ncidpl.eq.-99) rcod=nf90_open('hur.nc',nf90_nowrite,
     246     $           ncidpl)
    243247         endif
    244248c
    … …  
    251255          status=NF_INQ_DIMLEN(ncidpl,rid,nlev)
    252256         print *,'nlev guide', nlev
    253          call ncclos(ncidpl,rcod)
     257         rcod = nf90_close(ncidpl)
    254258c   Lecture du premier etat des reanalyses.
    255259         call read_reanalyse(1,ps

  • LMDZ4/branches/LMDZ4-dev/libf/dyn3d/read_reanalyse.F

    r1120 r1135  
    1313c   Declarations
    1414c -----------------------------------------------------------------
     15      use netcdf
     16
    1517      IMPLICIT NONE
    1618
    … …  
    7274            print *,'Vous êtes entrain de lire des données sur
    7375     .               niveaux modèle'
    74             ncidpl=NCOPN('apbp.nc',NCNOWRIT,rcode)
    75             varidap=NCVID(ncidpl,'AP',rcode)
    76             varidbp=NCVID(ncidpl,'BP',rcode)
     76            rcode=nf90_open('apbp.nc',nf90_nowrite,ncidpl)
     77            rcode = nf90_inq_varid(ncidpl, 'AP', varidap)
     78            rcode = nf90_inq_varid(ncidpl, 'BP', varidbp)
    7779            print*,'ncidpl,varidap',ncidpl,varidap
    7880            endif
    … …  
    8082c Vent zonal
    8183            if (guide_u) then
    82             ncidu=NCOPN('u.nc',NCNOWRIT,rcode)
    83             varidu=NCVID(ncidu,'UWND',rcode)
    84             print*,'ncidu,varidu',ncidu,varidu
    85             if (ncidpl.eq.-99) ncidpl=ncidu
     84               rcode=nf90_open('u.nc',nf90_nowrite,ncidu)
     85               rcode = nf90_inq_varid(ncidu, 'UWND', varidu)
     86               print*,'ncidu,varidu',ncidu,varidu
     87               if (ncidpl.eq.-99) ncidpl=ncidu
    8688            endif
    8789
    8890c Vent meridien
    8991            if (guide_v) then
    90             ncidv=NCOPN('v.nc',NCNOWRIT,rcode)
    91             varidv=NCVID(ncidv,'VWND',rcode)
     92            rcode=nf90_open('v.nc',nf90_nowrite,ncidv)
     93            rcode = nf90_inq_varid(ncidv, 'VWND', varidv)
    9294            print*,'ncidv,varidv',ncidv,varidv
    9395            if (ncidpl.eq.-99) ncidpl=ncidv
    … …  
    9698c Temperature
    9799            if (guide_T) then
    98             ncidt=NCOPN('T.nc',NCNOWRIT,rcode)
    99             varidt=NCVID(ncidt,'AIR',rcode)
     100            rcode=nf90_open('T.nc',nf90_nowrite,ncidt)
     101            rcode = nf90_inq_varid(ncidt, 'AIR', varidt)
    100102            print*,'ncidt,varidt',ncidt,varidt
    101103            if (ncidpl.eq.-99) ncidpl=ncidt
    … …  
    104106c Humidite
    105107            if (guide_Q) then
    106             ncidQ=NCOPN('hur.nc',NCNOWRIT,rcode)
    107             varidQ=NCVID(ncidQ,'RH',rcode)
     108            rcode=nf90_open('hur.nc',nf90_nowrite,ncidQ)
     109            rcode = nf90_inq_varid(ncidQ, 'RH', varidQ)
    108110            print*,'ncidQ,varidQ',ncidQ,varidQ
    109111            if (ncidpl.eq.-99) ncidpl=ncidQ
    … …  
    112114c Pression de surface
    113115            if ((guide_P).OR.(guide_modele)) then
    114             ncidps=NCOPN('ps.nc',NCNOWRIT,rcode)
    115             varidps=NCVID(ncidps,'SP',rcode)
     116            rcode=nf90_open('ps.nc',nf90_nowrite,ncidps)
     117            rcode = nf90_inq_varid(ncidps, 'SP', varidps)
    116118            print*,'ncidps,varidps',ncidps,varidps
    117119            endif
    … …  
    119121c Coordonnee verticale
    120122            if (.not.guide_modele) then
    121               if (ncep) then
    122                print*,'Vous etes entrain de lire des donnees NCEP'
    123                varidpl=NCVID(ncidpl,'LEVEL',rcode)
    124               else
    125                print*,'Vous etes entrain de lire des donnees ECMWF'
    126                varidpl=NCVID(ncidpl,'PRESSURE',rcode)
    127               endif
    128               print*,'ncidpl,varidpl',ncidpl,varidpl
     123               if (ncep) then
     124                  print*,'Vous etes entrain de lire des donnees NCEP'
     125                  rcode = nf90_inq_varid(ncidpl, 'LEVEL', varidpl)
     126               else
     127                  print*,'Vous etes entrain de lire des donnees ECMWF'
     128                  rcode = nf90_inq_varid(ncidpl, 'PRESSURE', varidpl)
     129               endif
     130               print*,'ncidpl,varidpl',ncidpl,varidpl
    129131            endif
    130132! endif (first)

  • LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/guide_p.F

    r1049 r1135  
    44      subroutine guide_pp(itau,ucov,vcov,teta,q,masse,ps)
    55      USE parallel
     6      use netcdf
    67
    78      IMPLICIT NONE
    … …  
    229230         IF (mpi_rank==0) THEN
    230231          if (guide_modele) then
    231             if (ncidpl.eq.-99) ncidpl=NCOPN('apbp.nc',NCNOWRIT,rcod)
    232           else 
    233            if (guide_u) then
    234             if (ncidpl.eq.-99) ncidpl=NCOPN('u.nc',NCNOWRIT,rcod)
    235           endif
    236 c
    237           if (guide_v) then
    238             if (ncidpl.eq.-99) ncidpl=NCOPN('v.nc',NCNOWRIT,rcod)
    239           endif
    240 c
    241           if (guide_T) then
    242             if (ncidpl.eq.-99) ncidpl=NCOPN('T.nc',NCNOWRIT,rcod)
    243           endif
    244 c
    245           if (guide_Q) then
    246             if (ncidpl.eq.-99) ncidpl=NCOPN('hur.nc',NCNOWRIT,rcod)
    247           endif
     232             if (ncidpl.eq.-99) rcod=nf90_open('apbp.nc',Nf90_NOWRITe,
     233     $            ncidpl)
     234          else
     235             if (guide_u) then
     236                if (ncidpl.eq.-99) rcod=nf90_open('u.nc',Nf90_NOWRITe,
     237     $               ncidpl)
     238             endif
     239c
     240             if (guide_v) then
     241                if (ncidpl.eq.-99) rcod=nf90_open('v.nc',nf90_nowrite,
     242     $               ncidpl)
     243             endif
     244c
     245             if (guide_T) then
     246                if (ncidpl.eq.-99) rcod=nf90_open('T.nc',nf90_nowrite,
     247     $               ncidpl)
     248             endif
     249c
     250             if (guide_Q) then
     251                if (ncidpl.eq.-99) rcod=nf90_open('hur.nc',nf90_nowrite,
     252     $               ncidpl)
     253             endif
    248254c
    249255          endif  !guide_modele
    … …  
    256262          status=NF_INQ_DIMLEN(ncidpl,rid,nlev)
    257263         print *,'nlev guide', nlev
    258          call ncclos(ncidpl,rcod)
     264         rcod = nf90_close(ncidpl)
    259265c   Lecture du premier etat des reanalyses.
    260266         call Gather_Field(ps,ip1jmp1,1,0)

  • LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/read_reanalyse.F

    r1121 r1135  
    1111
    1212       USE parallel
     13       use netcdf
    1314c -----------------------------------------------------------------
    1415c   Declarations
    … …  
    7273            print *,'Vous êtes entrain de lire des données sur
    7374     .               niveaux modèle'
    74             ncidpl=NCOPN('apbp.nc',NCNOWRIT,rcode)
    75             varidap=NCVID(ncidpl,'AP',rcode)
    76             varidbp=NCVID(ncidpl,'BP',rcode)
     75            rcode=nf90_open('apbp.nc',nf90_nowrite,ncidpl)
     76            rcode = nf90_inq_varid(ncidpl, 'AP', varidap)
     77            rcode = nf90_inq_varid(ncidpl, 'BP', varidbp)
    7778            print*,'ncidpl,varidap',ncidpl,varidap
    7879            endif
    … …  
    8081c Vent zonal
    8182            if (guide_u) then
    82                ncidu=NCOPN('u.nc',NCNOWRIT,rcode)
    83                varidu=NCVID(ncidu,'UWND',rcode)
     83               rcode=nf90_open('u.nc',nf90_nowrite,ncidu)
     84               rcode = nf90_inq_varid(ncidu, 'UWND', varidu)
    8485               print*,'ncidu,varidu',ncidu,varidu
    8586               if (ncidpl.eq.-99) ncidpl=ncidu
    … …  
    8889c Vent meridien
    8990            if (guide_v) then
    90                ncidv=NCOPN('v.nc',NCNOWRIT,rcode)
    91                varidv=NCVID(ncidv,'VWND',rcode)
     91               rcode=nf90_open('v.nc',nf90_nowrite,ncidv)
     92               rcode = nf90_inq_varid(ncidv, 'VWND', varidv)
    9293               print*,'ncidv,varidv',ncidv,varidv
    9394               if (ncidpl.eq.-99) ncidpl=ncidv
    … …  
    9697c Temperature
    9798            if (guide_T) then
    98                ncidt=NCOPN('T.nc',NCNOWRIT,rcode)
    99                varidt=NCVID(ncidt,'AIR',rcode)
     99               rcode=nf90_open('T.nc',nf90_nowrite,ncidt)
     100               rcode = nf90_inq_varid(ncidt, 'AIR', varidt)
    100101               print*,'ncidt,varidt',ncidt,varidt
    101102               if (ncidpl.eq.-99) ncidpl=ncidt
    … …  
    104105c Humidite
    105106            if (guide_Q) then
    106                ncidQ=NCOPN('hur.nc',NCNOWRIT,rcode)
    107                varidQ=NCVID(ncidQ,'RH',rcode)
     107               rcode=nf90_open('hur.nc',nf90_nowrite,ncidQ)
     108               rcode = nf90_inq_varid(ncidQ, 'RH', varidQ)
    108109               print*,'ncidQ,varidQ',ncidQ,varidQ
    109110               if (ncidpl.eq.-99) ncidpl=ncidQ
    … …  
    112113c Pression de surface
    113114            if ((guide_P).OR.(guide_modele)) then
    114                ncidps=NCOPN('ps.nc',NCNOWRIT,rcode)
    115                varidps=NCVID(ncidps,'SP',rcode)
     115               rcode=nf90_open('ps.nc',nf90_nowrite,ncidps)
     116               rcode = nf90_inq_varid(ncidps, 'SP', varidps)
    116117               print*,'ncidps,varidps',ncidps,varidps
    117118            endif
    … …  
    119120c Coordonnee verticale
    120121            if (.not.guide_modele) then
    121               if (ncep) then
    122                print*,'Vous etes entrain de lire des donnees NCEP'
    123                varidpl=NCVID(ncidpl,'LEVEL',rcode)
    124             else
    125                print*,'Vous etes entrain de lire des donnees ECMWF'
    126                varidpl=NCVID(ncidpl,'PRESSURE',rcode)
    127               endif
    128               print*,'ncidpl,varidpl',ncidpl,varidpl
     122               if (ncep) then
     123                  print*,'Vous etes entrain de lire des donnees NCEP'
     124                  rcode = nf90_inq_varid(ncidpl, 'LEVEL', varidpl)
     125               else
     126                  print*,'Vous etes entrain de lire des donnees ECMWF'
     127                  rcode = nf90_inq_varid(ncidpl, 'PRESSURE', varidpl)
     128               endif
     129               print*,'ncidpl,varidpl',ncidpl,varidpl
    129130            endif
    130131            print*,'ncidu,varidpl',ncidu,varidpl

  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/read_pstoke.F

    r940 r1135  
    1818C******************************************************************************
    1919
     20        use netcdf
    2021       USE dimphy
    2122       IMPLICIT NONE
    … …  
    116117      if (irec .eq. 0) then
    117118
    118             ncidp=NCOPN('phystoke.nc',NCNOWRIT,rcode)
    119 
    120             varidps=NCVID(ncidp,'phis',rcode)
     119            rcode=nf90_open('phystoke.nc',nf90_nowrite,ncidp)
     120
     121            rcode = nf90_inq_varid(ncidp, 'phis', varidps)
    121122            print*,'ncidp,varidps',ncidp,varidps
    122123
    123             varidpl=NCVID(ncidp,'sig_s',rcode)
     124            rcode = nf90_inq_varid(ncidp, 'sig_s', varidpl)
    124125            print*,'ncidp,varidpl',ncidp,varidpl
    125126
    126             varidai=NCVID(ncidp,'aire',rcode)
     127            rcode = nf90_inq_varid(ncidp, 'aire', varidai)
    127128            print*,'ncidp,varidai',ncidp,varidai
    128129
    129130c A FAIRE: Es-il necessaire de stocke t?
    130                 varidt=NCVID(ncidp,'t',rcode)
     131                rcode = nf90_inq_varid(ncidp, 't', varidt)
    131132                print*,'ncidp,varidt',ncidp,varidt
    132133
    133             varidmfu=NCVID(ncidp,'mfu',rcode)
     134            rcode = nf90_inq_varid(ncidp, 'mfu', varidmfu)
    134135            print*,'ncidp,varidmfu',ncidp,varidmfu
    135136
    136             varidmfd=NCVID(ncidp,'mfd',rcode)
     137            rcode = nf90_inq_varid(ncidp, 'mfd', varidmfd)
    137138            print*,'ncidp,varidmfd',ncidp,varidmfd
    138139
    139             varidenu=NCVID(ncidp,'en_u',rcode)
     140            rcode = nf90_inq_varid(ncidp, 'en_u', varidenu)
    140141            print*,'ncidp,varidenu',ncidp,varidenu
    141142
    142             variddeu=NCVID(ncidp,'de_u',rcode)
     143            rcode = nf90_inq_varid(ncidp, 'de_u', variddeu)
    143144            print*,'ncidp,variddeu',ncidp,variddeu
    144145
    145             varidend=NCVID(ncidp,'en_d',rcode)
     146            rcode = nf90_inq_varid(ncidp, 'en_d', varidend)
    146147            print*,'ncidp,varidend',ncidp,varidend
    147148       
    148             varidded=NCVID(ncidp,'de_d',rcode)
     149            rcode = nf90_inq_varid(ncidp, 'de_d', varidded)
    149150            print*,'ncidp,varidded',ncidp,varidded
    150151       
    151             varidch=NCVID(ncidp,'coefh',rcode)
     152            rcode = nf90_inq_varid(ncidp, 'coefh', varidch)
    152153            print*,'ncidp,varidch',ncidp,varidch
    153154       
    154155c abder (pour thermiques)
    155              varidfmth=NCVID(ncidp,'fm_th',rcode)
     156             rcode = nf90_inq_varid(ncidp, 'fm_th', varidfmth)
    156157             print*,'ncidp,varidfmth',ncidp,varidfmth
    157158
    158              varidenth=NCVID(ncidp,'en_th',rcode)
     159             rcode = nf90_inq_varid(ncidp, 'en_th', varidenth)
    159160             print*,'ncidp,varidenth',ncidp,varidenth
    160161
    161             varidfi=NCVID(ncidp,'frac_impa',rcode)
     162            rcode = nf90_inq_varid(ncidp, 'frac_impa', varidfi)
    162163            print*,'ncidp,varidfi',ncidp,varidfi
    163164       
    164             varidfn=NCVID(ncidp,'frac_nucl',rcode)
     165            rcode = nf90_inq_varid(ncidp, 'frac_nucl', varidfn)
    165166            print*,'ncidp,varidfn',ncidp,varidfn
    166167       
    167             varidyu1=NCVID(ncidp,'pyu1',rcode)
     168            rcode = nf90_inq_varid(ncidp, 'pyu1', varidyu1)
    168169            print*,'ncidp,varidyu1',ncidp,varidyu1
    169170       
    170             varidyv1=NCVID(ncidp,'pyv1',rcode)
     171            rcode = nf90_inq_varid(ncidp, 'pyv1', varidyv1)
    171172            print*,'ncidp,varidyv1',ncidp,varidyv1
    172173       
    173             varidfts1=NCVID(ncidp,'ftsol1',rcode)
     174            rcode = nf90_inq_varid(ncidp, 'ftsol1', varidfts1)
    174175            print*,'ncidp,varidfts1',ncidp,varidfts1
    175176       
    176             varidfts2=NCVID(ncidp,'ftsol2',rcode)
     177            rcode = nf90_inq_varid(ncidp, 'ftsol2', varidfts2)
    177178            print*,'ncidp,varidfts2',ncidp,varidfts2
    178179         
    179             varidfts3=NCVID(ncidp,'ftsol3',rcode)
     180            rcode = nf90_inq_varid(ncidp, 'ftsol3', varidfts3)
    180181            print*,'ncidp,varidfts3',ncidp,varidfts3
    181182 
    182             varidfts4=NCVID(ncidp,'ftsol4',rcode)
     183            rcode = nf90_inq_varid(ncidp, 'ftsol4', varidfts4)
    183184            print*,'ncidp,varidfts4',ncidp,varidfts4
    184185       
    185             varidpsr1=NCVID(ncidp,'psrf1',rcode)
     186            rcode = nf90_inq_varid(ncidp, 'psrf1', varidpsr1)
    186187            print*,'ncidp,varidpsr1',ncidp,varidpsr1
    187188       
    188             varidpsr2=NCVID(ncidp,'psrf2',rcode)
     189            rcode = nf90_inq_varid(ncidp, 'psrf2', varidpsr2)
    189190            print*,'ncidp,varidpsr2',ncidp,varidpsr2
    190191       
    191             varidpsr3=NCVID(ncidp,'psrf3',rcode)
     192            rcode = nf90_inq_varid(ncidp, 'psrf3', varidpsr3)
    192193            print*,'ncidp,varidpsr3',ncidp,varidpsr3
    193194
    194             varidpsr4=NCVID(ncidp,'psrf4',rcode)
     195            rcode = nf90_inq_varid(ncidp, 'psrf4', varidpsr4)
    195196            print*,'ncidp,varidpsr4',ncidp,varidpsr4
    196197       

  • 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.