New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10516 for NEMO/trunk/tests/ICE_AGRIF/EXPREF – NEMO

Ignore:
Timestamp:
2019-01-15T17:31:25+01:00 (5 years ago)
Author:
clem
Message:

change demo config name from ICEDYN to ICE_AGRIF

Location:
NEMO/trunk/tests/ICE_AGRIF
Files:
7 edited
1 moved

Legend:

Unmodified
Added
Removed

  • NEMO/trunk/tests/ICE_AGRIF/EXPREF/1_namelist_cfg

    r10073 r10516  
    22!! NEMO/OCE  Configuration namelist : overwrite defaults values defined in SHARED/namelist_ref 
    33!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    4 !!               ICEDYN (SAS BIPERIODIC) configuration                !! 
     4!!               ICE_AGRIF configuration                              !! 
    55!!====================================================================== 
    66!!              ***  Domain & Run management namelists  ***           !! 
    … …  
    1818! 
    1919!----------------------------------------------------------------------- 
    20 &namusr_def    !   ICEDYN user defined namelist 
     20&namusr_def    !   ICE_AGRIF user defined namelist 
    2121!----------------------------------------------------------------------- 
    2222   rn_dx       =    1000.   ! horizontal resolution in meters 
    … …  
    2929!----------------------------------------------------------------------- 
    3030   nn_no       =       0      !  job number (no more used...) 
    31    cn_exp      =  "ICEDYN"    !  experience name 
     31   cn_exp      =  "ICE_AGRIF" !  experience name 
    3232   nn_it000    =       1      !  first time step 
    3333   nn_itend    =    1500      !  last  time step (std 5475) 
    … …  
    4848   ln_read_cfg = .false.    !  (=T) read the domain configuration file 
    4949      !                     !  (=F) user defined configuration  ==>>>  see usrdef(_...) modules 
    50       cn_domcfg = "ICEDYN_domcfg"    ! domain configuration filename 
     50      cn_domcfg = "ICE_AGRIF_domcfg"    ! domain configuration filename 
    5151/ 
    5252!!====================================================================== 

  • NEMO/trunk/tests/ICE_AGRIF/EXPREF/1_namelist_ice_cfg

    r10432 r10516  
    2828&namitd         !   Ice discretization 
    2929!------------------------------------------------------------------------------ 
    30    rn_himin         =   0.1           !  minimum ice thickness (m) used in remapping 
    3130/ 
    3231!------------------------------------------------------------------------------ 

  • NEMO/trunk/tests/ICE_AGRIF/EXPREF/README

    r9789 r10516  
    22Purpose 
    33------- 
    4 This demonstration case can serve to different purposes: 
    5  
    6 1) The main one: study of the advection of a patch of ice in a bi-periodic channel with a slab ocean (i.e. 1 ocean layer) 
    7               ==> set ln_icethd=false in namelist_ice_cfg 
    8  
    9 2) Secondary: study of ice thermodynamics in the same basin 
    10               ==> set ln_icedyn=false in namelist_ice_cfg 
    11               ==> eventually change usrdef_sbc.F90 in MY_SRC to fit your needs 
     4This demonstration case can be used to study the advection of a patch of ice in a bi-periodic channel 
     5with a slab ocean (i.e. 1 ocean layer) and an agrif zoom in the center 
    126 
    137----------- 
    … …  
    1711 
    1812a) A simple channel at 3km horizontal resolution (slightly variable => +- 10%) 
    19               ==> remove key_agrif in cpp keys 
     13              ==> in MY_SRC/usrdef_hgr.F90, uncomment the part "variable scale factors" and comment part "constant scale factors" 
    2014 
    21 b) Same as a) but with an Agrif zoom in the center 1:3, and 3 ghost cells 
    22               ==> add key_agrif in cpp keys 
     15b) A simple channel at 3km horizontal resolution (constant scale factors) 
     16              ==> in MY_SRC/usrdef_hgr.F90, comment the part "variable scale factors" and uncomment part "constant scale factors" 
    2317 
    2418---------- 
    … …  
    2721 
    2822a) Compile and run the model once to get a mesh_mask.nc file with the following command: 
    29 ../../../makenemo -a TEST_CASES -n ICEDYN -m X64_ADA -j 4 
    30 poe ./opa -procs 1 
     23../../../makenemo -r ICE_AGRIF -n ICE_AGRIF -m X64_ADA -j 4 
     24mpirun ./nemo -np 1 
    3125 
    3226b) Create the initial condition file for sea-ice (initice.nc) by running this python script:  
    … …  
    3428 
    3529c) Run the model a second time 
    36 poe ./opa -procs 4 
     30mpirun ./nemo -np 4 
    3731 
    3832--------------- 
    3933What to look at 
    4034--------------- 
    41 In case of purpose 1, One can test 
    42 --- 
    43 a) the advection scheme: Ultimate-Macho (ln_adv_UMx=T) versus Prather (ln_adv_Pra=T) 
     35the advection through an agrif zoom 1:3 
    4436   for a square (ice concentration) or a gaussian (ice volume) 
    45       with either a constant velocity (ln_dynADV=T) 
     37      with either a constant velocity (ln_dynADV2D=T) 
    4638               or a constant ice-atm. stress, thus velocity is calculated by rheology (ln_dynRHGADV=T) 
    4739      with 1 or 5 ice categories (jpl=1 or 5) 
    48  
    49 b) the advection through an agrif zoom 1:3 
    50  
    51 In case of purpose 2, one can test conservation of properties: 
    52 --- 
    53 c) ice should not change at all if surface fluxes = 0 and SST = freezing temperature 
    54  
    55 -------------------------------------- 
    56 Interpretation of the results (remarks) 
    57 -------------------------------------- 
    58 - Prather conserves the max values but also creates side lobes 
    59 - UM does not conserve the max but does not create side lobes 
    60 - The "unmoving" ice (if any) is due to ice rheology which states that ice mass below a certain thresold (1kg/m2) is considered to move at the ocean velocity (thus 0 m/s) 
     40         (note that ln_dynADV2D=T only works with jpl=1) 

  • NEMO/trunk/tests/ICE_AGRIF/EXPREF/file_def_nemo-ice.xml

    r10413 r10516  
    1111     
    1212     
    13       <file_group id="1d" output_freq="1d"  output_level="10" enabled=".TRUE.">  <!-- 1 day files -->     
     13      <file_group id="5ts" output_freq="5ts"  output_level="10" enabled=".TRUE.">  <!-- 1 day files -->      
    1414 
    1515   <file id="file1" name_suffix="_icemod" description="ice variables" enabled=".true." > 

  • NEMO/trunk/tests/ICE_AGRIF/EXPREF/make_INITICE.py

    r9019 r10516  
    8989 
    9090# -------------------------------------- 
    91 # for basin=99x99km with dx=1km ; dy=3km 
    92 #sigx=-0.04 
    93 #sigy=-0.04*9. 
    94 #xshift=50.-1. 
    95 #yshift=17.-1. 
    96 #dlat=7 
    97 #dlon=21 
    98  
    99 # -------------------------------------- 
    10091# for basin=99x99km with dx=1km ; dy=1km 
    10192#sigx=-0.04 
    … …  
    114105#        cdfsmi[:,y,x] = rn_smi_ini 
    115106# 
    116 # --- Lipscomb 2004 experiment --- 
    117 #cdfhti[:,:,:]=1. 
    118 #cdfati[:,:,:]=0.1 
    119 #for y in np.arange(0,LAT1,1) : 
    120 #    for x in np.arange(0,LON1,1) : 
    121 #        if (x > ceil(0.5*xshift) and x < xshift): # and (y > ceil(0.5*yshift) and y < yshift): 
    122 #            cdfati[:,y,x] = rn_ati_ini / (ceil(0.5*xshift)*ceil(0.5*yshift)) * (x - ceil(0.5*xshift)) * (y - ceil(0.5*yshift)) 
    123 #            cdfsmi[:,y,x] = rn_smi_ini / (ceil(0.5*xshift)*ceil(0.5*yshift)) * (x - ceil(0.5*xshift)) * (y - ceil(0.5*yshift)) 
    124 #            cdfhti[:,y,x] = rn_hti_ini 
    125 #            cdfhts[:,y,x] = rn_hts_ini 
    126 #        if (x >= 25. and x < 50.): 
    127 #            cdfati[:,y,x] = 0.9 * (x - 25.) / 25. 
    128 #            cdfhti[:,y,x] = 1. 
    129 #        elif (x >= 50. and x <= 75.): 
    130 #            cdfati[:,y,x] = 0.9 
    131 #            cdfhti[:,y,x] = 1. 
    132 #        if (x > 30. and x < 70.): 
    133 #            cdfhti[:,y,x] = 0.2 
    134              
    135107             
    136108# ---------------------------------------------- 
    137 # for basin=99x99km with dx=1km ; dy=1km + AGRIF 
     109# for basin=300x300km with dx=3km ; dy=3km + AGRIF 
    138110sigx=-0.012 
    139111sigy=-0.012 

  • NEMO/trunk/tests/ICE_AGRIF/EXPREF/namelist_cfg

    r10073 r10516  
    22!! NEMO/OCE  Configuration namelist : overwrite default values defined in SHARED/namelist_ref 
    33!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    4 !!               ICEDYN (SAS BIPERIODIC) configuration                !! 
     4!!               ICE_AGRIF configuration                              !! 
    55!!====================================================================== 
    66!!              ***  Domain & Run management namelists  ***           !! 
    … …  
    1818! 
    1919!----------------------------------------------------------------------- 
    20 &namusr_def    !   ICEDYN user defined namelist 
     20&namusr_def    !   ICE_AGRIF user defined namelist 
    2121!----------------------------------------------------------------------- 
    2222   rn_dx       =    3000.   ! horizontal resolution in meters 
    … …  
    2929!----------------------------------------------------------------------- 
    3030   nn_no       =       0      !  job number (no more used...) 
    31    cn_exp      =  "ICEDYN"    !  experience name 
     31   cn_exp      =  "ICE_AGRIF" !  experience name 
    3232   nn_it000    =       1      !  first time step 
    3333   nn_itend    =     500      !  last  time step (std 5475) 
    … …  
    4848   ln_read_cfg = .false.    !  (=T) read the domain configuration file 
    4949      !                     !  (=F) user defined configuration  ==>>>  see usrdef(_...) modules 
    50       cn_domcfg = "ICEDYN_domcfg"    ! domain configuration filename 
     50      cn_domcfg = "ICE_AGRIF_domcfg"    ! domain configuration filename 
    5151/ 
    5252!!====================================================================== 

  • NEMO/trunk/tests/ICE_AGRIF/EXPREF/namelist_ice_cfg

    r10432 r10516  
    2828&namitd         !   Ice discretization 
    2929!------------------------------------------------------------------------------ 
    30    rn_himin         =   0.1           !  minimum ice thickness (m) used in remapping 
    3130/ 
    3231!------------------------------------------------------------------------------ 
Note: See TracChangeset for help on using the changeset viewer.