Ignore:
Timestamp:
2013-11-15T09:55:10+01:00 (8 years ago)
Author:
cbricaud
Message:

put dev_r4022_MERCATOR5_PAPA1D changes in dev_MERCATOR_2013

Location:
branches/2013/dev_MERCATOR_2013/NEMOGCM
Files:
5 edited
2 copied

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERCATOR_2013/NEMOGCM/NEMO/OPA_SRC/DOM/dom_oce.F90

    r3851 r4211  
    3030   !                                              !!* Namelist namdom : time & space domain * 
    3131   INTEGER , PUBLIC ::   nn_bathy     =    0       !: = 0/1 ,compute/read the bathymetry file 
     32   REAL(wp), PUBLIC ::   rn_bathy     =    0.0_wp  !: depth of flat bottom (active if nn_bathy=0; if =0 depth=jpkm1) 
    3233   REAL(wp), PUBLIC ::   rn_hmin      =   -3.0_wp  !: minimum ocean depth (>0) or minimum number of ocean levels (<0) 
    3334   REAL(wp), PUBLIC ::   rn_e3zps_min =    5.0_wp  !: miminum thickness for partial steps (meters) 
  • branches/2013/dev_MERCATOR_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

    r3764 r4211  
    124124         &             nn_it000, nn_itend  , nn_date0    , nn_leapy     , nn_istate , nn_stock ,   & 
    125125         &             nn_write, ln_dimgnnn, ln_mskland  , ln_clobber   , nn_chunksz 
    126       NAMELIST/namdom/ nn_bathy , rn_e3zps_min, rn_e3zps_rat, nn_msh    , rn_hmin,   & 
    127          &             nn_acc   , rn_atfp     , rn_rdt      , rn_rdtmin ,            & 
     126      NAMELIST/namdom/ nn_bathy, rn_bathy , rn_e3zps_min, rn_e3zps_rat, nn_msh, rn_hmin,   & 
     127         &             nn_acc   , rn_atfp     , rn_rdt      , rn_rdtmin ,                  & 
    128128         &             rn_rdtmax, rn_rdth     , nn_baro     , nn_closea 
    129129      NAMELIST/namcla/ nn_cla 
     
    207207         WRITE(numout,*) '   Namelist namdom : space & time domain' 
    208208         WRITE(numout,*) '      flag read/compute bathymetry      nn_bathy     = ', nn_bathy 
     209         WRITE(numout,*) '      Depth (if =0 bathy=jpkm1)         rn_bathy     = ', rn_bathy 
    209210         WRITE(numout,*) '      min depth of the ocean    (>0) or    rn_hmin   = ', rn_hmin 
    210211         WRITE(numout,*) '      min number of ocean level (<0)       ' 
  • branches/2013/dev_MERCATOR_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r4006 r4211  
    368368            IF(lwp) WRITE(numout,*) 
    369369            IF(lwp) WRITE(numout,*) '         bathymetry field: flat basin' 
    370             idta(:,:) = jpkm1                            ! before last level 
    371             zdta(:,:) = gdepw_0(jpk)                     ! last w-point depth 
    372             h_oce     = gdepw_0(jpk) 
     370            IF( rn_bathy > 0.01 ) THEN  
     371               IF(lwp) WRITE(numout,*) '         Depth = rn_bathy read in namelist' 
     372               zdta(:,:) = rn_bathy 
     373               IF( ln_sco ) THEN                                   ! s-coordinate (zsc       ): idta()=jpk 
     374                  idta(:,:) = jpkm1 
     375               ELSE                                                ! z-coordinate (zco or zps): step-like topography 
     376                  idta(:,:) = jpkm1 
     377                  DO jk = 1, jpkm1 
     378                     WHERE( gdept_0(jk) < zdta(:,:) .AND. zdta(:,:) <= gdept_0(jk+1) )   idta(:,:) = jk 
     379                  END DO 
     380               ENDIF 
     381            ELSE 
     382               IF(lwp) WRITE(numout,*) '         Depth = depthw(jpkm1)' 
     383               idta(:,:) = jpkm1                            ! before last level 
     384               zdta(:,:) = gdepw_0(jpk)                     ! last w-point depth 
     385               h_oce     = gdepw_0(jpk) 
     386            ENDIF 
    373387         ELSE                                         ! bump centered in the basin 
    374388            IF(lwp) WRITE(numout,*) 
  • branches/2013/dev_MERCATOR_2013/NEMOGCM/NEMO/OPA_SRC/DOM/istate.F90

    r3764 r4211  
    180180   END SUBROUTINE istate_t_s 
    181181 
    182  
    183182   SUBROUTINE istate_eel 
    184183      !!---------------------------------------------------------------------- 
     
    404403 
    405404   END SUBROUTINE istate_gyre 
    406  
    407405 
    408406   SUBROUTINE istate_uvg 
  • branches/2013/dev_MERCATOR_2013/NEMOGCM/NEMO/OPA_SRC/par_oce.F90

    r3294 r4211  
    8686   !!--------------------------------------------------------------------- 
    8787#             include "par_AMM_12km.h90" 
     88#elif defined key_c1d_papa 
     89   !!--------------------------------------------------------------------- 
     90   !!   'key_c1d_papa':                1D conf. based on PAPA station  
     91   !!--------------------------------------------------------------------- 
     92#             include "par_C1D_PAPA.h90" 
    8893#else 
    8994   !!--------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.