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

Changeset 3369


Ignore:
Timestamp:
2012-04-26T09:23:19+02:00 (12 years ago)
Author:
sga
Message:

NEMO branch dev_r3337_NOCS10_ICB: bugfix to calculation of iceberg domain and small text changes

Location:
branches/2012/dev_r3337_NOCS10_ICB
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_r3337_NOCS10_ICB/DOC/TexFiles/Chapters/Chap_SBC.tex

    r3368 r3369  
    925925A non-zero value represents how many timesteps between writes of information into the output file. 
    926926These output files are in NETCDF format. 
    927 When \key{mpp\_mpi} is defined, each output file contains only those icebergs in the corresponding processor, 
    928 so care is needed to recreate data for individual icebergs. 
     927When \key{mpp\_mpi} is defined, each output file contains only those icebergs in the corresponding processor. 
     928Trajectory points are written out in the order of their parent iceberg in the model's "linked list" of icebergs. 
     929So care is needed to recreate data for individual icebergs, since its trajectory data may be spread across 
     930multiple files. 
    929931 
    930932 
  • branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist

    r3344 r3369  
    324324      nn_verbose_write         = 15                   ! Timesteps between verbose messages 
    325325      nn_sample_rate           = 1                    ! Timesteps between sampling for trajectory storage 
     326                                                      ! Initial mass required for an iceberg of each class 
    326327      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   
    327329      rn_distribution          = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 
    328330                                                      ! Ratio between effective and real iceberg mass (non-dim) 
     331                                                      ! i.e. number of icebergs represented at a point          
    329332      rn_mass_scaling          = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 
    330                                                       ! Total thickness of newly calved bergs (m) 
     333                                                      ! thickness of newly calved bergs (m) 
    331334      rn_initial_thickness     = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. 
    332335      rn_rho_bergs             = 850.                 ! Density of icebergs 
     
    336339      rn_sicn_shift            = 0.                   ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1) 
    337340      ln_passive_mode          = .false.              ! iceberg - ocean decoupling    
    338       nn_test_icebergs         =  10                  ! Create icebergs in absence of a calving file (-1 = no) 
     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) 
    339343      rn_test_box              = 108.0,  116.0, -66.0, -58.0 
    340344      rn_speed_limit           = 0.                   ! CFL speed limit for a berg    
  • branches/2012/dev_r3337_NOCS10_ICB/NEMOGCM/NEMO/OPA_SRC/ICB/icbini.F90

    r3359 r3369  
    193193         i2 = INT( i3/nicbpack ) 
    194194         i1 = i3 - i2*nicbpack 
    195          IF( i1 == nimpp+ji-1 ) THEN 
     195         i3 = INT( src_calving_hflx(ji,jj) ) 
     196         IF( i1 == nimpp+ji-1 .AND. i3 == narea ) THEN 
    196197            IF( nicbdi < 0 ) THEN 
    197198               nicbdi = ji 
     
    210211         i2 = INT( i3/nicbpack ) 
    211212         i1 = i3 - i2*nicbpack 
    212          IF( i2 == njmpp+jj-1 ) THEN 
     213         i3 = INT( src_calving_hflx(ji,jj) ) 
     214         IF( i2 == njmpp+jj-1 .AND. i3 == narea ) THEN 
    213215            IF( nicbdj < 0 ) THEN 
    214216               nicbdj = jj 
Note: See TracChangeset for help on using the changeset viewer.