Changeset 1233 for LMDZ4/trunk/libf/phylmd/undefSTD.F

Timestamp:

Aug 28, 2009, 11:13:28 AM (15 years ago)

Author:

musat

Message:

Pour la coherence avec netcdf : on met les memes _Fillvalue pour les variables a un niveau de pression
IM

File:

• LMDZ4/trunk/libf/phylmd/undefSTD.F (modified) (5 diffs)

LMDZ4/trunk/libf/phylmd/undefSTD.F

@@ -5 +5 @@
      $           dtime,ecrit_hf,
      $           oknondef,tnondef,tsumSTD)
+      USE netcdf
       USE dimphy
       IMPLICIT none
@@ -13 +14 @@
 c
 c Calcul * du nombre de pas de temps (FLOAT(ecrit_XXX)-tnondef)) 
-c          ou la variable tlevSTD est bien definie (.NE.1.E+20), 
+c          ou la variable tlevSTD est bien definie (.NE.missing_val), 
 c et 
 c        * de la somme de tlevSTD => tsumSTD
@@ -45 +46 @@
       REAL tnondef(klon,klevSTD,nout)
 c
+      REAL missing_val
+c
+      missing_val=nf90_fill_real
+c
 c calcul variables tous les pas de temps de la physique 
 c
       DO k=1, nlevSTD
        DO i=1, klon
-        IF(tlevSTD(i,k).EQ.1E+20) THEN
+        IF(tlevSTD(i,k).EQ.missing_val) THEN
          IF(oknondef(i,k,1)) THEN          
           tnondef(i,k,1)=tnondef(i,k,1)+1.
@@ -58 +63 @@
          ENDIF !oknondef(i,k)
 c
-        ELSE IF(tlevSTD(i,k).NE.1E+20) THEN
+        ELSE IF(tlevSTD(i,k).NE.missing_val) THEN
          tsumSTD(i,k,1)=tsumSTD(i,k,1)+tlevSTD(i,k)
          tsumSTD(i,k,2)=tsumSTD(i,k,2)+tlevSTD(i,k)
@@ -71 +76 @@
        DO k=1, nlevSTD
         DO i=1, klon
-         IF(tlevSTD(i,k).EQ.1E+20) THEN
+         IF(tlevSTD(i,k).EQ.missing_val) THEN
           IF(oknondef(i,k,3)) THEN          
            tnondef(i,k,3)=tnondef(i,k,3)+1.
           ENDIF !oknondef(i,k)
 c
-         ELSE IF(tlevSTD(i,k).NE.1E+20) THEN
+         ELSE IF(tlevSTD(i,k).NE.missing_val) THEN
          tsumSTD(i,k,3)=tsumSTD(i,k,3)+tlevSTD(i,k)
          ENDIF