Changeset 2385

Timestamp:

Nov 6, 2015, 6:30:34 PM (9 years ago)

Author:

musat

Message:

Correction calcul sea level pressure (slp) cf. Arpege-IFS
Ajout initialisation cldfrarad

Location:

LMDZ5/trunk/libf/phylmd

Files:

• phys_local_var_mod.F90 (modified) (3 diffs)
• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_write_mod.F90 (modified) (3 diffs)
• physiq.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -239 +239 @@
 
 !Ajout de celles nécessaires au phys_output_write_mod
-      REAL, SAVE, ALLOCATABLE :: slp(:)
-!$OMP THREADPRIVATE(slp)
+      REAL, SAVE, ALLOCATABLE :: tal1(:), pal1(:), pab1(:), pab2(:)
+!$OMP THREADPRIVATE(tal1, pal1, pab1, pab2)
+      REAL, SAVE, ALLOCATABLE :: ptstar(:), pt0(:), slp(:)
+!$OMP THREADPRIVATE(ptstar, pt0, slp)
       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: sens, flwp, fiwp 
 !$OMP THREADPRIVATE(sens, flwp, fiwp)
@@ -507 +509 @@
 ! FH Ajout de celles nécessaires au phys_output_write_mod
 
-      ALLOCATE(slp(klon))
+      ALLOCATE(tal1(klon), pal1(klon), pab1(klon), pab2(klon))
+      ALLOCATE(ptstar(klon),pt0(klon),slp(klon))
       ALLOCATE(ale_wake(klon), alp_wake(klon), bils(klon))
       ALLOCATE(cdragm(klon), cdragh(klon), cldh(klon), cldl(klon))
@@ -709 +712 @@
 
 ! FH Ajout de celles nécessaires au phys_output_write_mod
-      DEALLOCATE(slp)
+      DEALLOCATE(tal1, pal1, pab1, pab2)
+      DEALLOCATE(ptstar, pt0, slp)
       DEALLOCATE(ale_wake, alp_wake, bils)
       DEALLOCATE(cdragm, cdragh, cldh, cldl)

LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -43 +43 @@
   TYPE(ctrl_out), SAVE :: o_flat = ctrl_out((/ 5, 1, 10, 10, 5, 10, 11, 11, 11 /), &
     'flat', 'Latent heat flux', 'W/m2', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_ptstar = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
+    'ptstar', 'Surface Temperature', 'K', (/ ('', i=1, 9) /))
+  TYPE(ctrl_out), SAVE :: o_pt0 = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
+    'pt0', 'Standard surface temperature', 'K', (/ ('', i=1, 9) /))
   TYPE(ctrl_out), SAVE :: o_slp = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
     'slp', 'Sea Level Pressure', 'Pa', (/ ('', i=1, 9) /))

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -29 +29 @@
     USE phys_output_ctrlout_mod, only: o_phis, o_aire, is_ter, is_lic, is_oce, &
          is_ave, is_sic, o_contfracATM, o_contfracOR, &
-         o_aireTER, o_flat, o_slp, o_tsol, &
+         o_aireTER, o_flat, o_slp, o_ptstar, o_pt0, o_tsol, &
          o_t2m, o_t2m_min, o_t2m_max, &
          o_t2m_min_mon, o_t2m_max_mon, &
@@ -183 +183 @@
          vphiSTD, wTSTD, u2STD, v2STD, T2STD, missing_val_nf90
 
-    USE phys_local_var_mod, only: zxfluxlat, slp, zxtsol, zt2m, &
+    USE phys_local_var_mod, only: zxfluxlat, slp, ptstar, pt0, zxtsol, zt2m, &
          t2m_min_mon, t2m_max_mon, evap, &
          zu10m, zv10m, zq2m, zustar, zxqsurf, &
@@ -350 +350 @@
 !!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
        CALL histwrite_phy(o_flat, zxfluxlat)
+       CALL histwrite_phy(o_ptstar, ptstar)
+       CALL histwrite_phy(o_pt0, pt0)
        CALL histwrite_phy(o_slp, slp)
        CALL histwrite_phy(o_tsol, zxtsol)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -3502 +3502 @@
   cldtaupirad = cldtaupi
   cldemirad   = cldemi
+  cldfrarad   = cldfra
 
   !
@@ -4186 +4187 @@
   include "calcul_STDlev.h"
   !
-  ! slp sea level pressure
-  slp(:) = paprs(:,1)*exp(pphis(:)/(RD*t_seri(:,1)))
+  ! slp sea level pressure derived from Arpege-IFS : CALL ctstar + CALL pppmer
+  CALL diag_slp(klon,t_seri,paprs,pplay,pphis,ptstar,pt0,slp)
   !
   !cc prw = eau precipitable