Gemclim 3.3.3.1 nudging : Différence entre versions

Un article de Informaticiens département des sciences de la Terre et l'atmosphère
Aller à: navigation, charcher
m (Controlling the profile)
m
Ligne 1: Ligne 1:
 
= Spectral nudging<br>  =
 
= Spectral nudging<br>  =
  
<br>  
+
Spectral nudging is a way to force (nudge) the larger scales of an RCM towards the ones of its driving GCM resp. RCM.<br>  
 +
 
 +
 
  
 
=== Grid size restrictions<br>  ===
 
=== Grid size restrictions<br>  ===
Ligne 11: Ligne 13:
 
Put the value of 'Pil_maxcfl' from your gemclim_settings.nml!<br> Your 'Grd_ni' and 'Grd_nj' will then have to be one of the values findfft spits out.<br> If your grid is larger than 500 points, just set '-gnimax' to a number large enough for you.  
 
Put the value of 'Pil_maxcfl' from your gemclim_settings.nml!<br> Your 'Grd_ni' and 'Grd_nj' will then have to be one of the values findfft spits out.<br> If your grid is larger than 500 points, just set '-gnimax' to a number large enough for you.  
  
== Activate spectral nudging in CRCM5 ==
 
  
To activate spectral nudging in CRCM5 one has to / can set several parameters in the file 'gemclim_settings.nml' into the namelist 'gem_cfgs'.
+
 
 +
== Activate spectral nudging in CRCM5  ==
 +
 
 +
To activate spectral nudging in CRCM5 one has to / can set several parameters in the file 'gemclim_settings.nml' into the namelist 'gem_cfgs'.  
 +
 
 +
 
  
 
=== Fields to nudge<br>  ===
 
=== Fields to nudge<br>  ===
  
In CRCM5 one has the option of nudging u-wind, v-wind, and/or temperature.<br>Tests have shown that it is sufficient to nudge u-wind and v-wind and that nudging the temperature as well does not change much the results. Therefore I&nbsp;suggest to set:
+
In CRCM5 one has the option of nudging u-wind, v-wind, and/or temperature.<br>Tests have shown that it is sufficient to nudge u-wind and v-wind and that nudging the temperature as well does not change much the results. Therefore I&nbsp;suggest to set:  
  
 
&nbsp; '''Lam_spn_nudging_S = 'UV'&nbsp;&nbsp;&nbsp;&nbsp; ,'''<br>  
 
&nbsp; '''Lam_spn_nudging_S = 'UV'&nbsp;&nbsp;&nbsp;&nbsp; ,'''<br>  
  
=== Half responds wavelength<br> ===
 
  
This is the wavelength at which equal weight is given to the driving data and RCM. For longer wavelength more weight is given to the driving data, for shorter wavelength more weight is given to the RCM. Therefore this is one way to control the strength of the nudging. Typical wavelengths are between 1000km and 2000km. For example set:<br>
 
  
'''&nbsp; Lam_spn_half_respon = 1000.&nbsp; ,'''
+
=== Half responds wavelength<br>  ===
 +
 
 +
This is the wavelength at which equal weight is given to the driving data and RCM. For longer wavelength more weight is given to the driving data, for shorter wavelength more weight is given to the RCM. Therefore this is one way to control the strength of the nudging. Typical wavelengths are between 1000km and 2000km. For example set:<br>
 +
 
 +
'''&nbsp; Lam_spn_half_respon = 1000.&nbsp; ,'''  
 +
 
  
 
=== Controlling the profile<br>  ===
 
=== Controlling the profile<br>  ===
Ligne 31: Ligne 40:
 
There are 3 parameters to control the nudging profile:  
 
There are 3 parameters to control the nudging profile:  
  
