Changeset 1188 for LMDZ4/branches/LMDZ4-dev/libf/dyn3d

Timestamp:

Jun 18, 2009, 5:50:11 PM (15 years ago)

Author:

lguez

Message:

Translated calls using NetCDF 2.4 interface to calls using NetCDF 3.6
Fortran 90 interface.
Corrected a few comments in "output.def".

File:

• LMDZ4/branches/LMDZ4-dev/libf/dyn3d/guide_mod.F90 (modified) (13 diffs)

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/guide_mod.F90

@@ -11 +11 @@
   USE getparam
   USE Write_Field
+  use netcdf, only: nf90_nowrite, nf90_open, nf90_inq_varid, nf90_close
 
   IMPLICIT NONE
@@ -140 +141 @@
     ncidpl=-99
     if (guide_modele) then
-       if (ncidpl.eq.-99) ncidpl=NCOPN('apbp.nc',NCNOWRIT,rcod)
+       if (ncidpl.eq.-99) rcod=nf90_open('apbp.nc',Nf90_NOWRITe, ncidpl)
     else
          if (guide_u) then
-           if (ncidpl.eq.-99) ncidpl=NCOPN('u.nc',NCNOWRIT,rcod)
+           if (ncidpl.eq.-99) rcod=nf90_open('u.nc',Nf90_NOWRITe,ncidpl)
          elseif (guide_v) then
-           if (ncidpl.eq.-99) ncidpl=NCOPN('v.nc',NCNOWRIT,rcod)
+           if (ncidpl.eq.-99) rcod=nf90_open('v.nc',nf90_nowrite,ncidpl)
          elseif (guide_T) then
-           if (ncidpl.eq.-99) ncidpl=NCOPN('T.nc',NCNOWRIT,rcod)
+           if (ncidpl.eq.-99) rcod=nf90_open('T.nc',nf90_nowrite,ncidpl)
          elseif (guide_Q) then
-           if (ncidpl.eq.-99) ncidpl=NCOPN('hur.nc',NCNOWRIT,rcod)
+           if (ncidpl.eq.-99) rcod=nf90_open('hur.nc',nf90_nowrite, ncidpl)
          endif
     endif 
@@ -160 +161 @@
     error=NF_INQ_DIMLEN(ncidpl,rid,nlevnc)
     print *,'Guide: nombre niveaux vert. nlevnc', nlevnc 
-    CALL NCCLOS(ncidpl,rcod)
+    rcod = nf90_close(ncidpl)
 
 ! ---------------------------------------------
@@ -992 +993 @@
          if (guide_modele) then
              print *,'Lecture du guidage sur niveaux mod�le'
-             ncidpl=NCOPN('apbp.nc',NCNOWRIT,rcode)
-             varidap=NCVID(ncidpl,'AP',rcode)
-             varidbp=NCVID(ncidpl,'BP',rcode)
+             rcode = nf90_open('apbp.nc', nf90_nowrite, ncidpl)
+             rcode = nf90_inq_varid(ncidpl, 'AP', varidap)
+             rcode = nf90_inq_varid(ncidpl, 'BP', varidbp)
              print*,'ncidpl,varidap',ncidpl,varidap
          endif
 ! Vent zonal
          if (guide_u) then
-             ncidu=NCOPN('u.nc',NCNOWRIT,rcode)
-             varidu=NCVID(ncidu,'UWND',rcode)
+             rcode = nf90_open('u.nc', nf90_nowrite, ncidu)
+             rcode = nf90_inq_varid(ncidu, 'UWND', varidu)
              print*,'ncidu,varidu',ncidu,varidu
              if (ncidpl.eq.-99) ncidpl=ncidu
@@ -1006 +1007 @@
 ! Vent meridien
          if (guide_v) then
-             ncidv=NCOPN('v.nc',NCNOWRIT,rcode)
-             varidv=NCVID(ncidv,'VWND',rcode)
+             rcode = nf90_open('v.nc', nf90_nowrite, ncidv)
+             rcode = nf90_inq_varid(ncidv, 'VWND', varidv)
              print*,'ncidv,varidv',ncidv,varidv
              if (ncidpl.eq.-99) ncidpl=ncidv
@@ -1013 +1014 @@
 ! Temperature
          if (guide_T) then
-             ncidt=NCOPN('T.nc',NCNOWRIT,rcode)
-             varidt=NCVID(ncidt,'AIR',rcode)
+             rcode = nf90_open('T.nc', nf90_nowrite, ncidt)
+             rcode = nf90_inq_varid(ncidt, 'AIR', varidt)
              print*,'ncidT,varidT',ncidt,varidt
              if (ncidpl.eq.-99) ncidpl=ncidt
@@ -1020 +1021 @@
 ! Humidite
          if (guide_Q) then
-             ncidQ=NCOPN('hur.nc',NCNOWRIT,rcode)
-             varidQ=NCVID(ncidQ,'RH',rcode)
+             rcode = nf90_open('hur.nc', nf90_nowrite, ncidQ)
+             rcode = nf90_inq_varid(ncidQ, 'RH', varidQ)
              print*,'ncidQ,varidQ',ncidQ,varidQ
              if (ncidpl.eq.-99) ncidpl=ncidQ
