Changeset 1135 for LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/guide_p.F

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.

File:

• LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/guide_p.F (modified) (3 diffs)

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

@@ -4 +4 @@
       subroutine guide_pp(itau,ucov,vcov,teta,q,masse,ps)
       USE parallel
+      use netcdf
 
       IMPLICIT NONE
@@ -229 +230 @@
          IF (mpi_rank==0) THEN
           if (guide_modele) then
-            if (ncidpl.eq.-99) ncidpl=NCOPN('apbp.nc',NCNOWRIT,rcod)
-          else  
-           if (guide_u) then
-            if (ncidpl.eq.-99) ncidpl=NCOPN('u.nc',NCNOWRIT,rcod)
-          endif
-c
-          if (guide_v) then
-            if (ncidpl.eq.-99) ncidpl=NCOPN('v.nc',NCNOWRIT,rcod)
-          endif
-c
-          if (guide_T) then
-            if (ncidpl.eq.-99) ncidpl=NCOPN('T.nc',NCNOWRIT,rcod)
-          endif
-c
-          if (guide_Q) then
-            if (ncidpl.eq.-99) ncidpl=NCOPN('hur.nc',NCNOWRIT,rcod)
-          endif
+             if (ncidpl.eq.-99) rcod=nf90_open('apbp.nc',Nf90_NOWRITe,
+     $            ncidpl)
+          else
+             if (guide_u) then
+                if (ncidpl.eq.-99) rcod=nf90_open('u.nc',Nf90_NOWRITe,
+     $               ncidpl)
+             endif
+c
+             if (guide_v) then
+                if (ncidpl.eq.-99) rcod=nf90_open('v.nc',nf90_nowrite,
+     $               ncidpl)
+             endif
+c
+             if (guide_T) then
+                if (ncidpl.eq.-99) rcod=nf90_open('T.nc',nf90_nowrite,
+     $               ncidpl)
+             endif
+c
+             if (guide_Q) then
+                if (ncidpl.eq.-99) rcod=nf90_open('hur.nc',nf90_nowrite,
+     $               ncidpl)
+             endif
 c
           endif  !guide_modele
@@ -256 +262 @@
           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 Gather_Field(ps,ip1jmp1,1,0)