Changeset 71 for LMDZ.3.3/trunk/libf/dyn3d

Timestamp:

Mar 7, 2000, 4:47:15 PM (25 years ago)

Author:

lmdzadmin

Message:

Les lectures netcdf en simple precision sont passees en double pour resoudre
un probleme netcdf sur VPP.
LF d'apres P.Braconnot

File:

• LMDZ.3.3/trunk/libf/dyn3d/create_limit.F (modified) (18 diffs)

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

@@ -137 +137 @@
         STOP
       ENDIF
-      ierr = NF_GET_VAR_REAL(ncid,dimid,dlon_msk)
+c     ierr = NF_GET_VAR_REAL(ncid,dimid,dlon_msk)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,dlon_msk)
+c$$$      zbidon=0.0
+c$$$      ierr = NF_GET_VAR_REAL(ncid,dimid,zbidon(1:imdep))
+c$$$      dlon_msk(1 : imdep) = dble(zbidon(1:imdep))
+c
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -153 +158 @@
         STOP
       ENDIF
-      ierr = NF_GET_VAR_REAL(ncid,dimid,dlat_msk)
-      if (ierr.ne.0) then
-        print *, NF_STRERROR(ierr)
-        STOP
-      ENDIF
-      ierr = NF_GET_VAR_REAL(ncid,varid,champ_msk)
+c      ierr = NF_GET_VAR_REAL(ncid,dimid,dlat_msk)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,dlat_msk)
+c$$$      zbidon=0.
+c$$$      ierr = NF_GET_VAR_REAL(ncid,dimid,zbidon(1:jmdep))
+c$$$      dlat_msk=dble(zbidon(1:jmdep))
+c
+      if (ierr.ne.0) then
+        print *, NF_STRERROR(ierr)
+        STOP
+      ENDIF
+c      ierr = NF_GET_VAR_REAL(ncid,varid,champ_msk)
+      ierr = NF_GET_VAR_DOUBLE(ncid,varid,champ_msk)
+c$$$      zbidon=0.
+c$$$      ierr = NF_GET_VAR_REAL(ncid,varid,zbidon(1:imdep*jmdep))
+c$$$      champ_msk(1: imdep*jmdep) = zbidon(1:imdep*jmdep)
+c
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -211 +226 @@
         STOP
       ENDIF 
-      ierr = NF_GET_VAR_REAL(ncid,dimid,dlon)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,dlon)
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -227 +242 @@
         STOP
       ENDIF 
-      ierr = NF_GET_VAR_REAL(ncid,dimid,dlat)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,dlat)
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -243 +258 @@
         STOP
       ENDIF
-      ierr = NF_GET_VAR_REAL(ncid,dimid,timecoord)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,timecoord)
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -260 +275 @@
          PRINT*,'Lecture temporelle et int. horizontale ',l,timecoord(l)
          print*,dimfirst,dimlast
-         ierr = NF_GET_VARA_REAL(ncid,varid,dimfirst,dimlast,champ)
+         ierr = NF_GET_VARA_DOUBLE(ncid,varid,dimfirst,dimlast,champ)
          if (ierr.ne.0) then
            print *, NF_STRERROR(ierr)
@@ -343 +358 @@
         STOP
       ENDIF 
-      ierr = NF_GET_VAR_REAL(ncid,dimid,dlon)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,dlon)
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -359 +374 @@
         STOP
       ENDIF 
-      ierr = NF_GET_VAR_REAL(ncid,dimid,dlat)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,dlat)
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -375 +390 @@
         STOP
       ENDIF 
-      ierr = NF_GET_VAR_REAL(ncid,dimid,timecoord)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,timecoord)
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -391 +406 @@
 c
          PRINT*,'Lecture temporelle et int. horizontale ',l,timecoord(l)
-         ierr = NF_GET_VARA_REAL(ncid,varid,dimfirst,dimlast,champ)
+         ierr = NF_GET_VARA_DOUBLE(ncid,varid,dimfirst,dimlast,champ)
          if (ierr.ne.0) then
            print *, NF_STRERROR(ierr)
@@ -480 +495 @@
         STOP
       ENDIF 
-      ierr = NF_GET_VAR_REAL(ncid,dimid,dlon)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,dlon)
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -496 +511 @@
         STOP
       ENDIF 
-      ierr = NF_GET_VAR_REAL(ncid,dimid,dlat)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,dlat)
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -512 +527 @@
         STOP
       ENDIF 
-      ierr = NF_GET_VAR_REAL(ncid,dimid,timecoord)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,timecoord)
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -528 +543 @@
 c
          PRINT*,'Lecture temporelle et int. horizontale ',l,timecoord(l)
-         ierr = NF_GET_VARA_REAL(ncid,varid,dimfirst,dimlast,champ)
+         ierr = NF_GET_VARA_DOUBLE(ncid,varid,dimfirst,dimlast,champ)
          if (ierr.ne.0) then
            print *, NF_STRERROR(ierr)
@@ -606 +621 @@
         STOP
       ENDIF 
-      ierr = NF_GET_VAR_REAL(ncid,dimid,dlon)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,dlon)
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -622 +637 @@
         STOP
       ENDIF 
-      ierr = NF_GET_VAR_REAL(ncid,dimid,dlat)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,dlat)
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -638 +653 @@
         STOP
       ENDIF 
-      ierr = NF_GET_VAR_REAL(ncid,dimid,timecoord)
+      ierr = NF_GET_VAR_DOUBLE(ncid,dimid,timecoord)
       if (ierr.ne.0) then
         print *, NF_STRERROR(ierr)
@@ -654 +669 @@
 c
          PRINT*,'Lecture temporelle et int. horizontale ',l,timecoord(l)
-         ierr = NF_GET_VARA_REAL(ncid,varid,dimfirst,dimlast,champ)
+         ierr = NF_GET_VARA_DOUBLE(ncid,varid,dimfirst,dimlast,champ)
          if (ierr.ne.0) then
            print *, NF_STRERROR(ierr)