Changeset 218 for LMDZ.3.3/trunk/libf/dyn3d/read_fstoke0.F

Timestamp:

May 14, 2001, 5:01:36 PM (23 years ago)

Author:

lmdz

Message:

Generalisation des if NC_DOUBLE pour passer sur VPP MAF
LF

File:

• LMDZ.3.3/trunk/libf/dyn3d/read_fstoke0.F (modified) (7 diffs)

LMDZ.3.3/trunk/libf/dyn3d/read_fstoke0.F

@@ -24 +24 @@
      . rlatu_dy(iip1,jjp1),rlatv_dy(iip1,jjm)
       integer*4 ncrec,ncim,ncjm,nclm
-      real zrec,zim,zjm,zlm
+      integer*4 zrec,zim,zjm,zlm
       integer*4 xid,yid,zid,tid
       real*4 zdtvr,ziadvtrac
@@ -122 +122 @@
 c niveaux de pression
         
+#ifdef NC_DOUBLE
+       status=NF_GET_VARA_DOUBLE(ncidf,varidpl,1,zlm,pl)
+#else
        status=NF_GET_VARA_REAL(ncidf,varidpl,1,zlm,pl)
+#endif
 
 c rlonu,rlonv,rlatu,rlatv
@@ -136 +140 @@
       count(4)=0
         
-        status=NF_GET_VARA_REAL(ncidf,varidnlo,start,count,rlonu_dy)
-        status=NF_GET_VARA_REAL(ncidf,varidnla,start,count,rlatu_dy)
-        
+#ifdef NC_DOUBLE
+      status=NF_GET_VARA_DOUBLE(ncidf,varidnlo,start,count,rlonu_dy)
+      status=NF_GET_VARA_DOUBLE(ncidf,varidnla,start,count,rlatu_dy)
+#else
+      status=NF_GET_VARA_REAL(ncidf,varidnlo,start,count,rlonu_dy)
+      status=NF_GET_VARA_REAL(ncidf,varidnla,start,count,rlatu_dy)
+#endif
+
       start(1)=1
       start(2)=1
@@ -149 +158 @@
       count(4)=0
 
-        status=NF_GET_VARA_REAL(ncidfv,varidnlov,start,count,rlonv_dy)
-        status=NF_GET_VARA_REAL(ncidfv,varidnlav,start,count,rlatv_dy)
+#ifdef NC_DOUBLE
+      status=NF_GET_VARA_DOUBLE(ncidfv,varidnlov,start,count,rlonv_dy)
+      status=NF_GET_VARA_DOUBLE(ncidfv,varidnlav,start,count,rlatv_dy)
+#else
+      status=NF_GET_VARA_REAL(ncidfv,varidnlov,start,count,rlonv_dy)
+      status=NF_GET_VARA_REAL(ncidfv,varidnlav,start,count,rlatv_dy)
+#endif
 
 c Lecture de phis et aire
@@ -165 +179 @@
 
 c phis
+#ifdef NC_DOUBLE
+      status=NF_GET_VARA_DOUBLE(ncidf,varidps,start,count,phis)
+#else
       status=NF_GET_VARA_REAL(ncidf,varidps,start,count,phis)
+#endif
 
 c aire
+#ifdef NC_DOUBLE
+       status=NF_GET_VARA_DOUBLE(ncidf,varidai,start,count,airedy)
+#else
        status=NF_GET_VARA_REAL(ncidf,varidai,start,count,airedy)
+#endif
 
        else 
@@ -193 +215 @@
 c masse 
         PRINT*,'LECTURE DE masse'
+#ifdef NC_DOUBLE
+      status=NF_GET_VARA_DOUBLE(ncidf,varidm,start,count,masse)
+#else
       status=NF_GET_VARA_REAL(ncidf,varidm,start,count,masse)
+#endif
 
 c pbaru 
+#ifdef NC_DOUBLE
+      status=NF_GET_VARA_DOUBLE(ncidf,varidpu,start,count,pbaru)
+#else
       status=NF_GET_VARA_REAL(ncidf,varidpu,start,count,pbaru)
+#endif
 
 c w 
+#ifdef NC_DOUBLE
+      status=NF_GET_VARA_DOUBLE(ncidf,varidw,start,count,w)
+#else
       status=NF_GET_VARA_REAL(ncidf,varidw,start,count,w)
+#endif
 
 c teta 
+#ifdef NC_DOUBLE
+      status=NF_GET_VARA_DOUBLE(ncidf,varidt,start,count,teta)
+#else
       status=NF_GET_VARA_REAL(ncidf,varidt,start,count,teta)
+#endif
 
 c phi 
+#ifdef NC_DOUBLE
+      status=NF_GET_VARA_DOUBLE(ncidf,varidph,start,count,phi)
+#else
       status=NF_GET_VARA_REAL(ncidf,varidph,start,count,phi)
+#endif
         
         count(2) = zjm-1 
 
 c  pbarv
+#ifdef NC_DOUBLE
+      status=NF_GET_VARA_DOUBLE(ncidfv,varidpv,start,count,pbarv)
+#else
       status=NF_GET_VARA_REAL(ncidfv,varidpv,start,count,pbarv)
+#endif
 
       start(3)=irec
@@ -222 +268 @@
       return
       end
-