Ignore:
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:
4 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

    r2333 r2385  
    239239
    240240!Ajout de celles nécessaires au phys_output_write_mod
    241       REAL, SAVE, ALLOCATABLE :: slp(:)
    242 !$OMP THREADPRIVATE(slp)
     241      REAL, SAVE, ALLOCATABLE :: tal1(:), pal1(:), pab1(:), pab2(:)
     242!$OMP THREADPRIVATE(tal1, pal1, pab1, pab2)
     243      REAL, SAVE, ALLOCATABLE :: ptstar(:), pt0(:), slp(:)
     244!$OMP THREADPRIVATE(ptstar, pt0, slp)
    243245      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: sens, flwp, fiwp
    244246!$OMP THREADPRIVATE(sens, flwp, fiwp)
     
    507509! FH Ajout de celles nécessaires au phys_output_write_mod
    508510
    509       ALLOCATE(slp(klon))
     511      ALLOCATE(tal1(klon), pal1(klon), pab1(klon), pab2(klon))
     512      ALLOCATE(ptstar(klon),pt0(klon),slp(klon))
    510513      ALLOCATE(ale_wake(klon), alp_wake(klon), bils(klon))
    511514      ALLOCATE(cdragm(klon), cdragh(klon), cldh(klon), cldl(klon))
     
    709712
    710713! FH Ajout de celles nécessaires au phys_output_write_mod
    711       DEALLOCATE(slp)
     714      DEALLOCATE(tal1, pal1, pab1, pab2)
     715      DEALLOCATE(ptstar, pt0, slp)
    712716      DEALLOCATE(ale_wake, alp_wake, bils)
    713717      DEALLOCATE(cdragm, cdragh, cldh, cldl)
  • LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

    r2333 r2385  
    4343  TYPE(ctrl_out), SAVE :: o_flat = ctrl_out((/ 5, 1, 10, 10, 5, 10, 11, 11, 11 /), &
    4444    'flat', 'Latent heat flux', 'W/m2', (/ ('', i=1, 9) /))
     45  TYPE(ctrl_out), SAVE :: o_ptstar = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
     46    'ptstar', 'Surface Temperature', 'K', (/ ('', i=1, 9) /))
     47  TYPE(ctrl_out), SAVE :: o_pt0 = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
     48    'pt0', 'Standard surface temperature', 'K', (/ ('', i=1, 9) /))
    4549  TYPE(ctrl_out), SAVE :: o_slp = ctrl_out((/ 1, 1, 1, 10, 10, 10, 11, 11, 11 /), &
    4650    'slp', 'Sea Level Pressure', 'Pa', (/ ('', i=1, 9) /))
  • LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

    r2377 r2385  
    2929    USE phys_output_ctrlout_mod, only: o_phis, o_aire, is_ter, is_lic, is_oce, &
    3030         is_ave, is_sic, o_contfracATM, o_contfracOR, &
    31          o_aireTER, o_flat, o_slp, o_tsol, &
     31         o_aireTER, o_flat, o_slp, o_ptstar, o_pt0, o_tsol, &
    3232         o_t2m, o_t2m_min, o_t2m_max, &
    3333         o_t2m_min_mon, o_t2m_max_mon, &
     
    183183         vphiSTD, wTSTD, u2STD, v2STD, T2STD, missing_val_nf90
    184184
    185     USE phys_local_var_mod, only: zxfluxlat, slp, zxtsol, zt2m, &
     185    USE phys_local_var_mod, only: zxfluxlat, slp, ptstar, pt0, zxtsol, zt2m, &
    186186         t2m_min_mon, t2m_max_mon, evap, &
    187187         zu10m, zv10m, zq2m, zustar, zxqsurf, &
     
    350350!!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    351351       CALL histwrite_phy(o_flat, zxfluxlat)
     352       CALL histwrite_phy(o_ptstar, ptstar)
     353       CALL histwrite_phy(o_pt0, pt0)
    352354       CALL histwrite_phy(o_slp, slp)
    353355       CALL histwrite_phy(o_tsol, zxtsol)
  • LMDZ5/trunk/libf/phylmd/physiq.F90

    r2372 r2385  
    35023502  cldtaupirad = cldtaupi
    35033503  cldemirad   = cldemi
     3504  cldfrarad   = cldfra
    35043505
    35053506  !
     
    41864187  include "calcul_STDlev.h"
    41874188  !
    4188   ! slp sea level pressure
    4189   slp(:) = paprs(:,1)*exp(pphis(:)/(RD*t_seri(:,1)))
     4189  ! slp sea level pressure derived from Arpege-IFS : CALL ctstar + CALL pppmer
     4190  CALL diag_slp(klon,t_seri,paprs,pplay,pphis,ptstar,pt0,slp)
    41904191  !
    41914192  !cc prw = eau precipitable
Note: See TracChangeset for help on using the changeset viewer.