Changeset 8706


Ignore:
Timestamp:
2017-11-15T12:30:49+01:00 (3 years ago)
Author:
flavoni
Message:

ENHANCE09: update SAS_BIPER

Location:
branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/CONFIG/SHARED/namelist_ref

    r8685 r8706  
    183183   ln_usr      = .false.   !  user defined formulation                  (T => check usrdef_sbc) 
    184184   ln_flx      = .false.   !  flux formulation                          (T => fill namsbc_flx ) 
    185    ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk ) 
     185   ln_blk      = .false.    !  Bulk formulation                          (T => fill namsbc_blk ) 
    186186                     ! Type of coupling (Ocean/Ice/Atmosphere) : 
    187187   ln_cpl      = .false.   !  atmosphere coupled   formulation          ( requires key_oasis3 ) 
  • branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/1_namelist_cfg

    r8637 r8706  
    1414   cn_exp      =  "SASBIPER"  !  experience name 
    1515   nn_it000    =       1      !  first time step 
    16    nn_itend    =    30300     !  last  time step (std 5475) 
     16   nn_itend    =    3030      !  last  time step (std 5475) 
     17   nn_istate   =       1      !  output the initial state (1) or not (0) 
    1718/ 
    1819!----------------------------------------------------------------------- 
     
    2627&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    2728!----------------------------------------------------------------------- 
    28    ln_linssh   = .false.   !  =T  linear free surface  ==>>  model level are fixed in time 
    29    ! 
    30    nn_msh      =    1      !  create (>0) a mesh file or not (=0) 
     29   ln_linssh   = .true.   !  =T  linear free surface  ==>>  model level are fixed in time 
     30   ! 
     31   nn_msh      =    1     !  create (>0) a mesh file or not (=0) 
    3132   ! 
    3233   rn_rdt      = 400.     !  time step for the dynamics (and tracer if nn_acc=0) 
     
    4445!----------------------------------------------------------------------- 
    4546   ln_usr      = .true.   !  user defined formulation                  (T => check usrdef_sbc) 
     47   ln_blk      = .false.  !  Bulk formulation                          (T => fill namsbc_blk ) 
    4648   nn_fsbc     =  1       !  frequency of surface boundary condition computation 
    4749/ 
     
    5153/ 
    5254!----------------------------------------------------------------------- 
     55&namsbc_sas    !   Stand Alone Surface boundary condition 
     56!----------------------------------------------------------------------- 
     57   l_sasread   = .false.   ! Read fields in a file if .TRUE. , or initialize to 0. in sbcssm.F90 if .FALSE. 
     58/ 
     59!----------------------------------------------------------------------- 
    5360&namtra_qsr    !   penetrative solar radiation 
    5461!----------------------------------------------------------------------- 
     
    6976&namberg       !   iceberg parameters 
    7077!----------------------------------------------------------------------- 
     78   ln_bergdia               = .false.               ! Calculate budgets 
    7179/ 
    7280!----------------------------------------------------------------------- 
     
    7785/ 
    7886!----------------------------------------------------------------------- 
    79 &nambfr        !   bottom friction 
    80 !----------------------------------------------------------------------- 
     87&namagrif      !  AGRIF zoom                                            ("key_agrif") 
     88!----------------------------------------------------------------------- 
     89   nn_cln_update =    1    !  baroclinic update frequency 
     90/ 
     91!----------------------------------------------------------------------- 
     92&namdrg            !   top/bottom drag coefficient                      (default: NO selection) 
     93!----------------------------------------------------------------------- 
     94   ln_NONE    = .true.     !  free-slip       : Cd = 0                   
    8195/ 
    8296!----------------------------------------------------------------------- 
     
    94108/ 
    95109!----------------------------------------------------------------------- 
    96 &namtra_adv    !   advection scheme for tracer 
     110&namtra_adv    !   advection scheme for tracer                          (default: NO selection) 
    97111!----------------------------------------------------------------------- 
    98112   ln_traadv_fct =  .true.    !  FCT scheme 
     
    138152&namtra_ldfeiv !   eddy induced velocity param. 
    139153!---------------------------------------------------------------------------------- 
     154   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     155   ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
     156   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     157   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     158   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     159   !                                !   =  0           constant  
     160   !                                !   = 10 F(k)      =ldf_c1d  
     161   !                                !   = 20 F(i,j)    =ldf_c2d  
     162   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     163   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    140164/ 
    141165!----------------------------------------------------------------------- 
     
    147171!----------------------------------------------------------------------- 
    148172   ln_dynadv_NONE= .false. !  linear dynamics (no momentum advection) 
    149    ln_dynadv_vec = .true.  !  vector form (T) or flux form (F) 
    150    nn_dynkeg     = 0       ! scheme for grad(KE): =0   C2  ;  =1   Hollingsworth correction 
     173   ln_dynadv_vec = .true. !  vector form - 2nd centered scheme 
     174     nn_dynkeg     = 0        ! grad(KE) scheme: =0   C2  ;  =1   Hollingsworth correction 
    151175   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
    152    ln_dynadv_ubs = .true. !  flux form - 3rd order UBS      scheme 
     176   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
    153177/ 
    154178!----------------------------------------------------------------------- 
     
    198222/ 
    199223!----------------------------------------------------------------------- 
    200 &namzdf        !   vertical physics 
    201 !----------------------------------------------------------------------- 
    202 / 
    203 !----------------------------------------------------------------------- 
    204 &nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi) 
     224&namzdf        !   vertical physics                                     (default: NO selection) 
     225!----------------------------------------------------------------------- 
     226   !                       ! type of vertical closure 
     227   ln_zdfcst   = .true.       !  constant mixing 
     228/ 
     229!----------------------------------------------------------------------- 
     230&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T) 
     231!----------------------------------------------------------------------- 
     232/ 
     233!----------------------------------------------------------------------- 
     234&namzdf_iwm    !   tidal mixing parameterization                        (ln_zdfiwm =T) 
     235!----------------------------------------------------------------------- 
     236/ 
     237!----------------------------------------------------------------------- 
     238&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi") 
    205239!----------------------------------------------------------------------- 
    206240/ 
  • branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/1_namelist_ice_cfg

    r8696 r8706  
    2222&nampar         !   Generic parameters 
    2323!------------------------------------------------------------------------------ 
    24    jpl              =   1             !  number of ice  categories 
    25    nlay_i           =   1             !  number of ice  layers 
    2624   ln_icedyn        = .true.          !  ice dynamics (T) or not (F) 
    2725   ln_icethd        = .false.         !  ice thermo   (T) or not (F) 
     
    8482&namini         !   Ice initialization 
    8583!------------------------------------------------------------------------------ 
     84   ln_iceini        = .false.          !  activate ice initialization (T) or not (F) 
    8685/ 
    8786!------------------------------------------------------------------------------ 
  • branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/file_def_nemo-lim.xml

    r7821 r8706  
    1717   <file id="file1" name_suffix="_icemod" description="ice variables" enabled=".true." > 
    1818 
    19         <field field_ref="snowthic_cea"     name="snthic"     long_name="surface_snow_thickness"   /> 
    2019        <field field_ref="snowvol"          name="snvolu" /> 
    21         <field field_ref="isnowhc"          name="snheco" /> 
    22  
    23         <field field_ref="icethic_cea"      name="sithic"     long_name="sea_ice_thickness"        /> 
     20        <field field_ref="icethic"          name="sithic" /> 
    2421        <field field_ref="icevolu"          name="sivolu" /> 
    25         <field field_ref="iceconc"          name="siconc"  /> 
     22        <field field_ref="iceconc"          name="siconc" /> 
    2623        <field field_ref="micesalt"         name="sisali" /> 
    2724        <field field_ref="micet"            name="sitemp" /> 
  • branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/namelist_cfg

    r8696 r8706  
    1515   nn_it000    =       1      !  first time step 
    1616   nn_itend    =    1010      !  last  time step (std 5475) 
     17   nn_istate   =       1      !  output the initial state (1) or not (0) 
    1718/ 
    1819!----------------------------------------------------------------------- 
     
    5253/ 
    5354!----------------------------------------------------------------------- 
     55&namsbc_sas    !   Stand Alone Surface boundary condition 
     56!----------------------------------------------------------------------- 
     57   l_sasread   = .false.   ! Read fields in a file if .TRUE. , or initialize to 0. in sbcssm.F90 if .FALSE. 
     58/ 
     59!----------------------------------------------------------------------- 
    5460&namtra_qsr    !   penetrative solar radiation 
    5561!----------------------------------------------------------------------- 
     
    6874/ 
    6975!----------------------------------------------------------------------- 
    70 &namsbc_sas    !   Stand Alone Surface boundary condition 
    71 !----------------------------------------------------------------------- 
    72    l_sasread   = .false.   ! Read fields in a file if .TRUE. , or initialize to 0. in sbcssm.F90 if .FALSE. 
    73 !----------------------------------------------------------------------- 
    7476&namberg       !   iceberg parameters 
    7577!----------------------------------------------------------------------- 
     78   ln_bergdia               = .false.               ! Calculate budgets 
    7679/ 
    7780!----------------------------------------------------------------------- 
     
    8083   rn_shlat    =    0.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    8184                           !  free slip  !   partial slip  !   no slip   ! strong slip 
     85/ 
     86!----------------------------------------------------------------------- 
     87&namagrif      !  AGRIF zoom                                            ("key_agrif") 
     88!----------------------------------------------------------------------- 
     89   nn_cln_update =    1    !  baroclinic update frequency 
    8290/ 
    8391!----------------------------------------------------------------------- 
  • branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/EXP00/namelist_ice_cfg

    r8637 r8706  
    2222&nampar         !   Generic parameters 
    2323!------------------------------------------------------------------------------ 
    24    jpl              =   1             !  number of ice  categories 
    25    nlay_i           =   1             !  number of ice  layers 
    2624   ln_icedyn        = .true.          !  ice dynamics (T) or not (F) 
    2725   ln_icethd        = .false.         !  ice thermo   (T) or not (F) 
     
    8482&namini         !   Ice initialization 
    8583!------------------------------------------------------------------------------ 
     84   ln_iceini        = .true.          !  activate ice initialization (T) or not (F) 
     85   ln_iceini_file   = .true.          !  netcdf file provided for initialization (T) or not (F) 
     86 
     87   sn_hti = 'initice'                 , -12 ,'hti'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     88   sn_hts = 'initice'                 , -12 ,'hts'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     89   sn_ati = 'initice'                 , -12 ,'ati'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     90   sn_tsu = 'initice'                 , -12 ,'tsu'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     91   sn_tmi = 'initice'                 , -12 ,'tmi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     92   sn_smi = 'initice'                 , -12 ,'smi'   ,  .false.  , .true., 'yearly'  , '' , '', '' 
     93   cn_dir='./' 
    8694/ 
    8795!------------------------------------------------------------------------------ 
  • branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/CONFIG/TEST_CASES/SAS_BIPER/MY_SRC/usrdef_zgr.F90

    r7355 r8706  
    6262      ! 
    6363      IF(lwp) WRITE(numout,*) 
    64       IF(lwp) WRITE(numout,*) 'usr_def_zgr : LOCK_EXCHANGE configuration (z-coordinate closed box ocean without cavities)' 
     64      IF(lwp) WRITE(numout,*) 'usr_def_zgr : SAS_BIPER configuration (slab ocean - advection of an ice patch in a biperiodic square box domain)' 
    6565      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 
    6666      ! 
     
    6868      ! type of vertical coordinate 
    6969      ! --------------------------- 
    70       ! set in usrdef_nam.F90 by reading the namusr_def namelist only ln_zco 
    71       ln_zco    = .TRUE.       ! z-partial-step coordinate 
    72       ln_zps    = .FALSE.      ! z-partial-step coordinate 
    73       ln_sco    = .FALSE.      ! s-coordinate 
     70      ld_zco    = .TRUE.       ! z-full-step coordinate 
     71      ld_zps    = .FALSE.      ! z-partial-step coordinate 
     72      ld_sco    = .FALSE.      ! s-coordinate 
    7473      ld_isfcav = .FALSE.      ! ISF Ice Shelves Flag 
    7574      ! 
  • branches/2017/dev_r7881_ENHANCE09_RK3/NEMOGCM/NEMO/LIM_SRC_3/icewri.F90

    r8637 r8706  
    330330      CALL histdef( kid, "vfxsum", "Ice surface melt"        , "m/s"    ,   & 
    331331      &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
     332      CALL histdef( kid, "vfxpnd", "Ice melt ponds flux"     , "m/s"    ,   & 
     333      &      jpi, jpj, kh_i, 1, 1, 1, -99, 32, "inst(x)", rdt, rdt ) 
    332334 
    333335      CALL histdef( kid, "sithicat", "Ice thickness"         , "m"      ,   & 
Note: See TracChangeset for help on using the changeset viewer.