Changeset 3110 for LMDZ6/trunk/libf


Ignore:
Timestamp:
Dec 4, 2017, 5:19:33 PM (7 years ago)
Author:
Laurent Fairhead
Message:

Added sza, solar zenithal angle, to outputs for CMIP6 dataresquest

Location:
LMDZ6/trunk/libf/phylmd
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

    r3106 r3110  
    5151
    5252!!! 2D
     53  TYPE(ctrl_out), SAVE :: o_sza = ctrl_out((/ 1, 1, 10, 10, 5, 10, 11, 11, 11, 11/), &
     54    'sza', 'Solar zenithal angle', 'degrees', (/ ('', i=1, 10) /))
    5355
    5456! Marine
  • LMDZ6/trunk/libf/phylmd/phys_output_var_mod.F90

    r3082 r3110  
    1414  REAL, SAVE, ALLOCATABLE :: snow_o(:), zfra_o(:)
    1515  !$OMP THREADPRIVATE(snow_o, zfra_o)
     16  REAL, SAVE, ALLOCATABLE :: sza_o(:) ! solar zenithal angle
     17  !$OMP THREADPRIVATE(sza_o)
    1618  INTEGER, SAVE, ALLOCATABLE ::  itau_con(:)       ! Nombre de pas ou rflag <= 1
    1719  !$OMP THREADPRIVATE(itau_con)
     
    144146
    145147    allocate(snow_o(klon), zfra_o(klon))
     148    allocate(sza_o(klon) )
    146149    allocate(itau_con(klon))
    147150    allocate(sens_prec_liq_o(klon,2))
     
    198201
    199202    deallocate(snow_o,zfra_o,itau_con)
     203    deallocate(sza_o)
    200204    deallocate (bils_ec,bils_ech,bils_tke,bils_diss,bils_kinetic,bils_enthalp,bils_latent)
    201205    deallocate (d_qw_col, d_ql_col, d_qs_col, d_qt_col, d_ek_col, d_h_dair_col &
  • LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

    r3106 r3110  
    178178         o_lat_prec_liq_oce, o_lat_prec_liq_sic, &
    179179         o_lat_prec_sol_oce, o_lat_prec_sol_sic, &
     180         o_sza, &
    180181! Marine
    181182         o_map_prop_hc, o_map_prop_hist, o_map_emis_hc, o_map_iwp_hc, &
     
    319320         itau_con, nfiles, clef_files, nid_files, dryaod_diag, &
    320321         zustr_gwd_hines, zvstr_gwd_hines,zustr_gwd_rando, zvstr_gwd_rando, &
    321          zustr_gwd_front, zvstr_gwd_front,     &
     322         zustr_gwd_front, zvstr_gwd_front, sza_o,    &
    322323         sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o, &
    323324         cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, &
     
    614615     ENDIF
    615616
     617       CALL histwrite_phy(o_sza, sza_o)
    616618       CALL histwrite_phy(o_flat, zxfluxlat)
    617619       CALL histwrite_phy(o_ptstar, ptstar)
  • LMDZ6/trunk/libf/phylmd/physiq_mod.F90

    r3107 r3110  
    11221122    logical, parameter :: mass_fixer=.false.
    11231123    real qql1(klon),qql2(klon),corrqql
     1124
     1125    REAL pi
     1126
     1127    pi = 4. * ATAN(1.)
    11241128
    11251129    ! Ehouarn: set value of jjmp1 since it is no longer a "fixed parameter"
     
    21902194       END SELECT
    21912195    ENDIF
     2196    sza_o = ACOS (rmu0) *180./pi
    21922197
    21932198    IF (mydebug) THEN
Note: See TracChangeset for help on using the changeset viewer.