&nbsp;&nbsp;&nbsp; Lam_spn_up_const_lev&nbsp;&nbsp;&nbsp;&nbsp;: above this level (in hybrid values) the profile is set to 1 (default 0.).<br> &nbsp;&nbsp;&nbsp; Lam_spn_start_lev &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;&nbsp; : below this level (in hybrid values) the profile is set to 0 (default 1.).<br>&nbsp;&nbsp;&nbsp; Lam_spn_trans_shape_S&nbsp;: Between these two levels the profile goes from 0 to 1. Either linearly (''''LINEAR'''') or with a cos2 (''''COS2'''') function. <br>  
+
&nbsp;&nbsp;&nbsp; Lam_spn_up_const_lev&nbsp;&nbsp;&nbsp;&nbsp;: above this level (in hybrid values) the profile is set to 1 (default 0.).<br> &nbsp;&nbsp;&nbsp; Lam_spn_start_lev &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;&nbsp;&nbsp;: below this level (in hybrid values) the profile is set to 0 (default 1.).<br>&nbsp;&nbsp;&nbsp; Lam_spn_trans_shape_S&nbsp;: Between these two levels the profile goes from 0 to 1. Either linearly (''''LINEAR'''') or with a cos2 (''''COS2'''') function. <br>  
  
 
For example set:  
 
For example set:  
  
 
'''&nbsp;&nbsp;&nbsp; Lam_spn_up_const_lev &nbsp; = 0.&nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp; &nbsp; , <br> &nbsp;&nbsp;&nbsp; Lam_spn_start_lev&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; = 0.5 &nbsp; &nbsp;&nbsp; &nbsp;&nbsp; ,<br> &nbsp;&nbsp;&nbsp; Lam_spn_trans_shape_S = 'LINEAR' ,<br>'''
 
'''&nbsp;&nbsp;&nbsp; Lam_spn_up_const_lev &nbsp; = 0.&nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp; &nbsp; , <br> &nbsp;&nbsp;&nbsp; Lam_spn_start_lev&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; = 0.5 &nbsp; &nbsp;&nbsp; &nbsp;&nbsp; ,<br> &nbsp;&nbsp;&nbsp; Lam_spn_trans_shape_S = 'LINEAR' ,<br>'''

Version depuis le 19 de septembre 2014 à 14:32

Spectral nudging

Spectral nudging is a way to force (nudge) the larger scales of an RCM towards the ones of its driving GCM resp. RCM.


Grid size restrictions

To be able to use spectral nudging in CRCM5, the core grid needs to match the FFT criteria in x- and y-direction. You can use 'findfft' to find get the possible values for 'Grd_ni' and 'Grd_nj' in your gemclim_settings.nml.
For example:

    findfft -gnimax 500 -cfl Pil_maxcfl

Put the value of 'Pil_maxcfl' from your gemclim_settings.nml!
Your 'Grd_ni' and 'Grd_nj' will then have to be one of the values findfft spits out.
If your grid is larger than 500 points, just set '-gnimax' to a number large enough for you.


Activate spectral nudging in CRCM5

To activate spectral nudging in CRCM5 one has to / can set several parameters in the file 'gemclim_settings.nml' into the namelist 'gem_cfgs'.


Fields to nudge

In CRCM5 one has the option of nudging u-wind, v-wind, and/or temperature.
Tests have shown that it is sufficient to nudge u-wind and v-wind and that nudging the temperature as well does not change much the results. Therefore I suggest to set:

  Lam_spn_nudging_S = 'UV'     ,


Half responds wavelength

This is the wavelength at which equal weight is given to the driving data and RCM. For longer wavelength more weight is given to the driving data, for shorter wavelength more weight is given to the RCM. Therefore this is one way to control the strength of the nudging. Typical wavelengths are between 1000km and 2000km. For example set:

  Lam_spn_half_respon = 1000.  ,


Controlling the profile

There are 3 parameters to control the nudging profile:

    Lam_spn_up_const_lev    : above this level (in hybrid values) the profile is set to 1 (default 0.).
    Lam_spn_start_lev           : below this level (in hybrid values) the profile is set to 0 (default 1.).
    Lam_spn_trans_shape_S : Between these two levels the profile goes from 0 to 1. Either linearly ('LINEAR') or with a cos2 ('COS2') function.

For example set:

    Lam_spn_up_const_lev   = 0.           ,
    Lam_spn_start_lev         = 0.5         ,
    Lam_spn_trans_shape_S = 'LINEAR' ,