Ignore:
Timestamp:
2013-09-09T15:43:56+02:00 (8 years ago)
Author:
clevy
Message:

Configuration setting, bugfixes for AGRIF, see ticket:#1074

Location:
branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/ARCH/CNRS/arch-macport_osx.fcm

    r3993 r4018  
    2727%NCDF_INC            -I/opt/local/include 
    2828%NCDF_LIB            -L/opt/local/lib -lnetcdf -lnetcdff 
    29 %XIOS_ROOT           /Users/clelod/Documents/WORK/XIOS 
     29%XIOS_ROOT           /Users/clairelevy/Documents/WORK/XIOS 
    3030%FC                mpif90  
    31 %FCFLAGS             -fdefault-real-8 -O3 -funroll-all-loops -fcray-pointer  
     31%FCFLAGS             -fdefault-real-8 -O3 -funroll-all-loops -fcray-pointer -ffree-line-length-none 
    3232%FFLAGS              %FCFLAGS 
    3333%LD                  %FC 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/1_namelist_cfg

    r3973 r4018  
    99   nn_stock    =   10950   !  frequency of creation of a restart file (modulo referenced to 1) 
    1010   nn_write    =   10950   !  frequency of write in the output file   (modulo referenced to nn_it000) 
     11   ln_clobber  = .true.    !  clobber (overwrite) an existing file 
    1112/ 
    1213!----------------------------------------------------------------------- 
    1314&namcfg        !   parameters of the configuration 
    1415!----------------------------------------------------------------------- 
    15    cp_cfg      =  "orca"               !  name of the configuration 
     16   cp_cfg      =  "default"               !  name of the configuration 
    1617   jp_cfg      =      -1               !  resolution of the configuration 
    1718   jpidta      =     182               !  1st lateral dimension ( >= jpi ) 
     
    2223   jpizoom     =       1               !  left bottom (i,j) indices of the zoom 
    2324   jpjzoom     =       1               !  in data domain indices 
    24    jperio      =       4               !  lateral cond. type (between 0 and 6) 
     25   jperio      =       0               !  lateral cond. type (between 0 and 6) 
    2526   jphgr_msh   =       0               !  type of horizontal mesh 
    2627   ppglam0     =  999999.0             !  longitude of first raw and column T-point (jphgr_msh = 1) 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/SHARED/1_namelist_ref

    r3973 r4018  
    4848&namcfg        !   parameters of the configuration 
    4949!----------------------------------------------------------------------- 
    50    cp_cfg      =  "default"            !  name of the configuration 
     50   cp_cfg      =  "Agu"                !  name of the configuration 
     51   cp_cfz      =  "Agulhas"            !  name of the zoom of configuration 
    5152   jp_cfg      =       0               !  resolution of the configuration 
    5253   jpidta      =      10               !  1st lateral dimension ( >= jpi ) 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/NST_SRC/agrif_user.F90

    r3973 r4018  
    3131      jpi     = ( jpiglo-2*jpreci + (jpni-1+0) ) / jpni + 2*jpreci 
    3232      jpj     = ( jpjglo-2*jprecj + (jpnj-1+0) ) / jpnj + 2*jprecj 
    33       jpk     = jpkdta 
    34       jpim1   = jpi-1 
    35       jpjm1   = jpj-1 
    36       jpkm1   = jpk-1                                         
    37       jpij    = jpi*jpj 
    3833      jpidta  = jpiglo 
    3934      jpjdta  = jpjglo 
     
    6358   IMPLICIT NONE 
    6459   !!---------------------------------------------------------------------- 
    65  
    6660   ! 0. Initializations 
    6761   !------------------- 
    68    IF( cp_cfg == 'orca' .AND. jp_cfg == 2 .OR. jp_cfg == 025 .OR. jp_cfg == 05 & 
     62   IF( cp_cfg == 'orca' ) then 
     63      IF ( jp_cfg == 2 .OR. jp_cfg == 025 .OR. jp_cfg == 05 & 
    6964  &                      .OR. jp_cfg == 4 ) THEN 
    70       jp_cfg = -1    ! set special value for jp_cfg on fine grids 
    71       cp_cfg = "default" 
     65         jp_cfg = -1    ! set special value for jp_cfg on fine grids 
     66         cp_cfg = "default" 
     67      ENDIF 
    7268   ENDIF 
    7369   ! Specific fine grid Initializations 
     
    7874   lk_bdy = .FALSE. 
    7975 
     76 
    8077   CALL nemo_init  ! Initializations of each fine grid 
     78 
    8179   CALL agrif_nemo_init 
    8280   CALL Agrif_InitValues_cont_dom 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90

    r3973 r4018  
    519519          
    520520#if defined key_agrif 
    521          IF cp_cfg == 'eel' .AND. jp_cfg == 6 ) THEN    ! for EEL6 configuration only 
     521         IF ( cp_cfg == 'eel' .AND. jp_cfg == 6 ) THEN    ! for EEL6 configuration only 
    522522            IF( .NOT. Agrif_Root() ) THEN 
    523523              zphi0 = ppgphi0 - FLOAT( Agrif_Parent(jpjglo)/2)*Agrif_Parent(ppe2_m) / (ra * rad) 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r3989 r4018  
    366366      IF(lwp) WRITE(numout,*) '    zgr_bat : defines level and meter bathymetry' 
    367367      IF(lwp) WRITE(numout,*) '    ~~~~~~~' 
    368  
    369368      !                                               ! ================== !  
    370369      IF( ntopo == 0 .OR. ntopo == -1 ) THEN          !   defined by hand  ! 
     
    613612      INTEGER ::   icompt, ibtest, ikmax         ! temporary integers 
    614613      REAL(wp), POINTER, DIMENSION(:,:) ::  zbathy 
     614 
    615615      !!---------------------------------------------------------------------- 
    616616      ! 
     
    622622      IF(lwp) WRITE(numout,*) '    zgr_bat_ctl : check the bathymetry' 
    623623      IF(lwp) WRITE(numout,*) '    ~~~~~~~~~~~' 
    624  
    625624      !                                          ! Suppress isolated ocean grid points 
    626625      IF(lwp) WRITE(numout,*) 
     
    656655         mbathy(:,:) = INT( zbathy(:,:) ) 
    657656      ENDIF 
    658  
    659657      !                                          ! East-west cyclic boundary conditions 
    660658      IF( nperio == 0 ) THEN 
     
    687685         !         STOP 'dom_mba' 
    688686      ENDIF 
    689  
    690687      !  Boundary condition on mbathy 
    691688      IF( .NOT.lk_mpp ) THEN  
     
    696693         mbathy(:,:) = INT( zbathy(:,:) ) 
    697694      ENDIF 
    698  
    699695      ! Number of ocean level inferior or equal to jpkm1 
    700696      ikmax = 0 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/ICB/icbini.F90

    r3954 r4018  
    370370         WRITE(numout,*) ' namelist namberg not read' 
    371371      ENDIF 
    372       ln_iceberg=.false.       
     372      ln_icebergs = .false.       
    373373#endif    
    374374      IF( .NOT. ln_icebergs ) THEN   ! no icebergs 
  • branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90

    r3993 r4018  
    229229         &             nn_bench, nn_timing 
    230230      NAMELIST/namcfg/ cp_cfg, cp_cfz, jp_cfg, jpidta, jpjdta, jpkdta, jpiglo, jpjglo, & 
    231          &             jpizoom, jpjzoom, jperio 
     231         &             jpizoom, jpjzoom, jperio, jphgr_msh, & 
     232         &             ppglam0, ppgphi0, ppe1_deg, ppe2_deg, ppe1_m, ppe2_m, & 
     233         &             ppsur, ppa0, ppa1, ppkth, ppacr, ppdzmin, pphmax, ldbletanh, & 
     234         &             ppa2, ppkth2, ppacr2 
    232235      !!---------------------------------------------------------------------- 
    233236      ! 
     
    257260904   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcfg in configuration namelist', lwp )    
    258261      WRITE( numond, namcfg ) 
     262! Force values for AGRIF zoom (cf. agrif_user.F90) 
     263#if defined key_agrif 
     264   IF( .NOT. Agrif_Root() ) THEN 
     265      jpiglo  = nbcellsx + 2 + 2*nbghostcells 
     266      jpjglo  = nbcellsy + 2 + 2*nbghostcells 
     267      jpi     = ( jpiglo-2*jpreci + (jpni-1+0) ) / jpni + 2*jpreci 
     268      jpj     = ( jpjglo-2*jprecj + (jpnj-1+0) ) / jpnj + 2*jprecj 
     269      jpidta  = jpiglo 
     270      jpjdta  = jpjglo 
     271      jpizoom = 1 
     272      jpjzoom = 1 
     273      nperio  = 0 
     274      jperio  = 0 
     275   ENDIF 
     276#endif 
    259277      ! 
    260278      !                             !--------------------------------------------! 
     
    309327         jpj = ( jpjglo-2*jprecj + (jpnj-1) ) / jpnj + 2*jprecj   ! second dim. 
    310328#endif 
     329      ENDIF 
    311330         jpk = jpkdta                                             ! third dim 
    312331         jpim1 = jpi-1                                            ! inner domain indices 
     
    314333         jpkm1 = jpk-1                                            !   "           " 
    315334         jpij  = jpi*jpj                                          !  jpi x j 
    316       ENDIF 
    317335 
    318336      IF(lwp) THEN                            ! open listing units 
     
    489507         WRITE(numout,*) '      left bottom j index of the zoom (in data domain) jpizoom = ', jpjzoom 
    490508         WRITE(numout,*) '      lateral cond. type (between 0 and 6) jperio = ', jperio    
    491  
     509         WRITE(numout,*) '      type of horizontal mesh jphgr_msh           = ', jphgr_msh 
     510         WRITE(numout,*) '      longitude of first raw and column T-point ppglam0 = ', ppglam0 
     511         WRITE(numout,*) '      latitude  of first raw and column T-point ppgphi0 = ', ppgphi0 
     512         WRITE(numout,*) '      zonal      grid-spacing (degrees) ppe1_deg        = ', ppe1_deg 
     513         WRITE(numout,*) '      meridional grid-spacing (degrees) ppe2_deg        = ', ppe2_deg 
     514         WRITE(numout,*) '      zonal      grid-spacing (degrees) ppe1_m          = ', ppe1_m 
     515         WRITE(numout,*) '      meridional grid-spacing (degrees) ppe2_m          = ', ppe2_m 
     516         WRITE(numout,*) '      ORCA r4, r2 and r05 coefficients  ppsur           = ', ppsur 
     517         WRITE(numout,*) '                                        ppa0            = ', ppa0 
     518         WRITE(numout,*) '                                        ppa1            = ', ppa1 
     519         WRITE(numout,*) '                                        ppkth           = ', ppkth 
     520         WRITE(numout,*) '                                        ppacr           = ', ppacr 
     521         WRITE(numout,*) '      Minimum vertical spacing ppdzmin                  = ', ppdzmin 
     522         WRITE(numout,*) '      Maximum depth pphmax                              = ', pphmax 
     523         WRITE(numout,*) '      Use double tanf function for vertical coordinates ldbletanh = ', ldbletanh 
     524         WRITE(numout,*) '      Double tanh function parameters ppa2              = ', ppa2 
     525         WRITE(numout,*) '                                      ppkth2            = ', ppkth2 
     526         WRITE(numout,*) '                                      ppacr2            = ', ppacr2 
    492527      ENDIF 
    493528      !                             ! Parameter control 
Note: See TracChangeset for help on using the changeset viewer.