Changeset 2271 for LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

May 12, 2015, 6:27:47 PM (9 years ago)

Author:

musat

Message:

Corrections for standard level pressure outputs hist*NMC.nc by XIOS and IOIPSL.
For XIOS one need to specify a default_value="1.e+20" in the field_def_lmdz.xml and
a detect_missing_value=".true." in the file_def_histins_lmdz.xml,
file_def_histday_lmdz.xml, etc files.

File:

• LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -174 +174 @@
          vqsumSTD, vTsumSTD, O3daysumSTD, wqsumSTD, &
          vphisumSTD, wTsumSTD, u2sumSTD, v2sumSTD, &
-         T2sumSTD, nlevSTD, du_gwd_rando, dv_gwd_rando
+         T2sumSTD, nlevSTD, du_gwd_rando, dv_gwd_rando, &
+         ulevSTD, vlevSTD, wlevSTD, philevSTD, qlevSTD, tlevSTD, &
+         missing_val_nf90
 
     USE phys_local_var_mod, only: zxfluxlat, slp, zxtsol, zt2m, &
@@ -244 +246 @@
     ! ug Pour les sorties XIOS
     USE xios, ONLY: xios_update_calendar
-    USE wxios, only: wxios_closedef
+    USE wxios, only: wxios_closedef, missing_val
 #endif
     USE phys_cal_mod, only : mth_len
@@ -290 +292 @@
     INTEGER, DIMENSION(iim*jjmp1*klev) :: ndex3d
     REAL, PARAMETER :: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
-    REAL, PARAMETER :: missing_val=nf90_fill_real
+!   REAL, PARAMETER :: missing_val=nf90_fill_real
+#ifndef CPP_XIOS
+    REAL :: missing_val
+#endif
     REAL, PARAMETER :: un_jour=86400.
 
@@ -658 +663 @@
        !       ENDIF
 
+  IF (.NOT.ok_all_xml) THEN
        ! ATTENTION, LES ANCIENS HISTWRITE ONT ETES CONSERVES EN ATTENDANT MIEUX:
        ! Champs interpolles sur des niveaux de pression
+#ifndef CPP_XIOS
+       missing_val=missing_val_nf90
+#endif
        DO iff=1, nfiles
           ll=0
@@ -681 +690 @@
           ENDDO
        ENDDO
-
+  ELSE
+!XIOS  CALL xios_get_field_attr("u850",default_value=missing_val)
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+       DO iff=1, nfiles
+          ll=0
+          DO k=1, nlevSTD
+             bb2=clevSTD(k) 
+             IF(bb2.EQ."850".OR.bb2.EQ."700".OR. &
+                bb2.EQ."500".OR.bb2.EQ."200".OR. &
+                bb2.EQ."100".OR. &
+                bb2.EQ."50".OR.bb2.EQ."10") THEN
+                ll=ll+1
+                CALL histwrite_phy(o_uSTDlevs(ll),ulevSTD(:,k), iff)
+                CALL histwrite_phy(o_vSTDlevs(ll),vlevSTD(:,k), iff)
+                CALL histwrite_phy(o_wSTDlevs(ll),wlevSTD(:,k), iff)
+                CALL histwrite_phy(o_zSTDlevs(ll),philevSTD(:,k), iff)
+                CALL histwrite_phy(o_qSTDlevs(ll),qlevSTD(:,k), iff)
+                CALL histwrite_phy(o_tSTDlevs(ll),tlevSTD(:,k), iff)
+             ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
+          ENDDO
+       ENDDO
+  ENDIF
        IF (vars_defined) THEN
           DO i=1, klon