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

Changeset 4667


Ignore:
Timestamp:
2014-06-12T10:47:52+02:00 (10 years ago)
Author:
cbricaud
Message:

correction pf 1D reference config. see ticket 1338

Location:
trunk/NEMOGCM
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/CONFIG/C1D_PAPA/EXP00/iodef.xml

    r4158 r4667  
    2121    --> 
    2222     
    23     <file_definition type="one_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="1d" min_digits="4"> 
     23    <file_definition type="multiple_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="1d" min_digits="4"> 
    2424     
    2525      <file_group id="1ts" output_freq="1ts"  output_level="10" enabled=".TRUE."/> <!-- 1 time step files --> 
  • trunk/NEMOGCM/CONFIG/C1D_PAPA/EXP00/namelist_cfg

    r4370 r4667  
    77   cn_exp      =  "PAPA"   !  experience name 
    88   nn_it000    =       1   !  first time step 
    9    nn_itend    =   78840   !  last  time step 
     9   nn_itend    =   87600   !  last  time step 
    1010   nn_date0    =  20100615 !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) 
    1111   nn_leapy    =       1   !  Leap year calendar (1) or not (0) 
    12    nn_stock    =   78840   !  frequency of creation of a restart file (modulo referenced to 1) 
    13    nn_write    =      60   !  frequency of write in the output file   (modulo referenced to nn_it000) 
     12   nn_stock    =   87600   !  frequency of creation of a restart file (modulo referenced to 1) 
     13   nn_write    =     240   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    1414/ 
    1515!----------------------------------------------------------------------- 
     
    3838&namdom        !   space and time domain (bathymetry, mesh, timestep) 
    3939!----------------------------------------------------------------------- 
     40   nn_bathy    =    0      !  compute (=0) or read (=1) the bathymetry file 
    4041   rn_bathy    = 4198.     !  value of the bathymetry. if (=0) bottom flat at jpkm1 
    41    nn_msh      =    3      !  create (=1) a mesh file or not (=0) 
    42    rn_rdt      =  400.     !  time step for the dynamics (and tracer if nn_acc=0) 
     42   nn_msh      =    0      !  create (=1) a mesh file or not (=0) 
     43   rn_rdt      =  360.     !  time step for the dynamics (and tracer if nn_acc=0) 
    4344   rn_rdtmin   =  360.           !  minimum time step on tracers (used if nn_acc=1) 
    4445   rn_rdtmax   =  360.           !  maximum time step on tracers (used if nn_acc=1) 
     
    7273/ 
    7374!----------------------------------------------------------------------- 
     75&namc1d        !   1D configuration options                             ("key_c1d") 
     76!----------------------------------------------------------------------- 
     77   ln_c1d_locpt= .false.   ! Localization of 1D config in a grid (T) or independant point (F) 
     78/ 
     79!----------------------------------------------------------------------- 
    7480&namtsd    !   data : Temperature  & Salinity 
    7581!----------------------------------------------------------------------- 
     
    8187   cn_dir        = './'     !  root directory for the location of the runoff files 
    8288   ln_tsd_init   = .true.   !  Initialisation of ocean T & S with T &S input data (T) or not (F) 
    83    ln_tsd_tradmp = .true.   !  damping of ocean T & S toward T &S input data (T) or not (F) 
     89   ln_tsd_tradmp = .false.   !  damping of ocean T & S toward T &S input data (T) or not (F) 
    8490/ 
    8591!----------------------------------------------------------------------- 
     
    112118!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    113119!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    114    sn_wndi     = 'forcing_PAPASTATION' ,         1         , 'wndwe'   ,   .true.     , .false. , 'yearly' , '' , '', '' 
    115    sn_wndj     = 'forcing_PAPASTATION' ,         1         , 'wndsn'   ,   .true.     , .false. , 'yearly' , '' , '', '' 
    116    sn_qsr      = 'forcing_PAPASTATION' ,         1         , 'qsr'     ,   .false.    , .false. , 'yearly' , '' , '', '' 
    117    sn_qlw      = 'forcing_PAPASTATION' ,         1         , 'qlw'     ,   .false.    , .false. , 'yearly' , '' , '', '' 
    118    sn_tair     = 'forcing_PAPASTATION' ,         1         , 'tair'    ,   .true.     , .false. , 'yearly' , '' , '', '' 
    119    sn_humi     = 'forcing_PAPASTATION' ,         1         , 'humi'    ,   .true.     , .false. , 'yearly' , '' , '', '' 
    120    sn_prec     = 'forcing_PAPASTATION' ,         1         , 'prec'    ,   .false.    , .false. , 'yearly' , '' , '', '' 
    121    sn_snow     = 'forcing_PAPASTATION' ,         1         , 'snow'    ,   .false.    , .false. , 'yearly' , '' , '', '' 
     120   sn_wndi     = 'forcing_PAPASTATION_1h' ,         1         , 'wndwe'   ,   .true.     , .false. , 'yearly' , '' , '', '' 
     121   sn_wndj     = 'forcing_PAPASTATION_1h' ,         1         , 'wndsn'   ,   .true.     , .false. , 'yearly' , '' , '', '' 
     122   sn_qsr      = 'forcing_PAPASTATION_1h' ,         1         , 'qsr'     ,   .false.    , .false. , 'yearly' , '' , '', '' 
     123   sn_qlw      = 'forcing_PAPASTATION_1h' ,         1         , 'qlw'     ,   .false.    , .false. , 'yearly' , '' , '', '' 
     124   sn_tair     = 'forcing_PAPASTATION_1h' ,         1         , 'tair'    ,   .true.     , .false. , 'yearly' , '' , '', '' 
     125   sn_humi     = 'forcing_PAPASTATION_1h' ,         1         , 'humi'    ,   .true.     , .false. , 'yearly' , '' , '', '' 
     126   sn_prec     = 'forcing_PAPASTATION_1h' ,         1         , 'prec'    ,   .false.    , .false. , 'yearly' , '' , '', '' 
     127   sn_snow     = 'forcing_PAPASTATION_1h' ,         1         , 'snow'    ,   .false.    , .false. , 'yearly' , '' , '', '' 
    122128   ln_2m         = .true.    !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
    123129/ 
     
    225231&namtra_ldf    !   lateral diffusion scheme for tracers 
    226232!---------------------------------------------------------------------------------- 
     233!---------------------------------------------------------------------------------- 
     234   ln_traldf_hor    =  .true.   !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
     235   ln_traldf_iso    =  .false.  !  iso-neutral                 (needs "key_ldfslp") 
    227236   rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    228237   rn_aht_0         =     0.    !  horizontal eddy diffusivity for tracers [m2/s] 
     
    244253/ 
    245254!----------------------------------------------------------------------- 
     255&namc1d_uvd    !   data: U & V currents                                 ("key_c1d") 
     256!----------------------------------------------------------------------- 
     257/ 
     258!----------------------------------------------------------------------- 
    246259&namdyn_hpg    !   Hydrostatic pressure gradient option 
    247260!----------------------------------------------------------------------- 
     
    258271&namzdf        !   vertical physics 
    259272!----------------------------------------------------------------------- 
    260    rn_avm0     =   5.0e-6  !rbb 1.2e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
    261    rn_avt0     =   5.0e-6  !rbb 1.2e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
     273!   rn_avm0     =   5.0e-6  !rbb 1.2e-4  !  vertical eddy viscosity   [m2/s]          (background Kz if not "key_zdfcst") 
     274!   rn_avt0     =   5.0e-6  !rbb 1.2e-5  !  vertical eddy diffusivity [m2/s]          (background Kz if not "key_zdfcst") 
    262275   ln_zdfevd   = .false. !rbb .true.    !  enhanced vertical diffusion (evd) (T) or not (F) 
    263276/ 
  • trunk/NEMOGCM/CONFIG/C1D_PAPA/cpp_C1D_PAPA.fcm

    r4247 r4667  
    1  bld::tool::fppkeys key_c1d key_dynspg_flt key_diahth key_zdfgls key_iomput key_mpp_mpi 
     1 bld::tool::fppkeys key_c1d key_zdfgls 
  • trunk/NEMOGCM/CONFIG/SHARED/namelist_ref

    r4384 r4667  
    183183/ 
    184184!----------------------------------------------------------------------- 
     185&namc1d        !   1D configuration options                             ("key_c1d") 
     186!----------------------------------------------------------------------- 
     187   rn_lat1d    =      50   !  Column latitude (default at PAPA station) 
     188   rn_lon1d    =    -145   !  Column longitude (default at PAPA station) 
     189   ln_c1d_locpt=  .true.   ! Localization of 1D config in a grid (T) or independant point (F) 
     190/ 
     191!----------------------------------------------------------------------- 
    185192&namtsd    !   data : Temperature  & Salinity 
    186193!----------------------------------------------------------------------- 
     
    957964                           !     (no physical validity of the results) 
    958965   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0) 
    959 / 
    960 !----------------------------------------------------------------------- 
    961 &namc1d        !   1D configuration options                             ("key_c1d") 
    962 !----------------------------------------------------------------------- 
    963    rn_lat      =    50     !  Column latitude 
    964    rn_lon      =    -145   !  Column longitude 
    965966/ 
    966967!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/NEMO/OPA_SRC/C1D/c1d.F90

    r4624 r4667  
    2727   REAL(wp), PUBLIC            ::  rn_lat1d     ! Column latitude 
    2828   REAL(wp), PUBLIC            ::  rn_lon1d     ! Column longitude 
     29   LOGICAL , PUBLIC            ::  ln_c1d_locpt ! Localization (or not) of 1D column in a grid 
    2930 
    3031   !!---------------------------------------------------------------------- 
     
    4445      !!---------------------------------------------------------------------- 
    4546      INTEGER ::   ios                 ! Local integer output status for namelist read 
    46       NAMELIST/namc1d/ rn_lat1d, rn_lon1d 
     47      NAMELIST/namc1d/ rn_lat1d, rn_lon1d , ln_c1d_locpt 
    4748      !!---------------------------------------------------------------------- 
    4849      ! 
     
    6364         WRITE(numout,*) '~~~~~~~~' 
    6465         WRITE(numout,*) '   Namelist namc1d : set options for the C1D model' 
    65          WRITE(numout,*) '      column latitude                 rn_lat1d = ', rn_lat1d 
    66          WRITE(numout,*) '      column longitude                rn_lon1d = ', rn_lon1d 
     66         WRITE(numout,*) '      column latitude                 rn_lat1d     = ', rn_lat1d 
     67         WRITE(numout,*) '      column longitude                rn_lon1d     = ', rn_lon1d 
     68         WRITE(numout,*) '      column localization in a grid   ln_c1d_locpt = ', ln_c1d_locpt 
    6769      ENDIF 
    6870      ! 
     
    7880   LOGICAL, PUBLIC, PARAMETER ::   lk_c1d = .FALSE.   !: 1D config. flag de-activated 
    7981   REAL(wp)                   ::   rn_lat1d, rn_lon1d 
     82   LOGICAL , PUBLIC           ::   ln_c1d_locpt = .FALSE.  
     83 
    8084CONTAINS 
    8185 
  • trunk/NEMOGCM/NEMO/OPA_SRC/C1D/domc1d.F90

    r4245 r4667  
    4343      !! ** Action  : Recalculate jpizoom, jpjzoom (indices of C1D zoom) 
    4444      !!---------------------------------------------------------------------- 
     45      NAMELIST/namdom/ nn_bathy, rn_bathy , rn_e3zps_min, rn_e3zps_rat, nn_msh, rn_hmin,   & 
     46         &             nn_acc   , rn_atfp     , rn_rdt      , rn_rdtmin ,                  & 
     47         &             rn_rdtmax, rn_rdth     , nn_closea , ln_crs,    & 
     48         &             jphgr_msh, & 
     49         &             ppglam0, ppgphi0, ppe1_deg, ppe2_deg, ppe1_m, ppe2_m, & 
     50         &             ppsur, ppa0, ppa1, ppkth, ppacr, ppdzmin, pphmax, ldbletanh, & 
     51         &             ppa2, ppkth2, ppacr2 
     52 
    4553      INTEGER  ::  ji, jj                          ! Dummy loop indices 
    4654      INTEGER  ::  inum                            ! Coordinate file handle (case 0) 
    4755      INTEGER  ::  ijeq                            ! Index of equator T point (case 4) 
     56      INTEGER  ::  ios                             ! Local integer output status for namelist read 
    4857 
    4958      INTEGER , DIMENSION(2) ::   iloc             ! Minloc returned indices 
     
    6372      IF( nn_timing == 1 )   CALL timing_start('dom_c1d') 
    6473 
     74      REWIND( numnam_ref )              ! Namelist namdom in reference namelist : space & time domain (bathymetry, mesh, timestep) 
     75      READ  ( numnam_ref, namdom, IOSTAT = ios, ERR = 901 ) 
     76901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdom in reference namelist', lwp ) 
     77   
     78      ! 
     79      REWIND( numnam_cfg )              ! Namelist namdom in configuration namelist : space & time domain (bathymetry, mesh, timestep) 
     80      READ  ( numnam_cfg, namdom, IOSTAT = ios, ERR = 902 ) 
     81902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdom in configuration namelist', lwp ) 
     82 
    6583      CALL wrk_alloc( jpidta, jpjdta, gphidta, glamdta, zdist ) 
    6684 
     
    8098         CALL iom_get( inum, jpdom_unknown, 'gphit', gphidta ) ! so use jpdom_unknown not jpdom_data 
    8199         CALL iom_close ( inum ) 
    82  
    83          PRINT *,'Check dom_c1d coordinates file data read in:' !!! 
    84          PRINT *,'Bottom-left most glamdta is ', glamdta(1,1)    !!! Need to check 
    85          PRINT *,'Bottom-left most gphidta is ', gphidta(1,1)    !!! field read 
    86          PRINT *,'We are using nimpp,njmpp = ' , nimpp,njmpp     !!! 
    87100 
    88101      CASE ( 1 )                 ! geographical mesh on the sphere with regular grid-spacing 
  • trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domcfg.F90

    r4245 r4667  
    8282      !!---------------------------------------------------------------------- 
    8383      !                              ! recalculate jpizoom/jpjzoom given lat/lon 
    84       IF( lk_c1d )  CALL dom_c1d( rn_lat1d, rn_lon1d ) 
     84      IF( lk_c1d .AND. ln_c1d_locpt )  CALL dom_c1d( rn_lat1d, rn_lon1d ) 
    8585      ! 
    8686      !                        ! ============== ! 
  • trunk/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90

    r4645 r4667  
    568568      ENDIF 
    569569      ! 
    570       IF( lk_c1d .AND. .NOT.lk_iomput )   CALL ctl_stop( 'nemo_ctl: The 1D configuration must be used ',   & 
    571          &                                               'with the IOM Input/Output manager. '         ,   & 
    572          &                                               'Compile with key_iomput enabled' ) 
    573       ! 
    574570      IF( 1_wp /= SIGN(1._wp,-0._wp)  )   CALL ctl_stop( 'nemo_ctl: The intrinsec SIGN function follows ',  & 
    575571         &                                               'f2003 standard. '                              ,  & 
Note: See TracChangeset for help on using the changeset viewer.