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

Timestamp:

Dec 5, 2013, 6:32:35 PM (10 years ago)

Author:

musat

Message:

1) Modifications pour faire des simulations par an avec un calendrier realiste (365 jours ou autre).

Il faut mettre une frequence de sortie de -1 (variable phys_out_filetimesteps dans config.def) pour
que IOIPSL calcule les moyennes mensuels en prenant en compte des longuers variables de chaque
mois. Par exemple, pour le fichier histmth (1er fichier) et histmthNMC (7eme).
phys_out_filetimesteps= -1 1day 6hr 6hr 6hr 1d -1 1day 6hr

2) Corrections titres variables niveaux de pression des fichiers histmth, histday, etc

File:

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

LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

@@ -44 +44 @@
     USE surface_data, ONLY : ok_snow
     USE phys_output_ctrlout_mod
+    USE mod_grid_phy_lmdz, only: klon_glo
 
 #ifdef CPP_XIOS
@@ -207 +208 @@
     ecrit_files(5) = ecrit_LES
     ecrit_files(6) = ecrit_ins
-    ecrit_files(7) = freq_outNMC(1)/86400.
-    ecrit_files(8) = freq_outNMC(2)/86400.
-    ecrit_files(9) = freq_outNMC(3)/86400.
+    ecrit_files(7) = freq_outNMC(1)
+    ecrit_files(8) = freq_outNMC(2)
+    ecrit_files(9) = freq_outNMC(3)
 
     !! Lectures des parametres de sorties dans physiq.def
@@ -248 +249 @@
     ! ug Réglage du calendrier xios
     !Temps julian => an, mois, jour, heure
-    CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
+    IF (klon_glo==1) THEN                                       
+        CALL ymds2ju(annee_ref, 1, day_ref, hour, zjulian)
+    ELSE
+        CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
+    END IF 
     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)
@@ -278 +283 @@
           ! ...)*86400.
           ecrit_files(iff)=ecrit_files(iff)*86400.
+      ELSE IF (chtimestep(iff).eq.'-1') then
+          PRINT*,'ecrit_files(',iff,') < 0 so IOIPSL work on different'
+          PRINT*,'months length'
+          ecrit_files(iff)=-1.
       else
        CALL convers_timesteps(chtimestep(iff),dtime,ecrit_files(iff)) 
@@ -304 +313 @@
 
           idayref = day_ref
-!          CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)       
-! correction pour l heure initiale                               !jyg 
-!                                                                !jyg 
-          CALL ymds2ju(annee_ref, 1, idayref, hour, zjulian)         !jyg
-! correction pour l heure initiale                               !jyg 
-!                                                                !jyg 
-!!!      CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)       !jyg 
-! correction pour l heure initiale                               !jyg 
-!                                                                !jyg 
-!      CALL ymds2ju(annee_ref, 1, idayref, hour, zjulian)         !jyg
-
+! A noter pour
+! l heure initiale - dans les fichiers histoire hist* - on met comme  
+! heure de debut soit la vraie heure (pour le 1D) soit 0h (pour le 3D) 
+! afin d avoir une seule sortie mensuelle par mois lorsque l on tourne 
+! par annee (IM).
+!
+     IF (klon_glo==1) THEN                                       
+         CALL ymds2ju(annee_ref, 1, idayref, hour, zjulian)     !jyg
+     ELSE
+         CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
+     END IF                                                      
 !!!!!!!!!!!!!!!!! Traitement dans le cas ou l'on veut stocker sur un domaine limite !!
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -479 +488 @@
     ecrit_LES = ecrit_files(5)
     ecrit_ins = ecrit_files(6)
-
+    freq_outNMC(1) = ecrit_files(7)
+    freq_outNMC(2) = ecrit_files(8)
+    freq_outNMC(3) = ecrit_files(9)
+
+    WRITE(lunout,*)'freq_outNMC(1)=',freq_outNMC(1)
+    WRITE(lunout,*)'freq_outNMC(2)=',freq_outNMC(2)
+    WRITE(lunout,*)'freq_outNMC(3)=',freq_outNMC(3)
     WRITE(lunout,*)'swaero_diag=',swaero_diag
     WRITE(lunout,*)'Fin phys_output_mod.F90'