Ignore:
Timestamp:
Jun 5, 2020, 1:59:19 PM (5 years ago)
Author:
dubos
Message:

simple_physics : output SW fluxes

Location:
dynamico_lmdz/simple_physics/config
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • dynamico_lmdz/simple_physics/config/DYNAMICO/TEST/field_def_dynamico.xml

    r4239 r4242  
    1313    <field id="preff" grid_ref="scalar"    unit="Pa"/>
    1414    <field id="ap" axis_ref="levp1"   long_name="hybrid A coefficient at layer interface" />
    15     <field id="bp" axis_ref="levp1"  long_name="hybrid B coefficient at layer interface" />
    16     <field id="mid_ap" axis_ref="lev"   long_name="hybrid A coefficient at midpoints" />
    17     <field id="mid_bp" axis_ref="lev"  long_name="hybrid B coefficient at midpoints" />
     15    <field id="bp" axis_ref="levp1"   long_name="hybrid B coefficient at layer interface" />
     16    <field id="mid_ap" axis_ref="lev" long_name="hybrid A coefficient at midpoints" />
     17    <field id="mid_bp" axis_ref="lev" long_name="hybrid B coefficient at midpoints" />
    1818   
    1919    <field_group domain_ref="i">
    20       <field id="ps" />
     20      <field id="ps"   standard_name="surface_pressure"     long_name="Surface pressure" unit="Pa"/>
     21      <field id="phis" standard_name="surface_geopotential" long_name="Surface geopotential" unit="m2/m2"/>
    2122      <field id="dps"  />
    2223      <field id="Ai" />
    23       <field id="phis" />
    2424      <field id="phi" />
    2525      <field id="precl" />
     
    2727      <field id="Q2_col_int" />     
    2828      <field id="ps_init" />
    29      
     29
    3030      <field_group axis_ref="lev">
    3131        <field id="mass" />
    32         <field id="p" />
    3332        <field id="dmass"/>
    34         <field id="theta" />
    3533        <field id="dyn_q" />
    3634        <field id="pk"/>
    37         <field id="ulon"/>
    38         <field id="ulat"/>
    39         <field id="uz"/>
    40         <field id="omega"/>
    41         <field id="temp"/>
     35        <field id="p"     standard_name="air_pressure"        long_name="Pressure"              unit="Pa"/>
     36        <field id="temp"  standard_name="air_temperature"     long_name="Temperature"           unit="K"/>
     37        <field id="theta" standard_name="pot_temp"            long_name="Potential temperature" unit="m/s"/>
     38        <field id="ulon"  standard_name="eastward_wind"       long_name="Zonal wind"            unit="m/s"/>
     39        <field id="ulat"  standard_name="northward_wind"      long_name="Meridional wind"       unit="m/s"/>
     40        <field id="uz"    standard_name="upward_air_velocity" long_name="Vertical velocity"     unit="m/s"/>
     41        <field id="omega" standard_name="lagrangian_tendency_of_air_pressure" long_name="Vertical pressure velocity"     unit="Pa/s"/>
    4242        <field id="kinetic_trisk"/>
    4343        <field id="kinetic"/>
     
    9393       
    9494      </field_group>
    95       <field id="u850" />
    96       <field id="v850" />
    97       <field id="w850" />
    98       <field id="t850" />
    99       <field id="omega850"/>                 
    100       <field id="u500" />
    101       <field id="v500" />
    102       <field id="w500" />
    103       <field id="t500" />
    104       <field id="omega500"/>                 
    105       <field id="SST"/>
    106      
     95
     96      <field id="u850" long_name="Zonal wind at 850 hPa"      unit="m/s"/>
     97      <field id="u500" long_name="Zonal wind at 500 hPa"      unit="m/s"/>
     98      <field id="v850" long_name="Meridional wind at 850 hPa" unit="m/s"/>
     99      <field id="v500" long_name="Meridional wind at 850 hPa" unit="m/s"/>
     100      <field id="w850" long_name="Vertical wind at 850 hPa"   unit="m/s"/>
     101      <field id="w500" long_name="Vertical wind at 500 hPa"   unit="m/s"/>
     102      <field id="t850" long_name="Temperature at 850 hPa"     unit="K"/>
     103      <field id="t500" long_name="Temperature at 500 hPa"     unit="K"/>
     104      <field id="omega850" long_name="Vertical pressure velocity at 850 hPa" unit="Pa/s"/>
     105      <field id="omega500" long_name="Vertical pressure velocity at 500 hPa" unit="Pa/s"/>
     106      <field id="SST"/>     
    107107      <field id="ulon_850_500" field_ref="p" axis_ref="lev_pressure" />
    108108     
    109109      <field_group axis_ref="levp1">
    110         <field id="geopot"/>
     110        <field id="geopot" standard_name="geopotential" long_name="geopotential" unit="m2/s2"/>
    111111        <field id="geopot_init"/>
    112112      </field_group>
     
    124124 
    125125</field_definition>
    126 
  • dynamico_lmdz/simple_physics/config/DYNAMICO/TEST/field_def_physics.xml

    r4239 r4242  
    11<!-- =========================================================================================================== -->
    2 <!-- field_def_phyparam.xml                                                                                      -->
     2<!-- field_def_physics.xml                                                                                      -->
    33<!-- Definition of all existing variables that can be output from phyparam                                       -->
    44<!-- =========================================================================================================== -->
     
    88  <field_group id="phyparam_output">     
    99    <field_group domain_ref="i">
    10       <field id="phyparam_mu0" />
     10      <field id="phyparam_ts"     name="Tsurf"  long_name="Surface temperature"     unit="K"/>
     11      <field id="phyparam_mu0"    name="mu0"    long_name="Cosine zenithal angle"/>
     12      <field id="phyparam_swsurf" name="swsurf" long_name="Upward SW flux at surface"/>
     13      <field id="phyparam_swtop"  name="swtop"  long_name="Downward SW fluw at TOA" unit="W/m2" />
     14      <field id="phyparam_lwsurf" name="lwsurf" long_name="LW flux at surface"      unit="W/m2" />
    1115      <field id="phyparam_alb" />
    12       <field id="phyparam_ts" />
    1316      <field id="phyparam_ps" />
    1417      <field id="phyparam_slp" />
    15       <field id="phyparam_swsurf" />
    16       <field id="phyparam_lwsurf" />
    1718
    1819      <field id="phyparam_coslon" />
     
    2021      <field id="phyparam_coslat" />
    2122      <field id="phyparam_sinlat" />
     23
     24      <field_group axis_ref="lev">
     25        <field id="phyparam_u"    name="u" long_name="Zonal velocity"      unit="m/s" />
     26        <field id="phyparam_v"    name="v" long_name="Meridional velocity" unit="m/s" />
     27        <field id="phyparam_temp" name="T" long_name="Temperature"         unit="K"   />
     28        <field id="phyparam_geop" />
     29        <field id="phyparam_plev" />
     30
     31        <field id="phyparam_du" name="du"    long_name="Zonal velocity tendency"      unit="m/s2" />
     32        <field id="phyparam_dv" name="dv"    long_name="Meridional velocity tendency" unit="m/s2" />
     33        <field id="phyparam_dt" name="dT"    long_name="Temperature tendency"         unit="K/s"  />
     34        <field id="phyparam_dtsw" name="dT_SW" long_name="Temperature tendency due to shortwave radiation" unit="K/s" />
     35        <field id="phyparam_dtlw" name="dT_LW" long_name="Temperature tendency due to longwave radiation"  unit="K/s" />
     36      </field_group>
     37
     38      <field_group axis_ref="levp1">
     39        <field id="phyparam_swflux_down" name="swflux_down" long_name="Downward SW radiative flux" unit="W/m2" />
     40      </field_group>
     41     
    2242    </field_group>
    2343  </field_group>
    2444 
    2545</field_definition>
    26 
    27 <!--       
    28        call writefield('u','Vent zonal moy','m/s',pu)
    29        call writefield('v','Vent meridien moy','m/s',pv)
    30        call writefield('temp','Temperature','K',pt)
    31        call writefield('geop','Geopotential','m2/s2',pphi)
    32        call writefield('plev','plev','Pa',pplev(:,1:nlayer))
    33 
    34        call writefield('du','du',' ',pdu)
    35        call writefield('dv','du',' ',pdv)
    36        call writefield('dt','du',' ',pdt)
    37        call writefield('dtsw','dtsw',' ',zdtsw)
    38        call writefield('dtlw','dtlw',' ',zdtlw)
    39 
    40        call writefield('ts','Surface temper','K',tsurf)
    41        call writefield('coslon','coslon',' ',coslon)
    42        call writefield('sinlon','sinlon',' ',sinlon)
    43        call writefield('coslat','coslat',' ',coslat)
    44        call writefield('sinlat','sinlat',' ',sinlat)
    45        call writefield('mu0','mu0',' ',mu0)
    46        call writefield('alb','alb',' ',albedo)
    47        call writefield('ps','Surface pressure','Pa',pplev(:,1))
    48        call writefield('slp','Sea level pressure','Pa',zpmer)
    49        call writefield('swsurf','SW surf','Pa',zfluxsw)
    50        call writefield('lwsurf','LW surf','Pa',zfluxlw)
    51 -->
  • dynamico_lmdz/simple_physics/config/DYNAMICO/TEST/file_def_dynamico.xml

    r4239 r4242  
    77  <file id="dynamico_native" convention="CF" enabled="false" output_freq="3h" sync_freq="3h" output_level="10" timeseries="none" description="Simple physics" >
    88   
    9     <field_group id="dcmip2016_output_field_once" operation="once" freq_offset="0ts" ts_enabled="true">
     9    <field_group id="output_field_once" operation="once" freq_offset="0ts" ts_enabled="true">
    1010      <field field_ref="timestep" name="mdt" enabled="false" />
    1111      <field field_ref="preff" name="P0" long_name="reference pressure" />
     
    1717   
    1818   
    19     <field_group id="dcmip2016_output_field" ts_enabled="true">
    20       <field field_ref="ps" name="PS"       standard_name="surface_pressure" long_name="Surface pressure"         unit="Pa" />
    21       <field field_ref="phis" operation="once" freq_offset="0ts" name="PHIS"   standard_name="surface_geopotential" long_name="Surface geopotential" unit="m2/m2"/>
     19    <field_group id="output_field" ts_enabled="true">
     20
     21      <field field_ref="phis" operation="once" freq_offset="0ts" />
     22      <field field_ref="ps" name="PS" />
     23      <field field_ref="ulon"  />
     24      <field field_ref="ulat"  />
     25      <field field_ref="temp"  />
     26      <field field_ref="theta" />
     27
     28      <field field_ref="u850" />
     29      <field field_ref="v850" />
     30      <field field_ref="t850" />
     31      <field field_ref="omega850" />
     32      <field field_ref="u500" />
     33      <field field_ref="v500" />
     34      <field field_ref="t500" />
     35      <field field_ref="omega500" />
    2236
    2337<!--     
    24       <field field_ref="geopot"  name="PHI"   standard_name="geopotential" long_name="geopotential" unit="m2/s2"/>
    25       <field field_ref="uz" name="W"        standard_name="upward_air_velocity" long_name="Vertical velocity"     unit="m/s"/>
    26       <field field_ref="omega" name="OMEGA"        standard_name="lagrangian_tendency_of_air_pressure" long_name="Vertical pressure velocity"     unit="Pa/s"/>
    27       <field field_ref="p"  name="P"        standard_name="air_pressure" long_name="Pressure"     unit="Pa"/>
    28       <field field_ref="w850" name="W850"       long_name="Vertical velocity at 850 hPa"     unit="m/s"/>
    29       <field field_ref="w500" name="W500"       long_name="Vertical velocity at 500 hPa"     unit="m/s"/>
     38      <field field_ref="geopot" />
     39      <field field_ref="uz" />
     40      <field field_ref="omega" />
     41      <field field_ref="p" />
     42      <field field_ref="w850" />
     43      <field field_ref="w500" />
    3044-->
    31       <field field_ref="ulon" name="U"      standard_name="eastward_wind" long_name="Zonal wind"                  unit="m/s"/>
    32       <field field_ref="ulat" name="V"      standard_name="northward_wind" long_name="Meridional wind"            unit="m/s"/>
    33       <field field_ref="temp"  name="T"        standard_name="air_temperature" long_name="Temperature"     unit="K"/>
    3445
    35       <field field_ref="u850" name="U850"       long_name="Zonal wind at 850 hPa"      unit="m/s"/>
    36       <field field_ref="v850" name="V850"       long_name="Meridional wind at 850 hPa"   unit="m/s"/>
    37       <field field_ref="t850" name="T850"       long_name="Temperature at 850 hPa"   unit="K"/>
    38       <field field_ref="omega850" name="OMEGA850"   long_name="Vertical pressure velocity at 850 hPa"     unit="Pa/s"/>
    39       <field field_ref="u500" name="U500"       long_name="Zonal wind at 500 hPa"      unit="m/s"/>
    40       <field field_ref="v500" name="V500"       long_name="Meridional wind at 850 hPa"   unit="m/s"/>
    41       <field field_ref="t500" name="T500"       long_name="Temperature at 850 hPa"   unit="K"/>
    42       <field field_ref="omega500" name="OMEGA500"   long_name="Vertical pressure velocity at 500 hPa"     unit="Pa/s"/>
    4346     
    4447    </field_group>
     
    5356  <file id="dynamico_regular" enabled="true" output_freq="24h" sync_freq="24h" output_level="10" description="Simple physics" timeseries="none" >
    5457   
    55     <field_group group_ref="dcmip2016_output_field_once"/>
    56    
    57     <field_group group_ref="dcmip2016_output_field" domain_ref="regular_one_degree" />
    58 
     58    <field_group group_ref="output_field_once"/>
     59    <field_group group_ref="output_field" domain_ref="regular_one_degree" />
    5960    <variable name="model" type="string" > dynamico </variable>
    6061    <variable name="grid" type="string" > hex </variable>
  • dynamico_lmdz/simple_physics/config/DYNAMICO/TEST/file_def_physics.xml

    r4239 r4242  
    44   
    55    <field_group id="phyparam_output_field" ts_enabled="true">
    6       <field field_ref="phyparam_ts" unit="K"/>
    7       <field field_ref="phyparam_swsurf" unit="W/m2"/>
    8       <field field_ref="phyparam_lwsurf" unit="W/m2"/>
     6
     7      <field field_ref="phyparam_mu0" />
     8      <field field_ref="phyparam_ts" />
     9      <field field_ref="phyparam_swsurf" />
     10      <field field_ref="phyparam_swtop" />
     11      <field field_ref="phyparam_lwsurf" />
     12
     13      <field field_ref="phyparam_u" />
     14      <field field_ref="phyparam_v" />
     15      <field field_ref="phyparam_temp" />
     16      <field field_ref="phyparam_du" />
     17      <field field_ref="phyparam_dv" />
     18      <field field_ref="phyparam_dtsw" />
     19      <field field_ref="phyparam_dtlw" />
     20      <field field_ref="phyparam_swflux_down" />
     21
    922    </field_group>
    1023
  • dynamico_lmdz/simple_physics/config/DYNAMICO/TEST/run.def

    r4239 r4242  
    3333
    3434#-------------- Physics -------------
    35 itau_physics=1
     35itau_physics=5
    3636physics=plugin
     37#physics_log_level=DBG
     38
    3739# 0.25 day = 6h
    3840period_sort=0.25
  • dynamico_lmdz/simple_physics/config/LMDZ/USEFUL

    r4230 r4242  
    33(cd TEST_PARAM/; time ./gcm.e | tee gcm.log | grep '*    pas ' && ./check.sh IRENE)
    44grep WRN TEST_PARAM/gcm.log | sort -u
     5
     6# IPSL / Camelot
     7module load gnu/7.2.0
Note: See TracChangeset for help on using the changeset viewer.