Changeset 2103 for LMDZ5/trunk/libf/phylmd

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

Location:

LMDZ5/trunk/libf/phylmd

Files:

• calcul_divers.h (modified) (1 diff)
• phys_local_var_mod.F90 (modified) (3 diffs)
• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_write_mod.F90 (modified) (3 diffs)

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
+!

LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -228 +228 @@
       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: wake_h, wbeff, zmax_th, zq2m, zt2m
 !$OMP THREADPRIVATE(wake_h, wbeff, zmax_th, zq2m, zt2m)
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: zt2m_min_mon, zt2m_max_mon
+!$OMP THREADPRIVATE(zt2m_min_mon, zt2m_max_mon)
+      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: t2m_min_mon, t2m_max_mon
+!$OMP THREADPRIVATE(t2m_min_mon, t2m_max_mon)
       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: weak_inversion
 !$OMP THREADPRIVATE(weak_inversion)
@@ -403 +407 @@
       ALLOCATE(rain_lsc(klon), wake_h(klon), wbeff(klon), zmax_th(klon))
       ALLOCATE(zq2m(klon), zt2m(klon), weak_inversion(klon))
+      ALLOCATE(zt2m_min_mon(klon), zt2m_max_mon(klon))
+      ALLOCATE(t2m_min_mon(klon), t2m_max_mon(klon))
       ALLOCATE(sens(klon), flwp(klon), fiwp(klon))
       ALLOCATE(ale_bl_stat(klon), alp_bl_conv(klon), alp_bl_det(klon))
@@ -554 +560 @@
       DEALLOCATE(rain_lsc, wake_h, wbeff, zmax_th)
       DEALLOCATE(zq2m, zt2m, weak_inversion)
+      DEALLOCATE(zt2m_min_mon, zt2m_max_mon)
+      DEALLOCATE(t2m_min_mon, t2m_max_mon)
       DEALLOCATE(sens, flwp, fiwp)
       DEALLOCATE(ale_bl_stat,alp_bl_conv,alp_bl_det)

LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -49 +49 @@
   TYPE(ctrl_out), SAVE :: o_t2m = ctrl_out((/ 1, 1, 1, 5, 10, 10, 11, 11, 11 /), &
     't2m', 'Temperature 2m', 'K', (/ ('', i=1, 9) /))
-  TYPE(ctrl_out), SAVE :: o_t2m_min = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
+  TYPE(ctrl_out), SAVE :: o_t2m_min = ctrl_out((/ 20, 1, 10, 10, 10, 10, 11, 11, 11 /), &
     't2m_min', 'Temp 2m min', 'K', &
       (/ "t_min(X)", "t_min(X)", "t_min(X)", "t_min(X)", "t_min(X)", "t_min(X)", "t_min(X)", "t_min(X)", "t_min(X)" /))
-  TYPE(ctrl_out), SAVE :: o_t2m_max = ctrl_out((/ 1, 1, 10, 10, 10, 10, 11, 11, 11 /), &
+  TYPE(ctrl_out), SAVE :: o_t2m_max = ctrl_out((/ 20, 1, 10, 10, 10, 10, 11, 11, 11 /), &
     't2m_max', 'Temp 2m max', 'K', &
       (/ "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)", &
          "t_max(X)", "t_max(X)", "t_max(X)", "t_max(X)" /))
+
+  TYPE(ctrl_out), SAVE :: o_t2m_min_mon = ctrl_out((/ 1, 20, 20, 20, 20, 20, 20, 20, 20 /), &
+    't2m_min_mon', 'Monthly average min 2m temperature', 'K', &
+      (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
+  TYPE(ctrl_out), SAVE :: o_t2m_max_mon = ctrl_out((/ 1, 20, 20, 20, 20, 20, 20, 20, 20 /), &
+    't2m_max_mon', 'Monthly average max 2m temperature', 'K', &
+      (/ "inst(X)", "inst(X)", "inst(X)", "inst(X)", "inst(X)", &
+         "inst(X)", "inst(X)", "inst(X)", "inst(X)" /))
 
   TYPE(ctrl_out), SAVE, DIMENSION(4) :: o_t2m_srf = (/ &

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -30 +30 @@
          o_aireTER, o_flat, o_slp, o_tsol, &
          o_t2m, o_t2m_min, o_t2m_max, &
+         o_t2m_min_mon, o_t2m_max_mon, &
          o_q2m, o_ustar, o_u10m, o_v10m, &
          o_wind10m, o_wind10max, o_sicf, &
@@ -172 +173 @@
 
     USE phys_local_var_mod, only: zxfluxlat, slp, zxtsol, zt2m, &
+         t2m_min_mon, t2m_max_mon, &
          zu10m, zv10m, zq2m, zustar, zxqsurf, qsol, &
          rain_lsc, snow_lsc, evap, bils, sens, fder, &
@@ -327 +329 @@
        CALL histwrite_phy(o_t2m_min, zt2m)
        CALL histwrite_phy(o_t2m_max, zt2m)
+       CALL histwrite_phy(o_t2m_max_mon, t2m_max_mon)
+       CALL histwrite_phy(o_t2m_min_mon, t2m_min_mon)
 
        IF (vars_defined) THEN