Changeset 1135 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/read_pstoke.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/phylmd/read_pstoke.F (modified) (2 diffs)

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

@@ -18 +18 @@
 C******************************************************************************
 
+        use netcdf
        USE dimphy
        IMPLICIT NONE
@@ -116 +117 @@
       if (irec .eq. 0) then
 
-            ncidp=NCOPN('phystoke.nc',NCNOWRIT,rcode)
-
-            varidps=NCVID(ncidp,'phis',rcode)
+            rcode=nf90_open('phystoke.nc',nf90_nowrite,ncidp)
+
+            rcode = nf90_inq_varid(ncidp, 'phis', varidps)
             print*,'ncidp,varidps',ncidp,varidps
 
-            varidpl=NCVID(ncidp,'sig_s',rcode)
+            rcode = nf90_inq_varid(ncidp, 'sig_s', varidpl)
             print*,'ncidp,varidpl',ncidp,varidpl
 
-            varidai=NCVID(ncidp,'aire',rcode)
+            rcode = nf90_inq_varid(ncidp, 'aire', varidai)
             print*,'ncidp,varidai',ncidp,varidai
 
 c A FAIRE: Es-il necessaire de stocke t?
-                varidt=NCVID(ncidp,'t',rcode)
+                rcode = nf90_inq_varid(ncidp, 't', varidt)
                 print*,'ncidp,varidt',ncidp,varidt
 
-            varidmfu=NCVID(ncidp,'mfu',rcode)
+            rcode = nf90_inq_varid(ncidp, 'mfu', varidmfu)
             print*,'ncidp,varidmfu',ncidp,varidmfu
 
-            varidmfd=NCVID(ncidp,'mfd',rcode)
+            rcode = nf90_inq_varid(ncidp, 'mfd', varidmfd)
             print*,'ncidp,varidmfd',ncidp,varidmfd
 
-            varidenu=NCVID(ncidp,'en_u',rcode)
+            rcode = nf90_inq_varid(ncidp, 'en_u', varidenu)
             print*,'ncidp,varidenu',ncidp,varidenu
 
-            variddeu=NCVID(ncidp,'de_u',rcode)
+            rcode = nf90_inq_varid(ncidp, 'de_u', variddeu)
             print*,'ncidp,variddeu',ncidp,variddeu
 
-            varidend=NCVID(ncidp,'en_d',rcode)
+            rcode = nf90_inq_varid(ncidp, 'en_d', varidend)
             print*,'ncidp,varidend',ncidp,varidend
         
-            varidded=NCVID(ncidp,'de_d',rcode)
+            rcode = nf90_inq_varid(ncidp, 'de_d', varidded)
             print*,'ncidp,varidded',ncidp,varidded
         
-            varidch=NCVID(ncidp,'coefh',rcode)
+            rcode = nf90_inq_varid(ncidp, 'coefh', varidch)
             print*,'ncidp,varidch',ncidp,varidch
         
 c abder (pour thermiques)
-             varidfmth=NCVID(ncidp,'fm_th',rcode)
+             rcode = nf90_inq_varid(ncidp, 'fm_th', varidfmth)
              print*,'ncidp,varidfmth',ncidp,varidfmth
 
-             varidenth=NCVID(ncidp,'en_th',rcode)
+             rcode = nf90_inq_varid(ncidp, 'en_th', varidenth)
              print*,'ncidp,varidenth',ncidp,varidenth
 
-            varidfi=NCVID(ncidp,'frac_impa',rcode)
+            rcode = nf90_inq_varid(ncidp, 'frac_impa', varidfi)
             print*,'ncidp,varidfi',ncidp,varidfi
         
-            varidfn=NCVID(ncidp,'frac_nucl',rcode)
+            rcode = nf90_inq_varid(ncidp, 'frac_nucl', varidfn)
             print*,'ncidp,varidfn',ncidp,varidfn
         
-            varidyu1=NCVID(ncidp,'pyu1',rcode)
+            rcode = nf90_inq_varid(ncidp, 'pyu1', varidyu1)
             print*,'ncidp,varidyu1',ncidp,varidyu1
         
-            varidyv1=NCVID(ncidp,'pyv1',rcode)
+            rcode = nf90_inq_varid(ncidp, 'pyv1', varidyv1)
             print*,'ncidp,varidyv1',ncidp,varidyv1
         
-            varidfts1=NCVID(ncidp,'ftsol1',rcode)
+            rcode = nf90_inq_varid(ncidp, 'ftsol1', varidfts1)
             print*,'ncidp,varidfts1',ncidp,varidfts1
         
-            varidfts2=NCVID(ncidp,'ftsol2',rcode)
+            rcode = nf90_inq_varid(ncidp, 'ftsol2', varidfts2)
             print*,'ncidp,varidfts2',ncidp,varidfts2
          
-            varidfts3=NCVID(ncidp,'ftsol3',rcode)
+            rcode = nf90_inq_varid(ncidp, 'ftsol3', varidfts3)
             print*,'ncidp,varidfts3',ncidp,varidfts3
   
-            varidfts4=NCVID(ncidp,'ftsol4',rcode)
+            rcode = nf90_inq_varid(ncidp, 'ftsol4', varidfts4)
             print*,'ncidp,varidfts4',ncidp,varidfts4
         
-            varidpsr1=NCVID(ncidp,'psrf1',rcode)
+            rcode = nf90_inq_varid(ncidp, 'psrf1', varidpsr1)
             print*,'ncidp,varidpsr1',ncidp,varidpsr1
         
-            varidpsr2=NCVID(ncidp,'psrf2',rcode)
+            rcode = nf90_inq_varid(ncidp, 'psrf2', varidpsr2)
             print*,'ncidp,varidpsr2',ncidp,varidpsr2
         
-            varidpsr3=NCVID(ncidp,'psrf3',rcode)
+            rcode = nf90_inq_varid(ncidp, 'psrf3', varidpsr3)
             print*,'ncidp,varidpsr3',ncidp,varidpsr3
 
-            varidpsr4=NCVID(ncidp,'psrf4',rcode)
+            rcode = nf90_inq_varid(ncidp, 'psrf4', varidpsr4)
             print*,'ncidp,varidpsr4',ncidp,varidpsr4