Changeset 2103 for LMDZ5/trunk/libf/phylmd/calcul_divers.h

Timestamp:

Jul 31, 2014, 3:46:23 PM (10 years ago)

Author:

musat

Message:

Ajout variables t2m_min_mon et t2m_max_mon dans les fichier
histmth. Ce sont les tempe©rature a minimales et
maximales journalie¨re) moyeneées sur le mo. Le calcul
est fait "manuelleÃment" dans calcul_divers.h.
Teste en sequentiel (en local) avec calendrier a 360 et 365 jours
et en parallele MPI avec calendrier 360 jours (sur adapp).
IM

• This line, and those below, will be ignored--

M phys_output_ctrlout_mod.F90
M phys_local_var_mod.F90
M phys_output_write_mod.F90
M calcul_divers.h

File:

• LMDZ5/trunk/libf/phylmd/calcul_divers.h (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/calcul_divers.h

@@ -33 +33 @@
 !       print*,'initialisation itapm1 ',itapm1
       ENDIF
+!
+! calcul temperatures minimale et maximale moyennees sur le mois 
+!
+!initialisation debut de mois ou de journee pour les fichiers mensuels
+  IF(itap.EQ.itapm1+1) THEN
+     t2m_min_mon=0.
+     t2m_max_mon=0.
+  ENDIF
+  IF(MOD(itap,NINT(un_jour/dtime)).EQ.1) THEN
+     zt2m_min_mon=zt2m
+     zt2m_max_mon=zt2m
+  ENDIF
+!calcul a chaque pas de temps pour les fichiers mensuels
+     DO i = 1, klon
+        zt2m_min_mon(i)=MIN(zt2m(i),zt2m_min_mon(i))
+        zt2m_max_mon(i)=MAX(zt2m(i),zt2m_max_mon(i))
+     ENDDO
+!fin de journee
+  IF(MOD(itap,NINT(un_jour/dtime)).EQ.0) THEN
+   t2m_min_mon=t2m_min_mon+zt2m_min_mon 
+   t2m_max_mon=t2m_max_mon+zt2m_max_mon 
+  ENDIF
+!fin mois
+  IF(itap==itapm1) THEN
+   t2m_min_mon=t2m_min_mon/mth_len
+   t2m_max_mon=t2m_max_mon/mth_len
+  ENDIF
+!