Changeset 997 for LMDZ4/trunk/libf/dyn3d

Timestamp:

Sep 9, 2008, 4:31:22 PM (16 years ago)

Author:

lsce

Message:

Remetais/ajoutais ecriture de limit.nc en double precision pour un compilation avec le cle NC_DOUBLE.

JG

File:

• LMDZ4/trunk/libf/dyn3d/limit_netcdf.F (modified) (4 diffs)

LMDZ4/trunk/libf/dyn3d/limit_netcdf.F

@@ -1177 +1177 @@
       ENDDO
 c
-c Ecriture de fichier limit toujours en simple precision, NF_FLOAT
       PRINT*, 'Ecriture du fichier limit'
 c
@@ -1190 +1189 @@
       dims(2) = ntim
 c
+#ifdef NC_DOUBLE
+      ierr = NF_DEF_VAR (nid, "TEMPS", NF_DOUBLE, 1,ntim, id_tim)
+#else
       ierr = NF_DEF_VAR (nid, "TEMPS", NF_FLOAT, 1,ntim, id_tim)
+#endif
       ierr = NF_PUT_ATT_TEXT (nid, id_tim, "title", 17,
      .                        "Jour dans l annee")
       IF (newlmt) THEN
 c
+#ifdef NC_DOUBLE
+        ierr = NF_DEF_VAR (nid, "FOCE", NF_DOUBLE, 2,dims, id_FOCE)
+#else
         ierr = NF_DEF_VAR (nid, "FOCE", NF_FLOAT, 2,dims, id_FOCE)
+#endif
         ierr = NF_PUT_ATT_TEXT (nid, id_FOCE, "title", 14,
      .                      "Fraction ocean")
 c
+#ifdef NC_DOUBLE
+        ierr = NF_DEF_VAR (nid, "FSIC", NF_DOUBLE, 2,dims, id_FSIC)
+#else
         ierr = NF_DEF_VAR (nid, "FSIC", NF_FLOAT, 2,dims, id_FSIC)
+#endif
         ierr = NF_PUT_ATT_TEXT (nid, id_FSIC, "title", 21,
      .                      "Fraction glace de mer")
 c
+#ifdef NC_DOUBLE
+        ierr = NF_DEF_VAR (nid, "FTER", NF_DOUBLE, 2,dims, id_FTER)
+#else
         ierr = NF_DEF_VAR (nid, "FTER", NF_FLOAT, 2,dims, id_FTER)
+#endif
         ierr = NF_PUT_ATT_TEXT (nid, id_FTER, "title", 14,
      .                      "Fraction terre")
 c
+#ifdef NC_DOUBLE
+        ierr = NF_DEF_VAR (nid, "FLIC", NF_DOUBLE, 2,dims, id_FLIC)
+#else
         ierr = NF_DEF_VAR (nid, "FLIC", NF_FLOAT, 2,dims, id_FLIC)
+#endif
         ierr = NF_PUT_ATT_TEXT (nid, id_FLIC, "title", 17,
      .                      "Fraction land ice")
 c
       ELSE 
+#ifdef NC_DOUBLE
+        ierr = NF_DEF_VAR (nid, "NAT", NF_DOUBLE, 2,dims, id_NAT)
+#else
         ierr = NF_DEF_VAR (nid, "NAT", NF_FLOAT, 2,dims, id_NAT)
+#endif
         ierr = NF_PUT_ATT_TEXT (nid, id_NAT, "title", 23,
      .                      "Nature du sol (0,1,2,3)")
       ENDIF 
+#ifdef NC_DOUBLE
+      ierr = NF_DEF_VAR (nid, "SST", NF_DOUBLE, 2,dims, id_SST)
+#else
       ierr = NF_DEF_VAR (nid, "SST", NF_FLOAT, 2,dims, id_SST)
+#endif
       ierr = NF_PUT_ATT_TEXT (nid, id_SST, "title", 35,
      .                      "Temperature superficielle de la mer")
+#ifdef NC_DOUBLE
+      ierr = NF_DEF_VAR (nid, "BILS", NF_DOUBLE, 2,dims, id_BILS)
+#else
       ierr = NF_DEF_VAR (nid, "BILS", NF_FLOAT, 2,dims, id_BILS)
+#endif
       ierr = NF_PUT_ATT_TEXT (nid, id_BILS, "title", 32,
      .                        "Reference flux de chaleur au sol")
+#ifdef NC_DOUBLE
+      ierr = NF_DEF_VAR (nid, "ALB", NF_DOUBLE, 2,dims, id_ALB)
+#else
       ierr = NF_DEF_VAR (nid, "ALB", NF_FLOAT, 2,dims, id_ALB)
+#endif
       ierr = NF_PUT_ATT_TEXT (nid, id_ALB, "title", 19,
      .                        "Albedo a la surface")
+#ifdef NC_DOUBLE
+      ierr = NF_DEF_VAR (nid, "RUG", NF_DOUBLE, 2,dims, id_RUG)
+#else
       ierr = NF_DEF_VAR (nid, "RUG", NF_FLOAT, 2,dims, id_RUG)
+#endif
       ierr = NF_PUT_ATT_TEXT (nid, id_RUG, "title", 8,
      .                        "Rugosite")
@@ -1238 +1277 @@
       epais(2) = 1
 c
+#ifdef NC_DOUBLE
+      ierr = NF_PUT_VAR1_DOUBLE (nid,id_tim,k,DBLE(k))
+c
+      IF (newlmt ) THEN
+          ierr = NF_PUT_VARA_DOUBLE (nid,id_FOCE,debut,epais
+     $        ,pctsrf_t(1,is_oce,k))
+          ierr = NF_PUT_VARA_DOUBLE (nid,id_FSIC,debut,epais
+     $        ,pctsrf_t(1,is_sic,k))
+          ierr = NF_PUT_VARA_DOUBLE (nid,id_FTER,debut,epais
+     $        ,pctsrf_t(1,is_ter,k))
+          ierr = NF_PUT_VARA_DOUBLE (nid,id_FLIC,debut,epais
+     $        ,pctsrf_t(1,is_lic,k))
+      ELSE 
+          ierr = NF_PUT_VARA_DOUBLE (nid,id_NAT,debut,epais
+     $        ,phy_nat(1,k))
+      ENDIF 
+c
+      ierr = NF_PUT_VARA_DOUBLE (nid,id_SST,debut,epais,phy_sst(1,k))
+      ierr = NF_PUT_VARA_DOUBLE (nid,id_BILS,debut,epais,phy_bil(1,k))
+      ierr = NF_PUT_VARA_DOUBLE (nid,id_ALB,debut,epais,phy_alb(1,k))
+      ierr = NF_PUT_VARA_DOUBLE (nid,id_RUG,debut,epais,phy_rug(1,k))
+#else
       ierr = NF_PUT_VAR1_REAL (nid,id_tim,k,FLOAT(k))
       IF (newlmt ) THEN
@@ -1256 +1317 @@
       ierr = NF_PUT_VARA_REAL (nid,id_ALB,debut,epais,phy_alb(1,k))
       ierr = NF_PUT_VARA_REAL (nid,id_RUG,debut,epais,phy_rug(1,k))
+#endif
 c
       ENDDO