Changeset 1234 for LMDZ4/branches

Timestamp:

Aug 28, 2009, 11:26:56 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/branches/LMDZ4-dev/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/branches/LMDZ4-dev/libf/phylmd/moy_undefSTD.F

@@ -5 +5 @@
      $           dtime,ecrit_day,ecrit_mth,ecrit_hf2mth,
      $           tnondef,tsumSTD)
+      USE netcdf
       USE dimphy
       IMPLICIT none
@@ -29 +30 @@
 cym#include "dimphy.h"
 c
+c
 c variables Input
       INTEGER nlevSTD, klevSTD, itap
@@ -42 +44 @@
       REAL tsumSTD(klon,klevSTD,nout)
 c
+      REAL missing_val
+c
+      missing_val=nf90_fill_real
+c
 c calcul 1 fois par jour
 c
@@ -51 +57 @@
      $    (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 +72 @@
      $    (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 +79 @@
      $    (ecrit_hf2mth-tnondef(i,k,3))
          ELSE
-          tsumSTD(i,k,3)=1.E+20
+          tsumSTD(i,k,3)=missing_val
          ENDIF !tnondef
 c

LMDZ4/branches/LMDZ4-dev/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/branches/LMDZ4-dev/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/branches/LMDZ4-dev/libf/phylmd/undefSTD.F

@@ -5 +5 @@
      $           dtime,ecrit_hf,
      $           oknondef,tnondef,tsumSTD)
+      USE netcdf
       USE dimphy
       IMPLICIT none
@@ -32 +33 @@
 cym#include "dimphy.h"
 c variables Input
+c
       INTEGER nlevSTD, klevSTD, itap
       PARAMETER(klevSTD=17)
@@ -45 +47 @@
       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 +64 @@
          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 +77 @@
        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