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 5067 for branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC – NEMO

Ignore:
Timestamp:
2015-02-06T19:12:57+01:00 (9 years ago)
Author:
clem
Message:

LIM3 change all namelist names to fit with NEMO convention

Location:
branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/DOM/dom_oce.F90

    r4990 r5067  
    162162   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::  gphit, gphiu   !: latitude  of t-, u-, v- and f-points (degre) 
    163163   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::  gphiv, gphif   !: 
    164    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, TARGET, DIMENSION(:,:) ::  e1t, e2t       !: horizontal scale factors at t-point (m) 
    165    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, TARGET, DIMENSION(:,:) ::  e1u, e2u       !: horizontal scale factors at u-point (m) 
    166    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, TARGET, DIMENSION(:,:) ::  e1v, e2v       !: horizontal scale factors at v-point (m) 
    167    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, TARGET, DIMENSION(:,:) ::  e1f, e2f       !: horizontal scale factors at f-point (m) 
     164   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, TARGET, DIMENSION(:,:) ::  e1t, e2t, r1_e1t, r1_e2t   !: horizontal scale factors and inverse at t-point (m) 
     165   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, TARGET, DIMENSION(:,:) ::  e1u, e2u, r1_e1u, r1_e2u   !: horizontal scale factors and inverse at u-point (m) 
     166   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, TARGET, DIMENSION(:,:) ::  e1v, e2v, r1_e1v, r1_e2v   !: horizontal scale factors and inverse at v-point (m) 
     167   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, TARGET, DIMENSION(:,:) ::  e1f, e2f, r1_e1f, r1_e2f   !: horizontal scale factors and inverse at f-point (m) 
    168168   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::  e1e2t          !: surface at t-point (m2) 
    169169   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::  ff             !: coriolis factor (2.*omega*sin(yphi) ) (s-1) 
    … …  
    345345         &      tpol(jpiglo)  , fpol(jpiglo)                               , STAT=ierr(2) ) 
    346346         ! 
    347       ALLOCATE( glamt(jpi,jpj) , gphit(jpi,jpj) , e1t(jpi,jpj) , e2t(jpi,jpj) ,                      &  
    348          &      glamu(jpi,jpj) , gphiu(jpi,jpj) , e1u(jpi,jpj) , e2u(jpi,jpj) ,                      &   
    349          &      glamv(jpi,jpj) , gphiv(jpi,jpj) , e1v(jpi,jpj) , e2v(jpi,jpj) , e1e2t(jpi,jpj) ,     &   
    350          &      glamf(jpi,jpj) , gphif(jpi,jpj) , e1f(jpi,jpj) , e2f(jpi,jpj) , ff   (jpi,jpj) , STAT=ierr(3) )      
     347      ALLOCATE( glamt(jpi,jpj) , gphit(jpi,jpj) , e1t(jpi,jpj) , e2t(jpi,jpj) , r1_e1t(jpi,jpj) , r1_e2t(jpi,jpj) ,   &  
     348         &      glamu(jpi,jpj) , gphiu(jpi,jpj) , e1u(jpi,jpj) , e2u(jpi,jpj) , r1_e1u(jpi,jpj) , r1_e2u(jpi,jpj) ,   &   
     349         &      glamv(jpi,jpj) , gphiv(jpi,jpj) , e1v(jpi,jpj) , e2v(jpi,jpj) , r1_e1v(jpi,jpj) , r1_e2v(jpi,jpj) ,   &   
     350         &      glamf(jpi,jpj) , gphif(jpi,jpj) , e1f(jpi,jpj) , e2f(jpi,jpj) , r1_e1f(jpi,jpj) , r1_e2f(jpi,jpj) ,   & 
     351         &      e1e2t(jpi,jpj) , ff   (jpi,jpj) , STAT=ierr(3) )      
    351352         ! 
    352353      ALLOCATE( gdep3w_0(jpi,jpj,jpk) , e3v_0(jpi,jpj,jpk) , e3f_0 (jpi,jpj,jpk) ,                         & 

  • branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90

    r4990 r5067  
    471471      re2u_e1u(:,:) = e2u(:,:) / e1u(:,:) 
    472472      re1v_e2v(:,:) = e1v(:,:) / e2v(:,:) 
     473      r1_e1t  (:,:) = 1._wp    / e1t(:,:) 
     474      r1_e1u  (:,:) = 1._wp    / e1u(:,:) 
     475      r1_e1v  (:,:) = 1._wp    / e1v(:,:) 
     476      r1_e1f  (:,:) = 1._wp    / e1f(:,:) 
     477      r1_e2t  (:,:) = 1._wp    / e2t(:,:) 
     478      r1_e2u  (:,:) = 1._wp    / e2u(:,:) 
     479      r1_e2v  (:,:) = 1._wp    / e2v(:,:) 
     480      r1_e2f  (:,:) = 1._wp    / e2f(:,:) 
    473481 
    474482      ! Control printing : Grid informations (if not restart) 

  • branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90

    r5059 r5067  
    236236         CALL lim_thd( kt )                         ! Ice thermodynamics  
    237237          
    238          CALL lim_update2( kt )                     ! Global variables update 
     238         CALL lim_update2( kt )                     ! Corrections 
    239239         ! 
    240240         CALL lim_sbc_flx( kt )                     ! Update surface ocean mass, heat and salt fluxes 
    … …  
    286286      IF(lwm) CALL ctl_opn( numoni, 'output.namelist.ice', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, 1 ) 
    287287 
    288       CALL ice_run                     ! set some ice run parameters including jpl, nlay_i and nlay_s 
     288      CALL ice_run                     ! set some ice run parameters 
     289      ! 
    289290      !                                ! Allocate the ice arrays 
    290291      ierr =        ice_alloc        ()      ! ice variables 
    … …  
    303304         &     'use more ocean levels or less ice/snow layers/categories.' ) 
    304305      ! 
     306      CALL lim_itd_init                ! ice thickness distribution initialization 
    305307      ! 
    306308      CALL lim_thd_init                ! set ice thermodynics parameters 
    … …  
    309311      ! 
    310312      CALL lim_msh                     ! ice mesh initialization 
    311       ! 
    312       CALL lim_itd_init                ! ice thickness distribution initialization 
    313313      ! 
    314314      CALL lim_itd_me_init             ! ice thickness distribution initialization for mecanical deformation 
    … …  
    352352      INTEGER  ::   ios                 ! Local integer output status for namelist read 
    353353      NAMELIST/namicerun/ jpl, nlay_i, nlay_s, cn_icerst_in, cn_icerst_out,   & 
    354          &                ln_limdyn, amax, ln_nicep, ln_limdiahsb, ln_limdiaout 
     354         &                ln_limdyn, rn_amax, ln_nicep, ln_limdiahsb, ln_limdiaout 
    355355      !!------------------------------------------------------------------- 
    356356      !                     
    … …  
    369369         WRITE(numout,*) 'ice_run : ice share parameters for dynamics/advection/thermo of sea-ice' 
    370370         WRITE(numout,*) ' ~~~~~~' 
    371          WRITE(numout,*) '   number of ice  categories                               = ', jpl 
    372          WRITE(numout,*) '   number of ice  layers                                   = ', nlay_i 
    373          WRITE(numout,*) '   number of snow layers                                   = ', nlay_s 
     371         WRITE(numout,*) '   number of ice  categories                                               = ', jpl 
     372         WRITE(numout,*) '   number of ice  layers                                                   = ', nlay_i 
     373         WRITE(numout,*) '   number of snow layers                                                   = ', nlay_s 
    374374         WRITE(numout,*) '   switch for ice dynamics (1) or not (0)      ln_limdyn   = ', ln_limdyn 
    375          WRITE(numout,*) '   maximum ice concentration                               = ', amax  
     375         WRITE(numout,*) '   maximum ice concentration                               = ', rn_amax  
    376376         WRITE(numout,*) '   Several ice points in the ice or not in ocean.output    = ', ln_nicep 
    377377         WRITE(numout,*) '   Diagnose heat/salt budget or not          ln_limdiahsb  = ', ln_limdiahsb 
    … …  
    404404      !!------------------------------------------------------------------- 
    405405      INTEGER  ::   ios                 ! Local integer output status for namelist read 
    406       NAMELIST/namiceitd/ nn_itdshp, rn_itmean 
     406      NAMELIST/namiceitd/ nn_catbnd, rn_himean 
    407407      ! 
    408408      INTEGER  ::   jl                   ! dummy loop index 
    … …  
    425425         WRITE(numout,*) 'ice_itd : ice cat distribution' 
    426426         WRITE(numout,*) ' ~~~~~~' 
    427          WRITE(numout,*) '   shape of ice categories distribution                          nn_itdshp = ', nn_itdshp 
    428          WRITE(numout,*) '   mean ice thickness in the domain (only active if nn_itdshp=2) rn_itmean = ', rn_itmean 
     427         WRITE(numout,*) '   shape of ice categories distribution                          nn_catbnd = ', nn_catbnd 
     428         WRITE(numout,*) '   mean ice thickness in the domain (only active if nn_catbnd=2) rn_himean = ', rn_himean 
    429429      ENDIF 
    430430 
    … …  
    438438      hi_max(:) = 0._wp 
    439439 
    440       SELECT CASE ( nn_itdshp  )        
     440      SELECT CASE ( nn_catbnd  )        
    441441                                   !---------------------- 
    442442         CASE (1)                  ! tanh function (CICE) 
    … …  
    456456         zalpha = 0.05             ! exponent of the transform function 
    457457 
    458          zhmax  = 3.*rn_itmean 
     458         zhmax  = 3.*rn_himean 
    459459 
    460460         DO jl = 1, jpl  
Note: See TracChangeset for help on using the changeset viewer.