Changeset 4147


Ignore:
Timestamp:
Sep 3, 2019, 1:54:17 PM (5 years ago)
Author:
jisesh
Message:

guided: implemented temperature and wind nudging from DYNAMICO data

Location:
dynamico_lmdz/guided/Experiments
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • dynamico_lmdz/guided/Experiments/filedef_dcmip2016.xml

    r4144 r4147  
    11<file_definition >
    22
    3   <file id="output_dcmip2016_native" convention="CF" enabled="true" output_freq="60s" sync_freq="60s" output_level="10" timeseries="none" description="DCMIP31" >
     3  <file id="output_dcmip2016_native" convention="CF" enabled="false" output_freq="1h" sync_freq="1h" output_level="10" timeseries="none" description="DCMIP41" >
    44   
    55    <field_group id="dcmip2016_output_field_once" operation="once" freq_offset="0ts" ts_enabled="true">
    … …  
    2222      <field field_ref="geopot_init" />
    2323    </field_group>
    24 
     24   
    2525    <field_group id="nudging_fields" ts_enabled="true">
    2626      <field field_ref="ulon_guided_out" name="U_guided"  unit="m/s"/>
    2727      <field field_ref="ulat_guided_out" name="V_guided"      unit="m/s"/>
     28      <field field_ref="T_guided_out" name="T_guided"      unit="m/s"/>
    2829    </field_group>
    29    
     30
    3031    <field_group id="dcmip2016_output_field" ts_enabled="true">
    3132      <field field_ref="ps" name="PS"       standard_name="surface_pressure" long_name="Surface pressure"         unit="Pa" />
    … …  
    7273
    7374    <variable name="model" type="string" > dynamico </variable>
    74     <variable name="test_case" type="string" > DCMIP31 </variable>
     75    <variable name="test_case" type="string" > DCMIP41 </variable>
    7576    <variable name="horizontal_resolution" type="string" > native (grid points : 80x80x10) </variable>
    7677    <variable name="levels" type="string" > L10 </variable>
    7778    <variable name="grid" type="string" > hex </variable>
    7879    <variable name="equation" type="string" > hydro </variable>
    79     <variable name="time_frequency" type="string" > 60s </variable>   
     80    <variable name="time_frequency" type="string" > 24h </variable>   
    8081
    8182  </file>
    8283
    8384
    84   <file id="output_dcmip2016_regular" enabled="true" output_freq="60s" sync_freq="60s" output_level="10" description="DCMIP31" timeseries="none" >
     85  <file id="output_dcmip2016_regular" enabled="true" output_freq="1d" sync_freq="1d" output_level="10" description="DCMIP41" timeseries="none" >
    8586   
    8687    <field_group group_ref="dcmip2016_output_field_once"/>
    … …  
    100101
    101102    <variable name="model" type="string" > dynamico </variable>
    102     <variable name="test_case" type="string" > DCMIP31 </variable>
     103    <variable name="test_case" type="string" > DCMIP41 </variable>
    103104    <variable name="horizontal_resolution" type="string" > two_degree interpolated (grid points : 80 x 80 x 10) </variable>
    104105    <variable name="levels" type="string" > L30 </variable>
    105106    <variable name="grid" type="string" > hex </variable>
    106107    <variable name="equation" type="string" > hydro </variable>
    107     <variable name="time_frequency" type="string" > 60s </variable>   
     108    <variable name="time_frequency" type="string" > 24h </variable>   
    108109
    109110  </file>

  • dynamico_lmdz/guided/Experiments/icosa.xml

    r4144 r4147  
    11<context id="icosagcm" >
    2   <calendar type="d360" start_date="2000-01-01 00:00:00" time_origin="2000-01-01 00:00:00"/>
     2  <calendar type="d360" start_date="2019-01-01 00:00:00" time_origin="2019-01-01 00:00:00"/>
    33 
    44  <field_definition src="./field_def.xml" />
    … …  
    1111      <field_group group_ref="standard_output"/>
    1212    </file>
     13
    1314<!---
    1415
    … …  
    2324
    2425    <file id="guided" name="ERA4" output_freq="1d" enabled="true" mode="read">
    25        <field_group enabled="true" domain_ref="domain_guided" freq_offset="1ts" operation="instant" freq_op="1d">
    26           <field id="T_guided"    name="t"   axis_ref="lev_guided" scale_factor="0.00211281711159895" add_offset="249.62701823744"/>
    27           <field id="ulon_guided" name="u"   axis_ref="lev_guided" scale_factor="0.00428658392668996" add_offset="25.9159551577437"/>
    28           <field id="ulat_guided" name="v"   axis_ref="lev_guided" scale_factor="0.00310116627116181" add_offset="15.7956601102238"/>
     26       <field_group enabled="true" domain_ref="domain_guided" freq_offset="1ts" operation="instant" freq_op="1d">
     27          <field id="T_guided"    name="t"   axis_ref="lev_guided" />
     28          <field id="ulon_guided" name="u"   axis_ref="lev_guided" />
     29          <field id="ulat_guided" name="v"   axis_ref="lev_guided" />
     30          <field id="ps_guided" name="PS" />
    2931       </field_group>
    3032    </file>
    3133
    32 
    33 
    3434    <file id="test_guided" output_freq="1ts" enabled="true" sync_freq="10ts" >
    35        <field_group enabled="true" domain_ref="i" operation="instant"> 
     35       <field_group enabled="true" domain_ref="i" operation="instant">
    3636          <field id="ulon_guided_out" name="U"   axis_ref="lev" />
    3737          <field id="ulat_guided_out" name="V"   axis_ref="lev" />
     38          <field id="T_guided_out" name="T"   axis_ref="lev" />
     39          <field id="PS_guided_out" name="PS" />
    3840       </field_group>
    3941    </file>
    … …  
    8587   
    8688  <axis_definition>
    87 <!--    <axis id="lev" name="lev" positive="down" standard_name="atmopshere_hybrid_sigma_pressure_coordinate" long_name="hybrid level at midpoints" axis_type="Z"/> -->
    8889    <axis id="lev" name="lev" positive="down" standard_name="atmopshere_hybrid_sigma_pressure_coordinate" long_name="hybrid level at midpoints" />
    8990    <axis id="lev_read" />
    9091    <axis id="levp1" name="ilev"  unit="level" positive="down" standard_name="atmopshere_hybride_sigma_pressure_coordinate" long_name="hybrid level at interface" />
    9192    <axis id="lev_pressure" positive="down"  n_glo="2" unit="Pa" value="(0,1)[85000 50000]" >
    92       <interpolate_axis order="1" coordinate="p" />
     93      <interpolate_axis order="2" coordinate="p" />
    9394    </axis>   
    9495    <axis id="nq" name="nq" n_glo="5" value="(0,4)[1 2 3 4 5]"/>
    9596
    96       <axis id="lev_guided" />
    97       <axis id="presnivs_mb" >
    98         <interpolate_axis order="1" />
    99       </axis>
     97    <axis id="lev_guided" />
     98    <axis id="presnivs_mb" >
     99      <interpolate_axis order="2" />
     100    </axis>
    100101
    101102  </axis_definition>   
    … …  
    103104  <domain_definition>
    104105   
    105        <domain_group id="i">
    106      
    107          <domain id="i" name="mesh"/>   
    108      
    109          <domain id="from_guided" >   
    110            <interpolate_domain order="1"/>
     106     <domain_group id="i">
     107       <domain id="i" name="mesh"/>   
     108
     109         <domain id="from_guided" >
     110           <interpolate_domain order="2"/>
    111111         </domain>
    112      
    113        </domain_group>
     112 
     113     </domain_group>
     114     
    114115     
    115116     <domain id="v" />
    … …  
    118119     <domain id="regular_half_degree" type="rectilinear" ni_glo="720" nj_glo="360" >
    119120       <generate_rectilinear_domain />
    120        <interpolate_domain order="1"/>
     121       <interpolate_domain order="2"/>
    121122     </domain>
    122123     
    123124     <domain id="regular_one_degree" type="rectilinear" ni_glo="360" nj_glo="180" >
    124125       <generate_rectilinear_domain />
    125        <interpolate_domain order="1"/>
     126       <interpolate_domain order="2"/>
    126127     </domain>
    127128     
    128129     <domain id="regular_two_degree" type="rectilinear" ni_glo="180" nj_glo="90" >
    129130       <generate_rectilinear_domain />
    130        <interpolate_domain order="1"/>
     131       <interpolate_domain order="2"/>
    131132     </domain>
    132133
    … …  
    134135         <generate_rectilinear_domain />
    135136       </domain>
    136  
     137     
    137138  </domain_definition>   
    138139 
    … …  
    161162  </grid_definition>   
    162163
    163    
    164     <field_definition>
    165       <field_group id="guided_read"  domain_ref="from_guided" axis_ref="presnivs_mb" read_access="true">
    166   <!--    <field_group id="guided_read"  domain_ref="from_guided" axis_ref="lev" read_access="true"> -->
    167         <field     id="T_guided_read"    field_ref="T_guided" />
    168         <field     id="ulon_guided_read" field_ref="ulon_guided" />
    169         <field     id="ulat_guided_read" field_ref="ulat_guided" />
     164  <field_definition>
     165     
     166      <field_group id="guided_read"  domain_ref="from_guided" read_access="true">
     167            <field_group axis_ref="presnivs_mb">
     168              <field     id="T_guided_read"    field_ref="T_guided" />
     169              <field     id="ulon_guided_read" field_ref="ulon_guided" />
     170              <field     id="ulat_guided_read" field_ref="ulat_guided" />
     171            </field_group>
     172         <field id="ps_guided_read" field_ref="ps_guided" />
    170173      </field_group>
    171    
    172          
    173     </field_definition>
    174174
     175  </field_definition>
     176 
    175177</context>

  • dynamico_lmdz/guided/Experiments/run.def

    r4144 r4147  
     1grid_type = icosahedral
     2
    13#------------- Planet ---------------
    24radius=6.37122e6
    3 omega=0.
     5omega=7.292e-5
    46g=9.80616
    57cpp=1004.5
    68kappa=0.2857143
    79preff=1e5
    8 scale_factor=125
    910#------------- Nudging  ------------
    1011guided_type=nudging
    1112guided_nudging_field=wind_temperature
    1213nudging_center_lat=90.
    13 nudging_radius=1.5e7
     14#nudging_radius=1.5e7
     15nudging_radius=0.
     16#nudging_time=86400
    1417nudging_time=86400
    1518#------------- Dynamics  ------------
    16 #hydrostatic=.FALSE.
    1719#-------------- Mesh ----------------
    18 nbp=41
    19 #optim_it=500
    20 optim_it=00
     20#nbp=80
     21nbp=40
    2122nsplit_i=2
    22 nsplit_j=3
    23 omp_level_size=2
    24 llm=10
    25 disvert=ncar
    26 ncar_disvert_c=1
    27 ncar_dz=1118
    28 ncar_T0=300
    29 ncar_p0=1e5
     23nsplit_j=2
     24optim_it=500
     25llm=30
     26disvert=ncarl30
    3027#------------ Numerics --------------
    31 dt=480
    32 dysl=.TRUE.
     28dt=300
     29nqtot=2
     30itau_adv=3
    3331#----------- Dissipation ------------
    3432nitergdiv=2
    35 tau_graddiv=90000
     33tau_graddiv=10000
    3634nitergrot=2
    37 tau_gradrot=90000
     35tau_gradrot=10000
    3836niterdivgrad=2
    39 tau_divgrad=90000
     37tau_divgrad=10000
    4038#-------------- Physics -------------
    4139#---------------- Run ---------------
    42 run_length=4.5e5
    43 write_period=7500
    44 etat0=dcmip3
     40run_length=2592000
     41#write_period=3600
     42write_period=300
     43etat0=dcmip4
     44dcmip4_testcase=1
    4545#------------ Diagnostics -----------
     46diagflux=.TRUE.
Note: See TracChangeset for help on using the changeset viewer.