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

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/dyn3d

Files:

• guide.F (modified) (3 diffs)
• read_reanalyse.F (modified) (7 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

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

@@ -13 +13 @@
 c   Declarations
 c -----------------------------------------------------------------
+      use netcdf
+
       IMPLICIT NONE
 
@@ -72 +74 @@
             print *,'Vous êtes entrain de lire des données 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
@@ -80 +82 @@
 c Vent zonal
             if (guide_u) then
-            ncidu=NCOPN('u.nc',NCNOWRIT,rcode)
-            varidu=NCVID(ncidu,'UWND',rcode)
-            print*,'ncidu,varidu',ncidu,varidu
-            if (ncidpl.eq.-99) ncidpl=ncidu
+               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
             endif
 
 c 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
@@ -96 +98 @@
 c 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
@@ -104 +106 @@
 c 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
@@ -112 +114 @@
 c 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
@@ -119 +121 @@
 c Coordonnee verticale
             if (.not.guide_modele) then
-              if (ncep) then
-               print*,'Vous etes entrain de lire des donnees NCEP'
-               varidpl=NCVID(ncidpl,'LEVEL',rcode)
-              else
-               print*,'Vous etes entrain de lire des donnees ECMWF'
-               varidpl=NCVID(ncidpl,'PRESSURE',rcode)
-              endif
-              print*,'ncidpl,varidpl',ncidpl,varidpl
+               if (ncep) then
+                  print*,'Vous etes entrain de lire des donnees NCEP'
+                  rcode = nf90_inq_varid(ncidpl, 'LEVEL', varidpl)
+               else
+                  print*,'Vous etes entrain de lire des donnees ECMWF'
+                  rcode = nf90_inq_varid(ncidpl, 'PRESSURE', varidpl)
+               endif
+               print*,'ncidpl,varidpl',ncidpl,varidpl
             endif
 ! endif (first)