Changeset 1135 for LMDZ4/branches/LMDZ4-dev/libf/dyn3d/guide.F

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:

• LMDZ4/branches/LMDZ4-dev/libf/dyn3d/guide.F (modified) (3 diffs)

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

@@ -3 +3 @@
 !
       subroutine guide(itau,ucov,vcov,teta,q,masse,ps)
+
+      use netcdf
 
       IMPLICIT NONE
@@ -225 +227 @@
 c lecture d'un fichier netcdf pour determiner le nombre de niveaux
          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)
          endif
 c
          if (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)
          endif
 c
          if (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)
          endif
 c
          if (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
 c
@@ -251 +255 @@
           status=NF_INQ_DIMLEN(ncidpl,rid,nlev)
          print *,'nlev guide', nlev 
-         call ncclos(ncidpl,rcod)
+         rcod = nf90_close(ncidpl)
 c   Lecture du premier etat des reanalyses.
          call read_reanalyse(1,ps