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

Timestamp:

Aug 6, 2013, 6:06:40 PM (11 years ago)

Author:

idelkadi

Message:

Gestion des fichiers de sorties sur les niveaux de pression NMC (hist*NMC.nc) via le meme mecanisme utilise pour les fichiers histoires (phys_output_mod, ...)
Ces sorties sont ainsi controles par les memes flags : phys_out_filekeys, phys_out_filenames, phys_out_filetimesteps, phys_out_filelevels, ...
(exemple dans config.def : phys_out_filenames=histmth histday histhf histins histLES histstn histmthNMC histdayNMC histhfNMC)

File:

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

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -20 +20 @@
    &                  pplay, lmax_th, aerosol_couple,         &
    &                  ok_ade, ok_aie, ivap, new_aod, ok_sync, &
-   &                  ptconv, read_climoz, clevSTD, ptconvth, &
+   &                  ptconv, read_climoz, clevSTD, freq_moyNMC, ptconvth, &
    &                  d_t, qx, d_qx, zmasse, flag_aerosol_strat)
 
@@ -40 +40 @@
     USE iophy
     USE mod_phys_lmdz_para
+    USE netcdf,  ONLY : nf90_fill_real
 
 #ifdef CPP_XIOS
@@ -70 +71 @@
     REAL, DIMENSION(klon, llm) :: zmasse
     LOGICAL :: flag_aerosol_strat
+    REAL, DIMENSION(3) :: freq_moyNMC
 
 ! Local
@@ -82 +84 @@
     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
 
      ! On calcul le nouveau tau:
@@ -974 +977 @@
       CALL histwrite_phy(o_rldcs4co2, lwdn0p)
       ENDIF
+!!!!!!!!!!!! Sorties niveaux de pression NMC !!!!!!!!!!!!!!!!!!!!
+DO iff=7, nfiles
+
+      CALL histwrite_phy(o_tnondef,tnondef(:,:,iff-6),iff)
+      CALL histwrite_phy(o_ta,twriteSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_zg,phiwriteSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_hus,qwriteSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_hur,rhwriteSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_ua,uwriteSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_va,vwriteSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_wap,wwriteSTD(:,:,iff-6),iff)
+IF(vars_defined) THEN
+       DO k=1, nlevSTD
+        DO i=1, klon
+         IF(tnondef(i,k,iff-6).NE.missing_val) THEN
+          zx_tmp_fi3d(i,k) = (100.*tnondef(i,k,iff-6))/freq_moyNMC(iff-6)
+         ELSE
+          zx_tmp_fi3d(i,k) = missing_val
+         ENDIF
+        ENDDO
+      ENDDO
+ENDIF
+      CALL histwrite_phy(o_psbg,zx_tmp_fi3d,iff)
+IF(vars_defined) THEN
+       DO k=1, nlevSTD
+        DO i=1, klon
+         IF(O3sumSTD(i,k,iff-6).NE.missing_val) THEN
+          zx_tmp_fi3d(i,k) = O3sumSTD(i,k,iff-6) * 1.e+9
+         ELSE
+          zx_tmp_fi3d(i,k) = missing_val
+         ENDIF
+        ENDDO
+       ENDDO !k=1, nlevSTD
+ENDIF
+      CALL histwrite_phy(o_tro3,zx_tmp_fi3d,iff)
+       if (read_climoz == 2) THEN
+IF(vars_defined) THEN
+       DO k=1, nlevSTD
+        DO i=1, klon
+         IF(O3daysumSTD(i,k,iff-6).NE.missing_val) THEN
+          zx_tmp_fi3d(i,k) = O3daysumSTD(i,k,iff-6) * 1.e+9
+         ELSE
+          zx_tmp_fi3d(i,k) = missing_val
+         ENDIF
+        ENDDO
+       ENDDO !k=1, nlevSTD
+ENDIF
+      CALL histwrite_phy(o_tro3_daylight,zx_tmp_fi3d,iff)
+      endif
+      CALL histwrite_phy(o_uxv,uvsumSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_vxq,vqsumSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_vxT,vTsumSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_wxq,wqsumSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_vxphi,vphisumSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_wxT,wTsumSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_uxu,u2sumSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_vxv,v2sumSTD(:,:,iff-6),iff)
+      CALL histwrite_phy(o_TxT,T2sumSTD(:,:,iff-6),iff)
+ENDDO !nfiles
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
         IF (nqtot.GE.3) THEN
          DO iq=3,nqtot