Changeset 1825 for LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

Timestamp:

Aug 2, 2013, 4:36:53 PM (11 years ago)

Author:

Ehouarn Millour

Message:

Première étape de l'implémentation de XIOS. Modifications isolées dans des flags CPP_XIOS. Sorties opérationnelles (sauf stations et régionalisation) en modes séquentiel et omp, pas mpi.
UG
...........................................
First step of the XIOS implementation. Modifications are confined into CPP_XIOS flags. Output is operationnal (except for stations and regionalization) in sequential and omp modes (not mpi).
UG

File:

• LMDZ5/trunk/libf/phylmd/phys_output_mod.F90 (modified) (7 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

@@ -42 +42 @@
     USE aero_mod, only : naero_spc,name_aero
     USE phys_output_ctrlout_mod
+
+#ifdef CPP_XIOS
+    ! ug Pour les sorties XIOS
+    USE wxios
+#endif
 
     IMPLICIT NONE
@@ -98 +103 @@
     CHARACTER(LEN=3)                      :: ctetaSTD(nbteta)
     REAL, DIMENSION(nfiles)               :: ecrit_files
-    CHARACTER(LEN=20), DIMENSION(nfiles)  :: phys_out_filenames
     INTEGER, DIMENSION(iim*jjmp1)         ::  ndex2d
     INTEGER, DIMENSION(iim*jjmp1*klev)    :: ndex3d
@@ -117 +121 @@
     REAL, DIMENSION(nfiles), SAVE     :: phys_out_latmin        = (/    -90.,    -90.,    -90.,     -90.,    -90.,    -90. /)
     REAL, DIMENSION(nfiles), SAVE     :: phys_out_latmax        = (/     90.,     90.,     90.,     90.,     90.,     90. /)
+
+#ifdef CPP_XIOS
+    ! ug Variables utilisées pour récupérer le calendrier pour xios
+    INTEGER :: x_an, x_mois, x_jour
+    REAL :: x_heure
+#endif
 
     WRITE(lunout,*) 'Debut phys_output_mod.F90'
@@ -209 +219 @@
     WRITE(lunout,*)'phys_out_filelevels=',lev_files
 
+#ifdef CPP_XIOS
+    ! ug Réglage du calendrier xios
+    !Temps julian => an, mois, jour, heure
+    CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
+    CALL ju2ymds(zjulian, x_an, x_mois, x_jour, x_heure)
+    CALL wxios_set_cal(dtime, calend, x_an, x_mois, x_jour, x_heure)
+#endif
+
 !!!!!!!!!!!!!!!!!!!!!!! Boucle sur les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     ! Appel de histbeg et histvert pour creer le fichier et les niveaux verticaux !!
@@ -257 +275 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
           IF (phys_out_regfkey(iff)) then
-
              imin_ins=1
              imax_ins=iim
@@ -297 +314 @@
                   itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
           else
-             CALL histbeg_phy(phys_out_filenames(iff),itau_phy,zjulian,dtime,nhorim(iff),nid_files(iff))
+             CALL histbeg_phy_all(phys_out_filenames(iff),itau_phy,zjulian,&
+                 dtime,chtimestep(iff),lev_files(iff),nhorim(iff),nid_files(iff))
           endif
 
@@ -329 +347 @@
           !                 1,preff,nvertp0(iff))
 
+#ifdef CPP_XIOS
+    ! ug déclaration des axes verticaux de chaque fichier:
+    CALL wxios_add_vaxis("presnivs", nid_files(iff), levmax(iff) - levmin(iff) + 1, presnivs(levmin(iff):levmax(iff)))
+    CALL wxios_add_vaxis("Ahyb", nid_files(iff), levmax(iff) - levmin(iff) + 1, Ahyb)
+    CALL wxios_add_vaxis("Bhyb", nid_files(iff), levmax(iff) - levmin(iff) + 1, Bhyb)
+    CALL wxios_add_vaxis("Ahyb", nid_files(iff), levmax(iff) - levmin(iff) + 1, Alt)
+#endif 
 
       IF (nqtot>=3) THEN