Changeset 3614


Ignore:
Timestamp:
2012-11-20T16:19:12+01:00 (8 years ago)
Author:
acc
Message:

Branch dev_NOC_2012_r3555. #1006. Step 6: Minor code changes and updated namelists to enable successful SETTE testing

Location:
branches/2012/dev_NOC_2012_rev3555/NEMOGCM
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_NOC_2012_rev3555/NEMOGCM/CONFIG/GYRE/EXP00/namelist

    r3306 r3614  
    315315   rn_alphc    =    0.65   !  compute albedo between two extremes values  
    316316   rn_alphdi   =    0.72   !  (Pyane, 1972) 
     317/ 
     318!----------------------------------------------------------------------- 
     319&namberg       !   iceberg parameters 
     320!----------------------------------------------------------------------- 
     321      ln_icebergs              = .false. 
     322      ln_bergdia               = .true.               ! Calculate budgets 
     323      nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0 
     324      nn_verbose_write         = 15                   ! Timesteps between verbose messages 
     325      nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage 
     326                                                      ! Initial mass required for an iceberg of each class 
     327      rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 
     328                                                      ! Proportion of calving mass to apportion to each class 
     329      rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
     330                                                      ! Ratio between effective and real iceberg mass (non-dim) 
     331                                                      ! i.e. number of icebergs represented at a point 
     332      rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
     333                                                      ! thickness of newly calved bergs (m) 
     334      rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 
     335      rn_rho_bergs             = 850.                 ! Density of icebergs 
     336      rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs 
     337      ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics 
     338      rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits 
     339      rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
     340      ln_passive_mode          = .false.              ! iceberg - ocean decoupling 
     341      nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no) 
     342                                                      ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 
     343      rn_test_box              = 108.0,  116.0, -66.0, -58.0 
     344      rn_speed_limit           = 0.                   ! CFL speed limit for a berg 
     345 
     346               ! filename ! freq (hours) ! variable ! time interp. ! clim  !'yearly' or ! weights  ! rotation ! 
     347               !          ! (<0  months) !   name   !  (logical)   ! (T/F) ! 'monthly'  ! filename ! pairing  ! 
     348      sn_icb =  'calving' ,     -1       , 'calvingmask',  .true.      , .true., 'yearly'   , ' '      , ' ' 
     349 
     350      cn_dir = './' 
    317351/ 
    318352 
  • branches/2012/dev_NOC_2012_rev3555/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist

    r3306 r3614  
    314314   rn_alphc    =    0.65   !  compute albedo between two extremes values  
    315315   rn_alphdi   =    0.72   !  (Pyane, 1972) 
     316/ 
     317!----------------------------------------------------------------------- 
     318&namberg       !   iceberg parameters 
     319!----------------------------------------------------------------------- 
     320      ln_icebergs              = .false. 
     321      ln_bergdia               = .true.               ! Calculate budgets 
     322      nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0 
     323      nn_verbose_write         = 15                   ! Timesteps between verbose messages 
     324      nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage 
     325                                                      ! Initial mass required for an iceberg of each class 
     326      rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 
     327                                                      ! Proportion of calving mass to apportion to each class 
     328      rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
     329                                                      ! Ratio between effective and real iceberg mass (non-dim) 
     330                                                      ! i.e. number of icebergs represented at a point 
     331      rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
     332                                                      ! thickness of newly calved bergs (m) 
     333      rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 
     334      rn_rho_bergs             = 850.                 ! Density of icebergs 
     335      rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs 
     336      ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics 
     337      rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits 
     338      rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
     339      ln_passive_mode          = .false.              ! iceberg - ocean decoupling 
     340      nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no) 
     341                                                      ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 
     342      rn_test_box              = 108.0,  116.0, -66.0, -58.0 
     343      rn_speed_limit           = 0.                   ! CFL speed limit for a berg 
     344 
     345               ! filename ! freq (hours) ! variable ! time interp. ! clim  !'yearly' or ! weights  ! rotation ! 
     346               !          ! (<0  months) !   name   !  (logical)   ! (T/F) ! 'monthly'  ! filename ! pairing  ! 
     347      sn_icb =  'calving' ,     -1       , 'calvingmask',  .true.      , .true., 'yearly'   , ' '      , ' ' 
     348 
     349      cn_dir = './' 
    316350/ 
    317351 
  • branches/2012/dev_NOC_2012_rev3555/NEMOGCM/CONFIG/ORCA2_SAS_LIM/EXP00/namelist

    r3362 r3614  
    330330   rn_alphc    =    0.65   !  compute albedo between two extremes values  
    331331   rn_alphdi   =    0.72   !  (Pyane, 1972) 
     332/ 
     333!----------------------------------------------------------------------- 
     334&namberg       !   iceberg parameters 
     335!----------------------------------------------------------------------- 
     336      ln_icebergs              = .false. 
     337      ln_bergdia               = .true.               ! Calculate budgets 
     338      nn_verbose_level         = 1                    ! Turn on more verbose output if level > 0 
     339      nn_verbose_write         = 15                   ! Timesteps between verbose messages 
     340      nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage 
     341                                                      ! Initial mass required for an iceberg of each class 
     342      rn_initial_mass          = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 
     343                                                      ! Proportion of calving mass to apportion to each class 
     344      rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
     345                                                      ! Ratio between effective and real iceberg mass (non-dim) 
     346                                                      ! i.e. number of icebergs represented at a point 
     347      rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
     348                                                      ! thickness of newly calved bergs (m) 
     349      rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 
     350      rn_rho_bergs             = 850.                 ! Density of icebergs 
     351      rn_LoW_ratio             = 1.5                  ! Initial ratio L/W for newly calved icebergs 
     352      ln_operator_splitting    = .true.               ! Use first order operator splitting for thermodynamics 
     353      rn_bits_erosion_fraction = 0.                   ! Fraction of erosion melt flux to divert to bergy bits 
     354      rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
     355      ln_passive_mode          = .false.              ! iceberg - ocean decoupling 
     356      nn_test_icebergs         =  10                  ! Create test icebergs of this class (-1 = no) 
     357                                                      ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2) 
     358      rn_test_box              = 108.0,  116.0, -66.0, -58.0 
     359      rn_speed_limit           = 0.                   ! CFL speed limit for a berg 
     360 
     361               ! filename ! freq (hours) ! variable ! time interp. ! clim  !'yearly' or ! weights  ! rotation ! 
     362               !          ! (<0  months) !   name   !  (logical)   ! (T/F) ! 'monthly'  ! filename ! pairing  ! 
     363      sn_icb =  'calving' ,     -1       , 'calvingmask',  .true.      , .true., 'yearly'   , ' '      , ' ' 
     364 
     365      cn_dir = './' 
    332366/ 
    333367 
  • branches/2012/dev_NOC_2012_rev3555/NEMOGCM/NEMO/OPA_SRC/ICB/icb_oce.F90

    r3383 r3614  
    9090   REAL(wp), PUBLIC, DIMENSION(:,:), ALLOCATABLE ::   ff_e 
    9191   REAL(wp), PUBLIC, DIMENSION(:,:), ALLOCATABLE ::   ua_e, va_e 
     92   REAL(wp), PUBLIC, DIMENSION(:,:), ALLOCATABLE ::   ssh_e 
    9293#if defined key_lim2 || defined key_lim3 
    9394   REAL(wp), PUBLIC, DIMENSION(:,:), ALLOCATABLE ::   ui_e, vi_e 
    94    REAL(wp), PUBLIC, DIMENSION(:,:), ALLOCATABLE ::   ssh_e 
    9595#endif 
    9696 
  • branches/2012/dev_NOC_2012_rev3555/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90

    r3607 r3614  
    2525   PRIVATE 
    2626 
    27    PUBLIC   sbc_ssm    ! routine called by step.F90 
     27   PUBLIC   sbc_ssm         ! routine called by step.F90 
     28   PUBLIC   sbc_ssm_init    ! routine called by sbcmod.F90 
    2829    
    2930   LOGICAL, SAVE  ::   l_ssm_mean = .FALSE.       ! keep track of whether means have been read 
     
    7172      ELSE 
    7273         !                                                ! ----------------------------------------------- ! 
    73          IF( kt == nit000 && .NOT. l_ssm_mean ) THEN      !   Initialisation: 1st time-step, no input means ! 
     74         IF( kt == nit000 .AND. .NOT. l_ssm_mean ) THEN   !   Initialisation: 1st time-step, no input means ! 
    7475            !                                             ! ----------------------------------------------- ! 
    7576            IF(lwp) WRITE(numout,*) 
     
    145146      !! ** Action  : - read parameters 
    146147      !!---------------------------------------------------------------------- 
    147       REAL(wp) ::   zcoef               ! local scalar 
     148      REAL(wp) ::   zcoef, zf_sbc       ! local scalar 
    148149      !!---------------------------------------------------------------------- 
    149150 
Note: See TracChangeset for help on using the changeset viewer.