Changeset 5282


Ignore:
Timestamp:
2015-05-18T17:19:50+02:00 (5 years ago)
Author:
diovino
Message:

Dev. branch CMCC4_simplification ticket #1456

Location:
branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM
Files:
3 added
10 deleted
111 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/ARCH/CMCC/arch-ifort_athena.fcm

    r4865 r5282  
    1515 
    1616 
    17 %NCDF_INC -I$NETCDF/include 
    18 %NCDF_LIB -L$NETCDF/lib -lnetcdf -lnetcdff 
     17%NCDF_INC -I${NETCDF}/include 
     18%NCDF_LIB -L${NETCDF}/lib -lnetcdf -lnetcdff 
    1919%CPP                 cpp 
    2020%FC mpiifort 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/ARCH/CMCC/arch-ifort_athena_xios.fcm

    r4997 r5282  
    3232 
    3333# required modules 
    34 # module load  INTEL/intel_xe_2013 NETCDF/netcdf-4.3_parallel NETCDF/parallel-netcdf-1.3.1 HDF5/hdf5-1.8.11_parallel 
     34# module load  INTEL/intel_xe_2013 NETCDF/netcdf-4.3_parallel_shared NETCDF/parallel-netcdf-1.3.1 HDF5/hdf5-1.8.11_parallel_shared 
    3535 
    3636# Environment variables set by user. Others should automatically define when loading modules. 
    37 # export XIOS=/users/home/models/nemo/xios 
     37#export XIOS=/users/home/models/nemo/xios 
     38#export HDF5=/users/home/opt/hdf5/hdf5-1.8.11_parallel_shared 
     39#export NETCDF=/users/home/opt/netcdf/netcdf-4.3_parallel_shared 
    3840 
    39 %NCDF_INC            -I$NETCDF/include -I$PNETCDF/include 
    40 %NCDF_LIB            -L$NETCDF/lib -lnetcdff -lnetcdf -L$PNETCDF/lib -lpnetcdf 
    41 %HDF5_INC            -I$PHDF5/include 
    42 %HDF5_LIB            -L$PHDF5/lib -lhdf5_hl -lhdf5 
     41%NCDF_INC            -I${NETCDF}/include  
     42%NCDF_LIB            -L${NETCDF}/lib -lnetcdff -lnetcdf 
     43%HDF5_INC            -I$HDF5/include 
     44%HDF5_LIB            -L$HDF5/lib -lhdf5_hl -lhdf5  
    4345%XIOS_INC            -I$XIOS/inc 
    4446%XIOS_LIB            -L$XIOS/lib -lxios 
    4547%CPP                 cpp 
    4648%FC                  mpiifort 
    47 %FCFLAGS             -r8 -O3 -xHost -fp-model source -traceback 
     49%FCFLAGS             -r8 -O3 -xHost -fp-model source -traceback ${CFLAGS}  
    4850%FFLAGS              %FCFLAGS 
    4951%LD                  mpiifort 
    5052%FPPFLAGS            -P -C -traditional 
    51 %LDFLAGS             -lstdc++ -lz -lgpfs -lcurl 
     53%LDFLAGS             -lstdc++ -lz -lgpfs -lcurl  ${LDFLAGS} 
    5254%AR                  ar  
    5355%ARFLAGS             -r 
    5456%MK                  gmake 
    55 %USER_INC            %XIOS_INC %NCDF_INC %HDF5_INC 
     57%USER_INC            %XIOS_INC %NCDF_INC %HDF5_INC  
    5658%USER_LIB            %XIOS_LIB %NCDF_LIB %HDF5_LIB 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/AMM12/EXP00/namelist_cfg

    r4839 r5282  
    5050!----------------------------------------------------------------------- 
    5151   rn_rdt      =   600.    !  time step for the dynamics (and tracer if nn_acc=0) 
    52    rn_rdtmin   =   600.          !  minimum time step on tracers (used if nn_acc=1) 
    53    rn_rdtmax   =   600.          !  maximum time step on tracers (used if nn_acc=1) 
    54    rn_rdth     =   600.          !  depth variation of tracer time step  (used if nn_acc=1) 
    5552   ppglam0     =  999999.0             !  longitude of first raw and column T-point (jphgr_msh = 1) 
    5653   ppgphi0     =  999999.0             ! latitude  of first raw and column T-point (jphgr_msh = 1) 
     
    422419!----------------------------------------------------------------------- 
    423420/ 
    424 !----------------------------------------------------------------------- 
    425 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    426 !----------------------------------------------------------------------- 
    427 / 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/C1D_PAPA/EXP00/namelist_cfg

    r4667 r5282  
    4141   rn_bathy    = 4198.     !  value of the bathymetry. if (=0) bottom flat at jpkm1 
    4242   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) 
    44    rn_rdtmin   =  360.           !  minimum time step on tracers (used if nn_acc=1) 
    45    rn_rdtmax   =  360.           !  maximum time step on tracers (used if nn_acc=1) 
     43   rn_rdt      =  360.     !  time step for the dynamics  
    4644   jphgr_msh   =       1                 !  type of horizontal mesh 
    4745   ppglam0     =    -150.0               !  longitude of first raw and column T-point (jphgr_msh = 1) 
     
    174172/ 
    175173!----------------------------------------------------------------------- 
    176 &namcla        !   cross land advection 
    177 !----------------------------------------------------------------------- 
    178 / 
    179 !----------------------------------------------------------------------- 
    180174&namobc        !   open boundaries parameters                           ("key_obc") 
    181175!----------------------------------------------------------------------- 
     
    356350!----------------------------------------------------------------------- 
    357351/ 
    358 !----------------------------------------------------------------------- 
    359 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    360 !----------------------------------------------------------------------- 
    361    ln_neptramp       = .false.  ! ramp down Neptune velocity in shallow water 
    362 / 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/GYRE/EXP00/namelist_cfg

    r4990 r5282  
    4343!----------------------------------------------------------------------- 
    4444   nn_bathy    =    0      !  compute (=0) or read (=1) the bathymetry file 
    45    rn_rdt      = 7200.     !  time step for the dynamics (and tracer if nn_acc=0) 
    46    rn_rdtmin   = 7200.           !  minimum time step on tracers (used if nn_acc=1) 
    47    rn_rdtmax   = 7200.           !  maximum time step on tracers (used if nn_acc=1) 
     45   rn_rdt      = 7200.     !  time step for the dynamics  
    4846   jphgr_msh   =       5                 !  type of horizontal mesh 
    4947   ppglam0     =       0.0               !  longitude of first raw and column T-point (jphgr_msh = 1) 
     
    156154/ 
    157155!----------------------------------------------------------------------- 
    158 &namcla        !   cross land advection 
    159 !----------------------------------------------------------------------- 
    160 / 
    161 !----------------------------------------------------------------------- 
    162156&namobc        !   open boundaries parameters                           ("key_obc") 
    163157!----------------------------------------------------------------------- 
     
    376370!----------------------------------------------------------------------- 
    377371/ 
    378 !----------------------------------------------------------------------- 
    379 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    380 !----------------------------------------------------------------------- 
    381    ln_neptramp       = .false.  ! ramp down Neptune velocity in shallow water 
    382 / 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_cfg

    r4990 r5282  
    4646!----------------------------------------------------------------------- 
    4747   nn_bathy    =    0      !  compute (=0) or read (=1) the bathymetry file 
    48    rn_rdt      = 7200.     !  time step for the dynamics (and tracer if nn_acc=0) 
    49    rn_rdtmin   = 7200.           !  minimum time step on tracers (used if nn_acc=1) 
    50    rn_rdtmax   = 7200.           !  maximum time step on tracers (used if nn_acc=1) 
     48   rn_rdt      = 7200.     !  time step for the dynamics  
    5149   jphgr_msh   =       5                 !  type of horizontal mesh 
    5250   ppglam0     =       0.0               !  longitude of first raw and column T-point (jphgr_msh = 1) 
     
    161159/ 
    162160!----------------------------------------------------------------------- 
    163 &namcla        !   cross land advection 
    164 !----------------------------------------------------------------------- 
    165 / 
    166 !----------------------------------------------------------------------- 
    167161&namobc        !   open boundaries parameters                           ("key_obc") 
    168162!----------------------------------------------------------------------- 
     
    358352!----------------------------------------------------------------------- 
    359353/ 
    360 !----------------------------------------------------------------------- 
    361 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    362 !----------------------------------------------------------------------- 
    363    ln_neptramp       = .false.  ! ramp down Neptune velocity in shallow water 
    364 / 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_cfg

    r4990 r5282  
    3636!----------------------------------------------------------------------- 
    3737   nn_bathy    =    0      !  compute (=0) or read (=1) the bathymetry file 
    38    rn_rdt      = 7200.     !  time step for the dynamics (and tracer if nn_acc=0) 
    39    rn_rdtmin   = 7200.           !  minimum time step on tracers (used if nn_acc=1) 
    40    rn_rdtmax   = 7200.           !  maximum time step on tracers (used if nn_acc=1) 
     38   rn_rdt      = 7200.     !  time step for the dynamics  
    4139   jphgr_msh   =       5                 !  type of horizontal mesh 
    4240   ppglam0     =       0.0               !  longitude of first raw and column T-point (jphgr_msh = 1) 
     
    104102/ 
    105103!----------------------------------------------------------------------- 
    106 &namcla        !   cross land advection 
    107 !----------------------------------------------------------------------- 
    108 / 
    109 !----------------------------------------------------------------------- 
    110104&nambfr        !   bottom friction 
    111105!----------------------------------------------------------------------- 
     
    213207!----------------------------------------------------------------------- 
    214208/ 
    215 !----------------------------------------------------------------------- 
    216 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    217 !----------------------------------------------------------------------- 
    218    ln_neptramp       = .false.  ! ramp down Neptune velocity in shallow water 
    219 / 
     209 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/GYRE_XIOS/EXP00/namelist_cfg

    r4373 r5282  
    4040!----------------------------------------------------------------------- 
    4141   nn_bathy    =    0      !  compute (=0) or read (=1) the bathymetry file 
    42    rn_rdt      = 7200.     !  time step for the dynamics (and tracer if nn_acc=0) 
     42   rn_rdt      = 7200.     !  time step for the dynamics  
    4343!   nn_baro     =   60      !  number of barotropic time step            ("key_dynspg_ts") 
    44    rn_rdtmin   = 7200.           !  minimum time step on tracers (used if nn_acc=1) 
    45    rn_rdtmax   = 7200.           !  maximum time step on tracers (used if nn_acc=1) 
    4644   jphgr_msh   =       5                 !  type of horizontal mesh 
    4745   ppglam0     =       0.0               !  longitude of first raw and column T-point (jphgr_msh = 1) 
     
    150148/ 
    151149!----------------------------------------------------------------------- 
    152 &namcla        !   cross land advection 
    153 !----------------------------------------------------------------------- 
    154 / 
    155 !----------------------------------------------------------------------- 
    156150&namobc        !   open boundaries parameters                           ("key_obc") 
    157151!----------------------------------------------------------------------- 
     
    195189&nameos        !   ocean physical parameters 
    196190!----------------------------------------------------------------------- 
    197    nn_eos      =   2       !  type of equation of state and Brunt-Vaisala frequency 
     191   nn_eos      =   0       !  type of equation of state and Brunt-Vaisala frequency 
    198192/ 
    199193!----------------------------------------------------------------------- 
     
    330324!----------------------------------------------------------------------- 
    331325/ 
    332 !----------------------------------------------------------------------- 
    333 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    334 !----------------------------------------------------------------------- 
    335    ln_neptramp       = .false.  ! ramp down Neptune velocity in shallow water 
    336 / 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/ISOMIP/EXP00/namelist_cfg

    r4924 r5282  
    55!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
    66!!                                    namsbc_apr, namsbc_ssr, namsbc_alb) 
    7 !!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 
     7!!              4 - lateral boundary (namlbc, namobc, namagrif, nambdy, nambdy_tide) 
    88!!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
    99!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_dmp) 
     
    4141   nn_stock    =  999999   !  frequency of creation of a restart file (modulo referenced to 1) 
    4242   nn_write    =  999999   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    43    ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    4443   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    4544   ln_clobber  = .false.   !  clobber (overwrite) an existing file 
     
    9998   rn_e3zps_rat=    0.1    !  rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 
    10099                           ! 
    101    rn_rdt      = 1800.     !  time step for the dynamics (and tracer if nn_acc=0) 
     100   rn_rdt      = 1800.     !  time step for the dynamics  
    102101   rn_atfp     =    0.1    !  asselin time filter parameter 
    103    nn_acc      =    0      !  acceleration of convergence : =1      used, rdt < rdttra(k) 
    104                                  !                          =0, not used, rdt = rdttra 
    105    rn_rdtmin   = 1800.           !  minimum time step on tracers (used if nn_acc=1) 
    106    rn_rdtmax   = 1800.           !  maximum time step on tracers (used if nn_acc=1) 
    107    rn_rdth     =  800.           !  depth variation of tracer time step  (used if nn_acc=1) 
    108102   ln_crs      = .false.      !  Logical switch for coarsening module 
    109103   jphgr_msh   =       1               !  type of horizontal mesh 
     
    290284!!====================================================================== 
    291285!!   namlbc        lateral momentum boundary condition 
    292 !!   namcla        cross land advection 
    293286!!   namobc        open boundaries parameters                           ("key_obc") 
    294287!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif") 
     
    302295   rn_shlat    =    0.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    303296                           !  free slip  !   partial slip  !   no slip   ! strong slip 
    304 / 
    305 !----------------------------------------------------------------------- 
    306 &namcla        !   cross land advection 
    307 !----------------------------------------------------------------------- 
    308297/ 
    309298!----------------------------------------------------------------------- 
     
    420409   !                                           (normally=0; not used with Griffies) 
    421410   rn_slpmax        =     0.01  !  slope limit 
    422    rn_chsmag        =     1.    !  multiplicative factor in Smagorinsky diffusivity 
    423    rn_smsh          =     1.    !  Smagorinsky diffusivity: = 0 - use only sheer 
    424411   rn_aht_m         =   100.    !  upper limit or stability criteria for lateral eddy diffusivity (m2/s) 
    425412/ 
     
    482469   rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s] 
    483470   rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s] 
    484    rn_cmsmag_1      =     3.    !  constant in laplacian Smagorinsky viscosity 
    485    rn_cmsmag_2      =     3     !  constant in bilaplacian Smagorinsky viscosity 
    486    rn_cmsh          =     1.    !  1 or 0 , if 0 -use only shear for Smagorinsky viscosity 
    487471   rn_ahm_m_blp     =    -1.e12 !  upper limit for bilap  abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 
    488472   rn_ahm_m_lap     =   600.    !  upper limit for lap  ahm < min(dx^2/16rdt, rn_ahm_m_lap) 
     
    659643!----------------------------------------------------------------------- 
    660644/ 
    661 !----------------------------------------------------------------------- 
    662 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    663 !----------------------------------------------------------------------- 
    664 / 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/1_namelist_cfg

    r4990 r5282  
    5050   ppkth2      =  999999.              ! 
    5151   ppacr2      =  999999.              ! 
    52    rn_rdt      = 2880.     !  time step for the dynamics (and tracer if nn_acc=0) 
    53    rn_rdtmin   = 14400.          !  minimum time step on tracers (used if nn_acc=1) 
    54    rn_rdtmax   = 14400.          !  maximum time step on tracers (used if nn_acc=1) 
     52   rn_rdt      = 2880.     !  time step for the dynamics  
    5553/ 
    5654!----------------------------------------------------------------------- 
     
    117115/ 
    118116!----------------------------------------------------------------------- 
    119 &namcla        !   cross land advection 
    120 !----------------------------------------------------------------------- 
    121 / 
    122 !----------------------------------------------------------------------- 
    123117&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    124118!----------------------------------------------------------------------- 
     
    218212/ 
    219213!----------------------------------------------------------------------- 
    220 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    221 !----------------------------------------------------------------------- 
    222 / 
    223 !----------------------------------------------------------------------- 
    224214&namobs       !  observation usage                                      ('key_diaobs') 
    225215!----------------------------------------------------------------------- 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist_cfg

    r4990 r5282  
    9393/ 
    9494!----------------------------------------------------------------------- 
    95 &namcla        !   cross land advection 
    96 !----------------------------------------------------------------------- 
    97 / 
    98 !----------------------------------------------------------------------- 
    9995&nambfr        !   bottom friction 
    10096!----------------------------------------------------------------------- 
     
    185181/ 
    186182!----------------------------------------------------------------------- 
    187 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    188 !----------------------------------------------------------------------- 
    189 / 
    190 !----------------------------------------------------------------------- 
    191183&namobs       !  observation usage                                      ('key_diaobs') 
    192184!----------------------------------------------------------------------- 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/1_namelist_cfg

    r4990 r5282  
    5050   ppkth2      =  999999.              ! 
    5151   ppacr2      =  999999.              ! 
    52    rn_rdt      = 2880.     !  time step for the dynamics (and tracer if nn_acc=0) 
    53    rn_rdtmin   = 14400.          !  minimum time step on tracers (used if nn_acc=1) 
    54    rn_rdtmax   = 14400.          !  maximum time step on tracers (used if nn_acc=1) 
     52   rn_rdt      = 2880.     !  time step for the dynamics 
    5553/ 
    5654!----------------------------------------------------------------------- 
     
    117115/ 
    118116!----------------------------------------------------------------------- 
    119 &namcla        !   cross land advection 
    120 !----------------------------------------------------------------------- 
    121 / 
    122 !----------------------------------------------------------------------- 
    123117&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    124118!----------------------------------------------------------------------- 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_cfg

    r4995 r5282  
    9393/ 
    9494!----------------------------------------------------------------------- 
    95 &namcla        !   cross land advection 
    96 !----------------------------------------------------------------------- 
    97 / 
    98 !----------------------------------------------------------------------- 
    9995&nambfr        !   bottom friction 
    10096!----------------------------------------------------------------------- 
     
    187183!----------------------------------------------------------------------- 
    188184/ 
    189 !----------------------------------------------------------------------- 
    190 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    191 !----------------------------------------------------------------------- 
    192 / 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/1_namelist_cfg

    r4147 r5282  
    55!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf,  
    66!!                                    namsbc_apr, namsbc_ssr, namsbc_alb) 
    7 !!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 
     7!!              4 - lateral boundary (namlbc, namobc, namagrif, nambdy, nambdy_tide) 
    88!!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
    99!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_dmp) 
     
    4040   nn_stock    =   10950   !  frequency of creation of a restart file (modulo referenced to 1) 
    4141   nn_write    =   10950   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    42    ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    4342   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    4443   ln_clobber  = .false.   !  clobber (overwrite) an existing file 
     
    8584   rn_e3zps_rat=    0.1    !  rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 
    8685                           ! 
    87    rn_rdt      = 2880.     !  time step for the dynamics (and tracer if nn_acc=0) 
     86   rn_rdt      = 2880.     !  time step for the dynamics  
    8887   nn_baro     =   64      !  number of barotropic time step            ("key_dynspg_ts") 
    8988   rn_atfp     =    0.1    !  asselin time filter parameter 
    90    nn_acc      =    0      !  acceleration of convergence : =1      used, rdt < rdttra(k) 
    91                                  !                          =0, not used, rdt = rdttra 
    92    rn_rdtmin   = 14400.          !  minimum time step on tracers (used if nn_acc=1) 
    93    rn_rdtmax   = 14400.          !  maximum time step on tracers (used if nn_acc=1) 
    94    rn_rdth     =  800.           !  depth variation of tracer time step  (used if nn_acc=1) 
    9589/ 
    9690!----------------------------------------------------------------------- 
     
    302296!!====================================================================== 
    303297!!   namlbc        lateral momentum boundary condition 
    304 !!   namcla        cross land advection 
    305298!!   namobc        open boundaries parameters                           ("key_obc") 
    306299!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")  
     
    314307   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    315308                           !  free slip  !   partial slip  !   no slip   ! strong slip 
    316 / 
    317 !----------------------------------------------------------------------- 
    318 &namcla        !   cross land advection 
    319 !----------------------------------------------------------------------- 
    320    nn_cla      =    0      !  advection between 2 ocean pts separates by land 
    321309/ 
    322310!----------------------------------------------------------------------- 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/namelist_cfg

    r4370 r5282  
    132132/ 
    133133!----------------------------------------------------------------------- 
    134 &namcla        !   cross land advection 
    135 !----------------------------------------------------------------------- 
    136 / 
    137 !----------------------------------------------------------------------- 
    138134&namobc        !   open boundaries parameters                           ("key_obc") 
    139135!----------------------------------------------------------------------- 
     
    291287!----------------------------------------------------------------------- 
    292288/ 
    293 !----------------------------------------------------------------------- 
    294 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    295 !----------------------------------------------------------------------- 
    296 / 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_cfg

    r4370 r5282  
    9393/ 
    9494!----------------------------------------------------------------------- 
    95 &namcla        !   cross land advection 
    96 !----------------------------------------------------------------------- 
    97 / 
    98 !----------------------------------------------------------------------- 
    9995&nambfr        !   bottom friction 
    10096!----------------------------------------------------------------------- 
     
    184180!----------------------------------------------------------------------- 
    185181/ 
    186 !----------------------------------------------------------------------- 
    187 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    188 !----------------------------------------------------------------------- 
    189 / 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_cfg

    r4370 r5282  
    3232!----------------------------------------------------------------------- 
    3333   nn_msh      =    1      !  create (=1) a mesh file or not (=0) 
    34    rn_rdt      = 21600.    !  time step for the dynamics (and tracer if nn_acc=0) 
    35    rn_rdtmin   = 21600.          !  minimum time step on tracers (used if nn_acc=1) 
    36    rn_rdtmax   = 21600.          !  maximum time step on tracers (used if nn_acc=1) 
     34   rn_rdt      = 21600.    !  time step for the dynamics  
    3735   jphgr_msh   =       0               !  type of horizontal mesh 
    3836   ppglam0     =  999999.0             !  longitude of first raw and column T-point (jphgr_msh = 1) 
     
    7775!----------------------------------------------------------------------- 
    7876/       
    79  
    80 !----------------------------------------------------------------------- 
    81 &namcla        !   cross land advection 
    82 !----------------------------------------------------------------------- 
    83 / 
    8477!----------------------------------------------------------------------- 
    8578&nambbl        !   bottom boundary layer scheme 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/SHARED/namelist_ref

    r5029 r5282  
    55!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
    66!!                                    namsbc_apr, namsbc_ssr, namsbc_alb) 
    7 !!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 
     7!!              4 - lateral boundary (namlbc, namobc, namagrif, nambdy, nambdy_tide) 
    88!!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
    99!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_dmp) 
     
    4141   nn_stock    =    5475   !  frequency of creation of a restart file (modulo referenced to 1) 
    4242   nn_write    =    5475   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    43    ln_dimgnnn  = .false.   !  DIMG file format: 1 file for all processors (F) or by processor (T) 
    4443   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%) 
    4544   ln_clobber  = .false.   !  clobber (overwrite) an existing file 
     
    121120   rn_e3zps_rat=    0.1    !  rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 
    122121                           ! 
    123    rn_rdt      = 5760.     !  time step for the dynamics (and tracer if nn_acc=0) 
     122   rn_rdt      = 5760.     !  time step for the dynamics  
    124123   rn_atfp     =    0.1    !  asselin time filter parameter 
    125    nn_acc      =    0      !  acceleration of convergence : =1      used, rdt < rdttra(k) 
    126                                  !                          =0, not used, rdt = rdttra 
    127    rn_rdtmin   = 28800.          !  minimum time step on tracers (used if nn_acc=1) 
    128    rn_rdtmax   = 28800.          !  maximum time step on tracers (used if nn_acc=1) 
    129    rn_rdth     =  800.           !  depth variation of tracer time step  (used if nn_acc=1) 
    130124   ln_crs      = .false.      !  Logical switch for coarsening module 
    131125   jphgr_msh   =       0               !  type of horizontal mesh 
     
    526520!!====================================================================== 
    527521!!   namlbc        lateral momentum boundary condition 
    528 !!   namcla        cross land advection 
    529522!!   namobc        open boundaries parameters                           ("key_obc") 
    530523!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif") 
     
    539532                           !  free slip  !   partial slip  !   no slip   ! strong slip 
    540533   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical eqs. 
    541 / 
    542 !----------------------------------------------------------------------- 
    543 &namcla        !   cross land advection 
    544 !----------------------------------------------------------------------- 
    545    nn_cla      =    0      !  advection between 2 ocean pts separates by land 
    546534/ 
    547535!----------------------------------------------------------------------- 
     
    783771   !                                           (normally=0; not used with Griffies) 
    784772   rn_slpmax        =     0.01  !  slope limit 
    785    rn_chsmag        =     1.    !  multiplicative factor in Smagorinsky diffusivity 
    786    rn_smsh          =     1.    !  Smagorinsky diffusivity: = 0 - use only sheer 
    787773   rn_aht_m         =  2000.    !  upper limit or stability criteria for lateral eddy diffusivity (m2/s) 
    788774/ 
     
    876862   rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s] 
    877863   rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s] 
    878    rn_cmsmag_1      =     3.    !  constant in laplacian Smagorinsky viscosity 
    879    rn_cmsmag_2      =     3     !  constant in bilaplacian Smagorinsky viscosity 
    880    rn_cmsh          =     1.    !  1 or 0 , if 0 -use only shear for Smagorinsky viscosity 
    881864   rn_ahm_m_blp     =    -1.e12 !  upper limit for bilap  abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 
    882865   rn_ahm_m_lap     = 40000.    !  upper limit for lap  ahm < min(dx^2/16rdt, rn_ahm_m_lap) 
     
    12591242   cn_dir_cdg  = './'  !  root directory for the location of drag coefficient files 
    12601243/ 
    1261 !----------------------------------------------------------------------- 
    1262 &namdyn_nept  !   Neptune effect (simplified: lateral and vertical diffusions removed) 
    1263 !----------------------------------------------------------------------- 
    1264    ! Suggested lengthscale values are those of Eby & Holloway (1994) for a coarse model 
    1265    ln_neptsimp       = .false.  ! yes/no use simplified neptune 
    1266  
    1267    ln_smooth_neptvel = .false.  ! yes/no smooth zunep, zvnep 
    1268    rn_tslse          =  1.2e4   ! value of lengthscale L at the equator 
    1269    rn_tslsp          =  3.0e3   ! value of lengthscale L at the pole 
    1270    ! Specify whether to ramp down the Neptune velocity in shallow 
    1271    ! water, and if so the depth range controlling such ramping down 
    1272    ln_neptramp       = .true.   ! ramp down Neptune velocity in shallow water 
    1273    rn_htrmin         =  100.0   ! min. depth of transition range 
    1274    rn_htrmax         =  200.0   ! max. depth of transition range 
    1275 / 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/CONFIG/cfg.txt

    r4990 r5282  
    1 GYRE_PISCES OPA_SRC TOP_SRC 
    21ORCA2_LIM_CFC_C14b OPA_SRC LIM_SRC_2 NST_SRC TOP_SRC 
    32GYRE_XIOS OPA_SRC 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/LIM_SRC_2/iceini_2.F90

    r4624 r5282  
    8888      CALL ice_run_2                   ! read in namelist some run parameters 
    8989      !           
    90       rdt_ice = nn_fsbc * rdttra(1)    ! sea-ice time step 
     90      rdt_ice = nn_fsbc * rdt           ! sea-ice time step 
    9191      numit   = nit000 - 1 
    9292      ! 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/LIM_SRC_2/limrst_2.F90

    r2528 r5282  
    6767            WRITE(numout,*) 
    6868            SELECT CASE ( jprstlib ) 
    69             CASE ( jprstdimg )   ;   WRITE(numout,*) '             open ice restart binary file: '//clname 
    7069            CASE DEFAULT         ;   WRITE(numout,*) '             open ice restart NetCDF file: '//clname 
    7170            END SELECT 
     
    185184      ENDIF 
    186185 
    187       IF ( jprstlib == jprstdimg ) THEN 
    188         ! eventually read netcdf file (monobloc)  for restarting on different number of processors 
    189         ! if {cn_icerst_in}.nc exists, then set jlibalt to jpnf90 
    190         INQUIRE( FILE = TRIM(cn_icerst_in)//'.nc', EXIST = llok ) 
    191         IF ( llok ) THEN ; jlibalt = jpnf90  ; ELSE ; jlibalt = jprstlib ; ENDIF 
    192       ENDIF 
    193  
    194186      CALL iom_open ( cn_icerst_in, numrir, kiolib = jlibalt ) 
    195187 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/LIM_SRC_2/limwri_2.F90

    r4696 r5282  
    8585 
    8686#if ! defined key_iomput 
    87 # if defined key_dimgout 
    88    !!---------------------------------------------------------------------- 
    89    !!   'key_dimgout'                                    Direct Access file 
    90    !!---------------------------------------------------------------------- 
    91 # include "limwri_dimg_2.h90" 
    92 # else 
    9387   SUBROUTINE lim_wri_2( kt ) 
    9488      !!------------------------------------------------------------------- 
     
    215209      ! 
    216210   END SUBROUTINE lim_wri_2 
    217       
    218 #endif      
    219211 
    220212   SUBROUTINE lim_wri_init_2 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/LIM_SRC_3/iceini.F90

    r4873 r5282  
    8282      CALL lim_thd_sal_init            ! set ice salinity parameters 
    8383      ! 
    84       rdt_ice   = nn_fsbc * rdttra(1)  ! sea-ice timestep 
     84      rdt_ice   = nn_fsbc * rdt        ! sea-ice timestep 
    8585      r1_rdtice = 1._wp / rdt_ice      ! sea-ice timestep inverse 
    8686      ! 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/LIM_SRC_3/limrst.F90

    r4990 r5282  
    7373            WRITE(numout,*) 
    7474            SELECT CASE ( jprstlib ) 
    75             CASE ( jprstdimg )   ;   WRITE(numout,*) '             open ice restart binary file: '//clname 
    7675            CASE DEFAULT         ;   WRITE(numout,*) '             open ice restart NetCDF file: '//clname 
    7776            END SELECT 
     
    326325      ENDIF 
    327326 
    328       IF ( jprstlib == jprstdimg ) THEN 
    329         ! eventually read netcdf file (monobloc)  for restarting on different number of processors 
    330         ! if {cn_icerst_in}.nc exists, then set jlibalt to jpnf90 
    331         INQUIRE( FILE = TRIM(cn_icerst_in)//'.nc', EXIST = llok ) 
    332         IF ( llok ) THEN ; jlibalt = jpnf90  ; ELSE ; jlibalt = jprstlib ; ENDIF 
    333       ENDIF 
    334  
    335327      CALL iom_open ( cn_icerst_in, numrir, kiolib = jprstlib ) 
    336328 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90

    r4990 r5282  
    4242CONTAINS 
    4343 
    44 #if defined key_dimgout 
    45 # include "limwri_dimg.h90" 
    46 #else 
    4744 
    4845   SUBROUTINE lim_wri( kindic ) 
     
    280277       
    281278   END SUBROUTINE lim_wri 
    282 #endif 
    283279 
    284280  
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OFF_SRC/domain.F90

    r4990 r5282  
    6969      ENDIF 
    7070 
    71       CALL dom_nam      ! read namelist ( namrun, namdom, namcla ) 
     71      CALL dom_nam      ! read namelist ( namrun, namdom ) 
    7272      CALL dom_zgr      ! Vertical mesh and bathymetry option 
    7373      CALL dom_rea      ! Create a domain file 
     
    112112      !! ** input   : - namrun namelist 
    113113      !!              - namdom namelist 
    114       !!              - namcla namelist 
    115114      !!---------------------------------------------------------------------- 
    116115      USE ioipsl 
     
    118117      NAMELIST/namrun/ nn_no   , cn_exp    , cn_ocerst_in, cn_ocerst_out, ln_rstart , nn_rstctl,   & 
    119118         &             nn_it000, nn_itend  , nn_date0    , nn_leapy     , nn_istate , nn_stock ,   & 
    120          &             nn_write, ln_dimgnnn, ln_mskland  , ln_clobber   , nn_chunksz, nn_euler 
     119         &             nn_write, ln_mskland  , ln_clobber   , nn_chunksz, nn_euler 
    121120      NAMELIST/namdom/ nn_bathy , rn_bathy, rn_e3zps_min, rn_e3zps_rat, nn_msh    , rn_hmin,   & 
    122          &             nn_acc   , rn_atfp     , rn_rdt      , rn_rdtmin ,            & 
    123          &             rn_rdtmax, rn_rdth     , nn_baro     , nn_closea , ln_crs, & 
    124          &             jphgr_msh, & 
     121         &             rn_atfp  , rn_rdt   , nn_baro     , nn_closea , ln_crs, jphgr_msh,  & 
    125122         &             ppglam0, ppgphi0, ppe1_deg, ppe2_deg, ppe1_m, ppe2_m, & 
    126123         &             ppsur, ppa0, ppa1, ppkth, ppacr, ppdzmin, pphmax, ldbletanh, & 
    127124         &             ppa2, ppkth2, ppacr2 
    128       NAMELIST/namcla/ nn_cla 
    129125#if defined key_netcdf4 
    130126      NAMELIST/namnc4/ nn_nchunks_i, nn_nchunks_j, nn_nchunks_k, ln_nc4zip 
     
    157153         WRITE(numout,*) '      frequency of restart file       nn_stock   = ', nn_stock 
    158154         WRITE(numout,*) '      frequency of output file        nn_write   = ', nn_write 
    159          WRITE(numout,*) '      multi file dimgout              ln_dimgnnn = ', ln_dimgnnn 
    160155         WRITE(numout,*) '      mask land points                ln_mskland = ', ln_mskland 
    161156         WRITE(numout,*) '      overwrite an existing file      ln_clobber = ', ln_clobber 
     
    188183      ! parameters correspondting to nit000 - 1 (as we start the step loop with a call to day) 
    189184      ndastp = ndate0 - 1        ! ndate0 read in the namelist in dom_nam, we assume that we start run at 00:00 
    190       adatrj = ( REAL( nit000-1, wp ) * rdttra(1) ) / rday 
     185      adatrj = ( REAL( nit000-1, wp ) * rdt ) / rday 
    191186 
    192187#if defined key_agrif 
     
    233228         WRITE(numout,*) '      asselin time filter parameter        rn_atfp   = ', rn_atfp 
    234229         WRITE(numout,*) '      time-splitting: nb of sub time-step  nn_baro   = ', nn_baro 
    235          WRITE(numout,*) '      acceleration of converge             nn_acc    = ', nn_acc 
    236          WRITE(numout,*) '        nn_acc=1: surface tracer rdt       rn_rdtmin = ', rn_rdtmin 
    237          WRITE(numout,*) '                  bottom  tracer rdt       rdtmax    = ', rn_rdtmax 
    238          WRITE(numout,*) '                  depth of transition      rn_rdth   = ', rn_rdth 
    239230         WRITE(numout,*) '      suppression of closed seas (=0)      nn_closea = ', nn_closea 
    240231         WRITE(numout,*) '      type of horizontal mesh jphgr_msh           = ', jphgr_msh 
     
    262253      e3zps_rat = rn_e3zps_rat 
    263254      nmsh      = nn_msh 
    264       nacc      = nn_acc 
    265255      atfp      = rn_atfp 
    266256      rdt       = rn_rdt 
    267       rdtmin    = rn_rdtmin 
    268       rdtmax    = rn_rdtmin 
    269       rdth      = rn_rdth 
    270  
    271       REWIND( numnam_ref )              ! Namelist namcla in reference namelist : Cross land advection 
    272       READ  ( numnam_ref, namcla, IOSTAT = ios, ERR = 905) 
    273 905   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in reference namelist', lwp ) 
    274  
    275       REWIND( numnam_cfg )              ! Namelist namcla in configuration namelist : Cross land advection 
    276       READ  ( numnam_cfg, namcla, IOSTAT = ios, ERR = 906 ) 
    277 906   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in configuration namelist', lwp ) 
    278       IF(lwm) WRITE( numond, namcla ) 
    279  
    280       IF(lwp) THEN 
    281          WRITE(numout,*) 
    282          WRITE(numout,*) '   Namelist namcla' 
    283          WRITE(numout,*) '      cross land advection                 nn_cla    = ', nn_cla 
    284       ENDIF 
    285  
    286257#if defined key_netcdf4 
    287258      !                             ! NetCDF 4 case   ("key_netcdf4" defined) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

    r4998 r5282  
    208208      ENDIF 
    209209 
    210       IF ( nacc /= 0 ) & 
    211          & CALL ctl_stop( ' nacc /= 0 and key_asminc :',  & 
    212          &                ' Assimilation increments have only been implemented', & 
    213          &                ' for synchronous time stepping' ) 
    214210 
    215211      IF ( ( ln_asmdin ).AND.( ln_asmiau ) )   & 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/BDY/bdytides.F90

    r4758 r5282  
    316316      ENDIF 
    317317 
    318       IF ( nsec_day == NINT(0.5_wp * rdttra(1)) .AND. zflag==1 ) THEN 
     318      IF ( nsec_day == NINT(0.5_wp * rdt) .AND. zflag==1 ) THEN 
    319319        ! 
    320320        kt_tide = kt 
     
    431431            ! We refresh nodal factors every day below 
    432432            ! This should be done somewhere else 
    433             IF ( nsec_day == NINT(0.5_wp * rdttra(1)) .AND. lk_first_btstp ) THEN 
    434                ! 
     433            IF ( nsec_day == NINT(0.5_wp * rdt) .AND. lk_first_btstp ) THEN 
     434                ! 
    435435               kt_tide = kt                
    436436               ! 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/C1D/domc1d.F90

    r4667 r5282  
    4444      !!---------------------------------------------------------------------- 
    4545      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, & 
     46         &             rn_atfp , rn_rdt,nn_closea , ln_crs,  jphgr_msh,        & 
    4947         &             ppglam0, ppgphi0, ppe1_deg, ppe2_deg, ppe1_m, ppe2_m, & 
    5048         &             ppsur, ppa0, ppa1, ppkth, ppacr, ppdzmin, pphmax, ldbletanh, & 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/CRS/crsfld.F90

    r4149 r5282  
    5252      !!      2. At time of output, rescale [1] by dimension and time 
    5353      !!         to yield the spatial and temporal average. 
    54       !!  See. diawri_dimg.h90, sbcmod.F90 
     54      !!  See. sbcmod.F90 
    5555      !! 
    5656      !! ** Method  :   
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/diaar5.F90

    r4990 r5282  
    2121   USE timing         ! preformance summary 
    2222   USE wrk_nemo       ! working arrays 
     23   USE fldread        ! type FLD_N 
     24   USE phycst         ! physical constant 
     25   USE in_out_manager  ! I/O manager 
    2326 
    2427   IMPLICIT NONE 
     
    197200      REAL(wp) ::   zztmp   
    198201      REAL(wp), POINTER, DIMENSION(:,:,:,:) ::   zsaldta   ! Jan/Dec levitus salinity 
     202      ! reading initial file 
     203      LOGICAL  ::   ln_tsd_init      !: T & S data flag 
     204      LOGICAL  ::   ln_tsd_tradmp    !: internal damping toward input data flag 
     205      CHARACTER(len=100)            ::   cn_dir 
     206      TYPE(FLD_N)                   ::  sn_tem,sn_sal 
     207      INTEGER  ::   ios=0 
     208 
     209      NAMELIST/namtsd/ ln_tsd_init,ln_tsd_tradmp,cn_dir,sn_tem,sn_sal 
     210      ! 
     211 
     212      REWIND( numnam_ref )              ! Namelist namtsd in reference namelist : 
     213      READ  ( numnam_ref, namtsd, IOSTAT = ios, ERR = 901) 
     214901   IF( ios /= 0 ) CALL ctl_nam ( ios , ' namtsd in reference namelist for dia_ar5', lwp ) 
     215      REWIND( numnam_cfg )              ! Namelist namtsd in configuration namelist : Parameters of the run 
     216      READ  ( numnam_cfg, namtsd, IOSTAT = ios, ERR = 902 ) 
     217902   IF( ios /= 0 ) CALL ctl_nam ( ios , ' namtsd in configuration namelist for dia_ar5', lwp ) 
     218      IF(lwm) WRITE ( numond, namtsd ) 
     219      ! 
    199220      !!---------------------------------------------------------------------- 
    200221      ! 
     
    216237      END DO 
    217238      IF( lk_mpp )   CALL mpp_sum( vol0 ) 
    218        
    219       CALL iom_open ( 'data_1m_salinity_nomask', inum ) 
    220       CALL iom_get  ( inum, jpdom_data, 'vosaline', zsaldta(:,:,:,1), 1  ) 
    221       CALL iom_get  ( inum, jpdom_data, 'vosaline', zsaldta(:,:,:,2), 12 ) 
     239 
     240      CALL iom_open ( TRIM( cn_dir )//TRIM(sn_sal%clname), inum ) 
     241      CALL iom_get  ( inum, jpdom_data, TRIM(sn_sal%clvar), zsaldta(:,:,:,1), 1  ) 
     242      CALL iom_get  ( inum, jpdom_data, TRIM(sn_sal%clvar), zsaldta(:,:,:,2), 12 ) 
    222243      CALL iom_close( inum ) 
    223244      sn0(:,:,:) = 0.5_wp * ( zsaldta(:,:,:,1) + zsaldta(:,:,:,2) )         
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/diafwb.F90

    r4990 r5282  
    115115          
    116116         ! Conversion in m3 
    117          a_fwf    = a_fwf * rdttra(1) * 1.e-3  
     117         a_fwf    = a_fwf * rdt * 1.e-3 
    118118          
    119119         ! fwf correction to bring back the mean ssh to zero 
     
    382382         WRITE(inum,*) 
    383383         WRITE(inum,*)    'Net freshwater budget ' 
    384          WRITE(inum,9010) '  fwf    = ',a_fwf,   ' m3 =', a_fwf   /(FLOAT(nitend-nit000+1)*rdttra(1)) * 1.e-6,' Sv' 
     384         WRITE(inum,9010) '  fwf    = ',a_fwf,   ' m3 =', a_fwf   /(FLOAT(nitend-nit000+1)*rdt) * 1.e-6,' Sv' 
    385385         WRITE(inum,*) 
    386386         WRITE(inum,9010) '  zarea =',zarea 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/diaharm.F90

    r4990 r5282  
    2323   USE ioipsl          ! NetCDF IPSL library 
    2424   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    25    USE diadimg         ! To write dimg 
    2625   USE timing          ! preformance summary 
    2726   USE wrk_nemo        ! working arrays 
     
    388387      !!---------------------------------------------------------------------- 
    389388 
    390 #if defined key_dimgout 
    391       cdfile_name_T = TRIM(cexper)//'_Tidal_harmonics_gridT.dimgproc' 
    392       cdfile_name_U = TRIM(cexper)//'_Tidal_harmonics_gridU.dimgproc' 
    393       cdfile_name_V = TRIM(cexper)//'_Tidal_harmonics_gridV.dimgproc' 
    394 #endif 
     389!#if  
     390!#endif 
    395391 
    396392      IF(lwp) WRITE(numout,*) '  ' 
    397393      IF(lwp) WRITE(numout,*) 'dia_wri_harm : Write harmonic analysis results' 
    398 #if defined key_dimgout 
    399       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~  Output files: ', TRIM(cdfile_name_T) 
    400       IF(lwp) WRITE(numout,*) '                             ', TRIM(cdfile_name_U) 
    401       IF(lwp) WRITE(numout,*) '                             ', TRIM(cdfile_name_V) 
    402 #endif 
     394!#if  
     395!#endif 
    403396      IF(lwp) WRITE(numout,*) '  ' 
    404397 
     
    406399      !///////////// 
    407400      ! 
    408 #if defined key_dimgout 
    409       cltext='Elevation amplitude and phase' 
    410       CALL dia_wri_dimg(TRIM(cdfile_name_T), TRIM(cltext), out_eta, 2*nb_ana, '2') 
    411 #else 
     401!#if  
     402!#else 
    412403      DO jh = 1, nb_ana 
    413404      CALL iom_put( TRIM(tname(jh))//'x', out_eta(:,:,jh) ) 
    414405      CALL iom_put( TRIM(tname(jh))//'y', out_eta(:,:,nb_ana+jh) ) 
    415406      END DO 
    416 #endif 
     407!#endif 
    417408 
    418409      ! B) ubar 
    419410      !///////// 
    420411      ! 
    421 #if defined key_dimgout 
    422       cltext='ubar amplitude and phase' 
    423       CALL dia_wri_dimg(TRIM(cdfile_name_U), TRIM(cltext), out_u, 2*nb_ana, '2') 
    424 #else 
     412!#if  
     413!#else 
    425414      DO jh = 1, nb_ana 
    426415      CALL iom_put( TRIM(tname(jh))//'x_u', out_u(:,:,jh) ) 
    427416      CALL iom_put( TRIM(tname(jh))//'y_u', out_u(:,:,nb_ana+jh) ) 
    428417      END DO 
    429 #endif 
     418!#endif 
    430419 
    431420      ! C) vbar 
    432421      !///////// 
    433422      ! 
    434 #if defined key_dimgout 
    435       cltext='vbar amplitude and phase' 
    436       CALL dia_wri_dimg(TRIM(cdfile_name_V), TRIM(cltext), out_v, 2*nb_ana, '2') 
    437 #else 
     423!#if  
     424!#else 
    438425      DO jh = 1, nb_ana 
    439426         CALL iom_put( TRIM(tname(jh))//'x_v', out_v(:,:,jh       ) ) 
    440427         CALL iom_put( TRIM(tname(jh))//'y_v', out_v(:,:,jh+nb_ana) ) 
    441428      END DO 
    442 #endif 
     429!#endif 
    443430      ! 
    444431   END SUBROUTINE dia_wri_harm 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/dianam.F90

    r2528 r5282  
    7272 
    7373      IF( llfsec .OR. kfreq < 0 ) THEN   ;   inbsec = kfreq                       ! output frequency already in seconds 
    74       ELSE                               ;   inbsec = kfreq * NINT( rdttra(1) )   ! from time-step to seconds 
     74      ELSE                               ;   inbsec = kfreq * NINT( rdt )   ! from time-step to seconds 
    7575      ENDIF 
    7676      iddss = NINT( rday          )                                         ! number of seconds in 1 day 
     
    116116      ! date of the beginning and the end of the run 
    117117 
    118       zdrun = rdttra(1) / rday * REAL( nitend - nit000, wp )                ! length of the run in days 
    119       zjul  = fjulday - rdttra(1) / rday 
     118      zdrun = rdt / rday * REAL( nitend - nit000, wp )                ! length of the run in days 
     119      zjul  = fjulday - rdt / rday 
    120120      CALL ju2ymds( zjul        , iyear1, imonth1, iday1, zsec1 )           ! year/month/day of the beginning of run 
    121121      CALL ju2ymds( zjul + zdrun, iyear2, imonth2, iday2, zsec2 )           ! year/month/day of the end       of run 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/diaptr.F90

    r4990 r5282  
    585585         niter = ( nit000 - 1 ) / nn_fptr 
    586586         zdt = rdt 
    587          IF( nacc == 1 )   zdt = rdtmin 
    588587         ! 
    589588         IF(lwp) THEN 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r4990 r5282  
    4343   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    4444   USE in_out_manager  ! I/O manager 
    45    USE diadimg         ! dimg direct access file format output 
    4645   USE iom 
    4746   USE ioipsl 
     
    9796  END FUNCTION dia_wri_alloc 
    9897 
    99 #if defined key_dimgout 
    100    !!---------------------------------------------------------------------- 
    101    !!   'key_dimgout'                                      DIMG output file 
    102    !!---------------------------------------------------------------------- 
    103 #   include "diawri_dimg.h90" 
    104  
    105 #else 
    10698   !!---------------------------------------------------------------------- 
    10799   !!   Default option                                   NetCDF output file 
    108100   !!---------------------------------------------------------------------- 
    109 # if defined key_iomput 
     101#if defined key_iomput 
    110102   !!---------------------------------------------------------------------- 
    111103   !!   'key_iomput'                                        use IOM library 
     
    418410      ! Define frequency of output and means 
    419411      zdt = rdt 
    420       IF( nacc == 1 ) zdt = rdtmin 
    421412      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
    422413      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     
    762753      IF( lk_vvl ) THEN 
    763754         zw3d(:,:,:) = ( ( fse3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 
    764          CALL histwrite( nid_T, "vovvle3t", it, fse3t_n (:,:,:) , ndim_T , ndex_T  )   ! level thickness 
     755         CALL histwrite( nid_T, "vovvle3t", it, fse3t_n(:,:,:) , ndim_T , ndex_T  )   ! level thickness 
    765756         CALL histwrite( nid_T, "vovvldep", it, fsdept_n(:,:,:) , ndim_T , ndex_T  )   ! t-point depth 
    766757         CALL histwrite( nid_T, "vovvldef", it, zw3d             , ndim_T , ndex_T  )   ! level thickness deformation 
     
    914905      ! 
    915906   END SUBROUTINE dia_wri 
    916 # endif 
    917  
    918 #endif 
     907#endif 
     908 
     909 
    919910 
    920911   SUBROUTINE dia_wri_state( cdfile_name, kt ) 
     
    10351026      ! ----------------- 
    10361027      CALL histclo( id_i ) 
    1037 #if ! defined key_iomput && ! defined key_dimgout 
     1028#if ! defined key_iomput 
    10381029      IF( ninist /= 1  ) THEN 
    10391030         CALL histclo( nid_T ) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DOM/daymod.F90

    r5002 r5282  
    2020   !! 
    2121   !!   we suppose that the time step is deviding the number of second of in a day 
    22    !!             ---> MOD( rday, rdttra(1) ) == 0 
     22   !!             ---> MOD( rday, rdt ) == 0 
    2323   !! 
    2424   !!           ----------- WARNING ----------- 
     
    7373      !!---------------------------------------------------------------------- 
    7474      ! 
    75       ! all calendar staff is based on the fact that MOD( rday, rdttra(1) ) == 0 
    76       IF( MOD( rday     , rdttra(1) ) /= 0. )   CALL ctl_stop( 'the time step must devide the number of second of in a day' ) 
     75      ! all calendar staff is based on the fact that MOD( rday, rdt ) == 0 
     76      IF( MOD( rday     , rdt ) /= 0. )   CALL ctl_stop( 'the time step must devide the number of second of in a day' ) 
    7777      IF( MOD( rday     , 2.        ) /= 0. )   CALL ctl_stop( 'the number of second of in a day must be an even number'    ) 
    78       IF( MOD( rdttra(1), 2.        ) /= 0. )   CALL ctl_stop( 'the time step (in second) must be an even number'           ) 
     78      IF( MOD( rdt, 2.        ) /= 0. )   CALL ctl_stop( 'the time step (in second) must be an even number'           ) 
    7979      nsecd   = NINT(rday           ) 
    8080      nsecd05 = NINT(0.5 * rday     ) 
    81       ndt     = NINT(      rdttra(1)) 
    82       ndt05   = NINT(0.5 * rdttra(1)) 
     81      ndt     = NINT(      rdt) 
     82      ndt05   = NINT(0.5 * rdt) 
    8383 
    8484      IF( .NOT. lk_offline ) CALL day_rst( nit000, 'READ' ) 
     
    218218      nsec_week  = nsec_week  + ndt 
    219219      nsec_day   = nsec_day   + ndt 
    220       adatrj  = adatrj  + rdttra(1) / rday 
    221       fjulday = fjulday + rdttra(1) / rday 
     220      adatrj  = adatrj  + rdt / rday 
     221      fjulday = fjulday + rdt / rday 
    222222      IF( ABS(fjulday - REAL(NINT(fjulday),wp)) < zprec )   fjulday = REAL(NINT(fjulday),wp)   ! avoid truncation error 
    223223      IF( ABS(adatrj  - REAL(NINT(adatrj ),wp)) < zprec )   adatrj  = REAL(NINT(adatrj ),wp)   ! avoid truncation error 
     
    334334               ! parameters correspondting to nit000 - 1 (as we start the step loop with a call to day) 
    335335               ndastp = ndate0 - 1     ! ndate0 read in the namelist in dom_nam, we assume that we start run at 00:00 
    336                adatrj = ( REAL( nit000-1, wp ) * rdttra(1) ) / rday 
     336               adatrj = ( REAL( nit000-1, wp ) * rdt ) / rday 
    337337               ! note this is wrong if time step has changed during run 
    338338            ENDIF 
     
    340340            ! parameters correspondting to nit000 - 1 (as we start the step loop with a call to day) 
    341341            ndastp = ndate0 - 1        ! ndate0 read in the namelist in dom_nam, we assume that we start run at 00:00 
    342             adatrj = ( REAL( nit000-1, wp ) * rdttra(1) ) / rday 
     342            adatrj = ( REAL( nit000-1, wp ) * rdt ) / rday 
    343343         ENDIF 
    344344         IF( ABS(adatrj  - REAL(NINT(adatrj),wp)) < 0.1 / rday )   adatrj = REAL(NINT(adatrj),wp)   ! avoid truncation error 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DOM/dom_oce.F90

    r4990 r5282  
    3535   REAL(wp), PUBLIC ::   rn_e3zps_rat    !: minimum thickness ration for partial steps 
    3636   INTEGER , PUBLIC ::   nn_msh          !: = 1 create a mesh-mask file 
    37    INTEGER , PUBLIC ::   nn_acc          !: = 0/1 use of the acceleration of convergence technique 
    3837   REAL(wp), PUBLIC ::   rn_atfp         !: asselin time filter parameter 
    39    REAL(wp), PUBLIC ::   rn_rdt          !: time step for the dynamics (and tracer if nacc=0) 
    40    REAL(wp), PUBLIC ::   rn_rdtmin       !: minimum time step on tracers 
    41    REAL(wp), PUBLIC ::   rn_rdtmax       !: maximum time step on tracers 
    42    REAL(wp), PUBLIC ::   rn_rdth         !: depth variation of tracer step 
     38   REAL(wp), PUBLIC ::   rn_rdt          !: time step for the dynamics and tracer 
    4339   INTEGER , PUBLIC ::   nn_closea       !: =0 suppress closed sea/lake from the ORCA domain or not (=1) 
    4440   INTEGER , PUBLIC ::   nn_euler        !: =0 start with forward time step or not (=1) 
     
    9490   REAL(wp), PUBLIC ::   e3zps_rat       !: minimum thickness ration for partial steps 
    9591   INTEGER , PUBLIC ::   nmsh            !: = 1 create a mesh-mask file 
    96    INTEGER , PUBLIC ::   nacc            !: = 0/1 use of the acceleration of convergence technique 
    9792   REAL(wp), PUBLIC ::   atfp            !: asselin time filter parameter 
    98    REAL(wp), PUBLIC ::   rdt             !: time step for the dynamics (and tracer if nacc=0) 
    99    REAL(wp), PUBLIC ::   rdtmin          !: minimum time step on tracers 
    100    REAL(wp), PUBLIC ::   rdtmax          !: maximum time step on tracers 
    101    REAL(wp), PUBLIC ::   rdth            !: depth variation of tracer step 
     93   REAL(wp), PUBLIC ::   rdt             !: time step for the dynamics and tracer 
    10294 
    10395   !                                                  !!! associated variables 
    10496   INTEGER , PUBLIC                 ::   neuler        !: restart euler forward option (0=Euler) 
    10597   REAL(wp), PUBLIC                 ::   atfp1         !: asselin time filter coeff. (atfp1= 1-2*atfp) 
    106    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   rdttra  !: vertical profile of tracer time step 
    107    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:) ::   r2dtra  !: = 2*rdttra except at nit000 (=rdttra) if neuler=0 
    108  
    109    !                                         !!* Namelist namcla : cross land advection 
    110    INTEGER, PUBLIC ::   nn_cla               !: =1 cross land advection for exchanges through some straits (ORCA2) 
     98   REAL(wp), PUBLIC ::   r2dt         !: = 2*rdt except at nit000 (=rdt) if neuler=0 
    11199 
    112100   !!---------------------------------------------------------------------- 
     
    336324      ierr(:) = 0 
    337325      ! 
    338       ALLOCATE( rdttra(jpk), r2dtra(jpk), mig(jpi), mjg(jpj), nfiimpp(jpni,jpnj),  & 
     326      ALLOCATE( mig(jpi), mjg(jpj), nfiimpp(jpni,jpnj),  & 
    339327         &      nfipproc(jpni,jpnj), nfilcit(jpni,jpnj), STAT=ierr(1) ) 
    340          ! 
     328        ! 
    341329      ALLOCATE( nimppt(jpnij) , ibonit(jpnij) , nlcit(jpnij) , nlcjt(jpnij) ,     & 
    342330         &      njmppt(jpnij) , ibonjt(jpnij) , nldit(jpnij) , nldjt(jpnij) ,     & 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

    r4990 r5282  
    8181      ENDIF 
    8282      ! 
    83                              CALL dom_nam      ! read namelist ( namrun, namdom, namcla ) 
     83                             CALL dom_nam      ! read namelist ( namrun, namdom ) 
    8484                             CALL dom_clo      ! Closed seas and lake 
    8585                             CALL dom_hgr      ! Horizontal mesh 
     
    131131      !! ** input   : - namrun namelist 
    132132      !!              - namdom namelist 
    133       !!              - namcla namelist 
    134133      !!              - namnc4 namelist   ! "key_netcdf4" only 
    135134      !!---------------------------------------------------------------------- 
     
    137136      NAMELIST/namrun/ nn_no   , cn_exp    , cn_ocerst_in, cn_ocerst_out, ln_rstart , nn_rstctl,   & 
    138137         &             nn_it000, nn_itend  , nn_date0    , nn_leapy     , nn_istate , nn_stock ,   & 
    139          &             nn_write, ln_dimgnnn, ln_mskland  , ln_clobber   , nn_chunksz, nn_euler 
     138         &             nn_write, ln_mskland  , ln_clobber   , nn_chunksz, nn_euler 
    140139      NAMELIST/namdom/ nn_bathy, rn_bathy , rn_e3zps_min, rn_e3zps_rat, nn_msh, rn_hmin,   & 
    141          &             nn_acc   , rn_atfp     , rn_rdt      , rn_rdtmin ,                  & 
    142          &             rn_rdtmax, rn_rdth     , nn_closea , ln_crs,    & 
    143          &             jphgr_msh, & 
     140         &             rn_atfp, rn_rdt ,  nn_closea , ln_crs,  jphgr_msh,       & 
    144141         &             ppglam0, ppgphi0, ppe1_deg, ppe2_deg, ppe1_m, ppe2_m, & 
    145142         &             ppsur, ppa0, ppa1, ppkth, ppacr, ppdzmin, pphmax, ldbletanh, & 
    146143         &             ppa2, ppkth2, ppacr2 
    147       NAMELIST/namcla/ nn_cla 
     144 
    148145#if defined key_netcdf4 
    149146      NAMELIST/namnc4/ nn_nchunks_i, nn_nchunks_j, nn_nchunks_k, ln_nc4zip 
     
    180177         WRITE(numout,*) '      frequency of restart file       nn_stock   = ', nn_stock 
    181178         WRITE(numout,*) '      frequency of output file        nn_write   = ', nn_write 
    182          WRITE(numout,*) '      multi file dimgout              ln_dimgnnn = ', ln_dimgnnn 
    183179         WRITE(numout,*) '      mask land points                ln_mskland = ', ln_mskland 
    184180         WRITE(numout,*) '      overwrite an existing file      ln_clobber = ', ln_clobber 
     
    259255         WRITE(numout,*) '      ocean time step                       rn_rdt    = ', rn_rdt 
    260256         WRITE(numout,*) '      asselin time filter parameter         rn_atfp   = ', rn_atfp 
    261          WRITE(numout,*) '      acceleration of converge              nn_acc    = ', nn_acc 
    262          WRITE(numout,*) '        nn_acc=1: surface tracer rdt        rn_rdtmin = ', rn_rdtmin 
    263          WRITE(numout,*) '                  bottom  tracer rdt        rdtmax    = ', rn_rdtmax 
    264          WRITE(numout,*) '                  depth of transition       rn_rdth   = ', rn_rdth 
    265257         WRITE(numout,*) '      suppression of closed seas (=0)       nn_closea = ', nn_closea 
    266258         WRITE(numout,*) '      online coarsening of dynamical fields ln_crs    = ', ln_crs 
     
    289281      e3zps_rat = rn_e3zps_rat 
    290282      nmsh      = nn_msh 
    291       nacc      = nn_acc 
    292283      atfp      = rn_atfp 
    293284      rdt       = rn_rdt 
    294       rdtmin    = rn_rdtmin 
    295       rdtmax    = rn_rdtmin 
    296       rdth      = rn_rdth 
    297  
    298       REWIND( numnam_ref )              ! Namelist namcla in reference namelist : Cross land advection 
    299       READ  ( numnam_ref, namcla, IOSTAT = ios, ERR = 905) 
    300 905   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in reference namelist', lwp ) 
    301  
    302       REWIND( numnam_cfg )              ! Namelist namcla in configuration namelist : Cross land advection 
    303       READ  ( numnam_cfg, namcla, IOSTAT = ios, ERR = 906 ) 
    304 906   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in configuration namelist', lwp ) 
    305       IF(lwm) WRITE( numond, namcla ) 
    306  
    307       IF(lwp) THEN 
    308          WRITE(numout,*) 
    309          WRITE(numout,*) '   Namelist namcla' 
    310          WRITE(numout,*) '      cross land advection                 nn_cla    = ', nn_cla 
    311       ENDIF 
    312       IF ( nn_cla .EQ. 1 ) THEN 
    313          IF  ( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN   ! ORCA R2  
    314             CONTINUE 
    315          ELSE 
    316             CALL ctl_stop( 'STOP', 'Cross land advation iplemented only for ORCA2 configuration: cp_cfg = "orca" and jp_cfg = 2 ' ) 
    317          ENDIF 
    318       ENDIF 
    319285 
    320286#if defined key_netcdf4 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90

    r4990 r5282  
    134134         IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN    ! ORCA R2 configuration 
    135135            !                                             ! ===================== 
    136             IF( nn_cla == 0 ) THEN 
    137                ! 
    138136               ii0 = 139   ;   ii1 = 140        ! Gibraltar Strait (e2u = 20 km) 
    139137               ij0 = 102   ;   ij1 = 102   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  20.e3 
     
    147145               IF(lwp) WRITE(numout,*) '             orca_r2: Bab el Mandeb: e2u reduced to 30 km' 
    148146               IF(lwp) WRITE(numout,*) '                                     e1v reduced to 18 km' 
    149             ENDIF 
    150147 
    151148            ii0 = 145   ;   ii1 = 146        ! Danish Straits (e2u = 10 km) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90

    r4990 r5282  
    366366      IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN   ! ORCA_R2 configuration 
    367367         !                                                 ! Increased lateral friction near of some straits 
    368          IF( nn_cla == 0 ) THEN 
    369368            !                                ! Gibraltar strait  : partial slip (fmask=0.5) 
    370369            ij0 = 101   ;   ij1 = 101 
     
    379378            ii0 = 159   ;   ii1 = 159   ;   fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) =  1._wp 
    380379            ! 
    381          ENDIF 
    382380         !                                ! Danish straits  : strong slip (fmask > 2) 
    383381! We keep this as an example but it is instable in this case  
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DOM/domstp.F90

    r4292 r5282  
    4141      !!      filter parameter read in namelist 
    4242      !!              - Model time step: 
    43       !!      nacc = 0 : synchronous time intergration.  
    44       !!      There is one time step only, defined by: rdt, rdttra(k)=rdt 
    45       !!      nacc = 1 : accelerating the convergence. There is 2 different 
    46       !!      time steps for dynamics and tracers: 
    47       !!        rdt      : dynamical part 
    48       !!        rdttra(k): temperature and salinity 
    49       !!      The tracer time step is a function of vertical level. the model 
    50       !!      reference time step ( i.e. for wind stress, surface heat and 
    51       !!      salt fluxes) is the surface tracer time step is rdttra(1). 
    52       !!         N.B. depth dependent acceleration of convergence is not im- 
    53       !!      plemented for s-coordinate. 
     43      !!                synchronous time intergration. 
     44      !!      There is one time step only, defined by: rdt for dynamics and 
     45      !!      tracer,wind stress, surface heat and salt fluxes 
    5446      !! 
    55       !! ** Action  : - rdttra   : vertical profile of tracer time step 
     47      !! ** Action  : [REMOVED - rdttra: vertical profile of tracer time step] 
    5648      !!              - atfp1    : = 1 - 2*atfp 
    5749      !! 
     
    7264      atfp1 = 1. - 2. * atfp 
    7365 
    74       SELECT CASE ( nacc ) 
     66      IF(lwp) WRITE(numout,*)'               synchronous time stepping' 
     67      IF(lwp) WRITE(numout,*)'               dynamics and tracer time step = ', rdt/3600., ' hours' 
    7568 
    76          CASE ( 0 )                ! Synchronous time stepping 
    77             IF(lwp) WRITE(numout,*)'               synchronous time stepping' 
    78             IF(lwp) WRITE(numout,*)'               dynamics and tracer time step = ', rdt/3600., ' hours' 
    79  
    80             rdttra(:) = rdt 
    81  
    82          CASE ( 1 )                ! Accelerating the convergence 
    83             IF(lwp) WRITE(numout,*) '              no tracer damping in the turbocline' 
    84             IF(lwp) WRITE(numout,*)'               accelerating the convergence' 
    85             IF(lwp) WRITE(numout,*)'               dynamics time step = ', rdt/3600., ' hours' 
    86             IF( ln_sco .AND. rdtmin /= rdtmax .AND. lk_vvl )   & 
    87                  & CALL ctl_stop ( ' depth dependent acceleration of convergence not implemented in s-coordinates & 
    88                  &                   nor in variable volume' ) 
    89             IF(lwp) WRITE(numout,*)'         tracers   time step :  dt (hours)  level' 
    90  
    91             DO jk = 1, jpk 
    92                IF( gdept_1d(jk) <= rdth ) rdttra(jk) = rdtmin 
    93                IF( gdept_1d(jk) >  rdth ) THEN 
    94                   rdttra(jk) = rdtmin + ( rdtmax - rdtmin )   & 
    95                                       * ( EXP( ( gdept_1d(jk ) - rdth ) / rdth ) - 1. )   & 
    96                                       / ( EXP( ( gdept_1d(jpk) - rdth ) / rdth ) - 1. ) 
    97                ENDIF 
    98                IF(lwp) WRITE(numout,"(36x,f5.2,5x,i3)") rdttra(jk)/3600., jk 
    99             END DO   
    100  
    101          CASE DEFAULT              ! E R R O R  
    102  
    103             WRITE(ctmp1,*) ' nacc value e r r o r, nacc= ',nacc 
    104             CALL ctl_stop( ctmp1 ) 
    105  
    106       END SELECT 
    10769 
    10870   END SUBROUTINE dom_stp 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

    r4998 r5282  
    156156      ! Vertical scale factor interpolations 
    157157      ! ------------------------------------ 
    158       CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3w_n (:,:,:), 'W'  ) 
     158      CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3w_n(:,:,:), 'W'  ) 
    159159      CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3uw_n(:,:,:), 'UW' ) 
    160160      CALL dom_vvl_interpol( fse3v_n(:,:,:), fse3vw_n(:,:,:), 'VW' ) 
    161       CALL dom_vvl_interpol( fse3t_b(:,:,:), fse3w_b (:,:,:), 'W'  ) 
     161      CALL dom_vvl_interpol( fse3t_b(:,:,:), fse3w_b(:,:,:), 'W'  ) 
    162162      CALL dom_vvl_interpol( fse3u_b(:,:,:), fse3uw_b(:,:,:), 'UW' ) 
    163163      CALL dom_vvl_interpol( fse3v_b(:,:,:), fse3vw_b(:,:,:), 'VW' ) 
     
    627627      ! - ML - fse3u_b and fse3v_b are allready computed in dynnxt 
    628628      ! - JC - hu_b, hv_b, hur_b, hvr_b also 
    629       CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3f_n (:,:,:), 'F'  ) 
     629      CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3f_n(:,:,:), 'F'  ) 
    630630      ! Vertical scale factor interpolations 
    631631      ! ------------------------------------ 
    632       CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3w_n (:,:,:), 'W'  ) 
     632      CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3w_n(:,:,:), 'W'  ) 
    633633      CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3uw_n(:,:,:), 'UW' ) 
    634634      CALL dom_vvl_interpol( fse3v_n(:,:,:), fse3vw_n(:,:,:), 'VW' ) 
    635       CALL dom_vvl_interpol( fse3t_b(:,:,:), fse3w_b (:,:,:), 'W'  ) 
     635      CALL dom_vvl_interpol( fse3t_b(:,:,:), fse3w_b(:,:,:), 'W'  ) 
    636636      CALL dom_vvl_interpol( fse3u_b(:,:,:), fse3uw_b(:,:,:), 'UW' ) 
    637637      CALL dom_vvl_interpol( fse3v_b(:,:,:), fse3vw_b(:,:,:), 'VW' ) 
     
    663663      ! Local depth and Inverse of the local depth of the water column at u- and v- points 
    664664      ! ---------------------------------------------------------------------------------- 
    665       hu (:,:) = hu_a (:,:) 
    666       hv (:,:) = hv_a (:,:) 
     665      hu(:,:) = hu_a(:,:) 
     666      hv(:,:) = hv_a(:,:) 
    667667 
    668668      ! Inverse of the local depth 
     
    680680      ! ============= 
    681681      z_e3t_def(:,:,:) = ( ( fse3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 
    682       CALL iom_put( "cellthc" , fse3t_n  (:,:,:) ) 
    683       CALL iom_put( "tpt_dep" , fsde3w_n (:,:,:) ) 
     682      CALL iom_put( "cellthc" , fse3t_n(:,:,:) ) 
     683      CALL iom_put( "tpt_dep" , fsde3w_n(:,:,:) ) 
    684684      CALL iom_put( "e3tdef"  , z_e3t_def(:,:,:) ) 
    685685 
     
    10571057         !                                             ! ===================== 
    10581058      !! acc 
    1059          IF( nn_cla == 0 ) THEN 
    1060             ! 
    10611059            ii0 = 139   ;   ii1 = 140        ! Gibraltar Strait (e2u was modified) 
    10621060            ij0 = 102   ;   ij1 = 102 
     
    11051103               END DO 
    11061104            END DO 
    1107          ENDIF 
    11081105 
    11091106         ii0 = 145   ;   ii1 = 146        ! Danish Straits (e2u was modified) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r5040 r5282  
    507507            IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN    ! ORCA R2 configuration 
    508508               !                                             ! ===================== 
    509                IF( nn_cla == 0 ) THEN 
    510                   ii0 = 140   ;   ii1 = 140                  ! Gibraltar Strait open  
     509                  ii0 = 140   ;   ii1 = 140                  ! Gibraltar Strait open 
    511510                  ij0 = 102   ;   ij1 = 102                  ! (Thomson, Ocean Modelling, 1995) 
    512511                  DO ji = mi0(ii0), mi1(ii1) 
     
    527526                  IF(lwp) WRITE(numout,*) 
    528527                  IF(lwp) WRITE(numout,*) '      orca_r2: Bab el Mandeb strait open at i=',ii0,' j=',ij0 
    529                ENDIF 
    530528               ! 
    531529            ENDIF 
     
    548546            IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN    ! ORCA R2 configuration 
    549547               ! 
    550               IF( nn_cla == 0 ) THEN 
    551                  ii0 = 140   ;   ii1 = 140                   ! Gibraltar Strait open  
     548                 ii0 = 140   ;   ii1 = 140                   ! Gibraltar Strait open 
    552549                 ij0 = 102   ;   ij1 = 102                   ! (Thomson, Ocean Modelling, 1995) 
    553550                 DO ji = mi0(ii0), mi1(ii1) 
     
    568565                 IF(lwp) WRITE(numout,*) 
    569566                 IF(lwp) WRITE(numout,*) '             orca_r2: Bab el Mandeb strait open at i=',ii0,' j=',ij0 
    570               ENDIF 
    571567              ! 
    572568           ENDIF 
     
    19951991 
    19961992      IF( nprint == 1 .AND. lwp )   THEN 
    1997          WRITE(numout,*) ' MAX val hif   t ', MAXVAL( hift (:,:) ), ' f ', MAXVAL( hiff (:,:) ),  & 
    1998             &                        ' u ',   MAXVAL( hifu (:,:) ), ' v ', MAXVAL( hifv (:,:) ) 
    1999          WRITE(numout,*) ' MIN val hif   t ', MINVAL( hift (:,:) ), ' f ', MINVAL( hiff (:,:) ),  & 
    2000             &                        ' u ',   MINVAL( hifu (:,:) ), ' v ', MINVAL( hifv (:,:) ) 
     1993         WRITE(numout,*) ' MAX val hif   t ', MAXVAL( hift(:,:) ), ' f ', MAXVAL( hiff(:,:) ),  & 
     1994            &                        ' u ',   MAXVAL( hifu(:,:) ), ' v ', MAXVAL( hifv(:,:) ) 
     1995         WRITE(numout,*) ' MIN val hif   t ', MINVAL( hift(:,:) ), ' f ', MINVAL( hiff(:,:) ),  & 
     1996            &                        ' u ',   MINVAL( hifu(:,:) ), ' v ', MINVAL( hifv(:,:) ) 
    20011997         WRITE(numout,*) ' MAX val hbat  t ', MAXVAL( hbatt(:,:) ), ' f ', MAXVAL( hbatf(:,:) ),  & 
    20021998            &                        ' u ',   MAXVAL( hbatu(:,:) ), ' v ', MAXVAL( hbatv(:,:) ) 
     
    20682064#endif 
    20692065 
    2070       fsdept(:,:,:) = gdept_0 (:,:,:) 
    2071       fsdepw(:,:,:) = gdepw_0 (:,:,:) 
     2066      fsdept(:,:,:) = gdept_0(:,:,:) 
     2067      fsdepw(:,:,:) = gdepw_0(:,:,:) 
    20722068      fsde3w(:,:,:) = gdep3w_0(:,:,:) 
    2073       fse3t (:,:,:) = e3t_0   (:,:,:) 
    2074       fse3u (:,:,:) = e3u_0   (:,:,:) 
    2075       fse3v (:,:,:) = e3v_0   (:,:,:) 
    2076       fse3f (:,:,:) = e3f_0   (:,:,:) 
    2077       fse3w (:,:,:) = e3w_0   (:,:,:) 
    2078       fse3uw(:,:,:) = e3uw_0  (:,:,:) 
    2079       fse3vw(:,:,:) = e3vw_0  (:,:,:) 
     2069      fse3t(:,:,:) = e3t_0(:,:,:) 
     2070      fse3u(:,:,:) = e3u_0(:,:,:) 
     2071      fse3v(:,:,:) = e3v_0(:,:,:) 
     2072      fse3f(:,:,:) = e3f_0(:,:,:) 
     2073      fse3w(:,:,:) = e3w_0(:,:,:) 
     2074      fse3uw(:,:,:) = e3uw_0(:,:,:) 
     2075      fse3vw(:,:,:) = e3vw_0(:,:,:) 
    20802076!! 
    20812077      ! HYBRID :  
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DOM/dtatsd.F90

    r4990 r5282  
    174174            END DO 
    175175         END DO 
    176          IF( nn_cla == 1 ) THEN                          ! Cross Land advection 
    177             il0 = 138   ;   il1 = 138                          ! set T & S profile at Gibraltar Strait 
    178             ij0 = 101   ;   ij1 = 102 
    179             ii0 = 139   ;   ii1 = 139 
    180             DO jl = mi0(il0), mi1(il1) 
    181                DO jj = mj0(ij0), mj1(ij1) 
    182                   DO ji = mi0(ii0), mi1(ii1) 
    183                      sf_tsd(jp_tem)%fnow(ji,jj,:) = sf_tsd(jp_tem)%fnow(jl,jj,:) 
    184                      sf_tsd(jp_sal)%fnow(ji,jj,:) = sf_tsd(jp_sal)%fnow(jl,jj,:) 
    185                   END DO 
    186                END DO 
    187             END DO 
    188             il0 = 164   ;   il1 = 164                          ! set T & S profile at Bab el Mandeb Strait 
    189             ij0 =  87   ;   ij1 =  88 
    190             ii0 = 161   ;   ii1 = 163 
    191             DO jl = mi0(il0), mi1(il1) 
    192                DO jj = mj0(ij0), mj1(ij1) 
    193                   DO ji = mi0(ii0), mi1(ii1) 
    194                      sf_tsd(jp_tem)%fnow(ji,jj,:) = sf_tsd(jp_tem)%fnow(jl,jj,:) 
    195                      sf_tsd(jp_sal)%fnow(ji,jj,:) = sf_tsd(jp_sal)%fnow(jl,jj,:) 
    196                   END DO 
    197                END DO 
    198             END DO 
    199          ELSE                                            ! No Cross Land advection 
    200176            ij0 =  87   ;   ij1 =  96                          ! Reduced temperature in Red Sea 
    201177            ii0 = 148   ;   ii1 = 160 
     
    203179            sf_tsd(jp_tem)%fnow( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 11:13 ) = 6.5_wp 
    204180            sf_tsd(jp_tem)%fnow( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 14:20 ) = 6.0_wp 
    205          ENDIF 
    206181      ENDIF 
    207182      ! 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DYN/divcur.F90

    r4990 r5282  
    1717   !!            3.3  ! 2010-09  (D.Storkey and E.O'Dea) bug fixes for BDY module 
    1818   !!             -   ! 2010-10  (R. Furner, G. Madec) runoff and cla added directly here 
     19   !!            3.6  ! 2014-15  cla removed 
    1920   !!---------------------------------------------------------------------- 
    2021 
     
    2829   USE sbcrnf          ! river runoff  
    2930   USE sbcisf          ! ice shelf  
    30    USE cla             ! cross land advection             (cla_div routine) 
    3131   USE in_out_manager  ! I/O manager 
    3232   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
     
    6868      !!         hdivn = 1/(e1t*e2t*e3t) ( di[e2u*e3u un] + dj[e1v*e3v vn] ) 
    6969      !!      correct hdiv with runoff inflow (div_rnf), ice shelf melting (div_isf) 
    70       !!      and cross land flow (div_cla)  
     70      !!      [REMOVED cross land flow (div_cla)] 
    7171      !!              II. vorticity : 
    7272      !!         - save the curl computed at the previous time-step 
     
    229229      IF( ln_rnf      )   CALL sbc_rnf_div( hdivn )          ! runoffs   (update hdivn field) 
    230230      IF( ln_divisf .AND. (nn_isf /= 0) )   CALL sbc_isf_div( hdivn )          ! ice shelf (update hdivn field) 
    231       IF( nn_cla == 1 )   CALL cla_div    ( kt )             ! Cross Land Advection (Update Hor. divergence) 
    232        
     231 
    233232      ! 4. Lateral boundary conditions on hdivn and rotn 
    234233      ! ---------------------------------=======---====== 
     
    259258      !!      - compute the now divergence given by : 
    260259      !!         hdivn = 1/(e1t*e2t*e3t) ( di[e2u*e3u un] + dj[e1v*e3v vn] ) 
    261       !!      correct hdiv with runoff inflow (div_rnf) and cross land flow (div_cla)  
     260      !!      correct hdiv with runoff inflow (div_rnf)  
     261      !!      [REMOVED cross land flow (div_cla)] 
    262262      !!              - Relavtive Vorticity : 
    263263      !!      - save the curl computed at the previous time-step (rotb = rotn) 
     
    328328      IF( ln_rnf      )   CALL sbc_rnf_div( hdivn )                            ! runoffs (update hdivn field) 
    329329      IF( ln_divisf .AND. (nn_isf .GT. 0) )   CALL sbc_isf_div( hdivn )          ! ice shelf (update hdivn field) 
    330       IF( nn_cla == 1 )   CALL cla_div    ( kt )             ! Cross Land Advection (update hdivn field) 
    331330      ! 
    332331      CALL lbc_lnk( hdivn, 'T', 1. )   ;   CALL lbc_lnk( rotn , 'F', 1. )     ! lateral boundary cond. (no sign change) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90

    r4990 r5282  
    753753               &                   * (  rhd(ji,jj,1)                                    & 
    754754               &                     + 0.5_wp * ( rhd(ji,jj,2) - rhd(ji,jj,1) )         & 
    755                &                              * ( fse3w (ji,jj,1) - fsde3w(ji,jj,1) )   & 
     755               &                              * ( fse3w(ji,jj,1) - fsde3w(ji,jj,1) )   & 
    756756               &                              / ( fsde3w(ji,jj,2) - fsde3w(ji,jj,1) )  ) 
    757757         END DO 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_bilap.F90

    r4990 r5282  
    140140          
    141141         ! Multiply by the eddy viscosity coef. (at u- and v-points) 
    142          zlu(:,:,jk) = zlu(:,:,jk) * ( fsahmu(:,:,jk) * (1-nkahm_smag) + nkahm_smag) 
    143  
    144          zlv(:,:,jk) = zlv(:,:,jk) * ( fsahmv(:,:,jk) * (1-nkahm_smag) + nkahm_smag) 
    145           
     142          zlu(:,:,jk) = zlu(:,:,jk) * ( fsahmu(:,:,jk) ) 
     143 
     144          zlv(:,:,jk) = zlv(:,:,jk) * ( fsahmv(:,:,jk) ) 
    146145         ! Contravariant "laplacian" 
    147146         zcu(:,:) = e1u(:,:) * zlu(:,:,jk) 
     
    197196                  &  + ( zut(ji,jj+1,jk) - zut(ji  ,jj,jk) ) / e2v(ji,jj) 
    198197               ! add it to the general momentum trends 
    199                ua(ji,jj,jk) = ua(ji,jj,jk) + zua * ( fsahmu(ji,jj,jk)*nkahm_smag +(1 -nkahm_smag )) 
    200                va(ji,jj,jk) = va(ji,jj,jk) + zva * ( fsahmv(ji,jj,jk)*nkahm_smag +(1 -nkahm_smag )) 
     198               ua(ji,jj,jk) = ua(ji,jj,jk) + zua 
     199               va(ji,jj,jk) = va(ji,jj,jk) + zva 
    201200            END DO 
    202201         END DO 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_bilapg.F90

    r4990 r5282  
    434434         ! II.3 Divergence of vertical fluxes added to the horizontal divergence 
    435435         ! --------------------------------------------------------------------- 
    436          IF( (kahm -nkahm_smag) ==1 ) THEN 
     436         IF( kahm ==1 ) THEN 
    437437            ! multiply the laplacian by the eddy viscosity coefficient 
    438438            DO jk = 1, jpkm1 
     
    449449               END DO 
    450450            END DO 
    451          ELSEIF( (kahm +nkahm_smag ) == 2 ) THEN 
     451         ELSEIF( kahm == 2 ) THEN 
    452452            ! second call, no multiplication 
    453453            DO jk = 1, jpkm1 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg.F90

    r4990 r5282  
    272272#endif 
    273273 
    274       !                        ! Control of timestep choice 
    275       IF( lk_dynspg_ts .OR. lk_dynspg_exp ) THEN 
    276          IF( nn_cla == 1 )   CALL ctl_stop( 'Crossland advection not implemented for this free surface formulation' ) 
    277       ENDIF 
    278  
    279274      !               ! Control of hydrostatic pressure choice 
    280275      IF( lk_dynspg_ts .AND. ln_dynhpg_imp ) THEN 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_flt.F90

    r4990 r5282  
    3636   USE bdydyn          ! ocean open boundary condition on dynamics 
    3737   USE bdyvol          ! ocean open boundary condition (bdy_vol routine) 
    38    USE cla             ! cross land advection 
    3938   USE trd_oce         ! trends: ocean variables 
    4039   USE trddyn          ! trend manager: dynamics 
     
    206205      CALL Agrif_dyn( kt )    ! Update velocities on each coarse/fine interfaces  
    207206#endif 
    208       IF( nn_cla == 1 .AND. cp_cfg == 'orca' .AND. jp_cfg == 2 )   CALL cla_dynspg( kt )      ! Cross Land Advection (update (ua,va)) 
    209207 
    210208      ! compute the next vertically averaged velocity (effect of the additional force not included) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf.F90

    r4990 r5282  
    3535 
    3636   INTEGER  ::   nzdf = 0   ! type vertical diffusion algorithm used, defined from ln_zdf... namlist logicals 
    37    REAL(wp) ::   r2dt       ! time-step, = 2 rdttra except at nit000 (=rdttra) if neuler=0 
    3837 
    3938   !! * Substitutions 
     
    6463      ! 
    6564      !                                          ! set time step 
    66       IF( neuler == 0 .AND. kt == nit000     ) THEN   ;   r2dt =      rdt   ! = rdtra (restart with Euler time stepping) 
    67       ELSEIF(               kt <= nit000 + 1 ) THEN   ;   r2dt = 2. * rdt   ! = 2 rdttra (leapfrog) 
     65      IF( neuler == 0 .AND. kt == nit000     ) THEN   ;   r2dt =      rdt   ! = rdt(restart with Euler time stepping) 
     66      ELSEIF(               kt <= nit000 + 1 ) THEN   ;   r2dt = 2. * rdt   ! = 2 rdt (leapfrog) 
    6867      ENDIF 
    6968 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/IOM/in_out_manager.F90

    r4990 r5282  
    3838   INTEGER       ::   nn_write         !: model standard output frequency 
    3939   INTEGER       ::   nn_stock         !: restart file frequency 
    40    LOGICAL       ::   ln_dimgnnn       !: type of dimgout. (F): 1 file for all proc 
    41                                                        !:                  (T): 1 file per proc 
    4240   LOGICAL       ::   ln_mskland       !: mask land points in NetCDF outputs (costly: + ~15%) 
    4341   LOGICAL       ::   ln_clobber       !: clobber (overwrite) an existing file 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/IOM/iom.F90

    r5043 r5282  
    88   !!            3.0  ! 2007-07  (D. Storkey) Changes to iom_gettime 
    99   !!            3.4  ! 2012-12  (R. Bourdalle-Badie and G. Reffray)  add C1D case   
     10   !!            3.6  ! 2014-15  DIMG format removed 
    1011   !!-------------------------------------------------------------------- 
    1112 
     
    2324   USE lbclnk          ! lateal boundary condition / mpp exchanges 
    2425   USE iom_def         ! iom variables definitions 
    25    USE iom_ioipsl      ! NetCDF format with IOIPSL library 
    2626   USE iom_nf90        ! NetCDF format with native NetCDF library 
    27    USE iom_rstdimg     ! restarts access direct format "dimg" style... 
    2827   USE in_out_manager  ! I/O manager 
    2928   USE lib_mpp           ! MPP library 
     
    205204      CHARACTER(LEN=256)    ::   clname    ! the name of the file based on cdname [[+clcpu]+clcpu] 
    206205      CHARACTER(LEN=256)    ::   cltmpn    ! tempory name to store clname (in writting mode) 
    207       CHARACTER(LEN=10)     ::   clsuffix  ! ".nc" or ".dimg" 
     206      CHARACTER(LEN=10)     ::   clsuffix  ! ".nc"  
    208207      CHARACTER(LEN=15)     ::   clcpu     ! the cpu number (max jpmax_digits digits) 
    209208      CHARACTER(LEN=256)    ::   clinfo    ! info character 
     
    268267      ! which suffix should we use? 
    269268      SELECT CASE (iolib) 
    270       CASE (jpioipsl ) ;   clsuffix = '.nc' 
    271269      CASE (jpnf90   ) ;   clsuffix = '.nc' 
    272       CASE (jprstdimg) ;   clsuffix = '.dimg' 
    273270      CASE DEFAULT     ;   clsuffix = '' 
    274          CALL ctl_stop( TRIM(clinfo), 'accepted IO library are only jpioipsl, jpnf90 and jprstdimg' ) 
     271         CALL ctl_stop( TRIM(clinfo), 'accepted IO library is only jpnf90 (jpioipsl option has been removed) ' ) 
    275272      END SELECT 
    276273      ! Add the suffix if needed 
     
    285282      IF( .NOT.llok ) THEN 
    286283         ! we try to add the cpu number to the name 
    287          IF( iolib == jprstdimg ) THEN   ;   WRITE(clcpu,*) narea 
    288          ELSE                            ;   WRITE(clcpu,*) narea-1 
    289          ENDIF 
     284         WRITE(clcpu,*) narea-1 
     285 
    290286         clcpu  = TRIM(ADJUSTL(clcpu)) 
    291287         iln = INDEX(clname,TRIM(clsuffix), back = .TRUE.) 
     
    334330         END SELECT 
    335331      ENDIF 
    336       ! Open the NetCDF or RSTDIMG file 
     332      ! Open the NetCDF file 
    337333      ! ============= 
    338334      ! do we have some free file identifier? 
     
    358354      IF( istop == nstop ) THEN   ! no error within this routine 
    359355         SELECT CASE (iolib) 
    360          CASE (jpioipsl )   ;   CALL iom_ioipsl_open(  clname, kiomid, llwrt, llok, idompar ) 
    361356         CASE (jpnf90   )   ;   CALL iom_nf90_open(    clname, kiomid, llwrt, llok, idompar ) 
    362          CASE (jprstdimg)   ;   CALL iom_rstdimg_open( clname, kiomid, llwrt, llok, idompar ) 
    363357         CASE DEFAULT 
    364             CALL ctl_stop( TRIM(clinfo)//' accepted IO library are only jpioipsl, jpnf90 and jprstdimg' ) 
     358            CALL ctl_stop( TRIM(clinfo)//' accepted IO library is only jpnf90 (jpioipsl option has been removed) ' ) 
    365359         END SELECT 
    366360      ENDIF 
     
    397391            IF( iom_file(jf)%nfid > 0 ) THEN 
    398392               SELECT CASE (iom_file(jf)%iolib) 
    399                CASE (jpioipsl )   ;   CALL iom_ioipsl_close(  jf ) 
    400393               CASE (jpnf90   )   ;   CALL iom_nf90_close(    jf ) 
    401                CASE (jprstdimg)   ;   CALL iom_rstdimg_close( jf ) 
    402394               CASE DEFAULT 
    403                   CALL ctl_stop( TRIM(clinfo)//' accepted IO library are only jpioipsl, jpnf90 and jprstdimg' ) 
     395                  CALL ctl_stop( TRIM(clinfo)//' accepted IO library is only jpnf90 (jpioipsl option has been removed)' ) 
    404396               END SELECT 
    405397               iom_file(jf)%nfid       = 0          ! free the id  
     
    456448               IF( iiv <= jpmax_vars ) THEN 
    457449                  SELECT CASE (iom_file(kiomid)%iolib) 
    458                   CASE (jpioipsl )   ;   iom_varid = iom_ioipsl_varid( kiomid, cdvar, iiv, kdimsz ) 
    459450                  CASE (jpnf90   )   ;   iom_varid = iom_nf90_varid  ( kiomid, cdvar, iiv, kdimsz, kndims ) 
    460                   CASE (jprstdimg)   ;   iom_varid = -1   ! all variables are listed in iom_file 
    461                   CASE DEFAULT    
    462                      CALL ctl_stop( TRIM(clinfo)//' accepted IO library are only jpioipsl, jpnf90 and jprstdimg' ) 
     451                  CASE DEFAULT 
     452                     CALL ctl_stop( TRIM(clinfo)//' accepted IO library is only jpnf90 (jpioipsl option has been removed)' ) 
    463453                  END SELECT 
    464454               ELSE 
     
    518508                                 &                         'Use ncwa -a to suppress the unnecessary dimensions' ) 
    519509            SELECT CASE (iom_file(kiomid)%iolib) 
    520             CASE (jpioipsl )   ;   CALL iom_ioipsl_get(  kiomid, idvar, pvar, itime ) 
    521510            CASE (jpnf90   )   ;   CALL iom_nf90_get(    kiomid, idvar, pvar, itime ) 
    522             CASE (jprstdimg)   ;   CALL iom_rstdimg_get( kiomid, idvar, pvar ) 
    523             CASE DEFAULT     
    524                CALL ctl_stop( 'iom_g0d: accepted IO library are only jpioipsl, jpnf90 and jprstdimg' ) 
     511            CASE DEFAULT 
     512               CALL ctl_stop( 'iom_g0d: accepted IO library is only jpnf90 (jpioipsl option has been removed)' ) 
    525513            END SELECT 
    526514         ENDIF 
     
    785773       
    786774         SELECT CASE (iom_file(kiomid)%iolib) 
    787          CASE (jpioipsl )   ;   CALL iom_ioipsl_get(  kiomid, idvar, inbdim, istart, icnt, ix1, ix2, iy1, iy2,   & 
    788             &                                         pv_r1d, pv_r2d, pv_r3d ) 
    789775         CASE (jpnf90   )   ;   CALL iom_nf90_get(    kiomid, idvar, inbdim, istart, icnt, ix1, ix2, iy1, iy2,   & 
    790776            &                                         pv_r1d, pv_r2d, pv_r3d ) 
    791          CASE (jprstdimg)   ;   CALL iom_rstdimg_get( kiomid, idom, idvar, ix1, ix2, iy1, iy2,   & 
    792             &                                         pv_r1d, pv_r2d, pv_r3d ) 
    793          CASE DEFAULT     
    794             CALL ctl_stop( TRIM(clinfo)//' accepted IO library are only jpioipsl, jpnf90 and jprstdimg' ) 
     777         CASE DEFAULT 
     778            CALL ctl_stop( TRIM(clinfo)//' accepted IO library is only jpnf90 (jpioipsl option has been removed)' ) 
    795779         END SELECT 
    796780 
     
    879863                  IF( iom_file(kiomid)%dimsz(1,idvar) <= size(ptime) ) THEN 
    880864                     SELECT CASE (iom_file(kiomid)%iolib) 
    881                      CASE (jpioipsl )   ;   CALL iom_ioipsl_gettime( kiomid, idvar, ptime, cdunits, cdcalendar ) 
    882865                     CASE (jpnf90   )   ;   CALL iom_nf90_gettime(   kiomid, idvar, ptime, cdunits, cdcalendar ) 
    883                      CASE (jprstdimg)   ;   CALL ctl_stop( TRIM(clinfo)//' case IO library == jprstdimg not coded...' ) 
    884                      CASE DEFAULT     
    885                         CALL ctl_stop( TRIM(clinfo)//' accepted IO library are only jpioipsl, jpnf90 and jprstdimg' ) 
     866                     CASE DEFAULT 
     867                        CALL ctl_stop( TRIM(clinfo)//' accepted IO library is only jpnf90 (jpioipsl option has been removed)' ) 
    886868                     END SELECT 
    887869                  ELSE 
     
    914896         IF( iom_file(kiomid)%nfid > 0 ) THEN 
    915897            SELECT CASE (iom_file(kiomid)%iolib) 
    916             CASE (jpioipsl )   ;   CALL ctl_stop('iom_getatt: only nf90 available') 
    917898            CASE (jpnf90   )   ;   CALL iom_nf90_getatt( kiomid, cdatt, pvar ) 
    918             CASE (jprstdimg)   ;   CALL ctl_stop('iom_getatt: only nf90 available') 
    919             CASE DEFAULT     
    920                CALL ctl_stop( 'iom_g0d_att: accepted IO library are only jpioipsl, jpnf90 and jprstdimg' ) 
     899            CASE DEFAULT 
     900               CALL ctl_stop( 'iom_g0d_att: accepted IO library is only jpnf90' ) 
    921901            END SELECT 
    922902         ENDIF 
     
    940920            ivid = iom_varid( kiomid, cdvar, ldstop = .FALSE. ) 
    941921            SELECT CASE (iom_file(kiomid)%iolib) 
    942             CASE (jpioipsl )   ;   CALL iom_ioipsl_rstput( kt, kwrite, kiomid, cdvar, ivid, ktype, pv_r0d = pvar ) 
    943922            CASE (jpnf90   )   ;   CALL iom_nf90_rstput(   kt, kwrite, kiomid, cdvar, ivid, ktype, pv_r0d = pvar ) 
    944             CASE (jprstdimg)   ;   IF( kt == kwrite )    CALL iom_rstdimg_rstput( kiomid, cdvar, ivid, pvar ) 
    945             CASE DEFAULT      
    946                CALL ctl_stop( 'iom_rp0d: accepted IO library are only jpioipsl, jpnf90 and jprstdimg' ) 
     923            CASE DEFAULT 
     924               CALL ctl_stop( 'iom_rp0d: accepted IO library is only jpnf90 (jpioipsl option has been removed)' ) 
    947925            END SELECT 
    948926         ENDIF 
     
    962940            ivid = iom_varid( kiomid, cdvar, ldstop = .FALSE. ) 
    963941            SELECT CASE (iom_file(kiomid)%iolib) 
    964             CASE (jpioipsl )   ;   CALL iom_ioipsl_rstput( kt, kwrite, kiomid, cdvar, ivid, ktype, pv_r1d = pvar ) 
    965942            CASE (jpnf90   )   ;   CALL iom_nf90_rstput(   kt, kwrite, kiomid, cdvar, ivid, ktype, pv_r1d = pvar ) 
    966             CASE (jprstdimg)   ;   IF( kt == kwrite )    CALL iom_rstdimg_rstput( kiomid, cdvar, ivid, pv_r1d = pvar ) 
    967             CASE DEFAULT      
    968                CALL ctl_stop( 'iom_rp1d: accepted IO library are only jpioipsl, jpnf90 and jprstdimg' ) 
     943            CASE DEFAULT 
     944               CALL ctl_stop( 'iom_rp1d: accepted IO library is only jpnf90 (jpioipsl option has been removed)' ) 
    969945            END SELECT 
    970946         ENDIF 
     
    984960            ivid = iom_varid( kiomid, cdvar, ldstop = .FALSE. ) 
    985961            SELECT CASE (iom_file(kiomid)%iolib) 
    986             CASE (jpioipsl )   ;   CALL iom_ioipsl_rstput( kt, kwrite, kiomid, cdvar, ivid, ktype, pv_r2d = pvar ) 
    987962            CASE (jpnf90   )   ;   CALL iom_nf90_rstput(   kt, kwrite, kiomid, cdvar, ivid, ktype, pv_r2d = pvar ) 
    988             CASE (jprstdimg)   ;   IF( kt == kwrite )   CALL iom_rstdimg_rstput( kiomid, cdvar, ivid, pv_r2d = pvar )  
    989             CASE DEFAULT      
    990                CALL ctl_stop( 'iom_rp2d: accepted IO library are only jpioipsl, jpnf90 and jprstdimg' ) 
     963            CASE DEFAULT 
     964               CALL ctl_stop( 'iom_rp2d: accepted IO library is only jpnf90 (jpioipsl option has been removed)' ) 
    991965            END SELECT 
    992966         ENDIF 
     
    1006980            ivid = iom_varid( kiomid, cdvar, ldstop = .FALSE. ) 
    1007981            SELECT CASE (iom_file(kiomid)%iolib) 
    1008             CASE (jpioipsl )   ;   CALL iom_ioipsl_rstput( kt, kwrite, kiomid, cdvar, ivid, ktype, pv_r3d = pvar ) 
    1009982            CASE (jpnf90   )   ;   CALL iom_nf90_rstput(   kt, kwrite, kiomid, cdvar, ivid, ktype, pv_r3d = pvar ) 
    1010             CASE (jprstdimg)   ;   IF( kt == kwrite )   CALL iom_rstdimg_rstput( kiomid, cdvar, ivid, pv_r3d = pvar ) 
    1011             CASE DEFAULT      
    1012                CALL ctl_stop( 'iom_rp3d: accepted IO library are only jpioipsl and jprstdimg' ) 
     983            CASE DEFAULT 
     984               CALL ctl_stop( 'iom_rp3d: accepted IO library is only jpnf90 (jpioipsl option has been removed)' ) 
    1013985            END SELECT 
    1014986         ENDIF 
     
    14171389            idx = INDEX(clname,'@startdate@') + INDEX(clname,'@STARTDATE@') 
    14181390            DO WHILE ( idx /= 0 )  
    1419                cldate = iom_sdate( fjulday - rdttra(1) / rday ) 
     1391               cldate = iom_sdate( fjulday - rdt / rday ) 
    14201392               clname = clname(1:idx-1)//TRIM(cldate)//clname(idx+11:LEN_TRIM(clname)) 
    14211393               idx = INDEX(clname,'@startdate@') + INDEX(clname,'@STARTDATE@') 
     
    14241396            idx = INDEX(clname,'@startdatefull@') + INDEX(clname,'@STARTDATEFULL@') 
    14251397            DO WHILE ( idx /= 0 )  
    1426                cldate = iom_sdate( fjulday - rdttra(1) / rday, ldfull = .TRUE. ) 
     1398               cldate = iom_sdate( fjulday - rdt / rday, ldfull = .TRUE. ) 
    14271399               clname = clname(1:idx-1)//TRIM(cldate)//clname(idx+15:LEN_TRIM(clname)) 
    14281400               idx = INDEX(clname,'@startdatefull@') + INDEX(clname,'@STARTDATEFULL@') 
     
    14311403            idx = INDEX(clname,'@enddate@') + INDEX(clname,'@ENDDATE@') 
    14321404            DO WHILE ( idx /= 0 )  
    1433                cldate = iom_sdate( fjulday + rdttra(1) / rday * REAL( nitend - nit000, wp ), ld24 = .TRUE. ) 
     1405               cldate = iom_sdate( fjulday + rdt / rday * REAL( nitend - nit000, wp ), ld24 = .TRUE. ) 
    14341406               clname = clname(1:idx-1)//TRIM(cldate)//clname(idx+9:LEN_TRIM(clname)) 
    14351407               idx = INDEX(clname,'@enddate@') + INDEX(clname,'@ENDDATE@') 
     
    14381410            idx = INDEX(clname,'@enddatefull@') + INDEX(clname,'@ENDDATEFULL@') 
    14391411            DO WHILE ( idx /= 0 )  
    1440                cldate = iom_sdate( fjulday + rdttra(1) / rday * REAL( nitend - nit000, wp ), ld24 = .TRUE., ldfull = .TRUE. ) 
     1412               cldate = iom_sdate( fjulday + rdt / rday * REAL( nitend - nit000, wp ), ld24 = .TRUE., ldfull = .TRUE. ) 
    14411413               clname = clname(1:idx-1)//TRIM(cldate)//clname(idx+13:LEN_TRIM(clname)) 
    14421414               idx = INDEX(clname,'@enddatefull@') + INDEX(clname,'@ENDDATEFULL@') 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/IOM/iom_def.F90

    r4205 r5282  
    2828   INTEGER, PARAMETER, PUBLIC ::   jpdom_autodta       = 9   !:  
    2929 
    30    INTEGER, PARAMETER, PUBLIC ::   jpioipsl    = 100      !: Use ioipsl (fliocom only) library 
    3130   INTEGER, PARAMETER, PUBLIC ::   jpnf90      = 101      !: Use nf90 library 
    32    INTEGER, PARAMETER, PUBLIC ::   jprstdimg   = 102      !: Use restart dimgs (fortran direct acces) library 
    33 #if defined key_dimgout 
    34    INTEGER, PARAMETER, PUBLIC ::   jprstlib  = jprstdimg  !: restarts io library 
    35 #else 
     31 
    3632   INTEGER, PARAMETER, PUBLIC ::   jprstlib  = jpnf90     !: restarts io library 
    37 #endif 
    3833 
    3934   INTEGER, PARAMETER, PUBLIC ::   jp_r8    = 200      !: write REAL(8) 
     
    5449      CHARACTER(LEN=240)                        ::   name     !: name of the file 
    5550      INTEGER                                   ::   nfid     !: identifier of the file (0 if closed) 
    56       INTEGER                                   ::   iolib    !: library used to read the file (jpioipsl, jpnf90 or jprstdimg) 
     51      INTEGER                                   ::   iolib    !: library used to read the file (jpnf90 or new formats, 
     52                                                              !: jpioipsl option has been removed) 
    5753      INTEGER                                   ::   nvars    !: number of identified varibles in the file 
    5854      INTEGER                                   ::   iduld    !: id of the unlimited dimension 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90

    r4990 r5282  
    8282            WRITE(numout,*) 
    8383            SELECT CASE ( jprstlib ) 
    84             CASE ( jprstdimg )   ;   WRITE(numout,*) '             open ocean restart binary file: '//clname 
    8584            CASE DEFAULT         ;   WRITE(numout,*) '             open ocean restart NetCDF file: '//clname 
    8685            END SELECT 
     
    110109      !!---------------------------------------------------------------------- 
    111110 
    112                      CALL iom_rstput( kt, nitrst, numrow, 'rdt'    , rdt       )   ! dynamics time step 
    113                      CALL iom_rstput( kt, nitrst, numrow, 'rdttra1', rdttra(1) )   ! surface tracer time step 
     111                     CALL iom_rstput( kt, nitrst, numrow, 'rdt'    , rdt       )   ! dynamics and tracer time step 
    114112 
    115113                     CALL iom_rstput( kt, nitrst, numrow, 'ub'     , ub        )     ! before fields 
     
    165163            SELECT CASE ( jprstlib ) 
    166164            CASE ( jpnf90    )   ;   WRITE(numout,*) 'rst_read : read oce NetCDF restart file' 
    167             CASE ( jprstdimg )   ;   WRITE(numout,*) 'rst_read : read oce binary restart file' 
    168165            END SELECT 
    169166            IF ( snc4set%luse )      WRITE(numout,*) 'rst_read : configured with NetCDF4 support' 
     
    171168         ENDIF 
    172169 
    173          IF ( jprstlib == jprstdimg ) THEN 
    174            ! eventually read netcdf file (monobloc)  for restarting on different number of processors 
    175            ! if {cn_ocerst_in}.nc exists, then set jlibalt to jpnf90 
    176            INQUIRE( FILE = TRIM(cn_ocerst_in)//'.nc', EXIST = llok ) 
    177            IF ( llok ) THEN ; jlibalt = jpnf90  ; ELSE ; jlibalt = jprstlib ; ENDIF 
    178          ENDIF 
    179170         CALL iom_open( cn_ocerst_in, numror, kiolib = jlibalt ) 
    180171      ENDIF 
     
    189180      !! ** Method  :   Read in restart.nc file fields which are necessary for restart 
    190181      !!---------------------------------------------------------------------- 
    191       REAL(wp) ::   zrdt, zrdttra1 
     182      REAL(wp) ::   zrdt 
    192183      INTEGER  ::   jk 
    193184      LOGICAL  ::   llok 
     
    201192         IF( zrdt /= rdt )   neuler = 0 
    202193      ENDIF 
    203       IF( iom_varid( numror, 'rdttra1', ldstop = .FALSE. ) > 0 )   THEN 
    204          CALL iom_get( numror, 'rdttra1', zrdttra1 ) 
    205          IF( zrdttra1 /= rdttra(1) )   neuler = 0 
    206       ENDIF 
    207       !  
     194      ! 
    208195      IF( iom_varid( numror, 'ub', ldstop = .FALSE. ) > 0 ) THEN 
    209196         CALL iom_get( numror, jpdom_autoglo, 'ub'     , ub      )   ! before fields 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/LBC/lib_mpp.F90

    r4990 r5282  
    297297 
    298298      IF( mynode == 0 ) THEN 
    299         CALL ctl_opn( kumond, 'output.namelist.dyn', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. , 1 ) 
     299        CALL ctl_opn( kumond, 'output.namelist.dyn', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -2, 6, .TRUE. , 1 ) 
    300300        WRITE(kumond, nammpp)       
    301301      ENDIF 
     
    31983198      IF( PRESENT( localComm ) .OR. .NOT.PRESENT( localComm ) )   function_value = 0 
    31993199      IF( .FALSE. )   ldtxt(:) = 'never done' 
    3200       CALL ctl_opn( kumond, 'output.namelist.dyn', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. , 1 ) 
     3200      CALL ctl_opn( kumond, 'output.namelist.dyn', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, 6, .TRUE. , 1 ) 
    32013201   END FUNCTION mynode 
    32023202 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/LBC/mppini.F90

    r4679 r5282  
    457457#  include "mppini_2.h90" 
    458458 
    459 # if defined key_dimgout 
    460    !!---------------------------------------------------------------------- 
    461    !!   'key_dimgout'                  NO use of NetCDF files 
    462    !!---------------------------------------------------------------------- 
    463    SUBROUTINE mpp_init_ioipsl       ! Dummy routine 
    464    END SUBROUTINE mpp_init_ioipsl   
    465 # else 
    466459   SUBROUTINE mpp_init_ioipsl 
    467460      !!---------------------------------------------------------------------- 
     
    509502   END SUBROUTINE mpp_init_ioipsl   
    510503 
    511 # endif 
    512504#endif 
    513505 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90

    r4624 r5282  
    6969         &                 ln_dynldf_level, ln_dynldf_hor  , ln_dynldf_iso,   & 
    7070         &                 rn_ahm_0_lap   , rn_ahmb_0      , rn_ahm_0_blp ,   & 
    71          &                 rn_cmsmag_1    , rn_cmsmag_2    , rn_cmsh,         & 
    7271         &                 rn_ahm_m_lap   , rn_ahm_m_blp 
    7372 
     
    155154      IF(lwp) WRITE(numout,*) '        ahm1 = ahm2 = ahm0 =  ',ahm0 
    156155#endif 
    157      nkahm_smag = 0 
    158 #if defined key_dynldf_smag 
    159      nkahm_smag = 1 
    160 #endif 
    161156 
    162157      ! 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_oce.F90

    r4147 r5282  
    2323   REAL(wp), PUBLIC ::   rn_ahm_0_blp       !: lateral bilaplacian eddy viscosity (m4/s) 
    2424   REAL(wp), PUBLIC ::   ahm0, ahmb0, ahm0_blp         !: OLD namelist names 
    25    REAL(wp), PUBLIC ::   rn_cmsmag_1        !: constant in laplacian Smagorinsky viscosity 
    26    REAL(wp), PUBLIC ::   rn_cmsmag_2        !: constant in bilaplacian Smagorinsky viscosity 
    27    REAL(wp), PUBLIC ::   rn_cmsh            !: 1 or 0 , if 0 -use only shear for Smagorinsky viscosity 
    2825   REAL(wp), PUBLIC ::   rn_ahm_m_blp       !: upper limit for bilap  abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 
    2926   REAL(wp), PUBLIC ::   rn_ahm_m_lap       !: upper limit for lap  ahm < min(dx^2/16rdt, rn_ahm_m_lap) 
    3027 
    31    INTEGER , PUBLIC ::   nkahm_smag      =  0          !:  
    3228 
    3329   !                                                                                  !!! eddy coeff. at U-,V-,W-pts [m2/s] 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90

    r4624 r5282  
    7070         &                 ln_triad_iso   , ln_botmix_grif ,                  & 
    7171         &                 rn_aht_0       , rn_ahtb_0      , rn_aeiv_0,       & 
    72          &                 rn_slpmax      , rn_chsmag      ,    rn_smsh,      & 
    73          &                 rn_aht_m 
     72         &                 rn_slpmax      , rn_aht_m 
    7473      !!---------------------------------------------------------------------- 
    7574 
     
    167166#endif 
    168167 
    169 #if defined key_traldf_smag && ! defined key_traldf_c3d 
    170         CALL ctl_stop( 'key_traldf_smag can only be used with key_traldf_c3d' ) 
    171 #endif 
    172 #if defined key_traldf_smag 
    173         IF(lwp) WRITE(numout,*)' SMAGORINSKY DIFFUSION' 
    174         IF(lwp .AND. rn_smsh < 1)  WRITE(numout,*)' only  shear is used ' 
    175         IF(lwp.and.ln_traldf_bilap) CALL ctl_stop(' SMAGORINSKY + BILAPLACIAN - UNSTABLE OR NON_CONSERVATIVE' ) 
    176 #endif 
    177  
    178168      ! 
    179169   END SUBROUTINE ldf_tra_init 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_oce.F90

    r4147 r5282  
    3030   REAL(wp), PUBLIC ::   rn_aeiv_0        !: eddy induced velocity coefficient (m2/s) 
    3131   REAL(wp), PUBLIC ::   rn_slpmax        !: slope limit 
    32    REAL(wp), PUBLIC ::   rn_chsmag        !:  multiplicative factor in Smagorinsky diffusivity 
    33    REAL(wp), PUBLIC ::   rn_smsh          !:  Smagorinsky diffusivity: = 0 - use only sheer 
    3432   REAL(wp), PUBLIC ::   rn_aht_m         !:  upper limit or stability criteria for lateral eddy diffusivity (m2/s) 
    3533 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/SBC/fldread.F90

    r4784 r5282  
    162162         isecsbc = nsec_year + nsec1jan000 + (kit+it_offset)*NINT( rdt/REAL(nn_baro,wp) ) 
    163163      ELSE                      ! middle of sbc time step 
    164          isecsbc = nsec_year + nsec1jan000 + NINT(0.5 * REAL(kn_fsbc - 1,wp) * rdttra(1)) + it_offset * NINT(rdttra(1)) 
     164        isecsbc = nsec_year + nsec1jan000 + NINT(0.5 * REAL(kn_fsbc - 1,wp) * rdt) + it_offset * NINT(rdt) 
    165165      ENDIF 
    166166      imf = SIZE( sd ) 
     
    189189               CALL fld_rec( kn_fsbc, sd(jf), kt_offset = it_offset, kit = kit )    ! update after record informations 
    190190 
    191                ! if kn_fsbc*rdttra is larger than nfreqh (which is kind of odd), 
     191               ! if kn_fsbc*rdt is larger than nfreqh (which is kind of odd), 
    192192               ! it is possible that the before value is no more the good one... we have to re-read it 
    193193               ! if before is not the last record of the file currently opened and after is the first record to be read 
     
    210210               IF( sd(jf)%ln_tint ) THEN 
    211211                   
    212                   ! if kn_fsbc*rdttra is larger than nfreqh (which is kind of odd), 
     212                  ! if kn_fsbc*rdt is larger than nfreqh (which is kind of odd), 
    213213                  ! it is possible that the before value is no more the good one... we have to re-read it 
    214214                  ! if before record is not just just before the after record... 
     
    241241                        ! year/month/week/day file to be not present. If the run continue further than the current 
    242242                        ! year/month/week/day, next year/month/week/day file must exist 
    243                         isecend = nsec_year + nsec1jan000 + (nitend - kt) * NINT(rdttra(1))   ! second at the end of the run  
     243                        isecend = nsec_year + nsec1jan000 + (nitend - kt) * NINT(rdt)   ! second at the end of the run 
    244244                        llstop = isecend > sd(jf)%nrec_a(2)                                   ! read more than 1 record of next year 
    245245                        ! we suppose that the date of next file is next day (should be ok even for weekly files...) 
     
    454454      IF( PRESENT(kt_offset) )   it_offset = kt_offset 
    455455      IF( PRESENT(kit) ) THEN   ;   it_offset = ( kit + it_offset ) * NINT( rdt/REAL(nn_baro,wp) ) 
    456       ELSE                      ;   it_offset =         it_offset   * NINT(       rdttra(1)      ) 
     456      ELSE                      ;   it_offset =         it_offset   * NINT(       rdt      ) 
    457457      ENDIF 
    458458      ! 
     
    531531         ELSE                                           ;   ztmp = REAL(nsec_year ,wp)  ! since 00h on Jan 1 of the current year 
    532532         ENDIF 
    533          ztmp = ztmp + 0.5 * REAL(kn_fsbc - 1, wp) * rdttra(1) + REAL( it_offset, wp )  ! centrered in the middle of sbc time step 
    534          ztmp = ztmp + 0.01 * rdttra(1)                                                 ! avoid truncation error  
     533         ztmp = ztmp + 0.5 * REAL(kn_fsbc - 1, wp) * rdt + REAL( it_offset, wp )  ! centrered in the middle of sbc time step 
     534         ztmp = ztmp + 0.01 * rdt                                                 ! avoid truncation error 
    535535         IF( sdjf%ln_tint ) THEN                ! time interpolation, shift by 1/2 record 
    536536            ! 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

    r5009 r5282  
    645645      LOGICAL ::    llnewtx, llnewtau      ! update wind stress components and module?? 
    646646      INTEGER  ::   ji, jj, jn             ! dummy loop indices 
    647       INTEGER  ::   isec                   ! number of seconds since nit000 (assuming rdttra did not change since nit000) 
     647      INTEGER  ::   isec                   ! number of seconds since nit000 (assuming rdt did not change since nit000) 
    648648      REAL(wp) ::   zcumulneg, zcumulpos   ! temporary scalars      
    649649      REAL(wp) ::   zcoef                  ! temporary scalar 
     
    658658      CALL wrk_alloc( jpi,jpj, ztx, zty ) 
    659659      !                                                 ! Receive all the atmos. fields (including ice information) 
    660       isec = ( kt - nit000 ) * NINT( rdttra(1) )             ! date of exchanges 
     660      isec = ( kt - nit000 ) * NINT( rdt )             ! date of exchanges 
    661661      DO jn = 1, jprcv                                       ! received fields sent by the atmosphere 
    662662         IF( srcv(jn)%laction )   CALL cpl_rcv( jn, isec, frcv(jn)%z3, xcplmask, nrcvinfo(jn) ) 
     
    13381338      CALL wrk_alloc( jpi,jpj,jpl, ztmp3, ztmp4 ) 
    13391339 
    1340       isec = ( kt - nit000 ) * NINT(rdttra(1))        ! date of exchanges 
     1340      isec = ( kt - nit000 ) * NINT(rdt)        ! date of exchanges 
    13411341 
    13421342      zfr_l(:,:) = 1.- fr_i(:,:) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/SBC/sbcdcy.F90

    r3764 r5282  
    9090 
    9191      ! When are we during the day (from 0 to 1) 
    92       zlo = ( REAL(nsec_day, wp) - 0.5_wp * rdttra(1) ) / rday 
    93       zup = zlo + ( REAL(nn_fsbc, wp)     * rdttra(1) ) / rday 
     92      zlo = ( REAL(nsec_day, wp) - 0.5_wp * rdt ) / rday 
     93      zup = zlo + ( REAL(nn_fsbc, wp)     * rdt ) / rday 
    9494      !                                           
    9595      IF( nday_qsr == -1 ) THEN       ! first time step only   
     
    189189         END DO   
    190190         ! 
    191          ztmp = rday / ( rdttra(1) * REAL(nn_fsbc, wp) ) 
     191         ztmp = rday / ( rdt * REAL(nn_fsbc, wp) ) 
    192192         rscal(:,:) = rscal(:,:) * ztmp 
    193193         ! 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/SBC/sbcfwb.F90

    r4990 r5282  
    125125         ENDIF    
    126126         !                                         ! Update fwfold if new year start 
    127          ikty = 365 * 86400 / rdttra(1)    !!bug  use of 365 days leap year or 360d year !!!!!!! 
     127         ikty = 365 * 86400 / rdt       !!bug  use of 365 days leap year or 360d year !!!!!!! 
    128128         IF( MOD( kt, ikty ) == 0 ) THEN 
    129129            a_fwb_b = a_fwb                           ! mean sea level taking into account the ice+snow 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

    r4990 r5282  
    255255            CALL dom_vvl_interpol( fse3t_b(:,:,:), fse3u_b(:,:,:), 'U' ) 
    256256            CALL dom_vvl_interpol( fse3t_b(:,:,:), fse3v_b(:,:,:), 'V' ) 
    257             CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3u_n(:,:,:), 'U' ) 
     257            CALL dom_vvl_interpol( fse333333333333333333333333333333333t_n(:,:,:), fse3u_n(:,:,:), 'U' ) 
    258258            CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3v_n(:,:,:), 'V' ) 
    259259            CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3f_n(:,:,:), 'F' ) 
    260260            ! Vertical scale factor interpolations 
    261261            ! ------------------------------------ 
    262             CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3w_n (:,:,:), 'W'  ) 
     262            CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3w_n(:,:,:), 'W'  ) 
    263263            CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3uw_n(:,:,:), 'UW' ) 
    264264            CALL dom_vvl_interpol( fse3v_n(:,:,:), fse3vw_n(:,:,:), 'VW' ) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/SBC/sbcrnf.F90

    r4990 r5282  
    461461      CALL iom_close( inum )                                      ! close file 
    462462      ! 
    463       IF( nn_closea == 1 )   CALL clo_rnf( rnfmsk )               ! closed sea inflow set as ruver mouth 
     463      IF( nn_closea == 1 )   CALL clo_rnf( rnfmsk )               ! closed sea inflow set as river mouth 
    464464      ! 
    465465      rnfmsk_z(:)   = 0._wp                                       ! vertical structure 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/SBC/sbctide.F90

    r4292 r5282  
    4949      !!---------------------------------------------------------------------- 
    5050 
    51       IF( nsec_day == NINT(0.5_wp * rdttra(1)) ) THEN      ! start a new day 
     51      IF( nsec_day == NINT(0.5_wp * rdt) ) THEN      ! start a new day 
    5252         ! 
    5353         IF( kt == nit000 ) THEN 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90

    r4990 r5282  
    2424   USE traadv_eiv      ! eddy induced velocity     (tra_adv_eiv    routine) 
    2525   USE traadv_mle      ! ML eddy induced velocity  (tra_adv_mle    routine) 
    26    USE cla             ! cross land advection      (cla_traadv     routine) 
    2726   USE ldftra_oce      ! lateral diffusion coefficient on tracers 
    2827   USE in_out_manager  ! I/O manager 
     
    8483      !                                          ! set time step 
    8584      IF( neuler == 0 .AND. kt == nit000 ) THEN     ! at nit000 
    86          r2dtra(:) =  rdttra(:)                          ! = rdtra (restarting with Euler time stepping) 
     85         r2dt =  rdt                          ! = rdt (restarting with Euler time stepping) 
    8786      ELSEIF( kt <= nit000 + 1) THEN                ! at nit000 or nit000+1 
    88          r2dtra(:) = 2._wp * rdttra(:)                   ! = 2 rdttra (leapfrog) 
    89       ENDIF 
    90       ! 
    91       IF( nn_cla == 1 .AND. cp_cfg == 'orca' .AND. jp_cfg == 2 )   CALL cla_traadv( kt )       !==  Cross Land Advection  ==! (hor. advection) 
     87         r2dt = 2._wp * rdt                   ! = 2 rdt (leapfrog) 
     88      ENDIF 
    9289      ! 
    9390      !                                               !==  effective transport  ==! 
     
    117114      SELECT CASE ( nadv )                            !==  compute advection trend and add it to general trend  ==! 
    118115      CASE ( 1 )   ;    CALL tra_adv_cen2  ( kt, nit000, 'TRA',         zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  2nd order centered 
    119       CASE ( 2 )   ;    CALL tra_adv_tvd   ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  TVD  
    120       CASE ( 3 )   ;    CALL tra_adv_muscl ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb,      tsa, jpts, ln_traadv_msc_ups )   !  MUSCL  
    121       CASE ( 4 )   ;    CALL tra_adv_muscl2( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  MUSCL2  
    122       CASE ( 5 )   ;    CALL tra_adv_ubs   ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  UBS  
    123       CASE ( 6 )   ;    CALL tra_adv_qck   ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  QUICKEST  
    124       CASE ( 7 )   ;   CALL tra_adv_tvd_zts( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  TVD ZTS 
     116      CASE ( 2 )   ;    CALL tra_adv_tvd   ( kt, nit000, 'TRA', r2dt, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  TVD 
     117      CASE ( 3 )   ;    CALL tra_adv_muscl ( kt, nit000, 'TRA', r2dt, zun, zvn, zwn, tsb,      tsa, jpts, ln_traadv_msc_ups )   !  MUSCL 
     118      CASE ( 4 )   ;    CALL tra_adv_muscl2( kt, nit000, 'TRA', r2dt, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  MUSCL2 
     119      CASE ( 5 )   ;    CALL tra_adv_ubs   ( kt, nit000, 'TRA', r2dt, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  UBS 
     120      CASE ( 6 )   ;    CALL tra_adv_qck   ( kt, nit000, 'TRA', r2dt, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  QUICKEST 
     121      CASE ( 7 )   ;   CALL tra_adv_tvd_zts( kt, nit000, 'TRA', r2dt, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  TVD ZTS 
    125122      ! 
    126123      CASE (-1 )                                      !==  esopa: test all possibility with control print  ==! 
     
    128125         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' adv0 - Ta: ', mask1=tmask,               & 
    129126            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    130          CALL tra_adv_tvd   ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )           
     127         CALL tra_adv_tvd   ( kt, nit000, 'TRA', r2dt, zun, zvn, zwn, tsb, tsn, tsa, jpts ) 
    131128         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' adv1 - Ta: ', mask1=tmask,               & 
    132129            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    133          CALL tra_adv_muscl ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb,      tsa, jpts, ln_traadv_msc_ups )           
     130         CALL tra_adv_muscl ( kt, nit000, 'TRA', r2dt, zun, zvn, zwn, tsb,      tsa, jpts, ln_traadv_msc_ups ) 
    134131         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' adv3 - Ta: ', mask1=tmask,               & 
    135132            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    136          CALL tra_adv_muscl2( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )           
     133         CALL tra_adv_muscl2( kt, nit000, 'TRA', r2dt, zun, zvn, zwn, tsb, tsn, tsa, jpts ) 
    137134         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' adv4 - Ta: ', mask1=tmask,               & 
    138135            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    139          CALL tra_adv_ubs   ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )           
     136         CALL tra_adv_ubs   ( kt, nit000, 'TRA', r2dt, zun, zvn, zwn, tsb, tsn, tsa, jpts ) 
    140137         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' adv5 - Ta: ', mask1=tmask,               & 
    141138            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    142          CALL tra_adv_qck   ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )           
     139         CALL tra_adv_qck   ( kt, nit000, 'TRA', r2dt, zun, zvn, zwn, tsb, tsn, tsa, jpts ) 
    143140         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' adv6 - Ta: ', mask1=tmask,               & 
    144141            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl.F90

    r4990 r5282  
    7272      INTEGER                              , INTENT(in   ) ::   kjpt            ! number of tracers 
    7373      LOGICAL                              , INTENT(in   ) ::   ld_msc_ups      ! use upstream scheme within muscl 
    74       REAL(wp), DIMENSION(        jpk     ), INTENT(in   ) ::   p2dt            ! vertical profile of tracer time-step 
     74      REAL(wp)                             , INTENT(in   ) ::   p2dt            ! tracer time-step 
    7575      REAL(wp), DIMENSION(jpi,jpj,jpk     ), INTENT(in   ) ::   pun, pvn, pwn   ! 3 ocean velocity components 
    7676      REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(in   ) ::   ptb             ! before tracer field 
     
    176176         !                                             !-- MUSCL horizontal advective fluxes 
    177177         DO jk = 1, jpkm1                                     ! interior values 
    178             zdt  = p2dt(jk) 
     178            zdt  = p2dt 
    179179            DO jj = 2, jpjm1 
    180180               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    259259         ! 
    260260         DO jk = 1, jpkm1                                     ! interior values 
    261             zdt  = p2dt(jk) 
     261            zdt  = p2dt 
    262262            DO jj = 2, jpjm1       
    263263               DO ji = fs_2, fs_jpim1   ! vector opt. 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl2.F90

    r4990 r5282  
    6363      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype          ! =TRA or TRC (tracer indicator) 
    6464      INTEGER                              , INTENT(in   ) ::   kjpt            ! number of tracers 
    65       REAL(wp), DIMENSION(        jpk     ), INTENT(in   ) ::   p2dt            ! vertical profile of tracer time-step 
     65      REAL(wp)                             , INTENT(in   ) ::   p2dt            ! tracer time-step 
    6666      REAL(wp), DIMENSION(jpi,jpj,jpk     ), INTENT(in   ) ::   pun, pvn, pwn   ! 3 ocean velocity components 
    6767      REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(in   ) ::   ptb, ptn        ! before & now tracer fields 
     
    133133        !                                             !-- MUSCL horizontal advective fluxes 
    134134         DO jk = 1, jpkm1                                     ! interior values 
    135             zdt  = p2dt(jk) 
     135            zdt  = p2dt 
    136136            DO jj = 2, jpjm1 
    137137               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    240240         ! 
    241241         DO jk = 1, jpkm1                                     ! interior values 
    242             zdt  = p2dt(jk) 
     242            zdt  = p2dt 
    243243            DO jj = 2, jpjm1 
    244244               DO ji = fs_2, fs_jpim1   ! vector opt. 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_qck.F90

    r4990 r5282  
    8888      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype          ! =TRA or TRC (tracer indicator) 
    8989      INTEGER                              , INTENT(in   ) ::   kjpt            ! number of tracers 
    90       REAL(wp), DIMENSION(        jpk     ), INTENT(in   ) ::   p2dt            ! vertical profile of tracer time-step 
     90      REAL(wp)                             , INTENT(in   ) ::   p2dt            ! tracer time-step 
    9191      REAL(wp), DIMENSION(jpi,jpj,jpk     ), INTENT(in   ) ::   pun, pvn, pwn   ! 3 ocean velocity components 
    9292      REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(in   ) ::   ptb, ptn        ! before and now tracer fields 
     
    125125      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype     ! =TRA or TRC (tracer indicator) 
    126126      INTEGER                              , INTENT(in   ) ::   kjpt       ! number of tracers 
    127       REAL(wp), DIMENSION(        jpk     ), INTENT(in   ) ::   p2dt       ! vertical profile of tracer time-step 
     127      REAL(wp)                             , INTENT(in   ) ::   p2dt       ! tracer time-step 
    128128      REAL(wp), DIMENSION(jpi,jpj,jpk     ), INTENT(in   ) ::   pun        ! i-velocity components 
    129129      REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(in   ) ::   ptb, ptn   ! before and now tracer fields 
     
    170170         ! 
    171171         DO jk = 1, jpkm1   
    172             zdt =  p2dt(jk) 
     172            zdt =  p2dt 
    173173            DO jj = 2, jpjm1 
    174174               DO ji = fs_2, fs_jpim1   ! vector opt.    
     
    246246      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype     ! =TRA or TRC (tracer indicator) 
    247247      INTEGER                              , INTENT(in   ) ::   kjpt       ! number of tracers 
    248       REAL(wp), DIMENSION(        jpk     ), INTENT(in   ) ::   p2dt       ! vertical profile of tracer time-step 
     248      REAL(wp)                             , INTENT(in   ) ::   p2dt       ! tracer time-step 
    249249      REAL(wp), DIMENSION(jpi,jpj,jpk     ), INTENT(in   ) ::   pvn        ! j-velocity components 
    250250      REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(in   ) ::   ptb, ptn   ! before and now tracer fields 
     
    293293         ! 
    294294         DO jk = 1, jpkm1   
    295             zdt =  p2dt(jk) 
     295            zdt =  p2dt 
    296296            DO jj = 2, jpjm1 
    297297               DO ji = fs_2, fs_jpim1   ! vector opt.    
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_tvd.F90

    r4990 r5282  
    7474      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype          ! =TRA or TRC (tracer indicator) 
    7575      INTEGER                              , INTENT(in   ) ::   kjpt            ! number of tracers 
    76       REAL(wp), DIMENSION(        jpk     ), INTENT(in   ) ::   p2dt            ! vertical profile of tracer time-step 
     76      REAL(wp)                             , INTENT(in   ) ::   p2dt            ! tracer time-step 
    7777      REAL(wp), DIMENSION(jpi,jpj,jpk     ), INTENT(in   ) ::   pun, pvn, pwn   ! 3 ocean velocity components 
    7878      REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(in   ) ::   ptb, ptn        ! before and now tracer fields 
     
    161161         ! total advective trend 
    162162         DO jk = 1, jpkm1 
    163             z2dtt = p2dt(jk) 
     163            z2dtt = p2dt 
    164164            DO jj = 2, jpjm1 
    165165               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    287287      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype          ! =TRA or TRC (tracer indicator) 
    288288      INTEGER                              , INTENT(in   ) ::   kjpt            ! number of tracers 
    289       REAL(wp), DIMENSION(        jpk     ), INTENT(in   ) ::   p2dt            ! vertical profile of tracer time-step 
     289      REAL(wp)                             , INTENT(in   ) ::   p2dt            ! tracer time-step 
    290290      REAL(wp), DIMENSION(jpi,jpj,jpk     ), INTENT(in   ) ::   pun, pvn, pwn   ! 3 ocean velocity components 
    291291      REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(in   ) ::   ptb, ptn        ! before and now tracer fields 
    292292      REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(inout) ::   pta             ! tracer trend  
    293293      ! 
    294       REAL(wp), DIMENSION( jpk )                           ::   zts             ! length of sub-timestep for vertical advection 
    295       REAL(wp), DIMENSION( jpk )                           ::   zr_p2dt         ! reciprocal of tracer timestep 
     294      REAL(wp)                                             ::   zts             ! length of sub-timestep for vertical advection 
     295      REAL(wp)                                             ::   zr_p2dt         ! reciprocal of tracer timestep 
    296296      INTEGER  ::   ji, jj, jk, jl, jn       ! dummy loop indices   
    297297      INTEGER  ::   jnzts = 5       ! number of sub-timesteps for vertical advection 
     
    330330      zwi(:,:,:) = 0._wp 
    331331      z_rzts = 1._wp / REAL( jnzts, wp ) 
    332       zr_p2dt(:) = 1._wp / p2dt(:) 
     332      zr_p2dt = 1._wp / p2dt 
    333333      ! 
    334334      !                                                          ! =========== 
     
    375375         ! total advective trend 
    376376         DO jk = 1, jpkm1 
    377             z2dtt = p2dt(jk) 
     377            z2dtt = p2dt 
    378378            DO jj = 2, jpjm1 
    379379               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    446446            IF( jl == 1 ) THEN              ! Euler forward to kick things off 
    447447              jtb = 1   ;   jtn = 1   ;   jta = 2 
    448               zts(:) = p2dt(:) * z_rzts 
     448              zts = p2dt * z_rzts 
    449449              jtaken = MOD( jnzts + 1 , 2)  ! Toggle to collect every second flux 
    450450                                            ! starting at jl =1 if jnzts is odd;  
     
    452452            ELSEIF( jl == 2 ) THEN          ! First leapfrog step 
    453453              jtb = 1   ;   jtn = 2   ;   jta = 3 
    454               zts(:) = 2._wp * p2dt(:) * z_rzts 
     454              zts = 2._wp * p2dt * z_rzts 
    455455            ELSE                            ! Shuffle pointers for subsequent leapfrog steps 
    456456              jtb = MOD(jtb,3) + 1 
     
    462462                  DO ji = fs_2, fs_jpim1 
    463463                     zwz(ji,jj,jk) = 0.5_wp * pwn(ji,jj,jk) * ( ztrs(ji,jj,jk,jtn) + ztrs(ji,jj,jk-1,jtn) ) 
    464                      IF( jtaken == 0 ) zwzts(ji,jj,jk) = zwzts(ji,jj,jk) + zwz(ji,jj,jk)*zts(jk)           ! Accumulate time-weighted vertcal flux 
     464                     IF( jtaken == 0 ) zwzts(ji,jj,jk) = zwzts(ji,jj,jk) + zwz(ji,jj,jk)*zts           ! Accumulate time-weighted vertcal flux 
    465465                  END DO 
    466466               END DO 
     
    475475                     ! total advective trends 
    476476                     ztra = - zbtr * (  zhdiv(ji,jj,jk) + zwz(ji,jj,jk) - zwz(ji  ,jj  ,jk+1) ) 
    477                      ztrs(ji,jj,jk,jta) = ztrs(ji,jj,jk,jtb) + zts(jk) * ztra 
     477                     ztrs(ji,jj,jk,jta) = ztrs(ji,jj,jk,jtb) + zts * ztra 
    478478                  END DO 
    479479               END DO 
     
    485485            DO jj = 2, jpjm1 
    486486               DO ji = fs_2, fs_jpim1 
    487                   zwz(ji,jj,jk) = zwzts(ji,jj,jk) * zr_p2dt(jk) - zwz_sav(ji,jj,jk) 
     487                  zwz(ji,jj,jk) = zwzts(ji,jj,jk) * zr_p2dt - zwz_sav(ji,jj,jk) 
    488488               END DO 
    489489            END DO 
     
    553553      !!       in-space based differencing for fluid 
    554554      !!---------------------------------------------------------------------- 
    555       REAL(wp), DIMENSION(jpk)         , INTENT(in   ) ::   p2dt            ! vertical profile of tracer time-step 
     555      REAL(wp)                         , INTENT(in   ) ::   p2dt            ! tracer time-step 
    556556      REAL(wp), DIMENSION (jpi,jpj,jpk), INTENT(in   ) ::   pbef, paft      ! before & after field 
    557557      REAL(wp), DIMENSION (jpi,jpj,jpk), INTENT(inout) ::   paa, pbb, pcc   ! monotonic fluxes in the 3 directions 
     
    584584            DO jk = mikt(ji,jj), jpkm1 
    585585               ikm1 = MAX(jk-1,mikt(ji,jj)) 
    586                z2dtt = p2dt(jk) 
     586               z2dtt = p2dt 
    587587                
    588588               ! search maximum in neighbourhood 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_ubs.F90

    r4990 r5282  
    8181      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype          ! =TRA or TRC (tracer indicator) 
    8282      INTEGER                              , INTENT(in   ) ::   kjpt            ! number of tracers 
    83       REAL(wp), DIMENSION(        jpk     ), INTENT(in   ) ::   p2dt            ! vertical profile of tracer time-step 
     83      REAL(wp)                             , INTENT(in   ) ::   p2dt            ! tracer time-step 
    8484      REAL(wp), DIMENSION(jpi,jpj,jpk     ), INTENT(in   ) ::   pun, pvn, pwn   ! 3 ocean transport components 
    8585      REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(in   ) ::   ptb, ptn        ! before and now tracer fields 
     
    207207         ! update and guess with monotonic sheme 
    208208         DO jk = 1, jpkm1 
    209             z2dtt = p2dt(jk) 
     209            z2dtt = p2dt 
    210210            DO jj = 2, jpjm1 
    211211               DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    281281      !!       in-space based differencing for fluid 
    282282      !!---------------------------------------------------------------------- 
    283       REAL(wp), INTENT(in   ), DIMENSION(jpk)          ::   p2dt   ! vertical profile of tracer time-step 
     283      REAL(wp), INTENT(in   )                          ::   p2dt   ! tracer time-step 
    284284      REAL(wp),                DIMENSION (jpi,jpj,jpk) ::   pbef   ! before field 
    285285      REAL(wp), INTENT(inout), DIMENSION (jpi,jpj,jpk) ::   paft   ! after field 
     
    340340 
    341341      DO jk = 1, jpkm1 
    342          z2dtt = p2dt(jk) 
     342         z2dtt = p2dt 
    343343         DO jj = 2, jpjm1 
    344344            DO ji = fs_2, fs_jpim1   ! vector opt. 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/trabbl.F90

    r4990 r5282  
    362362         DO ji = 1, jpi 
    363363            ik = mbkt(ji,jj)                             ! bottom T-level index 
    364             zts (ji,jj,jp_tem) = tsb(ji,jj,ik,jp_tem)    ! bottom before T and S 
    365             zts (ji,jj,jp_sal) = tsb(ji,jj,ik,jp_sal) 
     364            zts(ji,jj,jp_tem) = tsb(ji,jj,ik,jp_tem)    ! bottom before T and S 
     365            zts(ji,jj,jp_sal) = tsb(ji,jj,ik,jp_sal) 
    366366            ! 
    367367            zdep(ji,jj) = fsdept(ji,jj,ik)               ! bottom T-level reference depth 
    368             zub (ji,jj) = un(ji,jj,mbku(ji,jj))          ! bottom velocity 
    369             zvb (ji,jj) = vn(ji,jj,mbkv(ji,jj)) 
     368            zub(ji,jj) = un(ji,jj,mbku(ji,jj))          ! bottom velocity 
     369            zvb(ji,jj) = vn(ji,jj,mbkv(ji,jj)) 
    370370         END DO 
    371371      END DO 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/traldf.F90

    r4990 r5282  
    326326      ELSE 
    327327         DO jk = 1, jpkm1 
    328             t0_ldf(:,:,jk) = ( tsa(:,:,jk,jp_tem) - tsb(:,:,jk,jp_tem) ) / ( z12 *rdttra(jk) ) 
    329             s0_ldf(:,:,jk) = ( tsa(:,:,jk,jp_sal) - tsb(:,:,jk,jp_sal) ) / ( z12 *rdttra(jk) ) 
     328            t0_ldf(:,:,jk) = ( tsa(:,:,jk,jp_tem) - tsb(:,:,jk,jp_tem) ) / ( z12 *rdt ) 
     329            s0_ldf(:,:,jk) = ( tsa(:,:,jk,jp_sal) - tsb(:,:,jk,jp_sal) ) / ( z12 *rdt ) 
    330330         END DO 
    331331      ENDIF 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

    r4990 r5282  
    120120  
    121121      ! set time step size (Euler/Leapfrog) 
    122       IF( neuler == 0 .AND. kt == nit000 ) THEN   ;   r2dtra(:) =     rdttra(:)      ! at nit000             (Euler) 
    123       ELSEIF( kt <= nit000 + 1 )           THEN   ;   r2dtra(:) = 2._wp* rdttra(:)      ! at nit000 or nit000+1 (Leapfrog) 
     122IF( neuler == 0 .AND. kt == nit000 ) THEN   ;   r2dt =     rdt      ! at nit000             (Euler) 
     123ELSEIF( kt <= nit000 + 1 )           THEN   ;   r2dt = 2._wp* rdt      ! at nit000 or nit000+1 (Leapfrog) 
    124124      ENDIF 
    125125 
     
    155155      ! trends computation 
    156156      IF( l_trdtra ) THEN      ! trend of the Asselin filter (tb filtered - tb)/dt      
    157          DO jk = 1, jpkm1 
    158             zfact = 1._wp / r2dtra(jk)              
     157      zfact = 1._wp / r2dt 
     158        DO jk = 1, jpkm1 
     159!            zfact = 1._wp / r2dt 
    159160            ztrdt(:,:,jk) = ( tsb(:,:,jk,jp_tem) - ztrdt(:,:,jk) ) * zfact 
    160161            ztrds(:,:,jk) = ( tsb(:,:,jk,jp_sal) - ztrds(:,:,jk) ) * zfact 
     
    297298      DO jn = 1, kjpt       
    298299         DO jk = 1, jpkm1 
    299             zfact1 = atfp * rdttra(jk) 
     300            zfact1 = atfp * rdt 
    300301            zfact2 = zfact1 / rau0 
    301302            DO jj = 1, jpj 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90

    r4990 r5282  
    2828   USE sbcrnf          ! River runoff   
    2929   USE sbcisf          ! Ice shelf    
    30    USE sbcmod          ! ln_rnf   
    3130   USE iom 
    3231   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf.F90

    r4990 r5282  
    6666      ! 
    6767      IF( neuler == 0 .AND. kt == nit000 ) THEN     ! at nit000 
    68          r2dtra(:) =  rdttra(:)                          ! = rdtra (restarting with Euler time stepping) 
     68         r2dt =  rdt                          ! = rdt (restarting with Euler time stepping) 
    6969      ELSEIF( kt <= nit000 + 1) THEN                ! at nit000 or nit000+1 
    70          r2dtra(:) = 2. * rdttra(:)                      ! = 2 rdttra (leapfrog) 
     70         r2dt = 2. * rdt                      ! = 2 rdt (leapfrog) 
    7171      ENDIF 
    7272 
     
    7878 
    7979      SELECT CASE ( nzdf )                       ! compute lateral mixing trend and add it to the general trend 
    80       CASE ( 0 )    ;    CALL tra_zdf_exp( kt, nit000, 'TRA', r2dtra, nn_zdfexp, tsb, tsa, jpts )  !   explicit scheme  
    81       CASE ( 1 )    ;    CALL tra_zdf_imp( kt, nit000, 'TRA', r2dtra,            tsb, tsa, jpts )  !   implicit scheme  
     80      CASE ( 0 )    ;    CALL tra_zdf_exp( kt, nit000, 'TRA', r2dt, nn_zdfexp, tsb, tsa, jpts )  !   explicit scheme 
     81      CASE ( 1 )    ;    CALL tra_zdf_imp( kt, nit000, 'TRA', r2dt,            tsb, tsa, jpts )  !   implicit scheme 
    8282      CASE ( -1 )                                       ! esopa: test all possibility with control print 
    83          CALL tra_zdf_exp( kt, nit000, 'TRA', r2dtra, nn_zdfexp, tsb, tsa, jpts ) 
     83         CALL tra_zdf_exp( kt, nit000, 'TRA', r2dt, nn_zdfexp, tsb, tsa, jpts ) 
    8484         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' zdf0 - Ta: ', mask1=tmask,               & 
    8585         &             tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    86          CALL tra_zdf_imp( kt, nit000, 'TRA', r2dtra,            tsb, tsa, jpts )  
     86         CALL tra_zdf_imp( kt, nit000, 'TRA', r2dt,            tsb, tsa, jpts ) 
    8787         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' zdf1 - Ta: ', mask1=tmask,               & 
    8888         &             tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
     
    9191      IF( l_trdtra )   THEN                      ! save the vertical diffusive trends for further diagnostics 
    9292         DO jk = 1, jpkm1 
    93             ztrdt(:,:,jk) = ( ( tsa(:,:,jk,jp_tem) - tsb(:,:,jk,jp_tem) ) / r2dtra(jk) ) - ztrdt(:,:,jk) 
    94             ztrds(:,:,jk) = ( ( tsa(:,:,jk,jp_sal) - tsb(:,:,jk,jp_sal) ) / r2dtra(jk) ) - ztrds(:,:,jk) 
     93            ztrdt(:,:,jk) = ( ( tsa(:,:,jk,jp_tem) - tsb(:,:,jk,jp_tem) ) / r2dt ) - ztrdt(:,:,jk) 
     94            ztrds(:,:,jk) = ( ( tsa(:,:,jk,jp_sal) - tsb(:,:,jk,jp_sal) ) / r2dt ) - ztrds(:,:,jk) 
    9595         END DO 
    9696         CALL lbc_lnk( ztrdt, 'T', 1. ) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_exp.F90

    r3294 r5282  
    8181      INTEGER                              , INTENT(in   ) ::   kjpt        ! number of tracers 
    8282      INTEGER                              , INTENT(in   ) ::   kn_zdfexp   ! number of sub-time step 
    83       REAL(wp), DIMENSION(        jpk     ), INTENT(in   ) ::   p2dt        ! vertical profile of tracer time-step 
     83      REAL(wp)                             , INTENT(in   ) ::   p2dt        ! tracer time-step 
    8484      REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(in   ) ::   ptb         ! before and now tracer fields 
    8585      REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(inout) ::   pta         ! tracer trend  
     
    136136                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    137137                     ze3tr = zlavmr / fse3t_n(ji,jj,jk) 
    138                      zwx(ji,jj,jk) = zwx(ji,jj,jk) + p2dt(jk) * ( zwy(ji,jj,jk) - zwy(ji,jj,jk+1) ) * ze3tr 
     138                     zwx(ji,jj,jk) = zwx(ji,jj,jk) + p2dt * ( zwy(ji,jj,jk) - zwy(ji,jj,jk+1) ) * ze3tr 
    139139                  END DO 
    140140               END DO 
     
    150150                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    151151                     ze3tb = fse3t_b(ji,jj,jk) / fse3t(ji,jj,jk)                          ! before e3t 
    152                      ztra  = zwx(ji,jj,jk) - ptb(ji,jj,jk,jn) + p2dt(jk) * pta(ji,jj,jk,jn)       ! total trends * 2*rdt  
     152                     ztra  = zwx(ji,jj,jk) - ptb(ji,jj,jk,jn) + p2dt * pta(ji,jj,jk,jn)       ! total trends * 2*rdt 
    153153                     pta(ji,jj,jk,jn) = ( ze3tb * ptb(ji,jj,jk,jn) + ztra ) * tmask(ji,jj,jk) 
    154154                  END DO 
     
    159159               DO jj = 2, jpjm1  
    160160                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    161                      pta(ji,jj,jk,jn) = ( zwx(ji,jj,jk) + p2dt(jk) * pta(ji,jj,jk,jn) ) * tmask(ji,jj,jk) 
     161                     pta(ji,jj,jk,jn) = ( zwx(ji,jj,jk) + p2dt * pta(ji,jj,jk,jn) ) * tmask(ji,jj,jk) 
    162162                  END DO 
    163163               END DO 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90

    r4990 r5282  
    8282      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype   ! =TRA or TRC (tracer indicator) 
    8383      INTEGER                              , INTENT(in   ) ::   kjpt     ! number of tracers 
    84       REAL(wp), DIMENSION(        jpk     ), INTENT(in   ) ::   p2dt     ! vertical profile of tracer time-step 
     84      REAL(wp)                             , INTENT(in   ) ::   p2dt     ! tracer time-step 
    8585      REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(in   ) ::   ptb      ! before and now tracer fields 
    8686      REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(inout) ::   pta      ! tracer trend  
     
    154154                     ze3ta =  ( 1. - r_vvl ) +        r_vvl   * fse3t_a(ji,jj,jk)   ! after scale factor at T-point 
    155155                     ze3tn =         r_vvl   + ( 1. - r_vvl ) * fse3t_n(ji,jj,jk)   ! now   scale factor at T-point 
    156                      zwi(ji,jj,jk) = - p2dt(jk) * zwt(ji,jj,jk  ) / ( ze3tn * fse3w(ji,jj,jk  ) ) 
    157                      zws(ji,jj,jk) = - p2dt(jk) * zwt(ji,jj,jk+1) / ( ze3tn * fse3w(ji,jj,jk+1) ) 
     156                     zwi(ji,jj,jk) = - p2dt * zwt(ji,jj,jk  ) / ( ze3tn * fse3w(ji,jj,jk  ) ) 
     157                     zws(ji,jj,jk) = - p2dt * zwt(ji,jj,jk+1) / ( ze3tn * fse3w(ji,jj,jk+1) ) 
    158158                     zwd(ji,jj,jk) = ze3ta - zwi(ji,jj,jk) - zws(ji,jj,jk) 
    159159                 END DO 
     
    199199               ze3tn = ( 1. - r_vvl ) + r_vvl * fse3t(ji,jj,mikt(ji,jj)) 
    200200               pta(ji,jj,mikt(ji,jj),jn) = ze3tb * ptb(ji,jj,mikt(ji,jj),jn)                     & 
    201                   &                      + p2dt(mikt(ji,jj)) * ze3tn * pta(ji,jj,mikt(ji,jj),jn) 
     201                  &                      + p2dt * ze3tn * pta(ji,jj,mikt(ji,jj),jn) 
    202202               DO jk = mikt(ji,jj)+1, jpkm1 
    203203                  ze3tb = ( 1. - r_vvl ) + r_vvl * fse3t_b(ji,jj,jk) 
    204                   ze3tn = ( 1. - r_vvl ) + r_vvl * fse3t  (ji,jj,jk) 
    205                   zrhs = ze3tb * ptb(ji,jj,jk,jn) + p2dt(jk) * ze3tn * pta(ji,jj,jk,jn)   ! zrhs=right hand side  
     204                  ze3tn = ( 1. - r_vvl ) + r_vvl * fse3t(ji,jj,jk) 
     205                  zrhs = ze3tb * ptb(ji,jj,jk,jn) + p2dt * ze3tn * pta(ji,jj,jk,jn)   ! zrhs=right hand side 
    206206                  pta(ji,jj,jk,jn) = zrhs - zwi(ji,jj,jk) / zwt(ji,jj,jk-1) * pta(ji,jj,jk-1,jn) 
    207207               END DO 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRD/trdmxl.F90

    r4990 r5282  
    3232   USE ioipsl          ! NetCDF library 
    3333   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    34    USE diadimg         ! dimg direct access file format output 
    3534   USE trdmxl_rst      ! restart for diagnosing the ML trends 
    3635   USE prtctl          ! Print control 
     
    298297      !!  
    299298      !! ** Purpose :  Compute and cumulate the mixed layer trends over an analysis 
    300       !!               period, and write NetCDF (or dimg) outputs. 
     299      !!               period, and write NetCDF outputs. 
    301300      !! 
    302301      !! ** Method/usage : 
     
    349348      REAL(wp), POINTER, DIMENSION(:,:  ) ::  ztmltot2, zsmltot2, ztmlres2, zsmlres2, ztmlatf2, zsmlatf2, ztmltrdm2, zsmltrdm2   
    350349      REAL(wp), POINTER, DIMENSION(:,:,:) ::  ztmltrd2, zsmltrd2   ! only needed for mean diagnostics 
    351 #if defined key_dimgout 
    352       INTEGER ::  iyear,imon,iday 
    353       CHARACTER(LEN=80) :: cltext, clmode 
    354 #endif 
    355350      !!---------------------------------------------------------------------- 
    356351   
     
    801796      END IF 
    802797 
    803       IF( nn_cla == 1 )   CALL ctl_warn( '      You set n_cla = 1. Note that the Mixed-Layer diagnostics  ',   & 
    804          &                               '      are not exact along the corresponding straits.            ') 
    805  
    806798      !                                   ! allocate trdmxl arrays 
    807799      IF( trd_mxl_alloc()    /= 0 )   CALL ctl_stop( 'STOP', 'trd_mxl_init : unable to allocate trdmxl     arrays' ) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRD/trdmxl_rst.F90

    r4990 r5282  
    5959            WRITE(numout,*) 
    6060            SELECT CASE ( jprstlib ) 
    61             CASE ( jprstdimg )   ;   WRITE(numout,*) '             open ocean restart_mxl binary file: '//clname 
    6261            CASE DEFAULT         ;   WRITE(numout,*) '             open ocean restart_mxl NetCDF file: '//clname 
    6362            END SELECT 
     
    140139         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~' 
    141140      ENDIF 
    142       IF ( jprstlib == jprstdimg ) THEN 
    143          ! eventually read netcdf file (monobloc)  for restarting on different number of processors 
    144          ! if {cn_trdrst_in}.nc exists, then set jlibalt to jpnf90 
    145          INQUIRE( FILE = TRIM(cn_trdrst_in)//'.nc', EXIST = llok ) 
    146          IF ( llok ) THEN   ;   jlibalt = jpnf90    
    147          ELSE               ;   jlibalt = jprstlib    
    148          ENDIF 
    149       ENDIF 
    150141 
    151142      CALL iom_open( cn_trdrst_in, inum, kiolib = jlibalt )  
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRD/trdtra.F90

    r4990 r5282  
    3636   PUBLIC   trd_tra   ! called by all tra_... modules 
    3737 
    38    REAL(wp) ::   r2dt   ! time-step, = 2 rdttra except at nit000 (=rdttra) if neuler=0 
     38!   REAL(wp) ::   r2dt   ! time-step, = 2 rdt except at nit000 (=rdt) if neuler=0 
    3939 
    4040   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   trdtx, trdty, trdt   ! use to store the temperature trends 
     
    227227      !!---------------------------------------------------------------------- 
    228228 
    229       IF( neuler == 0 .AND. kt == nit000    ) THEN   ;   r2dt =      rdt      ! = rdtra (restart with Euler time stepping) 
    230       ELSEIF(               kt <= nit000 + 1) THEN   ;   r2dt = 2. * rdt      ! = 2 rdttra (leapfrog) 
     229      IF( neuler == 0 .AND. kt == nit000    ) THEN   ;   r2dt =      rdt      ! = rdt (restart with Euler time stepping) 
     230      ELSEIF(               kt <= nit000 + 1) THEN   ;   r2dt = 2. * rdt      ! = 2 rdt (leapfrog) 
    231231      ENDIF 
    232232 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/TRD/trdvor.F90

    r4990 r5282  
    8585      !!  
    8686      !! ** Purpose :  computation of cumulated trends over analysis period 
    87       !!               and make outputs (NetCDF or DIMG format) 
     87      !!               and make outputs (NetCDF format) 
    8888      !!---------------------------------------------------------------------- 
    8989      REAL(wp), DIMENSION(:,:,:), INTENT(inout) ::   putrd, pvtrd   ! U and V trends  
     
    318318      !!  
    319319      !! ** Purpose :  computation of cumulated trends over analysis period 
    320       !!               and make outputs (NetCDF or DIMG format) 
     320      !!               and make outputs (NetCDF format) 
    321321      !!---------------------------------------------------------------------- 
    322322      INTEGER                   , INTENT(in   ) ::   kt             ! time step 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90

    r4990 r5282  
    4242   !!---------------------------------------------------------------------- 
    4343   USE step_oce        ! module used in the ocean time stepping module 
    44    USE cla             ! cross land advection               (tra_cla routine) 
    4544   USE domcfg          ! domain configuration               (dom_cfg routine) 
    4645   USE mppini          ! shared/distributed memory setting (mpp_init routine) 
     
    388387         &                  CALL bdytide_init   ! Open boundaries initialisation of tidal harmonic forcing 
    389388 
    390                             CALL dyn_nept_init  ! simplified form of Neptune effect 
    391389 
    392390      !      
     
    429427 
    430428      !                                     ! Misc. options 
    431       IF( nn_cla == 1 .AND. cp_cfg == 'orca' .AND. jp_cfg == 2 )   CALL cla_init       ! Cross Land Advection 
    432429                            CALL icb_init( rdt, nit000)   ! initialise icebergs instance 
    433430      
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/step.F90

    r5012 r5282  
    158158      IF( lk_traldf_eiv )   CALL ldf_eiv( kstp )      ! eddy induced velocity coefficient 
    159159#endif 
    160 #if defined key_traldf_c3d && key_traldf_smag 
    161                           CALL ldf_tra_smag( kstp )      ! eddy induced velocity coefficient 
    162 #  endif 
    163 #if defined key_dynldf_c3d && key_dynldf_smag 
    164                           CALL ldf_dyn_smag( kstp )      ! eddy induced velocity coefficient 
    165 #  endif 
    166160 
    167161      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     
    185179          IF(  ln_asmiau .AND. & 
    186180             & ln_dyninc       )  CALL dyn_asm_inc  ( kstp )   ! apply dynamics assimilation increment 
    187           IF( ln_neptsimp )       CALL dyn_nept_cor ( kstp )   ! subtract Neptune velocities (simplified) 
    188181          IF( lk_bdy           )  CALL bdy_dyn3d_dmp( kstp )   ! bdy damping trends 
    189182                                  CALL dyn_adv      ( kstp )   ! advection (vector or flux form) 
    190183                                  CALL dyn_vor      ( kstp )   ! vorticity term including Coriolis 
    191184                                  CALL dyn_ldf      ( kstp )   ! lateral mixing 
    192           IF( ln_neptsimp )       CALL dyn_nept_cor ( kstp )   ! add Neptune velocities (simplified) 
    193185#if defined key_agrif 
    194186          IF(.NOT. Agrif_Root())  CALL Agrif_Sponge_dyn        ! momentum sponge 
     
    289281           & ln_dyninc      )  CALL dyn_asm_inc( kstp )     ! apply dynamics assimilation increment 
    290282        IF( ln_bkgwri )        CALL asm_bkg_wri( kstp )     ! output background fields 
    291         IF( ln_neptsimp )      CALL dyn_nept_cor( kstp )    ! subtract Neptune velocities (simplified) 
    292283        IF( lk_bdy          )  CALL bdy_dyn3d_dmp(kstp )    ! bdy damping trends 
    293284                               CALL dyn_adv( kstp )         ! advection (vector or flux form) 
    294285                               CALL dyn_vor( kstp )         ! vorticity term including Coriolis 
    295286                               CALL dyn_ldf( kstp )         ! lateral mixing 
    296         IF( ln_neptsimp )      CALL dyn_nept_cor( kstp )    ! add Neptune velocities (simplified) 
    297287#if defined key_agrif 
    298288        IF(.NOT. Agrif_Root()) CALL Agrif_Sponge_dyn        ! momemtum sponge 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/step_oce.F90

    r4990 r5282  
    5050   USE dynspg_oce       ! surface pressure gradient        (dyn_spg routine) 
    5151   USE dynspg           ! surface pressure gradient        (dyn_spg routine) 
    52    USE dynnept          ! simp. form of Neptune effect(dyn_nept_cor routine) 
    5352 
    5453   USE dynnxt           ! time-stepping                    (dyn_nxt routine) 
     
    6867   USE ldfslp           ! iso-neutral slopes               (ldf_slp routine) 
    6968   USE ldfeiv           ! eddy induced velocity coef.      (ldf_eiv routine) 
    70    USE ldftra_smag      ! Smagirinsky diffusion            (ldftra_smag routine) 
    71    USE ldfdyn_smag      ! Smagorinsky viscosity            (ldfdyn_smag routine)  
    7269 
    7370   USE zdftmx           ! tide-induced vertical mixing     (zdf_tmx routine) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/SAS_SRC/daymod.F90

    r4162 r5282  
    2020   !! 
    2121   !!   we suppose that the time step is deviding the number of second of in a day 
    22    !!             ---> MOD( rday, rdttra(1) ) == 0 
     22   !!             ---> MOD( rday, rdt ) == 0 
    2323   !! 
    2424   !!           ----------- WARNING ----------- 
     
    7171      !!---------------------------------------------------------------------- 
    7272      ! 
    73       ! all calendar staff is based on the fact that MOD( rday, rdttra(1) ) == 0 
    74       IF( MOD( rday     , rdttra(1) ) /= 0. )   CALL ctl_stop( 'the time step must devide the number of second of in a day' ) 
     73      ! all calendar staff is based on the fact that MOD( rday, rdt ) == 0 
     74      IF( MOD( rday     , rdt ) /= 0. )   CALL ctl_stop( 'the time step must devide the number of second of in a day' ) 
    7575      IF( MOD( rday     , 2.        ) /= 0. )   CALL ctl_stop( 'the number of second of in a day must be an even number'    ) 
    76       IF( MOD( rdttra(1), 2.        ) /= 0. )   CALL ctl_stop( 'the time step (in second) must be an even number'           ) 
     76      IF( MOD( rdt, 2.        ) /= 0. )   CALL ctl_stop( 'the time step (in second) must be an even number'           ) 
    7777      nsecd   = NINT(rday           ) 
    7878      nsecd05 = NINT(0.5 * rday     ) 
    79       ndt     = NINT(      rdttra(1)) 
    80       ndt05   = NINT(0.5 * rdttra(1)) 
     79      ndt     = NINT(      rdt) 
     80      ndt05   = NINT(0.5 * rdt) 
    8181 
    8282      ! parameters corresponding to nit000 - 1 (as we start the step loop with a call to day) 
    8383      ndastp = ndate0 - 1        ! ndate0 read in the namelist in dom_nam, we assume that we start run at 00:00 
    84       adatrj = ( REAL( nit000-1, wp ) * rdttra(1) ) / rday 
     84      adatrj = ( REAL( nit000-1, wp ) * rdt ) / rday 
    8585      IF( ABS(adatrj  - REAL(NINT(adatrj),wp)) < 0.1 / rday )   adatrj = REAL(NINT(adatrj),wp)   ! avoid truncation error 
    8686      ! 
     
    226226      nsec_week  = nsec_week  + ndt 
    227227      nsec_day   = nsec_day   + ndt                 
    228       adatrj  = adatrj  + rdttra(1) / rday 
    229       fjulday = fjulday + rdttra(1) / rday 
     228      adatrj  = adatrj  + rdt / rday 
     229      fjulday = fjulday + rdt / rday 
    230230      IF( ABS(fjulday - REAL(NINT(fjulday),wp)) < zprec )   fjulday = REAL(NINT(fjulday),wp)   ! avoid truncation error 
    231231      IF( ABS(adatrj  - REAL(NINT(adatrj ),wp)) < zprec )   adatrj  = REAL(NINT(adatrj ),wp)   ! avoid truncation error 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/SAS_SRC/diawri.F90

    r4292 r5282  
    4040   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    4141   USE in_out_manager  ! I/O manager 
    42    USE diadimg         ! dimg direct access file format output 
    4342   USE diaar5, ONLY :   lk_diaar5 
    4443   USE iom 
     
    8584  END FUNCTION dia_wri_alloc 
    8685 
    87 #if defined key_dimgout 
    88    !!---------------------------------------------------------------------- 
    89    !!   'key_dimgout'                                      DIMG output file 
    90    !!---------------------------------------------------------------------- 
    91 #   include "diawri_dimg.h90" 
    92  
    93 #else 
    9486   !!---------------------------------------------------------------------- 
    9587   !!   Default option                                   NetCDF output file 
    9688   !!---------------------------------------------------------------------- 
    97 # if defined key_iomput 
     89#if defined key_iomput 
    9890   !!---------------------------------------------------------------------- 
    9991   !!   'key_iomput'                                        use IOM library 
     
    164156      ! Define frequency of output and means 
    165157      zdt = rdt 
    166       IF( nacc == 1 ) zdt = rdtmin 
    167158      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
    168159      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
     
    338329      ! 
    339330   END SUBROUTINE dia_wri 
    340 # endif 
    341  
    342331#endif 
     332 
     333 
    343334 
    344335   SUBROUTINE dia_wri_state( cdfile_name, kt ) 
     
    435426      ! ----------------- 
    436427      CALL histclo( id_i ) 
    437 #if ! defined key_iomput && ! defined key_dimgout 
     428#if ! defined key_iomput 
    438429      IF( ninist /= 1  ) THEN 
    439430         CALL histclo( nid_T ) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/SAS_SRC/nemogcm.F90

    r4624 r5282  
    4242   USE step_oce        ! module used in the ocean time stepping module 
    4343   USE sbc_oce         ! surface boundary condition: ocean 
    44    USE cla             ! cross land advection               (tra_cla routine) 
    4544   USE domcfg          ! domain configuration               (dom_cfg routine) 
    4645   USE daymod          ! calendar 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/TOP_SRC/PISCES/SED/sedrst.F90

    r3443 r5282  
    5959 
    6060      ALLOCATE( zdta(jpi,jpj,jpksed,jptrased), zdta1(jpi,jpj,jpksed,2), zhipor(jpoce,jpksed) )  
    61  
    62       IF ( jprstlib == jprstdimg ) THEN 
    63         ! eventually read netcdf file (monobloc)  for restarting on different number of processors 
    64         ! if restart_sed.nc exists, then set jlibalt to jpnf90 
    65         INQUIRE( FILE = 'restart_sed.nc', EXIST = llok ) 
    66         IF ( llok ) THEN ; jlibalt = jpnf90  ; ELSE ; jlibalt = jprstlib ; ENDIF 
    67       ENDIF 
    6861 
    6962      CALL iom_open( 'restart_sed', numrsr, kiolib = jlibalt )      
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/TOP_SRC/PISCES/trcini_pisces.F90

    r4996 r5282  
    107107      CALL p4z_sms_init       !  Maint routine 
    108108      !                                            ! Time-step 
    109       rfact   = rdttrc(1)                          ! --------- 
     109      rfact   = rdttrc                          ! --------- 
    110110      rfactr  = 1. / rfact 
    111111      rfact2  = rfact / FLOAT( nrdttrc ) 
    112112      rfact2r = 1. / rfact2 
    113113 
    114       IF(lwp) WRITE(numout,*) '    Passive Tracer  time step    rfact  = ', rfact, ' rdt = ', rdttra(1) 
     114      IF(lwp) WRITE(numout,*) '    Passive Tracer  time step    rfact  = ', rfact, ' rdt = ', rdt 
    115115      IF(lwp) write(numout,*) '    PISCES  Biology time step    rfact2 = ', rfact2 
    116116 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/TOP_SRC/TRP/trcadv.F90

    r4610 r5282  
    3232 
    3333   PUBLIC   trc_adv          ! routine called by step module 
    34    PUBLIC   trc_adv_alloc    ! routine called by nemogcm module 
    35  
    36    INTEGER ::   nadv   ! choice of the type of advection scheme 
    37    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:) ::   r2dt  ! vertical profile time-step, = 2 rdttra 
    38    !                                                    ! except at nitrrc000 (=rdttra) if neuler=0 
     34 
     35   INTEGER ::   nadv    ! choice of the type of advection scheme 
     36   REAL(wp)::   r2dtrc  ! time-step, = 2 rdttrc 
     37!   !                   ! except at nitrrc000 (=rdttrc) if neuler=0 
    3938 
    4039   !! * Substitutions 
     
    4746   !!---------------------------------------------------------------------- 
    4847CONTAINS 
    49  
    50    INTEGER FUNCTION trc_adv_alloc() 
    51       !!---------------------------------------------------------------------- 
    52       !!                  ***  ROUTINE trc_adv_alloc  *** 
    53       !!---------------------------------------------------------------------- 
    54  
    55       ALLOCATE( r2dt(jpk), STAT=trc_adv_alloc ) 
    56  
    57       IF( trc_adv_alloc /= 0 ) CALL ctl_warn('trc_adv_alloc : failed to allocate array.') 
    58  
    59    END FUNCTION trc_adv_alloc 
    60  
    6148 
    6249   SUBROUTINE trc_adv( kt ) 
     
    8471 
    8572      IF( ln_top_euler) THEN 
    86          r2dt(:) =  rdttrc(:)              ! = rdttrc (use Euler time stepping) 
     73         r2dtrc =  rdttrc              ! = rdttrc (use Euler time stepping) 
    8774      ELSE 
    8875         IF( neuler == 0 .AND. kt == nittrc000 ) THEN     ! at nittrc000 
    89             r2dt(:) =  rdttrc(:)           ! = rdttrc (restarting with Euler time stepping) 
     76            r2dtrc =  rdttrc           ! = rdttrc (restarting with Euler time stepping) 
    9077         ELSEIF( kt <= nittrc000 + 1 ) THEN          ! at nittrc000 or nittrc000+1 
    91             r2dt(:) = 2. * rdttrc(:)       ! = 2 rdttrc (leapfrog) 
     78            r2dtrc = 2. * rdttrc       ! = 2 rdttrc (leapfrog) 
    9279         ENDIF 
    9380      ENDIF 
     
    118105      SELECT CASE ( nadv )                            !==  compute advection trend and add it to general trend  ==! 
    119106      CASE ( 1 )   ;    CALL tra_adv_cen2  ( kt, nittrc000, 'TRC',       zun, zvn, zwn, trb, trn, tra, jptra )   !  2nd order centered 
    120       CASE ( 2 )   ;    CALL tra_adv_tvd   ( kt, nittrc000, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )   !  TVD  
    121       CASE ( 3 )   ;    CALL tra_adv_muscl ( kt, nittrc000, 'TRC', r2dt, zun, zvn, zwn, trb,      tra, jptra, ln_trcadv_msc_ups )   !  MUSCL  
    122       CASE ( 4 )   ;    CALL tra_adv_muscl2( kt, nittrc000, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )   !  MUSCL2  
    123       CASE ( 5 )   ;    CALL tra_adv_ubs   ( kt, nittrc000, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )   !  UBS  
    124       CASE ( 6 )   ;    CALL tra_adv_qck   ( kt, nittrc000, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )   !  QUICKEST  
     107      CASE ( 2 )   ;    CALL tra_adv_tvd   ( kt, nittrc000, 'TRC', r2dtrc, zun, zvn, zwn, trb, trn, tra, jptra )   !  TVD 
     108      CASE ( 3 )   ;    CALL tra_adv_muscl ( kt, nittrc000, 'TRC', r2dtrc, zun, zvn, zwn, trb,      tra, jptra, ln_trcadv_msc_ups )   !  MUSCL 
     109      CASE ( 4 )   ;    CALL tra_adv_muscl2( kt, nittrc000, 'TRC', r2dtrc, zun, zvn, zwn, trb, trn, tra, jptra )   !  MUSCL2 
     110      CASE ( 5 )   ;    CALL tra_adv_ubs   ( kt, nittrc000, 'TRC', r2dtrc, zun, zvn, zwn, trb, trn, tra, jptra )   !  UBS 
     111      CASE ( 6 )   ;    CALL tra_adv_qck   ( kt, nittrc000, 'TRC', r2dtrc, zun, zvn, zwn, trb, trn, tra, jptra )   !  QUICKEST 
    125112      ! 
    126113      CASE (-1 )                                      !==  esopa: test all possibility with control print  ==! 
     
    128115         WRITE(charout, FMT="('adv1')")  ; CALL prt_ctl_trc_info(charout) 
    129116                                           CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 
    130          CALL tra_adv_tvd   ( kt, nittrc000, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )           
     117         CALL tra_adv_tvd   ( kt, nittrc000, 'TRC', r2dtrc, zun, zvn, zwn, trb, trn, tra, jptra ) 
    131118         WRITE(charout, FMT="('adv2')")  ; CALL prt_ctl_trc_info(charout) 
    132119                                           CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 
    133          CALL tra_adv_muscl ( kt, nittrc000, 'TRC', r2dt, zun, zvn, zwn, trb,      tra, jptra, ln_trcadv_msc_ups  )           
     120         CALL tra_adv_muscl ( kt, nittrc000, 'TRC', r2dtrc, zun, zvn, zwn, trb,      tra, jptra, ln_trcadv_msc_ups  ) 
    134121         WRITE(charout, FMT="('adv3')")  ; CALL prt_ctl_trc_info(charout) 
    135122                                           CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 
    136          CALL tra_adv_muscl2( kt, nittrc000, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )           
     123         CALL tra_adv_muscl2( kt, nittrc000, 'TRC', r2dtrc, zun, zvn, zwn, trb, trn, tra, jptra ) 
    137124         WRITE(charout, FMT="('adv4')")  ; CALL prt_ctl_trc_info(charout) 
    138125                                           CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 
    139          CALL tra_adv_ubs   ( kt, nittrc000, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )           
     126         CALL tra_adv_ubs   ( kt, nittrc000, 'TRC', r2dtrc, zun, zvn, zwn, trb, trn, tra, jptra ) 
    140127         WRITE(charout, FMT="('adv5')")  ; CALL prt_ctl_trc_info(charout) 
    141128                                           CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 
    142          CALL tra_adv_qck   ( kt, nittrc000, 'TRC', r2dt, zun, zvn, zwn, trb, trn, tra, jptra )           
     129         CALL tra_adv_qck   ( kt, nittrc000, 'TRC', r2dtrc, zun, zvn, zwn, trb, trn, tra, jptra ) 
    143130         WRITE(charout, FMT="('adv6')")  ; CALL prt_ctl_trc_info(charout) 
    144131                                           CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd') 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/TOP_SRC/TRP/trcnxt.F90

    r4990 r5282  
    4141 
    4242   PUBLIC   trc_nxt          ! routine called by step.F90 
    43    PUBLIC   trc_nxt_alloc    ! routine called by nemogcm.F90 
    4443 
    45    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:) ::   r2dt 
     44   REAL(wp) ::   r2dtrc 
    4645 
    4746   !!---------------------------------------------------------------------- 
     
    5150   !!---------------------------------------------------------------------- 
    5251CONTAINS 
    53  
    54    INTEGER FUNCTION trc_nxt_alloc() 
    55       !!---------------------------------------------------------------------- 
    56       !!                   ***  ROUTINE trc_nxt_alloc  *** 
    57       !!---------------------------------------------------------------------- 
    58       ALLOCATE( r2dt(jpk), STAT=trc_nxt_alloc ) 
    59       ! 
    60       IF( trc_nxt_alloc /= 0 )   CALL ctl_warn('trc_nxt_alloc : failed to allocate array') 
    61       ! 
    62    END FUNCTION trc_nxt_alloc 
    63  
    6452 
    6553   SUBROUTINE trc_nxt( kt ) 
     
    117105 
    118106      ! set time step size (Euler/Leapfrog) 
    119       IF( neuler == 0 .AND. kt ==  nittrc000 ) THEN  ;  r2dt(:) =     rdttrc(:)   ! at nittrc000             (Euler) 
    120       ELSEIF( kt <= nittrc000 + 1 )            THEN  ;  r2dt(:) = 2.* rdttrc(:)   ! at nit000 or nit000+1 (Leapfrog) 
     107      IF( neuler == 0 .AND. kt ==  nittrc000 ) THEN  ;  r2dtrc =     rdttrc   ! at nittrc000             (Euler) 
     108      ELSEIF( kt <= nittrc000 + 1 )            THEN  ;  r2dtrc = 2.* rdttrc   ! at nit000 or nit000+1 (Leapfrog) 
    121109      ENDIF 
    122110 
     
    143131 
    144132      ! trends computation 
    145       IF( l_trdtrc ) THEN                                      ! trends 
     133      IF( l_trdtrc ) THEN! trends 
    146134         DO jn = 1, jptra 
     135         zfact = 1.e0 / r2dtrc 
    147136            DO jk = 1, jpkm1 
    148                zfact = 1.e0 / r2dt(jk)   
    149137               ztrdt(:,:,jk,jn) = ( trb(:,:,jk,jn) - ztrdt(:,:,jk,jn) ) * zfact  
    150138               CALL trd_tra( kt, 'TRC', jn, jptra_atf, ztrdt ) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/TOP_SRC/TRP/trczdf.F90

    r4990 r5282  
    2727 
    2828   PUBLIC   trc_zdf          ! called by step.F90  
    29    PUBLIC   trc_zdf_alloc    ! called by nemogcm.F90  
    3029 
    3130   INTEGER ::   nzdf = 0               ! type vertical diffusion algorithm used 
    3231      !                                ! defined from ln_zdf...  namlist logicals) 
    33    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:) ::  r2dt   ! vertical profile time-step, = 2 rdttra 
    34       !                                                 ! except at nittrc000 (=rdttra) if neuler=0 
     32   REAL(wp) ::  r2dtrc   ! vertical profile time-step, = 2 rdt 
     33      !                  ! except at nittrc000 (=rdt) if neuler=0 
    3534 
    3635   !! * Substitutions 
     
    4443   !!---------------------------------------------------------------------- 
    4544CONTAINS 
    46     
    47    INTEGER FUNCTION trc_zdf_alloc() 
    48       !!---------------------------------------------------------------------- 
    49       !!                  ***  ROUTINE trc_zdf_alloc  *** 
    50       !!---------------------------------------------------------------------- 
    51       ALLOCATE( r2dt(jpk) , STAT=trc_zdf_alloc ) 
    52       ! 
    53       IF( trc_zdf_alloc /= 0 )   CALL ctl_warn('trc_zdf_alloc : failed to allocate array.') 
    54       ! 
    55    END FUNCTION trc_zdf_alloc 
    56  
    5745 
    5846   SUBROUTINE trc_zdf( kt ) 
     
    7462 
    7563      IF( ln_top_euler) THEN 
    76          r2dt(:) =  rdttrc(:)              ! = rdttrc (use Euler time stepping) 
     64         r2dtrc =  rdttrc              ! = rdttrc (use Euler time stepping) 
    7765      ELSE 
    7866         IF( neuler == 0 .AND. kt == nittrc000 ) THEN     ! at nittrc000 
    79             r2dt(:) =  rdttrc(:)           ! = rdttrc (restarting with Euler time stepping) 
     67            r2dtrc =  rdttrc           ! = rdttrc (restarting with Euler time stepping) 
    8068         ELSEIF( kt <= nittrc000 + 1 ) THEN          ! at nittrc000 or nittrc000+1 
    81             r2dt(:) = 2. * rdttrc(:)       ! = 2 rdttrc (leapfrog) 
     69            r2dtrc = 2. * rdttrc       ! = 2 rdttrc (leapfrog) 
    8270         ENDIF 
    8371      ENDIF 
     
    9078      SELECT CASE ( nzdf )                       ! compute lateral mixing trend and add it to the general trend 
    9179      CASE ( -1 )                                       ! esopa: test all possibility with control print 
    92          CALL tra_zdf_exp( kt, nittrc000, 'TRC', r2dt, nn_trczdf_exp, trb, tra, jptra )  
     80         CALL tra_zdf_exp( kt, nittrc000, 'TRC', r2dtrc, nn_trczdf_exp, trb, tra, jptra ) 
    9381         WRITE(charout, FMT="('zdf1 ')") ;  CALL prt_ctl_trc_info(charout) 
    9482                                            CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
    95          CALL tra_zdf_imp( kt, nittrc000, 'TRC', r2dt,                trb, tra, jptra )  
     83         CALL tra_zdf_imp( kt, nittrc000, 'TRC', r2dtrc,                trb, tra, jptra ) 
    9684         WRITE(charout, FMT="('zdf2 ')") ;  CALL prt_ctl_trc_info(charout) 
    9785                                            CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' ) 
    98       CASE ( 0 ) ;  CALL tra_zdf_exp( kt, nittrc000, 'TRC', r2dt, nn_trczdf_exp, trb, tra, jptra )    !   explicit scheme  
    99       CASE ( 1 ) ;  CALL tra_zdf_imp( kt, nittrc000, 'TRC', r2dt,                trb, tra, jptra )    !   implicit scheme           
     86      CASE ( 0 ) ;  CALL tra_zdf_exp( kt, nittrc000, 'TRC', r2dtrc, nn_trczdf_exp, trb, tra, jptra )    !   explicit scheme 
     87      CASE ( 1 ) ;  CALL tra_zdf_imp( kt, nittrc000, 'TRC', r2dtrc,                trb, tra, jptra )    !   implicit scheme 
    10088 
    10189      END SELECT 
     
    10492         DO jn = 1, jptra 
    10593            DO jk = 1, jpkm1 
    106                ztrtrd(:,:,jk,jn) = ( ( tra(:,:,jk,jn) - trb(:,:,jk,jn) ) / r2dt(jk) ) - ztrtrd(:,:,jk,jn) 
     94               ztrtrd(:,:,jk,jn) = ( ( tra(:,:,jk,jn) - trb(:,:,jk,jn) ) / r2dtrc ) - ztrtrd(:,:,jk,jn) 
    10795            END DO 
    10896            CALL trd_tra( kt, 'TRC', jn, jptra_zdf, ztrtrd(:,:,:,jn) ) 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/TOP_SRC/TRP/trdmxl_trc.F90

    r4990 r5282  
    331331      !!  
    332332      !! ** Purpose :  Compute and cumulate the mixed layer trends over an analysis 
    333       !!               period, and write NetCDF (or dimg) outputs. 
     333      !!               period, and write NetCDF outputs. 
    334334      !! 
    335335      !! ** Method/usage : 
     
    391391      ! 
    392392      CHARACTER (LEN=10) ::   clvar 
    393 #if defined key_dimgout 
    394       INTEGER ::   iyear,imon,iday 
    395       CHARACTER(LEN=80) ::   cltext, clmode 
    396 #endif 
    397393      !!---------------------------------------------------------------------- 
    398394 
     
    775771      ! ====================================================================== 
    776772 
    777       ! IV.1 Code for dimg mpp output 
    778       ! ----------------------------- 
    779  
    780 # if defined key_dimgout 
    781       STOP 'Not implemented' 
    782 # else 
     773!# if 
     774!# else 
    783775       
    784       ! IV.2 Code for IOIPSL/NetCDF output 
     776      ! IV.1 Code for IOIPSL/NetCDF output 
    785777      ! ---------------------------------- 
    786778 
     
    866858      icount = 1 
    867859 
    868 # endif /* key_dimgout */ 
     860!# endif 
    869861 
    870862      IF( MOD( itmod, nn_trd_trc ) == 0 ) THEN 
     
    897889      !! 
    898890      !! ** Purpose :  Compute and cumulate the mixed layer biological trends over an analysis 
    899       !!               period, and write NetCDF (or dimg) outputs. 
     891      !!               period, and write NetCDF outputs. 
    900892      !! 
    901893      !! ** Method/usage : 
     
    944936      LOGICAL  :: llwarn  = .TRUE., lldebug = .TRUE. 
    945937      REAL(wp) :: zfn, zfn2 
    946 #if defined key_dimgout 
    947       INTEGER ::  iyear,imon,iday 
    948       CHARACTER(LEN=80) :: cltext, clmode 
    949 #endif 
    950938      !!---------------------------------------------------------------------- 
    951939      ! ... Warnings 
     
    10561044      ! ====================================================================== 
    10571045 
    1058       ! IV.1 Code for dimg mpp output 
    1059       ! ----------------------------- 
    1060  
    1061 # if defined key_dimgout 
    1062       STOP 'Not implemented' 
    1063 # else 
    1064  
    1065       ! IV.2 Code for IOIPSL/NetCDF output 
     1046!# if 
     1047!# else 
     1048 
     1049      ! IV.1 Code for IOIPSL/NetCDF output 
    10661050      ! ---------------------------------- 
    10671051 
     
    11081092 
    11091093 
    1110 # endif /* key_dimgout */ 
     1094!# endif 
    11111095 
    11121096      IF( MOD( itmod, nn_trd_trc ) == 0 ) THEN 
     
    12591243      ! ====================================================================== 
    12601244 
    1261 #if defined key_dimgout  
    1262       ??? 
    1263 #else 
     1245!#if  
     1246!#else 
    12641247      ! clmxl = legend root for netCDF output 
    12651248      IF( nn_ctls_trc == 0 ) THEN                                   ! control surface = mixed-layer with density criterion 
     
    14041387#endif 
    14051388 
    1406 #endif        /* key_dimgout */ 
     1389!#endif 
    14071390   END SUBROUTINE trd_mxl_trc_init 
    14081391 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/TOP_SRC/TRP/trdmxl_trc_rst.F90

    r4990 r5282  
    141141      ENDIF 
    142142       
    143       IF ( jprstlib == jprstdimg ) THEN 
    144         ! eventually read netcdf file (monobloc)  for restarting on different number of processors 
    145         ! if {cn_trdrst_trc_in}.nc exists, then set jlibalt to jpnf90 
    146         INQUIRE( FILE = TRIM(cn_trdrst_trc_in)//'.nc', EXIST = llok ) 
    147         IF ( llok ) THEN ; jlibalt = jpnf90  ; ELSE ; jlibalt = jprstlib ; ENDIF 
    148       ENDIF 
    149  
    150143      CALL iom_open( cn_trdrst_trc_in, inum, kiolib = jlibalt )  
    151144       
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/TOP_SRC/trc.F90

    r4990 r5282  
    5555   CHARACTER(len = 80) , PUBLIC                                    ::  cn_trcrst_in   !: suffix of pass. tracer restart name (input) 
    5656   CHARACTER(len = 80) , PUBLIC                                    ::  cn_trcrst_out  !: suffix of pass. tracer restart name (output) 
    57    REAL(wp)            , PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:)   ::  rdttrc         !: vertical profile of passive tracer time step 
     57   REAL(wp)            , PUBLIC                                    ::  rdttrc         !: passive tracer time step 
    5858   LOGICAL             , PUBLIC                                    ::  ln_top_euler  !: boolean term for euler integration  
    5959   LOGICAL             , PUBLIC                                    ::  ln_trcdta      !: Read inputs data from files 
     
    187187         &      gtru (jpi,jpj,jptra)  , gtrv (jpi,jpj,jptra)                          ,       & 
    188188         &      gtrui(jpi,jpj,jptra)  , gtrvi(jpi,jpj,jptra)                          ,       & 
    189          &      cvol(jpi,jpj,jpk)     , rdttrc(jpk)           , trai(jptra)           ,       & 
     189         &      cvol(jpi,jpj,jpk)     , trai(jptra)           ,       & 
    190190         &      ctrcnm(jptra)         , ctrcln(jptra)         , ctrcun(jptra)         ,       &  
    191191         &      ln_trc_ini(jptra)     , ln_trc_wri(jptra)                             ,  STAT = trc_alloc  )   
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/TOP_SRC/trcini.F90

    r4990 r5282  
    7676 
    7777      IF( ltrcdm2dc )CALL ctl_warn( ' Diurnal cycle on physics but not in PISCES or LOBSTER ' ) 
    78  
    79       IF( nn_cla == 1 )   & 
    80          &  CALL ctl_stop( ' Cross Land Advection not yet implemented with passive tracer ; nn_cla must be 0' ) 
    8178 
    8279      CALL trc_nam      ! read passive tracers namelists 
     
    187184      !! ** Purpose :   Allocate all the dynamic arrays of the OPA modules 
    188185      !!---------------------------------------------------------------------- 
    189       USE trcadv        , ONLY:   trc_adv_alloc          ! TOP-related alloc routines... 
    190186      USE trc           , ONLY:   trc_alloc 
    191       USE trcnxt        , ONLY:   trc_nxt_alloc 
    192       USE trczdf        , ONLY:   trc_zdf_alloc 
    193187      USE trdtrc_oce    , ONLY:   trd_trc_oce_alloc 
    194188#if defined key_trdmxl_trc  
     
    199193      !!---------------------------------------------------------------------- 
    200194      ! 
    201       ierr =        trc_adv_alloc()          ! Start of TOP-related alloc routines... 
    202       ierr = ierr + trc_alloc    () 
    203       ierr = ierr + trc_nxt_alloc() 
    204       ierr = ierr + trc_zdf_alloc() 
     195      ierr =        trc_alloc() 
    205196      ierr = ierr + trd_trc_oce_alloc() 
    206197#if defined key_trdmxl_trc  
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/TOP_SRC/trcnam.F90

    r4990 r5282  
    109109      ENDIF 
    110110 
    111        
    112       rdttrc(:) = rdttra(:) * FLOAT( nn_dttrc )   ! vertical profile of passive tracer time-step 
     111      rdttrc = rdt * FLOAT( nn_dttrc )   !  passive tracer time-step 
    113112   
    114113      IF(lwp) THEN                   ! control print 
    115114        WRITE(numout,*)  
    116         WRITE(numout,*) '    Passive Tracer  time step    rdttrc  = ', rdttrc(1) 
     115        WRITE(numout,*) '    Passive Tracer  time step    rdttrc  = ', rdttrc 
    117116        WRITE(numout,*)  
    118117      ENDIF 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

    r4990 r5282  
    123123      !!---------------------------------------------------------------------- 
    124124      ! 
    125       CALL iom_rstput( kt, nitrst, numrtw, 'rdttrc1', rdttrc(1) )   ! surface passive tracer time step 
     125      CALL iom_rstput( kt, nitrst, numrtw, 'rdttrc1', rdttrc )   ! passive tracer time step 
    126126      ! prognostic variables  
    127127      ! --------------------  
     
    186186         IF(lwp) WRITE(numout,*) 'trc_rst_cal : read the TOP restart file for calendar' 
    187187         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 
    188  
    189          IF ( jprstlib == jprstdimg ) THEN 
    190            ! eventually read netcdf file (monobloc)  for restarting on different number of processors 
    191            ! if {cn_trcrst_in}.nc exists, then set jlibalt to jpnf90  
    192            INQUIRE( FILE = TRIM(cn_trcrst_in)//'.nc', EXIST = llok ) 
    193            IF ( llok ) THEN ; jlibalt = jpnf90  ; ELSE ; jlibalt = jprstlib ; ENDIF 
    194          ENDIF 
    195188 
    196189         CALL iom_open( cn_trcrst_in, numrtr, kiolib = jlibalt ) 
     
    226219            ELSE 
    227220               ndastp = ndate0 - 1     ! ndate0 read in the namelist in dom_nam 
    228                adatrj = ( REAL( nittrc000-1, wp ) * rdttra(1) ) / rday 
     221               adatrj = ( REAL( nittrc000-1, wp ) * rdt ) / rday 
    229222               ! note this is wrong if time step has changed during run 
    230223            ENDIF 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/SETTE/iodef_sette.xml

    r4773 r5282  
    3030      <file_group id="6h" output_freq="6h"  output_level="10" enabled=".TRUE."/> <!-- 6h files -->      
    3131      <file_group id="1d" output_freq="1d"  output_level="10" enabled=".TRUE."/> <!-- 1d files --> 
    32       <file_group id="3d" output_freq="3d"  output_level="10" enabled=".TRUE."/> <!-- 3d files -->     
     32      <file_group id="3d" output_freq="3d"  output_level="10" enabled=".TRUE."/> <!-- 3d files --> 
    3333      <file_group id="5d" output_freq="5d"  output_level="10" enabled=".TRUE."/>  <!-- 5d files -->   
    3434      <file_group id="1m" output_freq="1mo" output_level="10" enabled=".TRUE."/> <!-- real monthly files --> 
     
    102102      <file_group id="6h" output_freq="6h"  output_level="10" enabled=".TRUE."/> <!-- 6h files -->      
    103103      <file_group id="1d" output_freq="1d"  output_level="10" enabled=".TRUE."/> <!-- 1d files --> 
    104       <file_group id="3d" output_freq="3d"  output_level="10" enabled=".TRUE."/> <!-- 3d files -->     
     104      <file_group id="3d" output_freq="3d"  output_level="10" enabled=".TRUE."/> <!-- 3d files --> 
    105105      <file_group id="5d" output_freq="5d"  output_level="10" enabled=".TRUE."/>  <!-- 5d files -->   
    106106      <file_group id="1m" output_freq="1mo" output_level="10" enabled=".TRUE."/> <!-- real monthly files --> 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/SETTE/sette.sh

    r4990 r5282  
    8888# 
    8989# Compiler among those in NEMOGCM/ARCH 
    90 COMPILER=X64_ADA 
    91 export BATCH_COMMAND_PAR="llsubmit" 
     90COMPILER=ifort_athena 
     91export BATCH_COMMAND_PAR="bsub <" 
    9292export BATCH_COMMAND_SEQ=$BATCH_COMMAND_PAR 
    9393export INTERACT_FLAG="no" 
    9494export MPIRUN_FLAG="yes" 
    95 export USING_XIOS="yes" 
     95export USING_XIOS="no" 
    9696# 
    9797export DEL_KEYS="key_iomput" 
     
    142142# ORCA2_LIM_OBS:   15 
    143143# ORCA2_AGRIF_LIM :16 
    144 for config in  1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
    145  
     144for config in 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
    146145do 
    147146 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/SETTE/sette_xios.sh

    r4990 r5282  
    8484# 
    8585# Compiler among those in NEMOGCM/ARCH 
    86 COMPILER= 
    87 export BATCH_COMMAND_PAR="llsubmit" 
     86COMPILER=ifort_athena_xios 
     87export BATCH_COMMAND_PAR="bsub <" 
    8888export BATCH_COMMAND_SEQ=$BATCH_COMMAND_PAR 
    8989export INTERACT_FLAG="no" 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/TOOLS/COMPILE/tools.txt

    r2281 r5282  
    1 REBUILD  
     1REBUILD_NEMO  
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/TOOLS/NESTING/src/agrif_create_restart.f90

    r2455 r5282  
    411411        ENDIF 
    412412        ! 
    413      CASE('nfice','nfbulk','kt','ndastp','adatrj','rdt','rdttra1')  
     413     CASE('nfice','nfbulk','kt','ndastp','adatrj','rdt')  
    414414        IF(.NOT. dimg ) THEN 
    415415           WRITE(*,*) 'copy ',TRIM(Ncdf_varname(i)) 
     
    417417              CALL Read_Ncdf_var(TRIM(Ncdf_varname(i)),TRIM(restart_file),tabtemp0dreal)   
    418418              SELECT CASE (TRIM(Ncdf_varname(i))) 
    419               CASE('rdt','rdttra1') 
     419              CASE('rdt') 
    420420                 tabtemp0dreal = tabtemp0dreal/rhot 
    421421              CASE('kt') 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/TOOLS/OBSTOOLS/src/in_out_manager.F90

    r3000 r5282  
    3737   INTEGER       ::   nn_write      =   10        !: model standard output frequency 
    3838   INTEGER       ::   nn_stock      =   10        !: restart file frequency 
    39    LOGICAL       ::   ln_dimgnnn    = .FALSE.     !: type of dimgout. (F): 1 file for all proc 
    40                                                        !:                  (T): 1 file per proc 
    4139   LOGICAL       ::   ln_mskland    = .FALSE.     !: mask land points in NetCDF outputs (costly: + ~15%) 
    4240   LOGICAL       ::   ln_clobber    = .FALSE.     !: clobber (overwrite) an existing file 
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/TOOLS/SIREN/cfg/variable.cfg

    r5037 r5282  
    4141kt           |                |      |   |                                  |                    |                
    4242rdt          |                |      |   |                                  |                    |                
    43 rdttra1      |                |      |   |                                  |                    |                
    4443utau_b       |                | XYT  | U |surface_downward_eastward_stress  |                    |                
    4544vtau_b       |                | XYT  | V |surface_downward_northward_stress |                    |                
  • branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/TOOLS/SIREN/src/vgrid.f90

    r5037 r5282  
    229229            dd_gdepw(jk) = ( dl_zw - 1.0 ) * dl_za1 
    230230            dd_gdept(jk) = ( dl_zt - 1.0 ) * dl_za1 
    231             dd_e3w  (jk) =  dl_za1 
    232             dd_e3t  (jk) =  dl_za1 
     231            dd_e3w(jk) =  dl_za1 
     232            dd_e3t(jk) =  dl_za1 
    233233         END DO 
    234234 
     
    244244            &                dl_za1 * dl_zacr * LOG( COSH( (dl_zt-dl_zkth)/dl_zacr ) ) + & 
    245245            &                dl_za2 * dl_zacr2* LOG( COSH( (dl_zt-dl_zkth2)/dl_zacr2 ) ) ) 
    246             dd_e3w  (jk) =             dl_za0 + & 
     246            dd_e3w(jk) =             dl_za0 + & 
    247247            &                          dl_za1 * TANH(      (dl_zw-dl_zkth)/dl_zacr   ) + & 
    248248            &                          dl_za2 * TANH(      (dl_zw-dl_zkth2)/dl_zacr2 ) 
    249             dd_e3t  (jk) =             dl_za0 + & 
     249            dd_e3t(jk) =             dl_za0 + & 
    250250            &                          dl_za1 * TANH(      (dl_zt-dl_zkth)/dl_zacr   ) + & 
    251251            &                          dl_za2 * TANH(      (dl_zt-dl_zkth2)/dl_zacr2 ) 
Note: See TracChangeset for help on using the changeset viewer.