@@ -1027 +1028 @@
 ! Pression de surface
          if ((guide_P).OR.(guide_modele)) then
-             ncidps=NCOPN('ps.nc',NCNOWRIT,rcode)
-             varidps=NCVID(ncidps,'SP',rcode)
+             rcode = nf90_open('ps.nc', nf90_nowrite, ncidps)
+             rcode = nf90_inq_varid(ncidps, 'SP', varidps)
              print*,'ncidps,varidps',ncidps,varidps
          endif
 ! Coordonnee verticale
          if (.not.guide_modele) then
-              rcode=NF_INQ_VARID(ncidpl,'LEVEL',varidpl)
-              IF (rcode.NE.0) rcode=NF_INQ_VARID(ncidpl,'PRESSURE',varidpl)
+              rcode = nf90_inq_varid(ncidpl, 'LEVEL', varidpl)
+              IF (rcode.NE.0) rcode = nf90_inq_varid(ncidpl, 'PRESSURE', varidpl)
               print*,'ncidpl,varidpl',ncidpl,varidpl
          endif
@@ -1175 +1176 @@
          if (guide_modele) then
              print *,'Lecture du guidage sur niveaux mod�le'
-             ncidpl=NCOPN('apbp.nc',NCNOWRIT,rcode)
-             varidap=NCVID(ncidpl,'AP',rcode)
-             varidbp=NCVID(ncidpl,'BP',rcode)
+             rcode = nf90_open('apbp.nc', nf90_nowrite, ncidpl)
+             rcode = nf90_inq_varid(ncidpl, 'AP', varidap)
+             rcode = nf90_inq_varid(ncidpl, 'BP', varidbp)
              print*,'ncidpl,varidap',ncidpl,varidap
          endif
 ! Vent zonal
          if (guide_u) then
-             ncidu=NCOPN('u.nc',NCNOWRIT,rcode)
-             varidu=NCVID(ncidu,'UWND',rcode)
+             rcode = nf90_open('u.nc', nf90_nowrite, ncidu)
+             rcode = nf90_inq_varid(ncidu, 'UWND', varidu)
              print*,'ncidu,varidu',ncidu,varidu
              if (ncidpl.eq.-99) ncidpl=ncidu
@@ -1189 +1190 @@
 ! Vent meridien
          if (guide_v) then
-             ncidv=NCOPN('v.nc',NCNOWRIT,rcode)
-             varidv=NCVID(ncidv,'VWND',rcode)
+             rcode = nf90_open('v.nc', nf90_nowrite, ncidv)
+             rcode = nf90_inq_varid(ncidv, 'VWND', varidv)
              print*,'ncidv,varidv',ncidv,varidv
              if (ncidpl.eq.-99) ncidpl=ncidv
@@ -1196 +1197 @@
 ! Temperature
          if (guide_T) then
-             ncidt=NCOPN('T.nc',NCNOWRIT,rcode)
-             varidt=NCVID(ncidt,'AIR',rcode)
+             rcode = nf90_open('T.nc', nf90_nowrite, ncidt)
+             rcode = nf90_inq_varid(ncidt, 'AIR', varidt)
              print*,'ncidT,varidT',ncidt,varidt
              if (ncidpl.eq.-99) ncidpl=ncidt
@@ -1203 +1204 @@
 ! Humidite
          if (guide_Q) then
-             ncidQ=NCOPN('hur.nc',NCNOWRIT,rcode)
-             varidQ=NCVID(ncidQ,'RH',rcode)
+             rcode = nf90_open('hur.nc', nf90_nowrite, ncidQ)
+             rcode = nf90_inq_varid(ncidQ, 'RH', varidQ)
              print*,'ncidQ,varidQ',ncidQ,varidQ
              if (ncidpl.eq.-99) ncidpl=ncidQ
@@ -1210 +1211 @@
 ! Pression de surface
          if ((guide_P).OR.(guide_modele)) then
-             ncidps=NCOPN('ps.nc',NCNOWRIT,rcode)
-             varidps=NCVID(ncidps,'SP',rcode)
+             rcode = nf90_open('ps.nc', nf90_nowrite, ncidps)
+             rcode = nf90_inq_varid(ncidps, 'SP', varidps)
              print*,'ncidps,varidps',ncidps,varidps
          endif
 ! Coordonnee verticale
          if (.not.guide_modele) then
-              rcode=NF_INQ_VARID(ncidpl,'LEVEL',varidpl)
-              IF (rcode.NE.0) rcode=NF_INQ_VARID(ncidpl,'PRESSURE',varidpl)
+              rcode = nf90_inq_varid(ncidpl, 'LEVEL', varidpl)
+              IF (rcode.NE.0) rcode = nf90_inq_varid(ncidpl, 'PRESSURE', varidpl)
               print*,'ncidpl,varidpl',ncidpl,varidpl
          endif