Ignore:
Timestamp:
2011-02-02T15:10:08+01:00 (10 years ago)
Author:
cetlod
Message:

Offline:rename opa.F90 to nemogcm.F90 as it was done in OPA_SRC and update domain initialisation

Location:
trunk/NEMOGCM/NEMO/OFF_SRC
Files:
1 added
1 deleted
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OFF_SRC/domain.F90

    r2528 r2574  
    2727   PUBLIC dom_init       ! called by opa.F90 
    2828 
    29    !                                    !!! ** Namelist namzgr_sco ** 
    30    REAL(wp) ::   rn_sbot_min =  300.     ! minimum depth of s-bottom surface (>0) (m) 
    31    REAL(wp) ::   rn_sbot_max = 5250.     ! maximum depth of s-bottom surface (= ocean depth) (>0) (m) 
    32    REAL(wp) ::   rn_theta    =    6.0    ! surface control parameter (0<=rn_theta<=20) 
    33    REAL(wp) ::   rn_thetb    =    0.75   ! bottom control parameter  (0<=rn_thetb<= 1) 
    34    REAL(wp) ::   rn_rmax     =    0.15   ! maximum cut-off r-value allowed (0<rn_rmax<1) 
    35    LOGICAL  ::   ln_s_sigma  = .false.   ! use hybrid s-sigma -coordinate & stretching function fssig1 (ln_sco=T) 
    36    REAL(wp) ::   rn_bb       =    0.8    ! stretching parameter for song and haidvogel stretching 
    37    !                                     ! ( rn_bb=0; top only, rn_bb =1; top and bottom) 
    38    REAL(wp) ::   rn_hc       = 150.      ! Critical depth for s-sigma coordinates 
    39  
    4029   !! * Substitutions 
    4130#  include "domzgr_substitute.h90" 
     
    10291         &             nn_it000, nn_itend  , nn_date0    , nn_leapy     , nn_istate , nn_stock ,   & 
    10392         &             nn_write, ln_dimgnnn, ln_mskland  , ln_clobber   , nn_chunksz 
    104       NAMELIST/namdom/ nn_bathy , rn_e3zps_min, rn_e3zps_rat, nn_msh   ,   & 
    105          &             nn_acc   , rn_atfp     , rn_rdt      , rn_rdtmin,   & 
     93      NAMELIST/namdom/ nn_bathy , rn_e3zps_min, rn_e3zps_rat, nn_msh    , rn_hmin,     & 
     94         &             nn_acc   , rn_atfp     , rn_rdt      , rn_rdtmin ,              & 
    10695         &             rn_rdtmax, rn_rdth     , nn_baro     , nn_closea 
    10796      NAMELIST/namcla/ nn_cla 
     97#if defined key_netcdf4 
     98      NAMELIST/namnc4/ nn_nchunks_i, nn_nchunks_j, nn_nchunks_k, ln_nc4zip 
     99#endif 
    108100      !!---------------------------------------------------------------------- 
    109101 
     
    185177         WRITE(numout,*) '   Namelist namdom : space & time domain' 
    186178         WRITE(numout,*) '      flag read/compute bathymetry      nn_bathy     = ', nn_bathy 
     179         WRITE(numout,*) '      min depth of the ocean    (>0) or    rn_hmin   = ', rn_hmin 
    187180         WRITE(numout,*) '      minimum thickness of partial      rn_e3zps_min = ', rn_e3zps_min, ' (m)' 
    188181         WRITE(numout,*) '         step level                     rn_e3zps_rat = ', rn_e3zps_rat 
     
    222215      ENDIF 
    223216 
    224       IF( lk_mpp_rep .AND. nn_cla /= 0 )   CALL ctl_stop( ' Reproductibility tests (nbit_cmp=1) require n_cla = 0' ) 
     217#if defined key_netcdf4 
     218      !                             ! NetCDF 4 case   ("key_netcdf4" defined) 
     219      REWIND( numnam )                    ! Namelist namnc4 : netcdf4 chunking parameters 
     220      READ  ( numnam, namnc4 ) 
     221      IF(lwp) THEN                        ! control print 
     222         WRITE(numout,*) 
     223         WRITE(numout,*) '   Namelist namnc4 - Netcdf4 chunking parameters' 
     224         WRITE(numout,*) '      number of chunks in i-dimension      nn_nchunks_i   = ', nn_nchunks_i 
     225         WRITE(numout,*) '      number of chunks in j-dimension      nn_nchunks_j   = ', nn_nchunks_j 
     226         WRITE(numout,*) '      number of chunks in k-dimension      nn_nchunks_k   = ', nn_nchunks_k 
     227         WRITE(numout,*) '      apply netcdf4/hdf5 chunking & compression ln_nc4zip = ', ln_nc4zip 
     228      ENDIF 
     229 
     230      ! Put the netcdf4 settings into a simple structure (snc4set, defined in in_out_manager module) 
     231      ! Note the chunk size in the unlimited (time) dimension will be fixed at 1 
     232      snc4set%ni   = nn_nchunks_i 
     233      snc4set%nj   = nn_nchunks_j 
     234      snc4set%nk   = nn_nchunks_k 
     235      snc4set%luse = ln_nc4zip 
     236#else 
     237      snc4set%luse = .FALSE.        ! No NetCDF 4 case 
     238#endif 
    225239      ! 
    226240   END SUBROUTINE dom_nam 
Note: See TracChangeset for help on using the changeset viewer.