Changeset 1233

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

Location:

LMDZ4/trunk/libf/phylmd

Files:

• moy_undefSTD.F (modified) (6 diffs)
• plevel.F (modified) (3 diffs)
• plevel_new.F (modified) (3 diffs)
• undefSTD.F (modified) (5 diffs)

LMDZ4/trunk/libf/phylmd/moy_undefSTD.F

@@ -5 +5 @@
      $           dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
      $           tnondef,tsumSTD)
+      USE netcdf
       USE dimphy
       IMPLICIT none
@@ -13 +14 @@
 c
 c Moyenne - a des frequences differentes - des valeurs bien definies
-c         (.NE.1.E+20) des variables interpolees a un niveau de
+c         (.NE.missing_val) des variables interpolees a un niveau de
 c         pression.
 c 1) les variables de type "day" (nout=1) ou "mth" (nout=2) sont sommees
@@ -42 +43 @@
       REAL tsumSTD(klon,klevSTD,nout)
 c
+      REAL missing_val
+c
+      missing_val=nf90_fill_real
+c
 c calcul 1 fois par jour
 c
@@ -51 +56 @@
      $    (ecrit_day/dtime-tnondef(i,k,1))
          ELSE
-          tsumSTD(i,k,1)=1.E+20
+          tsumSTD(i,k,1)=missing_val
          ENDIF !tnondef
         ENDDO !i
@@ -66 +71 @@
      $    (ecrit_mth/dtime-tnondef(i,k,2))
          ELSE
-          tsumSTD(i,k,2)=1.E+20
+          tsumSTD(i,k,2)=missing_val
          ENDIF !tnondef
 c
@@ -73 +78 @@
      $    (ecrit_hf2mth-tnondef(i,k,3))
          ELSE
-          tsumSTD(i,k,3)=1.E+20
+          tsumSTD(i,k,3)=missing_val
          ENDIF !tnondef
 c

LMDZ4/trunk/libf/phylmd/plevel.F

@@ -7 +7 @@
 c================================================================
 c================================================================
+      USE netcdf
       USE dimphy
       IMPLICIT none
@@ -54 +55 @@
 c$OMP THREADPRIVATE(first)
       INTEGER i, k
+c
+      REAL missing_val
+c
+      missing_val=nf90_fill_real
 c
       if (first) then
@@ -119 +124 @@
       do i = 1, klon
          if (pgcm(i,1).LT.pres) THEN
-c           Qpres(i)=1e33
-            Qpres(i)=1e+20
+            Qpres(i)=missing_val
          endif
       enddo

LMDZ4/trunk/libf/phylmd/plevel_new.F

@@ -7 +7 @@
 c================================================================
 c================================================================
+      USE netcdf
       USE dimphy
       IMPLICIT none
@@ -55 +56 @@
 c$OMP THREADPRIVATE(first)
       INTEGER i, k
+c
+      REAL missing_val
+c
+      missing_val=nf90_fill_real
 c
       if (first) then
@@ -120 +125 @@
          DO i=1,klon
             IF (pgcm(i,1).LT.pres(nlev)) THEN
-c     Qpres(i)=1e33
-               Qpres(i,nlev) = 1e+20
+               Qpres(i,nlev) = missing_val
             ELSE
                Qpres(i,nlev) = 

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