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 – NEMO

Changeset 10516 for 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
Files:
1 deleted
13 edited
1 copied
1 moved

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/tests/ICE_ADV2D/EXPREF/file_def_nemo-ice.xml

    r10515 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." > 
     
    2121        <field field_ref="snwvolu"          name="snvolu" /> 
    2222        <field field_ref="icethic"          name="sithic" /> 
     23        <field field_ref="icethic"          name="sithic_max" operation="maximum" /> 
     24        <field field_ref="icethic"          name="sithic_min" operation="minimum" /> 
     25        <field field_ref="iceneg_pres"      name="sineg_pres" /> 
     26        <field field_ref="iceneg_volu"      name="sineg_volu" /> 
     27        <field field_ref="fasticepres"      name="fasticepres" /> 
    2328        <field field_ref="icevolu"          name="sivolu" /> 
    2429        <field field_ref="iceconc"          name="siconc" /> 
  • 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!------------------------------------------------------------------------------ 
  • NEMO/trunk/tests/ICE_AGRIF/MY_SRC/usrdef_hgr.F90

    r10074 r10516  
    33   !!                       ***  MODULE  usrdef_hgr  *** 
    44   !! 
    5    !!                      ===  ICEDYN configuration  === 
     5   !!                      ===  ICE_AGRIF configuration  === 
    66   !! 
    77   !! User defined :   mesh and Coriolis parameter of a user configuration 
     
    1111 
    1212   !!---------------------------------------------------------------------- 
    13    !!   usr_def_hgr    : initialize the horizontal mesh for ICEDYN configuration 
     13   !!   usr_def_hgr    : initialize the horizontal mesh for ICE_AGRIF configuration 
    1414   !!---------------------------------------------------------------------- 
    1515   USE dom_oce         ! ocean space and time domain 
     
    4545      !! 
    4646      !! ** Method  :   set all intent(out) argument to a proper value 
    47       !!                ICEDYN configuration : uniform grid spacing (rn_dx) 
     47      !!                ICE_AGRIF configuration : uniform grid spacing (rn_dx) 
    4848      !!                without Coriolis force (f=0) 
    4949      !! 
     
    6868      ! 
    6969      IF(lwp) WRITE(numout,*) 
    70       IF(lwp) WRITE(numout,*) 'usr_def_hgr : ICEDYN configuration bassin' 
     70      IF(lwp) WRITE(numout,*) 'usr_def_hgr : ICE_AGRIF configuration bassin' 
    7171      IF(lwp) WRITE(numout,*) 
    7272      IF(lwp) WRITE(numout,*) '          f-plane with irregular grid-spacing (+- 10%)' 
  • NEMO/trunk/tests/ICE_AGRIF/MY_SRC/usrdef_nam.F90

    r10161 r10516  
    33   !!                       ***  MODULE  usrdef_nam  *** 
    44   !! 
    5    !!                      ===  ICEDYN configuration  === 
     5   !!                      ===  ICE_AGRIF configuration  === 
    66   !! 
    77   !! User defined : set the domain characteristics of a user configuration 
     
    4848      !! ** Method  :   read in namusr_def containing all the user specific namelist parameter 
    4949      !! 
    50       !!                Here ICEDYN configuration 
     50      !!                Here ICE_AGRIF configuration 
    5151      !! 
    5252      !! ** input   : - namusr_def namelist found in namelist_cfg 
     
    8181      WRITE( ldnam(:), namusr_def ) 
    8282      ! 
    83       cd_cfg = 'ICEDYN'           ! name & resolution (not used) 
     83      cd_cfg = 'ICE_AGRIF'           ! name & resolution (not used) 
    8484      kk_cfg = INT( rn_dx ) 
    8585      ! 
    86       ! Global Domain size:  ICEDYN domain is  300 km x 300 Km x 10 m 
     86      ! Global Domain size:  ICE_AGRIF domain is  300 km x 300 Km x 10 m 
    8787      kpi = INT( 300.e3 / rn_dx ) -1 
    8888      kpj = INT( 300.e3 / rn_dy ) -1 
     
    103103      WRITE(ldtxt(ii),*) 'usr_def_nam  : read the user defined namelist (namusr_def) in namelist_cfg'   ;   ii = ii + 1 
    104104      WRITE(ldtxt(ii),*) '~~~~~~~~~~~ '                                                                 ;   ii = ii + 1 
    105       WRITE(ldtxt(ii),*) '   Namelist namusr_def : ICEDYN test case'                                    ;   ii = ii + 1 
     105      WRITE(ldtxt(ii),*) '   Namelist namusr_def : ICE_AGRIF test case'                                 ;   ii = ii + 1 
    106106      WRITE(ldtxt(ii),*) '      horizontal resolution                    rn_dx  = ', rn_dx, ' meters'   ;   ii = ii + 1 
    107107      WRITE(ldtxt(ii),*) '      horizontal resolution                    rn_dy  = ', rn_dy, ' meters'   ;   ii = ii + 1 
    108       WRITE(ldtxt(ii),*) '      ICEDYN domain = 300 km x 300Km x 1 grid-point '                         ;   ii = ii + 1 
     108      WRITE(ldtxt(ii),*) '      ICE_AGRIF domain = 300 km x 300Km x 1 grid-point '                      ;   ii = ii + 1 
    109109      WRITE(ldtxt(ii),*) '         LX [km]: ', zlx                                                      ;   ii = ii + 1 
    110110      WRITE(ldtxt(ii),*) '         LY [km]: ', zly                                                      ;   ii = ii + 1 
     
    115115      ! 
    116116      !                             ! Set the lateral boundary condition of the global domain 
    117       kperio = 7                    ! ICEDYN configuration : bi-periodic basin 
     117      kperio = 7                    ! ICE_AGRIF configuration : bi-periodic basin 
    118118#if defined key_agrif 
    119119      IF( .NOT. Agrif_Root() ) THEN 
     
    124124      WRITE(ldtxt(ii),*) '   '                                                                          ;   ii = ii + 1 
    125125      WRITE(ldtxt(ii),*) '   Lateral boundary condition of the global domain'                           ;   ii = ii + 1 
    126       WRITE(ldtxt(ii),*) '      ICEDYN : bi-periodic basin               jperio = ', kperio             ;   ii = ii + 1 
     126      WRITE(ldtxt(ii),*) '      ICE_AGRIF : bi-periodic basin               jperio = ', kperio          ;   ii = ii + 1 
    127127      ! 
    128128   END SUBROUTINE usr_def_nam 
  • NEMO/trunk/tests/ICE_AGRIF/MY_SRC/usrdef_sbc.F90

    r10074 r10516  
    33   !!                       ***  MODULE  usrdef_sbc  *** 
    44   !!  
    5    !!                      ===  ICEDYN configuration  === 
     5   !!                      ===  ICE_AGRIF configuration  === 
    66   !! 
    77   !! User defined :   surface forcing of a user configuration 
     
    1111 
    1212   !!---------------------------------------------------------------------- 
    13    !!   usr_def_sbc    : user defined surface bounday conditions in ICEDYN case 
     13   !!   usr_def_sbc    : user defined surface bounday conditions in ICE_AGRIF case 
    1414   !!---------------------------------------------------------------------- 
    1515   USE oce             ! ocean dynamics and tracers 
     
    4949      !!              condition, i.e. the momentum, heat and freshwater fluxes. 
    5050      !! 
    51       !! ** Method  :   all 0 fields, for ICEDYN case 
     51      !! ** Method  :   all 0 fields, for ICE_AGRIF case 
    5252      !!                CAUTION : never mask the surface stress field ! 
    5353      !! 
     
    6161      IF( kt == nit000 ) THEN 
    6262         ! 
    63          IF(lwp)   WRITE(numout,*)' usrdef_sbc_oce : ICEDYN case: NO surface forcing' 
     63         IF(lwp)   WRITE(numout,*)' usrdef_sbc_oce : ICE_AGRIF case: NO surface forcing' 
    6464         ! --- oce variables --- ! 
    6565         utau(:,:) = 0._wp  
     
    9191      INTEGER, INTENT(in) ::   kt   ! ocean time step 
    9292      !!--------------------------------------------------------------------- 
    93       IF( kt==nit000 .AND. lwp)   WRITE(numout,*)' usrdef_sbc_ice : ICEDYN case: constant stress forcing' 
     93      IF( kt==nit000 .AND. lwp)   WRITE(numout,*)' usrdef_sbc_ice : ICE_AGRIF case: constant stress forcing' 
    9494      ! 
    9595      utau_ice(:,:) = 1.3_wp   ! <=> 0.5 m/s 
     
    112112      !!--------------------------------------------------------------------- 
    113113      ! 
    114       IF( kt==nit000 .AND. lwp)   WRITE(numout,*)' usrdef_sbc_ice : ICEDYN case: NO flux forcing' 
     114      IF( kt==nit000 .AND. lwp)   WRITE(numout,*)' usrdef_sbc_ice : ICE_AGRIF case: NO flux forcing' 
    115115      ! 
    116116      ! ocean variables (renaming) 
  • NEMO/trunk/tests/ICE_AGRIF/MY_SRC/usrdef_zgr.F90

    r10074 r10516  
    33   !!                   ***  MODULE  usrdef_zgr  *** 
    44   !! 
    5    !!                   ===      ICEDYN case     === 
     5   !!                   ===      ICE_AGRIF case     === 
    66   !! 
    77   !! Ocean domain : user defined vertical coordinate system  
     
    5959      ! 
    6060      IF(lwp) WRITE(numout,*) 
    61       IF(lwp) WRITE(numout,*) 'usr_def_zgr : ICEDYN configuration (slab ocean - advection of an ice patch in a biperiodic square box domain)' 
    62       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
     61      IF(lwp) WRITE(numout,*) 'usr_def_zgr : ICE_AGRIF configuration ' 
     62      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   (slab ocean - advection of an ice patch in a biperiodic square box domain)' 
    6363      ! 
    6464      ! 
    65       ! type of vertical coordinate  ==>>>   here ICEDYN : slab ocean always 
     65      ! type of vertical coordinate  ==>>>   here ICE_AGRIF : slab ocean always 
    6666      ! --------------------------- 
    6767      ld_zco    = .TRUE.       ! z-full-step coordinate 
  • NEMO/trunk/tests/demo_cfgs.txt

    r10515 r10516  
    33LOCK_EXCHANGE OCE 
    44OVERFLOW OCE 
    5 ICEDYN OCE NST SAS ICE 
     5ICE_AGRIF OCE NST SAS ICE 
    66ICE_ADV1D OCE SAS ICE 
    77ICE_ADV2D OCE SAS ICE 
Note: See TracChangeset for help on using the changeset viewer.