Ignore:
Timestamp:
Feb 6, 2015, 6:53:57 PM (10 years ago)
Author:
jescribano
Message:

Dust emission scheme changes: (1) Included possibility of use previous dust emission scheme (with 2 dust bins). (2) Parameter of Marticorena and Bergametti 1995 set to it's original value (2.61) for testing purposes with pdtphys=15min. (3) Cleaning ustar calculation.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/spla_output_write.h

    r2176 r2196  
    11!! Output SPLA
    22! JE20141223
    3       CALL histwrite_phy(o_m1dflux,m1dflux)
    4       CALL histwrite_phy(o_m2dflux,m2dflux)
    5       CALL histwrite_phy(o_m3dflux,m3dflux)
     3      IF(ok_chimeredust) CALL histwrite_phy(o_m1dflux,m1dflux)
     4      IF(ok_chimeredust) CALL histwrite_phy(o_m2dflux,m2dflux)
     5      IF(ok_chimeredust) CALL histwrite_phy(o_m3dflux,m3dflux)
    66! traceur
    77
     
    99      CALL histwrite_phy( o_taue670          , diag_aod670_tot   )
    1010      CALL histwrite_phy( o_taue865          , diag_aod865_tot  )
    11       CALL histwrite_phy( o_taue550_tr2      , diff_aod550_tr2  )
    12       CALL histwrite_phy( o_taue670_tr2      , diag_aod670_tr2   )
    13       CALL histwrite_phy( o_taue865_tr2      , diag_aod865_tr2  )
    14       CALL histwrite_phy( o_taue550_ss       , diag_aod550_ss       )
    15       CALL histwrite_phy( o_taue670_ss       , diag_aod670_ss   )
    16       CALL histwrite_phy( o_taue865_ss       , diag_aod865_ss    )
    17       CALL histwrite_phy( o_taue550_dust     , diag_aod550_dust   )
    18       CALL histwrite_phy( o_taue670_dust     , diag_aod670_dust )
    19       CALL histwrite_phy( o_taue865_dust     , diag_aod865_dust )
    20       CALL histwrite_phy( o_taue550_dustsco  , diag_aod550_dustsco )
    21       CALL histwrite_phy( o_taue670_dustsco  , diag_aod670_dustsco )
    22       CALL histwrite_phy( o_taue865_dustsco  , diag_aod865_dustsco )
     11      IF(id_fine>0) CALL histwrite_phy( o_taue550_tr2      , diff_aod550_tr2  )
     12      IF(id_fine>0) CALL histwrite_phy( o_taue670_tr2      , diag_aod670_tr2   )
     13      IF(id_fine>0) CALL histwrite_phy( o_taue865_tr2      , diag_aod865_tr2  )
     14      IF(id_coss>0) CALL histwrite_phy( o_taue550_ss       , diag_aod550_ss       )
     15      IF(id_coss>0) CALL histwrite_phy( o_taue670_ss       , diag_aod670_ss   )
     16      IF(id_coss>0) CALL histwrite_phy( o_taue865_ss       , diag_aod865_ss    )
     17      IF(id_codu>0) CALL histwrite_phy( o_taue550_dust     , diag_aod550_dust   )
     18      IF(id_codu>0) CALL histwrite_phy( o_taue670_dust     , diag_aod670_dust )
     19      IF(id_codu>0) CALL histwrite_phy( o_taue865_dust     , diag_aod865_dust )
     20      IF(id_scdu>0) CALL histwrite_phy( o_taue550_dustsco  , diag_aod550_dustsco )
     21      IF(id_scdu>0) CALL histwrite_phy( o_taue670_dustsco  , diag_aod670_dustsco )
     22      IF(id_scdu>0) CALL histwrite_phy( o_taue865_dustsco  , diag_aod865_dustsco )
    2323      CALL histwrite_phy( o_taue550_aqua     , aod550_aqua  )
    2424      CALL histwrite_phy( o_taue550_terra    , aod550_terra    )
     
    2727      CALL histwrite_phy( o_taue865_aqua     , aod865_aqua    )
    2828      CALL histwrite_phy( o_taue865_terra    , aod865_terra    )
    29       CALL histwrite_phy( o_trm01 , trm01 )
    30       CALL histwrite_phy( o_trm02 , trm02 )
    31       CALL histwrite_phy( o_trm03 , trm03 )
    32       CALL histwrite_phy( o_trm04 , trm04 )
    33       CALL histwrite_phy( o_trm05 , trm05 )
    34       CALL histwrite_phy( o_sconc01 , sconc01 )
    35       CALL histwrite_phy( o_sconc02 , sconc02 )
    36       CALL histwrite_phy( o_sconc03 , sconc03 )
    37       CALL histwrite_phy( o_sconc04 , sconc04 )
    38       CALL histwrite_phy( o_sconc05 , sconc05 )
     29      IF(id_prec>0) CALL histwrite_phy( o_trm01 , trm01 )
     30      IF(id_fine>0) CALL histwrite_phy( o_trm02 , trm02 )
     31      IF(id_coss>0) CALL histwrite_phy( o_trm03 , trm03 )
     32      IF(id_codu>0) CALL histwrite_phy( o_trm04 , trm04 )
     33      IF(id_scdu>0) CALL histwrite_phy( o_trm05 , trm05 )
     34      IF(id_prec>0) CALL histwrite_phy( o_sconc01 , sconc01 )
     35      IF(id_fine>0) CALL histwrite_phy( o_sconc02 , sconc02 )
     36      IF(id_coss>0) CALL histwrite_phy( o_sconc03 , sconc03 )
     37      IF(id_codu>0) CALL histwrite_phy( o_sconc04 , sconc04 )
     38      IF(id_scdu>0) CALL histwrite_phy( o_sconc05 , sconc05 )
    3939
    4040! Lessivage
    4141
    42       CALL histwrite_phy( o_flux01       , flux01 )
    43       CALL histwrite_phy( o_flux02       , flux02 )
    44       CALL histwrite_phy( o_flux03       , flux03 )
    45       CALL histwrite_phy( o_flux04       , flux04 )
    46       CALL histwrite_phy( o_flux05       , flux05 )
    47       CALL histwrite_phy( o_ds01         , ds01 )
    48       CALL histwrite_phy( o_ds02         , ds02 )
    49       CALL histwrite_phy( o_ds03         , ds03 )
    50       CALL histwrite_phy( o_ds04         , ds04 )
    51       CALL histwrite_phy( o_ds05         , ds05 )
    52       CALL histwrite_phy( o_dh01         , dh01 )
    53       CALL histwrite_phy( o_dh02         , dh02 )
    54       CALL histwrite_phy( o_dh03         , dh03 )
    55       CALL histwrite_phy( o_dh04         , dh04 )
    56       CALL histwrite_phy( o_dh05         , dh05 )
    57       CALL histwrite_phy( o_dtrconv01    , dtrconv01 )
    58       CALL histwrite_phy( o_dtrconv02    , dtrconv02 )
    59       CALL histwrite_phy( o_dtrconv03    , dtrconv03 )
    60       CALL histwrite_phy( o_dtrconv04    , dtrconv04 )
    61       CALL histwrite_phy( o_dtrconv05    , dtrconv05 )
    62       CALL histwrite_phy( o_dtherm01     , dtherm01 )
    63       CALL histwrite_phy( o_dtherm02     , dtherm02 )
    64       CALL histwrite_phy( o_dtherm03     , dtherm03 )
    65       CALL histwrite_phy( o_dtherm04     , dtherm04 )
    66       CALL histwrite_phy( o_dtherm05     , dtherm05 )
    67       CALL histwrite_phy( o_dhkecv01     , dhkecv01 )
    68       CALL histwrite_phy( o_dhkecv02     , dhkecv02 )
    69       CALL histwrite_phy( o_dhkecv03     , dhkecv03 )
    70       CALL histwrite_phy( o_dhkecv04     , dhkecv04 )
    71       CALL histwrite_phy( o_dhkecv05     , dhkecv05 )
    72       CALL histwrite_phy( o_dhkelsc01    , dhkelsc01 )
    73       CALL histwrite_phy( o_dhkelsc02    , dhkelsc02 )
    74       CALL histwrite_phy( o_dhkelsc03    , dhkelsc03 )
    75       CALL histwrite_phy( o_dhkelsc04    , dhkelsc04 )
    76       CALL histwrite_phy( o_dhkelsc05    , dhkelsc05 )
    77       CALL histwrite_phy( o_d_tr_cv01    , d_tr_cv01 )
    78       CALL histwrite_phy( o_d_tr_cv02    , d_tr_cv02 )
    79       CALL histwrite_phy( o_d_tr_cv03    , d_tr_cv03 )
    80       CALL histwrite_phy( o_d_tr_cv04    , d_tr_cv04 )
    81       CALL histwrite_phy( o_d_tr_cv05    , d_tr_cv05 )
    82       CALL histwrite_phy( o_d_tr_trsp01  , d_tr_trsp01 )
    83       CALL histwrite_phy( o_d_tr_trsp02  , d_tr_trsp02 )
    84       CALL histwrite_phy( o_d_tr_trsp03  , d_tr_trsp03 )
    85       CALL histwrite_phy( o_d_tr_trsp04  , d_tr_trsp04 )
    86       CALL histwrite_phy( o_d_tr_trsp05  , d_tr_trsp05 )
    87       CALL histwrite_phy( o_d_tr_sscav01 , d_tr_sscav01 )
    88       CALL histwrite_phy( o_d_tr_sscav02 , d_tr_sscav02 )
    89       CALL histwrite_phy( o_d_tr_sscav03 , d_tr_sscav03 )
    90       CALL histwrite_phy( o_d_tr_sscav04 , d_tr_sscav04 )
    91       CALL histwrite_phy( o_d_tr_sscav05 , d_tr_sscav05 )
    92       CALL histwrite_phy( o_d_tr_sat01   , d_tr_sat01 )
    93       CALL histwrite_phy( o_d_tr_sat02   , d_tr_sat02 )
    94       CALL histwrite_phy( o_d_tr_sat03   , d_tr_sat03 )
    95       CALL histwrite_phy( o_d_tr_sat04   , d_tr_sat04 )
    96       CALL histwrite_phy( o_d_tr_sat05   , d_tr_sat05 )
    97       CALL histwrite_phy( o_d_tr_uscav01 , d_tr_uscav01 )
    98       CALL histwrite_phy( o_d_tr_uscav02 , d_tr_uscav02 )
    99       CALL histwrite_phy( o_d_tr_uscav03 , d_tr_uscav03 )
    100       CALL histwrite_phy( o_d_tr_uscav04 , d_tr_uscav04 )
    101       CALL histwrite_phy( o_d_tr_uscav05 , d_tr_uscav05 )
     42      IF(id_prec>0) CALL histwrite_phy( o_flux01       , flux01 )
     43      IF(id_fine>0) CALL histwrite_phy( o_flux02       , flux02 )
     44      IF(id_coss>0) CALL histwrite_phy( o_flux03       , flux03 )
     45      IF(id_codu>0) CALL histwrite_phy( o_flux04       , flux04 )
     46      IF(id_scdu>0) CALL histwrite_phy( o_flux05       , flux05 )
     47      IF(id_prec>0) CALL histwrite_phy( o_ds01         , ds01 )
     48      IF(id_fine>0) CALL histwrite_phy( o_ds02         , ds02 )
     49      IF(id_coss>0) CALL histwrite_phy( o_ds03         , ds03 )
     50      IF(id_codu>0) CALL histwrite_phy( o_ds04         , ds04 )
     51      IF(id_scdu>0) CALL histwrite_phy( o_ds05         , ds05 )
     52      IF(id_prec>0) CALL histwrite_phy( o_dh01         , dh01 )
     53      IF(id_fine>0) CALL histwrite_phy( o_dh02         , dh02 )
     54      IF(id_coss>0) CALL histwrite_phy( o_dh03         , dh03 )
     55      IF(id_codu>0) CALL histwrite_phy( o_dh04         , dh04 )
     56      IF(id_scdu>0) CALL histwrite_phy( o_dh05         , dh05 )
     57      IF(id_prec>0) CALL histwrite_phy( o_dtrconv01    , dtrconv01 )
     58      IF(id_fine>0) CALL histwrite_phy( o_dtrconv02    , dtrconv02 )
     59      IF(id_coss>0) CALL histwrite_phy( o_dtrconv03    , dtrconv03 )
     60      IF(id_codu>0) CALL histwrite_phy( o_dtrconv04    , dtrconv04 )
     61      IF(id_scdu>0) CALL histwrite_phy( o_dtrconv05    , dtrconv05 )
     62      IF(id_prec>0) CALL histwrite_phy( o_dtherm01     , dtherm01 )
     63      IF(id_fine>0) CALL histwrite_phy( o_dtherm02     , dtherm02 )
     64      IF(id_coss>0) CALL histwrite_phy( o_dtherm03     , dtherm03 )
     65      IF(id_codu>0) CALL histwrite_phy( o_dtherm04     , dtherm04 )
     66      IF(id_scdu>0) CALL histwrite_phy( o_dtherm05     , dtherm05 )
     67      IF(id_prec>0) CALL histwrite_phy( o_dhkecv01     , dhkecv01 )
     68      IF(id_fine>0) CALL histwrite_phy( o_dhkecv02     , dhkecv02 )
     69      IF(id_coss>0) CALL histwrite_phy( o_dhkecv03     , dhkecv03 )
     70      IF(id_codu>0) CALL histwrite_phy( o_dhkecv04     , dhkecv04 )
     71      IF(id_scdu>0) CALL histwrite_phy( o_dhkecv05     , dhkecv05 )
     72      IF(id_prec>0) CALL histwrite_phy( o_dhkelsc01    , dhkelsc01 )
     73      IF(id_fine>0) CALL histwrite_phy( o_dhkelsc02    , dhkelsc02 )
     74      IF(id_coss>0) CALL histwrite_phy( o_dhkelsc03    , dhkelsc03 )
     75      IF(id_codu>0) CALL histwrite_phy( o_dhkelsc04    , dhkelsc04 )
     76      IF(id_scdu>0) CALL histwrite_phy( o_dhkelsc05    , dhkelsc05 )
     77      IF(id_prec>0) CALL histwrite_phy( o_d_tr_cv01    , d_tr_cv01 )
     78      IF(id_fine>0) CALL histwrite_phy( o_d_tr_cv02    , d_tr_cv02 )
     79      IF(id_coss>0) CALL histwrite_phy( o_d_tr_cv03    , d_tr_cv03 )
     80      IF(id_codu>0) CALL histwrite_phy( o_d_tr_cv04    , d_tr_cv04 )
     81      IF(id_scdu>0) CALL histwrite_phy( o_d_tr_cv05    , d_tr_cv05 )
     82      IF(id_prec>0) CALL histwrite_phy( o_d_tr_trsp01  , d_tr_trsp01 )
     83      IF(id_fine>0) CALL histwrite_phy( o_d_tr_trsp02  , d_tr_trsp02 )
     84      IF(id_coss>0) CALL histwrite_phy( o_d_tr_trsp03  , d_tr_trsp03 )
     85      IF(id_codu>0) CALL histwrite_phy( o_d_tr_trsp04  , d_tr_trsp04 )
     86      IF(id_scdu>0) CALL histwrite_phy( o_d_tr_trsp05  , d_tr_trsp05 )
     87      IF(id_prec>0) CALL histwrite_phy( o_d_tr_sscav01 , d_tr_sscav01 )
     88      IF(id_fine>0) CALL histwrite_phy( o_d_tr_sscav02 , d_tr_sscav02 )
     89      IF(id_coss>0) CALL histwrite_phy( o_d_tr_sscav03 , d_tr_sscav03 )
     90      IF(id_codu>0) CALL histwrite_phy( o_d_tr_sscav04 , d_tr_sscav04 )
     91      IF(id_scdu>0) CALL histwrite_phy( o_d_tr_sscav05 , d_tr_sscav05 )
     92      IF(id_prec>0) CALL histwrite_phy( o_d_tr_sat01   , d_tr_sat01 )
     93      IF(id_fine>0) CALL histwrite_phy( o_d_tr_sat02   , d_tr_sat02 )
     94      IF(id_coss>0) CALL histwrite_phy( o_d_tr_sat03   , d_tr_sat03 )
     95      IF(id_codu>0) CALL histwrite_phy( o_d_tr_sat04   , d_tr_sat04 )
     96      IF(id_scdu>0) CALL histwrite_phy( o_d_tr_sat05   , d_tr_sat05 )
     97      IF(id_prec>0) CALL histwrite_phy( o_d_tr_uscav01 , d_tr_uscav01 )
     98      IF(id_fine>0) CALL histwrite_phy( o_d_tr_uscav02 , d_tr_uscav02 )
     99      IF(id_coss>0) CALL histwrite_phy( o_d_tr_uscav03 , d_tr_uscav03 )
     100      IF(id_codu>0) CALL histwrite_phy( o_d_tr_uscav04 , d_tr_uscav04 )
     101      IF(id_scdu>0) CALL histwrite_phy( o_d_tr_uscav05 , d_tr_uscav05 )
    102102
    103103
    104       CALL histwrite_phy( o_sed_ss , sed_ss)
    105       CALL histwrite_phy( o_sed_dust , sed_dust)
    106       CALL histwrite_phy( o_sed_dustsco , sed_dustsco)
    107       CALL histwrite_phy( o_g2p_gas , his_g2pgas )
    108       CALL histwrite_phy( o_g2p_aer , his_g2paer)
     104      IF(id_coss>0) CALL histwrite_phy( o_sed_ss , sed_ss)
     105      IF(id_codu>0) CALL histwrite_phy( o_sed_dust , sed_dust)
     106      IF(id_scdu>0) CALL histwrite_phy( o_sed_dustsco , sed_dustsco)
     107      IF(id_prec>0) CALL histwrite_phy( o_g2p_gas , his_g2pgas )
     108      IF(id_fine>0) CALL histwrite_phy( o_g2p_aer , his_g2paer)
    109109
    110110! histrac_spl
    111111!
    112 
    113112      CALL histwrite_phy( o_fluxbb              , fluxbb               )
    114113      CALL histwrite_phy( o_fluxff              , fluxff               )
     
    140139      CALL histwrite_phy( o_fluxh2sbio          , fluxh2sbio           )
    141140      CALL histwrite_phy( o_fluxdustec          , fluxdustec           )
    142      CALL histwrite_phy( o_fluxddfine          , fluxddfine           )
    143      CALL histwrite_phy( o_fluxddcoa           , fluxddcoa            )
    144      CALL histwrite_phy( o_fluxddsco           , fluxddsco            )
    145      CALL histwrite_phy( o_fluxdd              , fluxdd               )
     141      CALL histwrite_phy( o_fluxddfine          , fluxddfine           )
     142      CALL histwrite_phy( o_fluxddcoa           , fluxddcoa            )
     143      CALL histwrite_phy( o_fluxddsco           , fluxddsco            )
     144      CALL histwrite_phy( o_fluxdd              , fluxdd               )
    146145      CALL histwrite_phy( o_fluxssfine          , fluxssfine           )
    147146      CALL histwrite_phy( o_fluxsscoa           , fluxsscoa            )
     
    158157      CALL histwrite_phy( o_v10m_ss             , v10m_ss              )
    159158
    160 
    161159!      CALL histwrite_phy(  ,)
    162160!
Note: See TracChangeset for help on using the changeset viewer.