New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7403 – NEMO

Changeset 7403


Ignore:
Timestamp:
2016-11-30T17:56:53+01:00 (6 years ago)
Author:
timgraham
Message:

Merge dev_INGV_METO_merge_2016 into branch

Location:
branches/2016/dev_merge_2016/NEMOGCM
Files:
7 deleted
118 edited
14 copied

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_merge_2016/NEMOGCM/ARCH/INGV/arch-IBM_EKMAN_INGV.fcm

    r5656 r7403  
    3434%USER_INC            -I%XIOS_ROOT/inc %NCDF_INC %MPI_INTEL -I/srv/lib/zlib-last/include 
    3535%USER_LIB            -L%XIOS_ROOT/lib -lxios %NCDF_LIB -L/srv/lib/zlib-last/lib -lz 
     36%CC                  icc 
     37%CFLAGS              -O0 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_pisces_cfg

    r4147 r7403  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! GYRE_PISCES: Configuration namelist used to overwrite SHARED/namelist_pisces_ref 
     2!! PISCES  (key_pisces) reference namelist (see below for key_pisces_reduced) 
     3!!              1  - air-sea exchange                         (nampisext) 
     4!!              2  - biological parameters                    (nampisbio) 
     5!!              3  - parameters for nutrient limitations      (nampislim)     
     6!!              4  - parameters for phytoplankton             (nampisprod,nampismort) 
     7!!              5  - parameters for zooplankton               (nampismes,nampiszoo) 
     8!!              6  - parameters for remineralization          (nampisrem) 
     9!!              7  - parameters for calcite chemistry         (nampiscal) 
     10!!              8  - parameters for inputs deposition         (nampissed) 
     11!!              11 - Damping                                  (nampisdmp) 
     12!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     13!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     14&nampismod     !  Model used  
     15!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     16    ln_p4z = .false. 
     17    ln_p2z = .true. 
     18/ 
     19!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     20&nampisext     !   air-sea exchange 
     21!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     22/ 
     23!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     24&nampisatm     !  Atmospheric prrssure  
     25!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     26/ 
     27!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     28&nampisbio     !   biological parameters 
     29!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     30/ 
     31!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     32&nampislim     !   parameters for nutrient limitations 
     33!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     34/ 
     35!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     36&nampisopt     !   parameters for optics 
     37!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     38/  
     39!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     40&nampisprod     !   parameters for phytoplankton growth 
     41!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     42/ 
     43!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     44&nampismort     !   parameters for phytoplankton sinks 
     45!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     46/ 
     47!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     48&nampismes     !   parameters for mesozooplankton 
     49!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     50/ 
     51!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     52&nampiszoo     !   parameters for microzooplankton 
     53!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     54/ 
     55!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     56&nampisfer     !   parameters for iron chemistry 
     57!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     58 
     59!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     60&nampisrem     !   parameters for remineralization 
     61!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     62/ 
     63!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     64&nampiscal     !   parameters for Calcite chemistry 
     65!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     66/ 
     67!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     68&nampissbc     !   parameters for inputs deposition 
     69!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     70/ 
     71!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     72&nampisice      !  Prescribed sea ice tracers 
     73!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     74/ 
     75!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     76&nampisdmp     !  Damping  
     77!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     78/ 
     79!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     80&nampismass     !  Mass conservation 
     81!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     82/ 
     83!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     84!! PISCES reduced (key_pisces_reduced, ex LOBSTER) : namelists 
     85!!              1  - biological parameters for phytoplankton    (namlobphy) 
     86!!              2  - biological parameters for nutrients        (namlobnut) 
     87!!              3  - biological parameters for zooplankton      (namlobzoo)     
     88!!              4  - biological parameters for detritus         (namlobdet) 
     89!!              5  - biological parameters for DOM              (namlobdom) 
     90!!              6  - parameters from aphotic layers to sediment (namlobsed) 
     91!!              7  - general coefficients                       (namlobrat) 
     92!!              8  - optical parameters                         (namlobopt) 
     93!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    394!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    495&namlobphy     !   biological parameters for phytoplankton 
     
    15106!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    16107&namlobdet     !   biological parameters for detritus 
    17 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,        
     108!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    18109/ 
    19110!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_top_cfg

    r5836 r7403  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! NEMO-TOP Configuration namelist for GYRE_PISCES configuration used to overwrite SHARED/namelist_top_ref 
     2!! NEMO/TOP1 :  Configuration namelist : used to overwrite defaults values defined in SHARED/namelist_top_ref 
     3!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    34!----------------------------------------------------------------------- 
    45&namtrc_run     !   run information 
    56!----------------------------------------------------------------------- 
    6    nn_writetrc   =  60     !  time step frequency for sn_tracer outputs 
    77/ 
    88!----------------------------------------------------------------------- 
    99&namtrc     !   tracers definition 
    1010!----------------------------------------------------------------------- 
    11    ln_trcdta     =   .false.    !  Initialisation from data input file (T) or not (F) 
     11   jp_bgc        =  6 
    1212! 
    13 !                ! name  !     title of the field          !   units       ! initial data ! save   ! 
    14 !                !       !                                 !               ! from file    ! or not ! 
    15 !                !       !                                 !               ! or not       !        ! 
    16    sn_tracer(1)   = 'DET'   , 'Detritus                   ',  'mmole-N/m3' ,  .false.     ,  .false. 
    17    sn_tracer(2)   = 'ZOO'   , 'Zooplankton concentration  ',  'mmole-N/m3' ,  .false.     ,  .false. 
    18    sn_tracer(3)   = 'PHY'   , 'Phytoplankton concentration',  'mmole-N/m3' ,  .false.     ,  .false. 
    19    sn_tracer(4)   = 'NO3'   , 'Nitrate concentration      ',  'mmole-N/m3' ,  .false.     ,  .true. 
    20    sn_tracer(5)   = 'NH4'   , 'Ammonium concentration     ',  'mmole-N/m3' ,  .false.     ,  .false. 
    21    sn_tracer(6)   = 'DOM'   , 'Dissolved organic matter   ',  'mmole-N/m3' ,  .false.     ,  .false. 
     13   ln_pisces     =  .true. 
     14   ln_age        =  .false. 
     15   ln_cfc11      =  .false. 
     16   ln_cfc12      =  .false. 
     17   ln_c14        =  .false. 
     18   ln_my_trc     =  .false. 
     19! 
     20!                 !           !                             !               ! 
     21!                 !    name   !   title of the field        !   units       ! initial data from file or not ! 
     22!                 !           !                                             !                               ! 
     23   sn_tracer(1)   = 'DET'     , 'Detritus                   ',  'mmole-N/m3' ,  .false.     
     24   sn_tracer(2)   = 'ZOO'     , 'Zooplankton concentration  ',  'mmole-N/m3' ,  .false.    
     25   sn_tracer(3)   = 'PHY'     , 'Phytoplankton concentration',  'mmole-N/m3' ,  .false.   
     26   sn_tracer(4)   = 'NO3'     , 'Nitrate concentration      ',  'mmole-N/m3' ,  .false.  
     27   sn_tracer(5)   = 'NH4'     , 'Ammonium concentration     ',  'mmole-N/m3' ,  .false.    
     28   sn_tracer(6)   = 'DOM'     , 'Dissolved organic matter   ',  'mmole-N/m3' ,  .false.   
    2229/ 
    2330!----------------------------------------------------------------------- 
    24 &namtrc_adv    !   advection scheme for passive tracer  
    25 !-----------------------------------------------------------------------    
     31&namtrc_dta      !    Initialisation from data input file 
     32!----------------------------------------------------------------------- 
     33/ 
    2634   ln_trcadv_fct =  .true.   !  FCT scheme 
    27       nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
    28       nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     35      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order 
     36      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order 
    2937      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
    3038      !                             !        (number of sub-timestep = nn_fct_zts) 
    3139/ 
    3240!----------------------------------------------------------------------- 
    33 &namtrc_ldf    !   lateral diffusion scheme for passive tracer  
     41&namtrc_ldf    !   lateral diffusion scheme for passive tracer 
    3442!----------------------------------------------------------------------- 
    3543   rn_ahtrc_0       =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     
    4048/ 
    4149!----------------------------------------------------------------------- 
    42 &namtrc_rad        !  treatment of negative concentrations  
     50&namtrc_rad        !  treatment of negative concentrations 
    4351!----------------------------------------------------------------------- 
    4452   ln_trcrad   =  .false.  !  artificially correct negative concentrations (T) or not (F) 
    4553/ 
    4654!----------------------------------------------------------------------- 
    47 &namtrc_dia       !   parameters for passive tracer additional diagnostics 
     55&namtrc_dmp      !   passive tracer newtonian damping 
     56!----------------------------------------------------------------------- 
     57/ 
     58!----------------------------------------------------------------------- 
     59&namtrc_ice      !    Representation of sea ice growth & melt effects 
     60!----------------------------------------------------------------------- 
     61/ 
     62!----------------------------------------------------------------------- 
     63&namtrc_trd      !   diagnostics on tracer trends                       ('key_trdtrc') 
     64!                          or mixed-layer trends                        ('key_trdmld_trc') 
    4865!---------------------------------------------------------------------- 
    49    nn_writedia   =  60     !  time step frequency for diagnostics 
    5066/ 
    5167!---------------------------------------------------------------------- 
    52 &namtrc_bc        !   data for boundary conditions 
     68&namtrc_bc       !   data for boundary conditions 
    5369!----------------------------------------------------------------------- 
    5470/ 
     71!---------------------------------------------------------------------- 
     72&namtrc_bdy      !   Setup of tracer boundary conditions 
     73!----------------------------------------------------------------------- 
     74/ 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/GYRE_PISCES/cpp_GYRE_PISCES.fcm

    r5930 r7403  
    1 bld::tool::fppkeys key_zdftke key_top key_pisces_reduced key_mpp_mpi 
     1bld::tool::fppkeys key_zdftke key_top key_mpp_mpi 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/iodef.xml

    r5385 r7403  
    176176          <field field_ref="CO3sat"   /> 
    177177          <field field_ref="PAR"      /> 
    178           <field field_ref="PPPHY"    /> 
    179           <field field_ref="PPPHY2"   /> 
     178          <field field_ref="PPPHYN"    /> 
     179          <field field_ref="PPPHYD"   /> 
    180180          <field field_ref="PPNEWN"   /> 
    181181          <field field_ref="PPNEWD"   /> 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_pisces_cfg

    r4147 r7403  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! PISCES  :   Configuration namelist : used to overwrite defaults values defined in SHARED/namelist_pis_ref 
    3 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!! PISCES reference namelist  
     3!!              1  - air-sea exchange                         (nampisext) 
     4!!              2  - biological parameters                    (nampisbio) 
     5!!              3  - parameters for nutrient limitations      (nampislim)     
     6!!              4  - parameters for phytoplankton             (nampisprod,nampismort) 
     7!!              5  - parameters for zooplankton               (nampismes,nampiszoo) 
     8!!              6  - parameters for remineralization          (nampisrem) 
     9!!              7  - parameters for calcite chemistry         (nampiscal) 
     10!!              8  - parameters for inputs deposition         (nampissed) 
     11!!              11 - Damping                                  (nampisdmp) 
     12!----------------------------------------------------------------------- 
     13&nampismod     !  Model used  
     14!----------------------------------------------------------------------- 
     15/ 
     16!----------------------------------------------------------------------- 
     17&nampisext     !   air-sea exchange 
     18!----------------------------------------------------------------------- 
     19/ 
     20!----------------------------------------------------------------------- 
     21&nampisatm     !  Atmospheric prrssure  
     22!----------------------------------------------------------------------- 
     23/ 
     24!----------------------------------------------------------------------- 
     25&nampisbio     !   biological parameters 
     26!----------------------------------------------------------------------- 
     27/ 
     28!----------------------------------------------------------------------- 
     29&namp4zlim     !   parameters for nutrient limitations for PISCES std - ln_p4z 
     30!----------------------------------------------------------------------- 
     31/ 
     32!----------------------------------------------------------------------- 
     33&namp5zlim     !   parameters for nutrient limitations PISCES QUOTA - ln_p5z 
     34!----------------------------------------------------------------------- 
     35/ 
     36!----------------------------------------------------------------------- 
     37&namp5zquota    !   parameters for nutrient limitations PISCES quota - ln_p5z 
     38!----------------------------------------------------------------------- 
     39/ 
     40!----------------------------------------------------------------------- 
     41&nampisopt     !   parameters for optics 
     42!----------------------------------------------------------------------- 
     43/  
     44!----------------------------------------------------------------------- 
     45&namp4zprod     !   parameters for phytoplankton growth for PISCES std - ln_p4z 
     46!----------------------------------------------------------------------- 
     47/ 
    448!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    5 &nampisext     !   air-sea exchange 
     49&namp5zprod     !   parameters for phytoplankton growth for PISCES quota - ln_p5z 
    650!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    751/ 
    8 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    9 &nampisatm     !  Atmospheric prrssure  
    10 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     52!----------------------------------------------------------------------- 
     53&namp4zmort     !   parameters for phytoplankton sinks for PISCES std - ln_p4z 
     54!----------------------------------------------------------------------- 
    1155/ 
    12 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    13 &nampisbio     !   biological parameters 
    14 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     56!----------------------------------------------------------------------- 
     57&namp5zmort     !   parameters for phytoplankton sinks for PISCES quota - ln_p5z 
     58!----------------------------------------------------------------------- 
    1559/ 
    16 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    17 &nampislim     !   parameters for nutrient limitations 
    18 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     60!----------------------------------------------------------------------- 
     61&namp4zmes     !   parameters for mesozooplankton for PISCES std - ln_p4z 
     62!----------------------------------------------------------------------- 
    1963/ 
    20 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    21 &nampisopt     !   parameters for optics 
    22 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    23 /  
    24 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    25 &nampisprod     !   parameters for phytoplankton growth 
    26 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     64!----------------------------------------------------------------------- 
     65&namp5zmes     !   parameters for mesozooplankton 
     66!----------------------------------------------------------------------- 
    2767/ 
    28 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    29 &nampismort     !   parameters for phytoplankton sinks 
    30 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     68!----------------------------------------------------------------------- 
     69&namp4zzoo     !   parameters for microzooplankton for PISCES std - ln_p4z 
     70!----------------------------------------------------------------------- 
    3171/ 
    32 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    33 &nampismes     !   parameters for mesozooplankton 
    34 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     72!----------------------------------------------------------------------- 
     73&namp5zzoo     !   parameters for microzooplankton 
     74!----------------------------------------------------------------------- 
    3575/ 
    36 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    37 &nampiszoo     !   parameters for microzooplankton 
    38 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     76!----------------------------------------------------------------------- 
     77&nampisfer     !   parameters for iron chemistry 
     78!----------------------------------------------------------------------- 
    3979/ 
    40 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    41 &nampisfer     !   parameters for iron chemistry 
    42 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    43  
    44 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     80!-----------------------------------------------------------------------   
    4581&nampisrem     !   parameters for remineralization 
    46 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     82!----------------------------------------------------------------------- 
    4783/ 
    48 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     84!----------------------------------------------------------------------- 
     85&nampispoc     !   parameters for organic particles 
     86!----------------------------------------------------------------------- 
     87/ 
     88!----------------------------------------------------------------------- 
    4989&nampiscal     !   parameters for Calcite chemistry 
    50 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     90!----------------------------------------------------------------------- 
    5191/ 
    52 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     92!----------------------------------------------------------------------- 
    5393&nampissbc     !   parameters for inputs deposition 
    54 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     94!----------------------------------------------------------------------- 
    5595/ 
    56 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     96!----------------------------------------------------------------------- 
     97&nampislig      !  Namelist parameters for ligands, nampislig 
     98!----------------------------------------------------------------------- 
     99/ 
     100!----------------------------------------------------------------------- 
     101&nampisice      !  Prescribed sea ice tracers 
     102!----------------------------------------------------------------------- 
     103/ 
     104!----------------------------------------------------------------------- 
    57105&nampisdmp     !  Damping  
    58 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     106!----------------------------------------------------------------------- 
    59107/ 
    60 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     108!----------------------------------------------------------------------- 
    61109&nampismass     !  Mass conservation 
    62 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     110!----------------------------------------------------------------------- 
    63111/ 
     112!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     113!! PISCES reduced (key_pisces_reduced, ex LOBSTER) : namelists 
     114!!              1  - biological parameters for phytoplankton    (namlobphy) 
     115!!              2  - biological parameters for nutrients        (namlobnut) 
     116!!              3  - biological parameters for zooplankton      (namlobzoo)     
     117!!              4  - biological parameters for detritus         (namlobdet) 
     118!!              5  - biological parameters for DOM              (namlobdom) 
     119!!              6  - parameters from aphotic layers to sediment (namlobsed) 
     120!!              7  - general coefficients                       (namlobrat) 
     121!!              8  - optical parameters                         (namlobopt) 
     122!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     123!----------------------------------------------------------------------- 
     124&namlobphy     !   biological parameters for phytoplankton 
     125!----------------------------------------------------------------------- 
     126/ 
     127!----------------------------------------------------------------------- 
     128&namlobnut     !   biological parameters for nutrients 
     129!----------------------------------------------------------------------- 
     130/ 
     131!----------------------------------------------------------------------- 
     132&namlobzoo     !   biological parameters for zooplankton 
     133!----------------------------------------------------------------------- 
     134/ 
     135!----------------------------------------------------------------------- 
     136&namlobdet     !   biological parameters for detritus 
     137!----------------------------------------------------------------------- 
     138/ 
     139!----------------------------------------------------------------------- 
     140&namlobdom     !   biological parameters for DOM 
     141!----------------------------------------------------------------------- 
     142/ 
     143!----------------------------------------------------------------------- 
     144&namlobsed     !   parameters from aphotic layers to sediment 
     145!----------------------------------------------------------------------- 
     146/ 
     147!----------------------------------------------------------------------- 
     148&namlobrat     !   general coefficients 
     149!----------------------------------------------------------------------- 
     150/ 
     151!----------------------------------------------------------------------- 
     152&namlobopt     !   optical parameters 
     153!----------------------------------------------------------------------- 
     154/ 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_top_cfg

    r6140 r7403  
    55&namtrc_run     !   run information 
    66!----------------------------------------------------------------------- 
     7   ln_top_euler  = .true. 
    78/ 
    89!----------------------------------------------------------------------- 
    910&namtrc     !   tracers definition 
    1011!----------------------------------------------------------------------- 
    11 !                !    name   !           title of the field              ! initial data ! initial data ! save   ! 
    12 !                !           !                                           !  units       ! from file    ! or not !  
    13 !                !           !                                           !              ! or not       !        ! 
    14    sn_tracer(1)   = 'DIC     ' , 'Dissolved inorganic Concentration      ',  'mol-C/L' ,  .true.     ,  .true. 
    15    sn_tracer(2)   = 'Alkalini' , 'Total Alkalinity Concentration         ',  'eq/L '   ,  .true.     ,  .true. 
    16    sn_tracer(3)   = 'O2      ' , 'Dissolved Oxygen Concentration         ',  'mol-C/L' ,  .true.     ,  .true. 
    17    sn_tracer(4)   = 'CaCO3   ' , 'Calcite Concentration                  ',  'mol-C/L' ,  .false.    ,  .true. 
    18    sn_tracer(5)   = 'PO4     ' , 'Phosphate Concentration                ',  'mol-C/L' ,  .true.     ,  .true. 
    19    sn_tracer(6)   = 'POC     ' , 'Small organic carbon Concentration     ',  'mol-C/L' ,  .false.    ,  .true. 
    20    sn_tracer(7)   = 'Si      ' , 'Silicate Concentration                 ',  'mol-C/L' ,  .true.     ,  .true. 
    21    sn_tracer(8)   = 'PHY     ' , 'Nanophytoplankton Concentration        ',  'mol-C/L' ,  .false.    ,  .true. 
    22    sn_tracer(9)   = 'ZOO     ' , 'Microzooplankton Concentration         ',  'mol-C/L' ,  .false.    ,  .true. 
    23    sn_tracer(10)  = 'DOC     ' , 'Dissolved organic Concentration        ',  'mol-C/L' ,  .true.     ,  .true. 
    24    sn_tracer(11)  = 'PHY2    ' , 'Diatoms Concentration                  ',  'mol-C/L' ,  .false.    ,  .true. 
    25    sn_tracer(12)  = 'ZOO2    ' , 'Mesozooplankton Concentration          ',  'mol-C/L' ,  .false.    ,  .true. 
    26    sn_tracer(13)  = 'DSi     ' , 'Diatoms Silicate Concentration         ',  'mol-C/L' ,  .false.    ,  .true. 
    27    sn_tracer(14)  = 'Fer     ' , 'Dissolved Iron Concentration           ',  'mol-C/L' ,  .true.     ,  .true. 
    28    sn_tracer(15)  = 'BFe     ' , 'Big iron particles Concentration       ',  'mol-C/L' ,  .false.    ,  .true. 
    29    sn_tracer(16)  = 'GOC     ' , 'Big organic carbon Concentration       ',  'mol-C/L' ,  .false.    ,  .true. 
    30    sn_tracer(17)  = 'SFe     ' , 'Small iron particles Concentration     ',  'mol-C/L' ,  .false.    ,  .true. 
    31    sn_tracer(18)  = 'DFe     ' , 'Diatoms iron  Concentration            ',  'mol-C/L' ,  .false.    ,  .true. 
    32    sn_tracer(19)  = 'GSi     ' , 'Sinking biogenic Silicate Concentration',  'mol-C/L' ,  .false.    ,  .true. 
    33    sn_tracer(20)  = 'NFe     ' , 'Nano iron Concentration                ',  'mol-C/L' ,  .false.    ,  .true. 
    34    sn_tracer(21)  = 'NCHL    ' , 'Nano chlorophyl Concentration          ',  'mol-C/L' ,  .false.    ,  .true. 
    35    sn_tracer(22)  = 'DCHL    ' , 'Diatoms chlorophyl Concentration       ',  'mol-C/L' ,  .false.    ,  .true. 
    36    sn_tracer(23)  = 'NO3     ' , 'Nitrates Concentration                 ',  'mol-C/L' ,  .true.     ,  .true. 
    37    sn_tracer(24)  = 'NH4     ' , 'Ammonium Concentration                 ',  'mol-C/L' ,  .false.    ,  .true. 
     12   jp_bgc        =  24 
     13! 
     14   ln_pisces     =  .true. 
     15   ln_age        =  .false. 
     16   ln_cfc11      =  .false. 
     17   ln_cfc12      =  .false. 
     18   ln_c14        =  .false. 
     19   ln_my_trc     =  .false. 
     20! 
     21   ln_trcdta     =  .true.  !  Initialisation from data input file (T) or not (F) 
     22!                 !           !                                          !             ! 
     23!                 !    name   !           title of the field             !   units     ! initial data from file or not ! 
     24!                 !           !                                          !             ! 
     25   sn_tracer(1)   = 'DIC     ' , 'Dissolved inorganic Concentration      ',  'mol-C/L' ,  .true. 
     26   sn_tracer(2)   = 'Alkalini' , 'Total Alkalinity Concentration         ',  'eq/L '   ,  .true. 
     27   sn_tracer(3)   = 'O2      ' , 'Dissolved Oxygen Concentration         ',  'mol-C/L' ,  .true. 
     28   sn_tracer(4)   = 'CaCO3   ' , 'Calcite Concentration                  ',  'mol-C/L' ,  .false. 
     29   sn_tracer(5)   = 'PO4     ' , 'Phosphate Concentration                ',  'mol-C/L' ,  .true. 
     30   sn_tracer(6)   = 'POC     ' , 'Small organic carbon Concentration     ',  'mol-C/L' ,  .false. 
     31   sn_tracer(7)   = 'Si      ' , 'Silicate Concentration                 ',  'mol-C/L' ,  .true. 
     32   sn_tracer(8)   = 'PHY     ' , 'Nanophytoplankton Concentration        ',  'mol-C/L' ,  .false. 
     33   sn_tracer(9)   = 'ZOO     ' , 'Microzooplankton Concentration         ',  'mol-C/L' ,  .false. 
     34   sn_tracer(10)  = 'DOC     ' , 'Dissolved organic Concentration        ',  'mol-C/L' ,  .true. 
     35   sn_tracer(11)  = 'PHY2    ' , 'Diatoms Concentration                  ',  'mol-C/L' ,  .false. 
     36   sn_tracer(12)  = 'ZOO2    ' , 'Mesozooplankton Concentration          ',  'mol-C/L' ,  .false. 
     37   sn_tracer(13)  = 'DSi     ' , 'Diatoms Silicate Concentration         ',  'mol-C/L' ,  .false. 
     38   sn_tracer(14)  = 'Fer     ' , 'Dissolved Iron Concentration           ',  'mol-C/L' ,  .true. 
     39   sn_tracer(15)  = 'BFe     ' , 'Big iron particles Concentration       ',  'mol-C/L' ,  .false. 
     40   sn_tracer(16)  = 'GOC     ' , 'Big organic carbon Concentration       ',  'mol-C/L' ,  .false. 
     41   sn_tracer(17)  = 'SFe     ' , 'Small iron particles Concentration     ',  'mol-C/L' ,  .false. 
     42   sn_tracer(18)  = 'DFe     ' , 'Diatoms iron  Concentration            ',  'mol-C/L' ,  .false. 
     43   sn_tracer(19)  = 'GSi     ' , 'Sinking biogenic Silicate Concentration',  'mol-C/L' ,  .false. 
     44   sn_tracer(20)  = 'NFe     ' , 'Nano iron Concentration                ',  'mol-C/L' ,  .false. 
     45   sn_tracer(21)  = 'NCHL    ' , 'Nano chlorophyl Concentration          ',  'mol-C/L' ,  .false. 
     46   sn_tracer(22)  = 'DCHL    ' , 'Diatoms chlorophyl Concentration       ',  'mol-C/L' ,  .false. 
     47   sn_tracer(23)  = 'NO3     ' , 'Nitrates Concentration                 ',  'mol-C/L' ,  .true. 
     48   sn_tracer(24)  = 'NH4     ' , 'Ammonium Concentration                 ',  'mol-C/L' ,  .false. 
    3849/ 
    3950!----------------------------------------------------------------------- 
     
    7889/ 
    7990!----------------------------------------------------------------------- 
    80 &namtrc_dia       !   parameters for passive tracer additional diagnostics 
     91&namtrc_dmp      !   passive tracer newtonian damping 
     92!----------------------------------------------------------------------- 
     93/ 
     94!----------------------------------------------------------------------- 
     95&namtrc_ice      !    Representation of sea ice growth & melt effects 
     96!----------------------------------------------------------------------- 
     97/ 
     98!----------------------------------------------------------------------- 
     99&namtrc_trd      !   diagnostics on tracer trends                       ('key_trdtrc') 
     100!                          or mixed-layer trends                        ('key_trdmld_trc') 
    81101!---------------------------------------------------------------------- 
    82102/ 
    83103!---------------------------------------------------------------------- 
    84 &namtrc_bc        !   data for boundary conditions 
     104&namtrc_bc       !   data for boundary conditions 
    85105!----------------------------------------------------------------------- 
    86106/ 
     107!---------------------------------------------------------------------- 
     108&namtrc_bdy      !   Setup of tracer boundary conditions 
     109!----------------------------------------------------------------------- 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/cpp_ORCA2_LIM_PISCES.fcm

    r5930 r7403  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_zdftke key_zdfddm key_zdftmx key_top key_pisces key_mpp_mpi key_iomput 
     1bld::tool::fppkeys key_trabbl key_lim2 key_zdftke key_zdfddm key_zdftmx key_top key_mpp_mpi key_iomput 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/iodef.xml

    r5385 r7403  
    108108          <field field_ref="CO3sat"   /> 
    109109          <field field_ref="PAR"      /> 
    110           <field field_ref="PPPHY"    /> 
    111           <field field_ref="PPPHY2"   /> 
     110          <field field_ref="PPPHYN"    /> 
     111          <field field_ref="PPPHYD"   /> 
    112112          <field field_ref="PPNEWN"   /> 
    113113          <field field_ref="PPNEWD"   /> 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_pisces_cfg

    r4147 r7403  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! PISCES  :   ORCA2_OFF_PISCES configuration namelsit used to overwrite SHARED/namelist_pisces_ref 
    3 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!! PISCES reference namelist  
     3!!              1  - air-sea exchange                         (nampisext) 
     4!!              2  - biological parameters                    (nampisbio) 
     5!!              3  - parameters for nutrient limitations      (nampislim)     
     6!!              4  - parameters for phytoplankton             (nampisprod,nampismort) 
     7!!              5  - parameters for zooplankton               (nampismes,nampiszoo) 
     8!!              6  - parameters for remineralization          (nampisrem) 
     9!!              7  - parameters for calcite chemistry         (nampiscal) 
     10!!              8  - parameters for inputs deposition         (nampissed) 
     11!!              11 - Damping                                  (nampisdmp) 
     12!----------------------------------------------------------------------- 
     13&nampismod     !  Model used  
     14!----------------------------------------------------------------------- 
     15/ 
     16!----------------------------------------------------------------------- 
     17&nampisext     !   air-sea exchange 
     18!----------------------------------------------------------------------- 
     19/ 
     20!----------------------------------------------------------------------- 
     21&nampisatm     !  Atmospheric prrssure  
     22!----------------------------------------------------------------------- 
     23/ 
     24!----------------------------------------------------------------------- 
     25&nampisbio     !   biological parameters 
     26!----------------------------------------------------------------------- 
     27   nrdttrc    =  4        ! time step frequency for biology 
     28/ 
     29!----------------------------------------------------------------------- 
     30&namp4zlim     !   parameters for nutrient limitations for PISCES std - ln_p4z 
     31!----------------------------------------------------------------------- 
     32/ 
     33!----------------------------------------------------------------------- 
     34&namp5zlim     !   parameters for nutrient limitations PISCES QUOTA - ln_p5z 
     35!----------------------------------------------------------------------- 
     36/ 
     37!----------------------------------------------------------------------- 
     38&namp5zquota    !   parameters for nutrient limitations PISCES quota - ln_p5z 
     39!----------------------------------------------------------------------- 
     40/ 
     41!----------------------------------------------------------------------- 
     42&nampisopt     !   parameters for optics 
     43!----------------------------------------------------------------------- 
     44/  
     45!----------------------------------------------------------------------- 
     46&namp4zprod     !   parameters for phytoplankton growth for PISCES std - ln_p4z 
     47!----------------------------------------------------------------------- 
     48/ 
    449!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    5 &nampisext     !   air-sea exchange 
     50&namp5zprod     !   parameters for phytoplankton growth for PISCES quota - ln_p5z 
    651!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    752/ 
    8 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    9 &nampisatm     !  Atmospheric prrssure  
    10 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     53!----------------------------------------------------------------------- 
     54&namp4zmort     !   parameters for phytoplankton sinks for PISCES std - ln_p4z 
     55!----------------------------------------------------------------------- 
    1156/ 
    12 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    13 &nampisbio     !   biological parameters 
    14 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    15    nrdttrc    =  4        ! time step frequency for biology 
     57!----------------------------------------------------------------------- 
     58&namp5zmort     !   parameters for phytoplankton sinks for PISCES quota - ln_p5z 
     59!----------------------------------------------------------------------- 
    1660/ 
    17 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    18 &nampislim     !   parameters for nutrient limitations 
    19 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     61!----------------------------------------------------------------------- 
     62&namp4zmes     !   parameters for mesozooplankton for PISCES std - ln_p4z 
     63!----------------------------------------------------------------------- 
    2064/ 
    21 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    22 &nampisopt     !   parameters for optics 
    23 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    24 /  
    25 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    26 &nampisprod     !   parameters for phytoplankton growth 
    27 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     65!----------------------------------------------------------------------- 
     66&namp5zmes     !   parameters for mesozooplankton 
     67!----------------------------------------------------------------------- 
    2868/ 
    29 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    30 &nampismort     !   parameters for phytoplankton sinks 
    31 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     69!----------------------------------------------------------------------- 
     70&namp4zzoo     !   parameters for microzooplankton for PISCES std - ln_p4z 
     71!----------------------------------------------------------------------- 
    3272/ 
    33 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    34 &nampismes     !   parameters for mesozooplankton 
    35 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     73!----------------------------------------------------------------------- 
     74&namp5zzoo     !   parameters for microzooplankton 
     75!----------------------------------------------------------------------- 
    3676/ 
    37 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    38 &nampiszoo     !   parameters for microzooplankton 
    39 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     77!----------------------------------------------------------------------- 
     78&nampisfer     !   parameters for iron chemistry 
     79!----------------------------------------------------------------------- 
    4080/ 
    41 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    42 &nampisfer     !   parameters for iron chemistry 
    43 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    44  
    45 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     81!-----------------------------------------------------------------------   
    4682&nampisrem     !   parameters for remineralization 
    47 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     83!----------------------------------------------------------------------- 
    4884/ 
    49 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     85!----------------------------------------------------------------------- 
     86&nampispoc     !   parameters for organic particles 
     87!----------------------------------------------------------------------- 
     88/ 
     89!----------------------------------------------------------------------- 
    5090&nampiscal     !   parameters for Calcite chemistry 
    51 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     91!----------------------------------------------------------------------- 
    5292/ 
    53 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     93!----------------------------------------------------------------------- 
    5494&nampissbc     !   parameters for inputs deposition 
    55 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     95!----------------------------------------------------------------------- 
    5696/ 
    57 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     97!----------------------------------------------------------------------- 
     98&nampislig      !  Namelist parameters for ligands, nampislig 
     99!----------------------------------------------------------------------- 
     100/ 
     101!----------------------------------------------------------------------- 
     102&nampisice      !  Prescribed sea ice tracers 
     103!----------------------------------------------------------------------- 
     104/ 
     105!----------------------------------------------------------------------- 
    58106&nampisdmp     !  Damping  
    59 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    60    nn_pisdmp    =  1460       !  Frequency of Relaxation  
     107!----------------------------------------------------------------------- 
     108   nn_pisdmp    =  1460       !  Frequency of Relaxation 
    61109/ 
    62 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     110!----------------------------------------------------------------------- 
    63111&nampismass     !  Mass conservation 
    64 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     112!----------------------------------------------------------------------- 
    65113/ 
     114!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     115!! PISCES reduced (key_pisces_reduced, ex LOBSTER) : namelists 
     116!!              1  - biological parameters for phytoplankton    (namlobphy) 
     117!!              2  - biological parameters for nutrients        (namlobnut) 
     118!!              3  - biological parameters for zooplankton      (namlobzoo)     
     119!!              4  - biological parameters for detritus         (namlobdet) 
     120!!              5  - biological parameters for DOM              (namlobdom) 
     121!!              6  - parameters from aphotic layers to sediment (namlobsed) 
     122!!              7  - general coefficients                       (namlobrat) 
     123!!              8  - optical parameters                         (namlobopt) 
     124!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     125!----------------------------------------------------------------------- 
     126&namlobphy     !   biological parameters for phytoplankton 
     127!----------------------------------------------------------------------- 
     128/ 
     129!----------------------------------------------------------------------- 
     130&namlobnut     !   biological parameters for nutrients 
     131!----------------------------------------------------------------------- 
     132/ 
     133!----------------------------------------------------------------------- 
     134&namlobzoo     !   biological parameters for zooplankton 
     135!----------------------------------------------------------------------- 
     136/ 
     137!----------------------------------------------------------------------- 
     138&namlobdet     !   biological parameters for detritus 
     139!----------------------------------------------------------------------- 
     140/ 
     141!----------------------------------------------------------------------- 
     142&namlobdom     !   biological parameters for DOM 
     143!----------------------------------------------------------------------- 
     144/ 
     145!----------------------------------------------------------------------- 
     146&namlobsed     !   parameters from aphotic layers to sediment 
     147!----------------------------------------------------------------------- 
     148/ 
     149!----------------------------------------------------------------------- 
     150&namlobrat     !   general coefficients 
     151!----------------------------------------------------------------------- 
     152/ 
     153!----------------------------------------------------------------------- 
     154&namlobopt     !   optical parameters 
     155!----------------------------------------------------------------------- 
     156/ 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_top_cfg

    r6140 r7403  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! NEMO/TOP1 : ORCA2_OFF_PISCES configuration namelist used to overwrite SHARED/namelist_top 
     2!! NEMO/TOP1 :  Configuration namelist : used to overwrite defaults values defined in SHARED/namelist_top_ref 
    33!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    44!----------------------------------------------------------------------- 
    55&namtrc_run     !   run information 
    66!----------------------------------------------------------------------- 
    7    nn_writetrc   =  1460     !  time step frequency for sn_tracer outputs 
     7   ln_top_euler  = .true. 
    88/ 
    99!----------------------------------------------------------------------- 
    1010&namtrc     !   tracers definition 
    1111!----------------------------------------------------------------------- 
     12   jp_bgc        =  24 
    1213! 
    13 !              !    name   !           title of the field              !   units    ! initial data ! save   ! 
    14 !              !           !                                           !            ! from file    ! or not !  
    15 !              !           !                                           !            ! or not       !        ! 
    16    sn_tracer(1)   = 'DIC     ' , 'Dissolved inorganic Concentration      ',  'mol-C/L' ,  .true.     ,  .true. 
    17    sn_tracer(2)   = 'Alkalini' , 'Total Alkalinity Concentration         ',  'eq/L '   ,  .true.     ,  .true. 
    18    sn_tracer(3)   = 'O2      ' , 'Dissolved Oxygen Concentration         ',  'mol-C/L' ,  .true.     ,  .true. 
    19    sn_tracer(4)   = 'CaCO3   ' , 'Calcite Concentration                  ',  'mol-C/L' ,  .false.    ,  .true. 
    20    sn_tracer(5)   = 'PO4     ' , 'Phosphate Concentration                ',  'mol-C/L' ,  .true.     ,  .true. 
    21    sn_tracer(6)   = 'POC     ' , 'Small organic carbon Concentration     ',  'mol-C/L' ,  .false.    ,  .true. 
    22    sn_tracer(7)   = 'Si      ' , 'Silicate Concentration                 ',  'mol-C/L' ,  .true.     ,  .true. 
    23    sn_tracer(8)   = 'PHY     ' , 'Nanophytoplankton Concentration        ',  'mol-C/L' ,  .false.    ,  .true. 
    24    sn_tracer(9)   = 'ZOO     ' , 'Microzooplankton Concentration         ',  'mol-C/L' ,  .false.    ,  .true. 
    25    sn_tracer(10)  = 'DOC     ' , 'Dissolved organic Concentration        ',  'mol-C/L' ,  .false.    ,  .true. 
    26    sn_tracer(11)  = 'PHY2    ' , 'Diatoms Concentration                  ',  'mol-C/L' ,  .false.    ,  .true. 
    27    sn_tracer(12)  = 'ZOO2    ' , 'Mesozooplankton Concentration          ',  'mol-C/L' ,  .false.    ,  .true. 
    28    sn_tracer(13)  = 'DSi     ' , 'Diatoms Silicate Concentration         ',  'mol-C/L' ,  .false.    ,  .true. 
    29    sn_tracer(14)  = 'Fer     ' , 'Dissolved Iron Concentration           ',  'mol-C/L' ,  .true.     ,  .true. 
    30    sn_tracer(15)  = 'BFe     ' , 'Big iron particles Concentration       ',  'mol-C/L' ,  .false.    ,  .true. 
    31    sn_tracer(16)  = 'GOC     ' , 'Big organic carbon Concentration       ',  'mol-C/L' ,  .false.    ,  .true. 
    32    sn_tracer(17)  = 'SFe     ' , 'Small iron particles Concentration     ',  'mol-C/L' ,  .false.    ,  .true. 
    33    sn_tracer(18)  = 'DFe     ' , 'Diatoms iron  Concentration            ',  'mol-C/L' ,  .false.    ,  .true. 
    34    sn_tracer(19)  = 'GSi     ' , 'Sinking biogenic Silicate Concentration',  'mol-C/L' ,  .false.    ,  .true. 
    35    sn_tracer(20)  = 'NFe     ' , 'Nano iron Concentration                ',  'mol-C/L' ,  .false.    ,  .true. 
    36    sn_tracer(21)  = 'NCHL    ' , 'Nano chlorophyl Concentration          ',  'mol-C/L' ,  .false.    ,  .true. 
    37    sn_tracer(22)  = 'DCHL    ' , 'Diatoms chlorophyl Concentration       ',  'mol-C/L' ,  .false.    ,  .true. 
    38    sn_tracer(23)  = 'NO3     ' , 'Nitrates Concentration                 ',  'mol-C/L' ,  .true.     ,  .true. 
    39    sn_tracer(24)  = 'NH4     ' , 'Ammonium Concentration                 ',  'mol-C/L' ,  .false.    ,  .true. 
     14   ln_pisces     =  .true. 
     15   ln_my_trc     =  .false. 
     16   ln_age        =  .false. 
     17   ln_cfc11      =  .false. 
     18   ln_cfc12      =  .false. 
     19   ln_c14        =  .false. 
     20! 
     21   ln_trcdta     =  .true.  !  Initialisation from data input file (T) or not (F) 
     22!                 !           !                                          !             ! 
     23!                 !    name   !           title of the field             !   units     ! initial data from file or not ! 
     24!                 !           !                                          !             ! 
     25   sn_tracer(1)   = 'DIC     ' , 'Dissolved inorganic Concentration      ',  'mol-C/L' ,  .true. 
     26   sn_tracer(2)   = 'Alkalini' , 'Total Alkalinity Concentration         ',  'eq/L '   ,  .true. 
     27   sn_tracer(3)   = 'O2      ' , 'Dissolved Oxygen Concentration         ',  'mol-C/L' ,  .true. 
     28   sn_tracer(4)   = 'CaCO3   ' , 'Calcite Concentration                  ',  'mol-C/L' ,  .false. 
     29   sn_tracer(5)   = 'PO4     ' , 'Phosphate Concentration                ',  'mol-C/L' ,  .true. 
     30   sn_tracer(6)   = 'POC     ' , 'Small organic carbon Concentration     ',  'mol-C/L' ,  .false. 
     31   sn_tracer(7)   = 'Si      ' , 'Silicate Concentration                 ',  'mol-C/L' ,  .true. 
     32   sn_tracer(8)   = 'PHY     ' , 'Nanophytoplankton Concentration        ',  'mol-C/L' ,  .false. 
     33   sn_tracer(9)   = 'ZOO     ' , 'Microzooplankton Concentration         ',  'mol-C/L' ,  .false. 
     34   sn_tracer(10)  = 'DOC     ' , 'Dissolved organic Concentration        ',  'mol-C/L' ,  .true. 
     35   sn_tracer(11)  = 'PHY2    ' , 'Diatoms Concentration                  ',  'mol-C/L' ,  .false. 
     36   sn_tracer(12)  = 'ZOO2    ' , 'Mesozooplankton Concentration          ',  'mol-C/L' ,  .false. 
     37   sn_tracer(13)  = 'DSi     ' , 'Diatoms Silicate Concentration         ',  'mol-C/L' ,  .false. 
     38   sn_tracer(14)  = 'Fer     ' , 'Dissolved Iron Concentration           ',  'mol-C/L' ,  .true. 
     39   sn_tracer(15)  = 'BFe     ' , 'Big iron particles Concentration       ',  'mol-C/L' ,  .false. 
     40   sn_tracer(16)  = 'GOC     ' , 'Big organic carbon Concentration       ',  'mol-C/L' ,  .false. 
     41   sn_tracer(17)  = 'SFe     ' , 'Small iron particles Concentration     ',  'mol-C/L' ,  .false. 
     42   sn_tracer(18)  = 'DFe     ' , 'Diatoms iron  Concentration            ',  'mol-C/L' ,  .false. 
     43   sn_tracer(19)  = 'GSi     ' , 'Sinking biogenic Silicate Concentration',  'mol-C/L' ,  .false. 
     44   sn_tracer(20)  = 'NFe     ' , 'Nano iron Concentration                ',  'mol-C/L' ,  .false. 
     45   sn_tracer(21)  = 'NCHL    ' , 'Nano chlorophyl Concentration          ',  'mol-C/L' ,  .false. 
     46   sn_tracer(22)  = 'DCHL    ' , 'Diatoms chlorophyl Concentration       ',  'mol-C/L' ,  .false. 
     47   sn_tracer(23)  = 'NO3     ' , 'Nitrates Concentration                 ',  'mol-C/L' ,  .true. 
     48   sn_tracer(24)  = 'NH4     ' , 'Ammonium Concentration                 ',  'mol-C/L' ,  .false. 
    4049/ 
    4150!----------------------------------------------------------------------- 
    4251&namtrc_dta      !    Initialisation from data input file 
    4352!----------------------------------------------------------------------- 
    44 ! 
    4553!          !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    4654!          !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    47    sn_trcdta(1)  = 'data_DIC_nomask'        ,        -12        ,  'DIC'     ,    .false.   , .true. , 'yearly'  , ''       , ''  , '' 
    48    sn_trcdta(2)  = 'data_Alkalini_nomask'   ,        -12        ,  'Alkalini',    .false.   , .true. , 'yearly'  , ''       , ''  , '' 
    49    sn_trcdta(3)  = 'data_O2_nomask'      ,        -1         ,  'O2'      ,    .true.    , .true. , 'yearly'  , ''       , ''  , '' 
    50    sn_trcdta(5)  = 'data_PO4_nomask'     ,        -1         ,  'PO4'     ,    .true.    , .true. , 'yearly'  , ''       , ''  , '' 
    51    sn_trcdta(7)  = 'data_Si_nomask'      ,        -1         ,  'Si'      ,    .true.    , .true. , 'yearly'  , ''       , ''  , '' 
    52    sn_trcdta(14) = 'data_Fer_nomask'        ,        -12        ,  'Fer'     ,    .false.   , .true. , 'yearly'  , ''       , ''  , '' 
    53    sn_trcdta(23) = 'data_NO3_nomask'     ,        -1         ,  'NO3'     ,    .true.    , .true. , 'yearly'  , ''       , ''  , '' 
    54 ! 
    55    cn_dir        =  './'      !  root directory for the location of the data files 
     55   sn_trcdta(1)  = 'data_DIC_nomask'        ,        -12        ,  'DIC'     ,    .false.   , .true. , 'yearly'  , ''       , ''   , '' 
     56   sn_trcdta(2)  = 'data_Alkalini_nomask'   ,        -12        ,  'Alkalini',    .false.   , .true. , 'yearly'  , ''       , ''   , '' 
     57   sn_trcdta(3)  = 'data_O2_nomask'         ,        -1         ,  'O2'      ,    .true.    , .true. , 'yearly'  , ''       , ''   , '' 
     58   sn_trcdta(5)  = 'data_PO4_nomask'        ,        -1         ,  'PO4'     ,    .true.    , .true. , 'yearly'  , ''       , ''   , '' 
     59   sn_trcdta(7)  = 'data_Si_nomask'         ,        -1         ,  'Si'      ,    .true.    , .true. , 'yearly'  , ''       , ''   , '' 
     60   sn_trcdta(10) = 'data_DOC_nomask'        ,        -12        ,  'DOC'     ,    .false.   , .true. , 'yearly'  , ''       , ''   , '' 
     61   sn_trcdta(14) = 'data_Fer_nomask'        ,        -12        ,  'Fer'     ,    .false.   , .true. , 'yearly'  , ''       , ''   , '' 
     62   sn_trcdta(23) = 'data_NO3_nomask'        ,        -1         ,  'NO3'     ,    .true.    , .true. , 'yearly'  , ''       , ''   , '' 
    5663   rn_trfac(1)   =   1.0e-06  !  multiplicative factor 
    5764   rn_trfac(2)   =   1.0e-06  !  -      -      -     - 
     
    5966   rn_trfac(5)   = 122.0e-06  !  -      -      -     - 
    6067   rn_trfac(7)   =   1.0e-06  !  -      -      -     - 
     68   rn_trfac(10)  =   1.0      !  -      -      -     - 
    6169   rn_trfac(14)  =   1.0      !  -      -      -     - 
    6270   rn_trfac(23)  =   7.6e-06  !  -      -      -     - 
     
    8189/ 
    8290!----------------------------------------------------------------------- 
    83 &namtrc_dia       !   parameters for passive tracer additional diagnostics 
     91&namtrc_dmp      !   passive tracer newtonian damping 
     92!----------------------------------------------------------------------- 
     93/ 
     94!----------------------------------------------------------------------- 
     95&namtrc_ice      !    Representation of sea ice growth & melt effects 
     96!----------------------------------------------------------------------- 
     97/ 
     98!----------------------------------------------------------------------- 
     99&namtrc_trd      !   diagnostics on tracer trends                       ('key_trdtrc') 
     100!                          or mixed-layer trends                        ('key_trdmld_trc') 
    84101!---------------------------------------------------------------------- 
    85    nn_writedia   =  1460     !  time step frequency for diagnostics 
    86102/ 
    87103!---------------------------------------------------------------------- 
    88 &namtrc_bc        !   data for boundary conditions 
     104&namtrc_bc       !   data for boundary conditions 
    89105!----------------------------------------------------------------------- 
    90106/ 
     107!---------------------------------------------------------------------- 
     108&namtrc_bdy      !   Setup of tracer boundary conditions 
     109!----------------------------------------------------------------------- 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/cpp_ORCA2_OFF_PISCES.fcm

    r5836 r7403  
    1 bld::tool::fppkeys key_trabbl key_top key_offline key_pisces key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_trabbl key_top key_offline key_iomput key_mpp_mpi 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/SHARED/field_def.xml

    r6472 r7403  
    4141         <field id="sstgrad2"     long_name="square of module of sst gradient"                                                        unit="degC2/m2" /> 
    4242         <field id="sbt"          long_name="sea bottom temperature"                                                                  unit="degC"     /> 
     43         <field id="tosmint"      long_name="vertical integral of temperature times density"   standard_name="integral_wrt_depth_of_product_of_density_and_potential_temperature"  unit="(kg m2) degree_C" /> 
    4344         <field id="sst_wl"       long_name="Delta SST of warm layer"                                                                 unit="degC"     /> 
    4445         <field id="sst_cs"       long_name="Delta SST of cool skin"                                                                  unit="degC"     /> 
     
    4950         <field id="sssmax"       long_name="max of sea surface salinity"   field_ref="sss"   operation="maximum"                 /> 
    5051         <field id="sssmin"       long_name="min of sea surface salinity"   field_ref="sss"   operation="minimum"                 /> 
    51          <field id="sbs"          long_name="sea bottom salinity"                                                     unit="1e-3" /> 
     52         <field id="sbs"          long_name="sea bottom salinity"                                                     unit="0.001" /> 
     53         <field id="somint"       long_name="vertical integral of salinity times density"   standard_name="integral_wrt_depth_of_product_of_density_and_salinity"  unit="(kg m2) x (1e-3)" />  
    5254 
    5355         <field id="taubot"       long_name="bottom stress module"                                                    unit="N/m2" />  
     
    9395 
    9496         <!-- variables available with key_diaar5 --> 
    95          <field id="botpres"      long_name="Pressure at sea floor"   standard_name="sea_water_pressure_at_sea_floor"   unit="dbar" /> 
     97         <field id="botpres"      long_name="Sea Water Pressure at Sea Floor"   standard_name="sea_water_pressure_at_sea_floor"   unit="dbar" /> 
     98         <field id="sshdyn"       long_name="dynamic sea surface height"     standard_name="dynamic_sea_surface_height_above_geoid"     unit="m" /> 
     99         <field id="sshdyn2"      long_name="square of dynamic sea surface height"     standard_name="dynamic_sea_surface_height_above_geoid_squared"     unit="m2" > sshdyn * sshdyn </field> 
     100         <field id="tnpeo"      long_name="Tendency of ocean potential energy content"          unit="W/m2"                           /> 
    96101 
    97102         <!-- variables available with key_vvl --> 
     
    205210         <!-- * variable related to ice shelf forcing * --> 
    206211         <field id="fwfisf"       long_name="Ice shelf melting"                                            unit="Kg/m2/s"  /> 
    207          <field id="qisf"         long_name="Ice Shelf Heat Flux"                                          unit="W/m2"     /> 
     212         <field id="fwfisf3d"     long_name="Ice shelf melting"                             unit="kg/m2/s"  grid_ref="grid_T_3D" /> 
     213         <field id="qlatisf"      long_name="Ice shelf latent heat flux"                    unit="W/m2"     /> 
     214         <field id="qlatisf3d"    long_name="Ice shelf latent heat flux"                    unit="W/m2"     grid_ref="grid_T_3D" /> 
     215         <field id="qhcisf"       long_name="Ice shelf heat content flux"                   unit="W/m2"     /> 
     216         <field id="qhcisf3d"     long_name="Ice shelf heat content flux"                   unit="W/m2"     grid_ref="grid_T_3D" /> 
    208217         <field id="isfgammat"    long_name="transfert coefficient for isf (temperature) "                 unit="m/s"      /> 
    209218         <field id="isfgammas"    long_name="transfert coefficient for isf (salinity)    "                 unit="m/s"      /> 
     
    418427         <field id="utbl"         long_name="zonal current in the Losh tbl"     unit="m/s" /> 
    419428 
    420          <!-- variables available with key_diaar5 --> 
    421          <field id="u_masstr"     long_name="ocean eulerian mass transport along i-axis"    standard_name="ocean_mass_x_transport"                          unit="kg/s"        grid_ref="grid_U_3D" /> 
     429         <field id="u_masstr"     long_name="Ocean Mass X Transport"    standard_name="ocean_mass_x_transport"                          unit="kg/s"        grid_ref="grid_U_3D" /> 
     430         <field id="u_masstr_vint" long_name="vertical integral of ocean eulerian mass transport along i-axis"    standard_name="vertical_integral_of_ocean_mass_x_transport"  unit="kg/s" /> 
    422431         <field id="u_heattr"     long_name="ocean eulerian heat transport along i-axis"    standard_name="ocean_heat_x_transport"                          unit="W"                                /> 
    423432         <field id="u_salttr"     long_name="ocean eulerian salt transport along i-axis"    standard_name="ocean_salt_x_transport"                          unit="1e-3*kg/s"                        /> 
     433         <field id="uadv_heattr"  long_name="ocean advective heat transport along i-axis"    standard_name="advectice_ocean_heat_x_transport"               unit="W"                                /> 
     434         <field id="uadv_salttr"  long_name="ocean advective salt transport along i-axis"    standard_name="advectice_ocean_salt_x_transport"               unit="1e-3*kg/s"                      /> 
    424435         <field id="ueiv_heattr"  long_name="ocean bolus heat transport along i-axis"       standard_name="ocean_heat_x_transport_due_to_bolus_advection"   unit="W"                                /> 
     436         <field id="ueiv_salttr"  long_name="ocean bolus salt transport along i-axis"       standard_name="ocean_salt_x_transport_due_to_bolus_advection"   unit="Kg"                                /> 
     437         <field id="ueiv_heattr3d" long_name="ocean bolus heat transport along i-axis"    standard_name="ocean_heat_x_transport_due_to_bolus_advection"   unit="W"    grid_ref="grid_U_3D" /> 
     438         <field id="ueiv_salttr3d" long_name="ocean bolus salt transport along i-axis"    standard_name="ocean_salt_x_transport_due_to_bolus_advection"   unit="kg"   grid_ref="grid_U_3D" /> 
    425439         <field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis"   standard_name="ocean_heat_x_transport_due_to_diffusion"         unit="W"                                /> 
     440         <field id="udiff_salttr" long_name="ocean diffusion salt transport along i-axis"   standard_name="ocean_salt_x_transport_due_to_diffusion"         unit="1e-3*kg/s"                                /> 
    426441      </field_group> 
    427442       
     
    464479         <field id="v_heattr"     long_name="ocean eulerian heat transport along j-axis"    standard_name="ocean_heat_y_transport"                          unit="W"                                /> 
    465480         <field id="v_salttr"     long_name="ocean eulerian salt transport along i-axis"    standard_name="ocean_salt_y_transport"                          unit="1e-3*kg/s"                        /> 
     481         <field id="vadv_heattr"  long_name="ocean advective heat transport along j-axis"   standard_name="advectice_ocean_heat_y_transport"                unit="W"                      /> 
     482         <field id="vadv_salttr"  long_name="ocean advective salt transport along j-axis"   standard_name="advectice_ocean_salt_y_transport"                unit="1e-3*kg/s"              /> 
    466483         <field id="veiv_heattr"  long_name="ocean bolus heat transport along j-axis"       standard_name="ocean_heat_y_transport_due_to_bolus_advection"   unit="W"                                /> 
     484         <field id="veiv_salttr"  long_name="ocean bolus salt transport along j-axis"       standard_name="ocean_salt_x_transport_due_to_bolus_advection"   unit="Kg"                                /> 
     485         <field id="veiv_heattr3d" long_name="ocean bolus heat transport along j-axis"    standard_name="ocean_heat_y_transport_due_to_bolus_advection"   unit="W"    grid_ref="grid_V_3D" /> 
     486         <field id="veiv_salttr3d" long_name="ocean bolus salt transport along j-axis"    standard_name="ocean_salt_y_transport_due_to_bolus_advection"   unit="kg"   grid_ref="grid_V_3D" /> 
    467487         <field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis"   standard_name="ocean_heat_y_transport_due_to_diffusion"         unit="W"                                /> 
     488         <field id="vdiff_salttr" long_name="ocean diffusion salt transport along j-axis"   standard_name="ocean_salt_y_transport_due_to_diffusion"         unit="1e-3*kg/s"                        /> 
    468489      </field_group> 
    469490       
     
    625646 
    626647      <!-- Poleward transport : ptr -->      
    627       <field_group id="diaptr" domain_ref="ptr" >  
     648      <field_group id="diaptr" domain_ref="ptr" >   
    628649        <field id="zomsfglo"          long_name="Meridional Stream-Function: Global"           unit="Sv"       grid_ref="gznl_W_3D" /> 
    629650        <field id="zomsfatl"          long_name="Meridional Stream-Function: Atlantic"         unit="Sv"       grid_ref="gznl_W_3D" /> 
     
    631652        <field id="zomsfind"          long_name="Meridional Stream-Function: Indian"           unit="Sv"       grid_ref="gznl_W_3D" /> 
    632653        <field id="zomsfipc"          long_name="Meridional Stream-Function: Pacific+Indian"   unit="Sv"       grid_ref="gznl_W_3D" /> 
    633         <field id="zotemglo"          long_name="Zonal Mean Temperature : Global"              unit="degC"     grid_ref="gznl_T_3D" /> 
    634         <field id="zotematl"          long_name="Zonal Mean Temperature : Atlantic"            unit="degC"     grid_ref="gznl_T_3D" /> 
    635         <field id="zotempac"          long_name="Zonal Mean Temperature : Pacific"             unit="degC"     grid_ref="gznl_T_3D" /> 
    636         <field id="zotemind"          long_name="Zonal Mean Temperature : Indian"              unit="degC"     grid_ref="gznl_T_3D" /> 
    637         <field id="zotemipc"          long_name="Zonal Mean Temperature : Pacific+Indian"      unit="degC"     grid_ref="gznl_T_3D" /> 
    638         <field id="zosalglo"          long_name="Zonal Mean Salinity : Global"                 unit="1e-3"     grid_ref="gznl_T_3D" /> 
    639         <field id="zosalatl"          long_name="Zonal Mean Salinity : Atlantic"               unit="1e-3"     grid_ref="gznl_T_3D" /> 
    640         <field id="zosalpac"          long_name="Zonal Mean Salinity : Pacific"                unit="1e-3"     grid_ref="gznl_T_3D" /> 
    641         <field id="zosalind"          long_name="Zonal Mean Salinity : Indian"                 unit="1e-3"     grid_ref="gznl_T_3D" /> 
    642         <field id="zosalipc"          long_name="Zonal Mean Salinity : Pacific+Indian"         unit="1e-3"     grid_ref="gznl_T_3D" /> 
     654        <field id="zotemglo"          long_name="Zonal Mean Temperature : Global"              unit="degree_C"     grid_ref="gznl_T_3D" /> 
     655        <field id="zotematl"          long_name="Zonal Mean Temperature : Atlantic"            unit="degree_C"     grid_ref="gznl_T_3D" /> 
     656        <field id="zotempac"          long_name="Zonal Mean Temperature : Pacific"             unit="degree_C"     grid_ref="gznl_T_3D" /> 
     657        <field id="zotemind"          long_name="Zonal Mean Temperature : Indian"              unit="degree_C"     grid_ref="gznl_T_3D" /> 
     658        <field id="zotemipc"          long_name="Zonal Mean Temperature : Pacific+Indian"      unit="degree_C"     grid_ref="gznl_T_3D" /> 
     659        <field id="zosalglo"          long_name="Zonal Mean Salinity : Global"                 unit="0.001"     grid_ref="gznl_T_3D" /> 
     660        <field id="zosalatl"          long_name="Zonal Mean Salinity : Atlantic"               unit="0.001"     grid_ref="gznl_T_3D" /> 
     661        <field id="zosalpac"          long_name="Zonal Mean Salinity : Pacific"                unit="0.001"     grid_ref="gznl_T_3D" /> 
     662        <field id="zosalind"          long_name="Zonal Mean Salinity : Indian"                 unit="0.001"     grid_ref="gznl_T_3D" /> 
     663        <field id="zosalipc"          long_name="Zonal Mean Salinity : Pacific+Indian"         unit="0.001"     grid_ref="gznl_T_3D" /> 
    643664        <field id="zosrfglo"          long_name="Zonal Mean Surface"                           unit="m2"       grid_ref="gznl_T_3D" /> 
    644665        <field id="zosrfatl"          long_name="Zonal Mean Surface : Atlantic"                unit="m2"       grid_ref="gznl_T_3D" /> 
     
    647668        <field id="zosrfipc"          long_name="Zonal Mean Surface : Pacific+Indian"          unit="m2"       grid_ref="gznl_T_3D" /> 
    648669        <field id="sophtadv"          long_name="Advective Heat Transport"                     unit="PW"       grid_ref="gznl_T_2D" /> 
     670        <field id="sophtadv_atl"      long_name="Advective Heat Transport: Atlantic"           unit="PW"       grid_ref="gznl_T_2D" /> 
     671        <field id="sophtadv_pac"      long_name="Advective Heat Transport: Pacific"            unit="PW"       grid_ref="gznl_T_2D" /> 
     672        <field id="sophtadv_ind"      long_name="Advective Heat Transport: Indian"             unit="PW"       grid_ref="gznl_T_2D" /> 
     673        <field id="sophtadv_ipc"      long_name="Advective Heat Transport: Pacific+Indian"     unit="PW"       grid_ref="gznl_T_2D" /> 
    649674        <field id="sophtldf"          long_name="Diffusive Heat Transport"                     unit="PW"       grid_ref="gznl_T_2D" /> 
     675        <field id="sophtldf_atl"      long_name="Diffusive Heat Transport: Atlantic"           unit="PW"       grid_ref="gznl_T_2D" /> 
     676        <field id="sophtldf_pac"      long_name="Diffusive Heat Transport: Pacific"            unit="PW"       grid_ref="gznl_T_2D" /> 
     677        <field id="sophtldf_ind"      long_name="Diffusive Heat Transport: Indian"             unit="PW"       grid_ref="gznl_T_2D" /> 
     678        <field id="sophtldf_ipc"      long_name="Diffusive Heat Transport: Pacific+Indian"     unit="PW"       grid_ref="gznl_T_2D" /> 
     679        <field id="sophtove"          long_name="Overturning Heat Transport"                     unit="PW"       grid_ref="gznl_T_2D" /> 
     680        <field id="sophtove_atl"      long_name="Overturning Heat Transport: Atlantic"           unit="PW"       grid_ref="gznl_T_2D" /> 
     681        <field id="sophtove_pac"      long_name="Overturning Heat Transport: Pacific"            unit="PW"       grid_ref="gznl_T_2D" /> 
     682        <field id="sophtove_ind"      long_name="Overturning Heat Transport: Indian"             unit="PW"       grid_ref="gznl_T_2D" /> 
     683        <field id="sophtove_ipc"      long_name="Overturning Heat Transport: Pacific+Indian"     unit="PW"       grid_ref="gznl_T_2D" /> 
     684        <field id="sophtbtr"          long_name="Barotropic Heat Transport"                     unit="PW"       grid_ref="gznl_T_2D" /> 
     685        <field id="sophtbtr_atl"      long_name="Barotropic Heat Transport: Atlantic"           unit="PW"       grid_ref="gznl_T_2D" /> 
     686        <field id="sophtbtr_pac"      long_name="Barotropic Heat Transport: Pacific"            unit="PW"       grid_ref="gznl_T_2D" /> 
     687        <field id="sophtbtr_ind"      long_name="Barotropic Heat Transport: Indian"             unit="PW"       grid_ref="gznl_T_2D" /> 
     688        <field id="sophtbtr_ipc"      long_name="Barotropic Heat Transport: Pacific+Indian"     unit="PW"       grid_ref="gznl_T_2D" /> 
     689        <field id="sophteiv"          long_name="Heat Transport from mesoscale eddy advection"                     unit="PW"       grid_ref="gznl_T_2D" /> 
     690        <field id="sophteiv_atl"      long_name="Heat Transport from mesoscale eddy advection: Atlantic"           unit="PW"       grid_ref="gznl_T_2D" /> 
     691        <field id="sophteiv_pac"      long_name="Heat Transport from mesoscale eddy advection: Pacific"            unit="PW"       grid_ref="gznl_T_2D" /> 
     692        <field id="sophteiv_ind"      long_name="Heat Transport from mesoscale eddy advection: Indian"             unit="PW"       grid_ref="gznl_T_2D" /> 
     693        <field id="sophteiv_ipc"      long_name="Heat Transport from mesoscale eddy advection: Pacific+Indian"     unit="PW"       grid_ref="gznl_T_2D" /> 
    650694        <field id="sopstadv"          long_name="Advective Salt Transport"                     unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     695        <field id="sopstadv_atl"      long_name="Advective Salt Transport: Atlantic"           unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     696        <field id="sopstadv_pac"      long_name="Advective Salt Transport: Pacific"            unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     697        <field id="sopstadv_ind"      long_name="Advective Salt Transport: Indian"             unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     698        <field id="sopstadv_ipc"      long_name="Advective Salt Transport: Pacific+Indian"     unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     699        <field id="sopstove"          long_name="Overturning Salt Transport"                     unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     700        <field id="sopstove_atl"      long_name="Overturning Salt Transport: Atlantic"           unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     701        <field id="sopstove_pac"      long_name="Overturning Salt Transport: Pacific"            unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     702        <field id="sopstove_ind"      long_name="Overturning Salt Transport: Indian"             unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     703        <field id="sopstove_ipc"      long_name="Overturning Salt Transport: Pacific+Indian"     unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     704        <field id="sopstbtr"          long_name="Barotropic Salt Transport"                     unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     705        <field id="sopstbtr_atl"      long_name="Barotropic Salt Transport: Atlantic"           unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     706        <field id="sopstbtr_pac"      long_name="Barotropic Salt Transport: Pacific"            unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     707        <field id="sopstbtr_ind"      long_name="Barotropic Salt Transport: Indian"             unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     708        <field id="sopstbtr_ipc"      long_name="Barotropic Salt Transport: Pacific+Indian"     unit="Giga g/s" grid_ref="gznl_T_2D" /> 
    651709        <field id="sopstldf"          long_name="Diffusive Salt Transport"                     unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     710        <field id="sopstldf_atl"      long_name="Diffusive Salt Transport: Atlantic"           unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     711        <field id="sopstldf_pac"      long_name="Diffusive Salt Transport: Pacific"            unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     712        <field id="sopstldf_ind"      long_name="Diffusive Salt Transport: Indian"             unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     713        <field id="sopstldf_ipc"      long_name="Diffusive Salt Transport: Pacific+Indian"     unit="Giga g/s" grid_ref="gznl_T_2D" /> 
     714        <field id="sopsteiv"          long_name="Salt Transport from mesoscale eddy advection"                     unit="Giga g/s"       grid_ref="gznl_T_2D" /> 
     715        <field id="sopsteiv_atl"      long_name="Salt Transport from mesoscale eddy advection: Atlantic"           unit="Giga g/s"       grid_ref="gznl_T_2D" /> 
     716        <field id="sopsteiv_pac"      long_name="Salt Transport from mesoscale eddy advection: Pacific"            unit="Giga g/s"       grid_ref="gznl_T_2D" /> 
     717        <field id="sopsteiv_ind"      long_name="Salt Transport from mesoscale eddy advection: Indian"             unit="Giga g/s"       grid_ref="gznl_T_2D" /> 
     718        <field id="sopsteiv_ipc"      long_name="Salt Transport from mesoscale eddy advection: Pacific+Indian"     unit="Giga g/s"       grid_ref="gznl_T_2D" />        
    652719      </field_group> 
    653720 
     
    668735      <field id="ttrd_ad"       long_name="temperature-trend: advection"               standard_name="tendency_of_sea_water_temperature_due_to_advection"           unit="degC/s"                         > sqrt( ttrd_xad^2 + ttrd_yad^2 + ttrd_zad^2 ) </field> 
    669736      <field id="strd_ad"       long_name="salinity   -trend: advection"               standard_name="tendency_of_sea_water_salinity_due_to_advection"              unit="1e-3/s"                         > sqrt( strd_xad^2 + strd_yad^2 + strd_zad^2 ) </field> 
     737      <field id="ttrd_totad"    long_name="temperature-trend: total advection"         standard_name="tendency_of_sea_water_salinity_due_to_advection"              unit="degC/s"                        /> 
     738      <field id="strd_totad"    long_name="salinity   -trend: total advection"         standard_name="tendency_of_sea_water_salinity_due_to_advection"              unit="1e-3/s"                        /> 
    670739      <field id="ttrd_sad"      long_name="temperature-trend: surface adv. (no-vvl)"                                                                                unit="degC/s"   grid_ref="grid_T_2D" /> 
    671740      <field id="strd_sad"      long_name="salinity   -trend: surface adv. (no-vvl)"                                                                                unit="1e-3/s"   grid_ref="grid_T_2D" /> 
     
    674743      <field id="ttrd_zdf"      long_name="temperature-trend: vertical diffusion"      standard_name="tendency_of_sea_water_temperature_due_to_vertical_mixing"     unit="degC/s"                        /> 
    675744      <field id="strd_zdf"      long_name="salinity   -trend: vertical diffusion"      standard_name="tendency_of_sea_water_salinity_due_to_vertical_mixing"        unit="1e-3/s"                        /> 
     745      <field id="ttrd_evd"      long_name="temperature-trend: EVD convection"                                                                                       unit="degC/s"                        /> 
     746      <field id="strd_evd"      long_name="salinity   -trend: EVD convection"                                                                                       unit="1e-3/s"                        /> 
    676747 
    677748      <!-- ln_traldf_iso=T only (iso-neutral diffusion) --> 
     749      <field id="ttrd_iso"      long_name="temperature-trend: isopycnal diffusion"                             unit="degC/s" > ttrd_ldf + ttrd_zdf - ttrd_zdfp </field> 
     750      <field id="strd_iso"      long_name="salinity   -trend: isopycnal diffusion"                             unit="1e-3/s" > strd_ldf + strd_zdf - strd_zdfp </field> 
    678751      <field id="ttrd_zdfp"     long_name="temperature-trend: pure vert. diffusion"   unit="degC/s" /> 
    679752      <field id="strd_zdfp"     long_name="salinity   -trend: pure vert. diffusion"   unit="1e-3/s" /> 
     
    692765      <field id="ttrd_atf"      long_name="temperature-trend: asselin time filter"       unit="degC/s" /> 
    693766      <field id="strd_atf"      long_name="salinity   -trend: asselin time filter"       unit="1e-3/s" /> 
     767      <field id="ttrd_tot"      long_name="temperature-trend: total model trend"         unit="degC/s" /> 
     768      <field id="strd_tot"      long_name="salinity   -trend: total model trend"         unit="1e-3/s" /> 
     769 
     770      <!-- Thickness weighted versions: --> 
     771      <field id="ttrd_xad_e3t"      unit="degC/s * m" >  ttrd_xad * e3t </field> 
     772      <field id="strd_xad_e3t"      unit="1e-3/s * m" >  strd_xad * e3t </field> 
     773      <field id="ttrd_yad_e3t"      unit="degC/s * m" >  ttrd_yad * e3t </field> 
     774      <field id="strd_yad_e3t"      unit="1e-3/s * m" >  strd_yad * e3t </field> 
     775      <field id="ttrd_zad_e3t"      unit="degC/s * m" >  ttrd_zad * e3t </field> 
     776      <field id="strd_zad_e3t"      unit="1e-3/s * m" >  strd_zad * e3t </field> 
     777      <field id="ttrd_ad_e3t"       unit="degC/s * m" >  ttrd_ad  * e3t </field> 
     778      <field id="strd_ad_e3t"       unit="1e-3/s * m" >  strd_ad  * e3t </field> 
     779      <field id="ttrd_totad_e3t"    unit="degC/s * m" >  ttrd_totad  * e3t </field> 
     780      <field id="strd_totad_e3t"    unit="1e-3/s * m" >  strd_totad  * e3t </field> 
     781      <field id="ttrd_ldf_e3t"      unit="degC/s * m" >  ttrd_ldf * e3t </field> 
     782      <field id="strd_ldf_e3t"      unit="1e-3/s * m" >  strd_ldf * e3t </field> 
     783      <field id="ttrd_zdf_e3t"      unit="degC/s * m" >  ttrd_zdf * e3t </field> 
     784      <field id="strd_zdf_e3t"      unit="1e-3/s * m" >  strd_zdf * e3t </field> 
     785      <field id="ttrd_evd_e3t"      unit="degC/s * m" >  ttrd_evd * e3t </field> 
     786      <field id="strd_evd_e3t"      unit="1e-3/s * m" >  strd_evd * e3t </field> 
     787 
     788      <!-- ln_traldf_iso=T only (iso-neutral diffusion) --> 
     789      <field id="ttrd_iso_e3t"      unit="degC/s * m"  >  ttrd_iso * e3t </field> 
     790      <field id="strd_iso_e3t"      unit="1e-3/s * m"  >  strd_iso * e3t </field> 
     791      <field id="ttrd_zdfp_e3t"     unit="degC/s * m"  >  ttrd_zdfp * e3t </field> 
     792      <field id="strd_zdfp_e3t"     unit="1e-3/s * m"  >  strd_zdfp * e3t </field> 
     793 
     794      <!-- --> 
     795      <field id="ttrd_dmp_e3t"      unit="degC/s * m"  >  ttrd_dmp * e3t </field> 
     796      <field id="strd_dmp_e3t"      unit="1e-3/s * m"  >  strd_dmp * e3t </field> 
     797      <field id="ttrd_bbl_e3t"      unit="degC/s * m"  >  ttrd_bbl * e3t </field> 
     798      <field id="strd_bbl_e3t"      unit="1e-3/s * m"  >  strd_bbl * e3t </field> 
     799      <field id="ttrd_npc_e3t"      unit="degC/s * m"  >  ttrd_npc * e3t </field> 
     800      <field id="strd_npc_e3t"      unit="1e-3/s * m"  >  strd_npc * e3t </field> 
     801      <field id="ttrd_qns_e3t"      unit="degC/s * m"  >  ttrd_qns * e3t_surf </field> 
     802      <field id="strd_cdt_e3t"      unit="degC/s * m"  >  strd_cdt * e3t_surf </field> 
     803      <field id="ttrd_qsr_e3t"      unit="degC/s * m"  >  ttrd_qsr * e3t </field> 
     804      <field id="ttrd_bbc_e3t"      unit="degC/s * m"  >  ttrd_bbc * e3t </field> 
     805      <field id="ttrd_atf_e3t"      unit="degC/s * m"  >  ttrd_atf * e3t </field> 
     806      <field id="strd_atf_e3t"      unit="1e-3/s * m"  >  strd_atf * e3t </field> 
     807      <field id="ttrd_tot_e3t"      unit="degC/s * m"  >  ttrd_tot * e3t </field> 
     808      <field id="strd_tot_e3t"      unit="1e-3/s * m"  >  strd_tot * e3t </field> 
     809 
    694810 
    695811      <!-- variables available with ln_KE_trd --> 
     
    785901 
    786902     <field_group id="ptrc_T" grid_ref="grid_T_3D"> 
     903       <!-- PISCES standard : variables available with ln_p4z  --> 
    787904       <field id="DIC"          long_name="Dissolved inorganic Concentration"        unit="mmol/m3" /> 
    788        <field id="DIC_E3T"      long_name="DIC * E3T"                                unit="mmol/m2" > DIC * e3t </field > 
     905       <field id="DIC_e3t"      long_name="DIC * e3t"                                unit="mmol/m2" > DIC * e3t </field > 
    789906       <field id="Alkalini"     long_name="Total Alkalinity Concentration"           unit="mmol/m3" /> 
    790        <field id="Alkalini_E3T" long_name="Alkalini * E3T"                           unit="mmol/m2"  > Alkalini * e3t </field > 
     907       <field id="Alkalini_e3t" long_name="Alkalini * e3t"                           unit="mmol/m2"  > Alkalini * e3t </field > 
    791908       <field id="O2"           long_name="Oxygen Concentration"                     unit="mmol/m3" /> 
    792        <field id="O2_E3T"       long_name="O2 * E3T"                                 unit="mmol/m2"  > O2 * e3t </field > 
     909       <field id="O2_e3t"       long_name="O2 * e3t"                                 unit="mmol/m2"  > O2 * e3t </field > 
    793910       <field id="CaCO3"        long_name="Calcite Concentration"                    unit="mmol/m3" /> 
    794        <field id="CaCO3_E3T"    long_name="CaCO3 * E3T"                              unit="mmol/m2"  > CaCO3 * e3t </field > 
     911       <field id="CaCO3_e3t"    long_name="CaCO3 * e3t"                              unit="mmol/m2"  > CaCO3 * e3t </field > 
    795912       <field id="PO4"          long_name="Phosphate Concentration"                  unit="mmol/m3" /> 
    796        <field id="PO4_E3T"      long_name="PO4 * E3T"                                unit="mmol/m2"  > PO4 * e3t </field > 
     913       <field id="PO4_e3t"      long_name="PO4 * e3t"                                unit="mmol/m2"  > PO4 * e3t </field > 
    797914       <field id="POC"          long_name="Small organic carbon Concentration"       unit="mmol/m3" /> 
    798        <field id="POC_E3T"      long_name="POC * E3T"                                unit="mmol/m2"  > POC * e3t </field > 
     915       <field id="POC_e3t"      long_name="POC * e3t"                                unit="mmol/m2"  > POC * e3t </field > 
    799916       <field id="Si"           long_name="Silicate Concentration"                   unit="mmol/m3" /> 
    800        <field id="Si_E3T"       long_name="Si * E3T"                                 unit="mmol/m2"  > Si * e3t </field > 
     917       <field id="Si_e3t"       long_name="Si * e3t"                                 unit="mmol/m2"  > Si * e3t </field > 
    801918       <field id="PHY"          long_name="(Nano)Phytoplankton Concentration"        unit="mmol/m3" /> 
    802        <field id="PHY_E3T"      long_name="PHY * E3T"                                unit="mmol/m2"  > PHY * e3t </field > 
     919       <field id="PHY_e3t"      long_name="PHY * e3t"                                unit="mmol/m2"  > PHY * e3t </field > 
    803920       <field id="ZOO"          long_name="(Micro)Zooplankton Concentration"         unit="mmol/m3" /> 
    804        <field id="ZOO_E3T"      long_name="ZOO2 * E3T"                               unit="mmol/m2"  > ZOO * e3t </field > 
     921       <field id="ZOO_e3t"      long_name="ZOO2 * e3t"                               unit="mmol/m2"  > ZOO * e3t </field > 
    805922       <field id="DOC"          long_name="Dissolved organic Concentration"          unit="mmol/m3" /> 
    806        <field id="DOC_E3T"      long_name="DOC * E3T"                                unit="mmol/m2"  > DOC * e3t </field > 
     923       <field id="DOC_e3t"      long_name="DOC * e3t"                                unit="mmol/m2"  > DOC * e3t </field > 
    807924       <field id="PHY2"         long_name="Diatoms Concentration"                    unit="mmol/m3" /> 
    808        <field id="PHY2_E3T"     long_name="PHY2 * E3T"                               unit="mmol/m2"  > PHY2 * e3t </field > 
     925       <field id="PHY2_e3t"     long_name="PHY2 * e3t"                               unit="mmol/m2"  > PHY2 * e3t </field > 
    809926       <field id="ZOO2"         long_name="Mesozooplankton Concentration"            unit="mmol/m3" /> 
    810        <field id="ZOO2_E3T"     long_name="ZOO2 * E3T"                               unit="mmol/m2"  > ZOO2 * e3t </field > 
     927       <field id="ZOO2_e3t"     long_name="ZOO2 * e3t"                               unit="mmol/m2"  > ZOO2 * e3t </field > 
    811928       <field id="DSi"          long_name="Diatoms Silicate Concentration"           unit="mmol/m3" /> 
    812        <field id="DSi_E3T"      long_name="Dsi * E3T"                                unit="mmol/m2"  > DSi * e3t </field > 
     929       <field id="DSi_e3t"      long_name="Dsi * e3t"                                unit="mmol/m2"  > DSi * e3t </field > 
    813930       <field id="Fer"          long_name="Dissolved Iron Concentration"             unit="mmol/m3" /> 
    814        <field id="Fer_E3T"      long_name="Fer * E3T"                                unit="mmol/m2"  > Fer * e3t </field > 
     931       <field id="Fer_e3t"      long_name="Fer * e3t"                                unit="mmol/m2"  > Fer * e3t </field > 
    815932       <field id="BFe"          long_name="Big iron particles Concentration"         unit="mmol/m3" /> 
    816        <field id="BFe_E3T"      long_name="BFe * E3T"                                unit="mmol/m2"  > BFe * e3t </field > 
     933       <field id="BFe_e3t"      long_name="BFe * e3t"                                unit="mmol/m2"  > BFe * e3t </field > 
    817934       <field id="GOC"          long_name="Big organic carbon Concentration"         unit="mmol/m3" /> 
    818        <field id="GOC_E3T"      long_name="GOC * E3T"                                unit="mmol/m2"  > GOC * e3t </field > 
     935       <field id="GOC_e3t"      long_name="GOC * e3t"                                unit="mmol/m2"  > GOC * e3t </field > 
    819936       <field id="SFe"          long_name="Small iron particles Concentration"       unit="mmol/m3" /> 
    820        <field id="SFe_E3T"      long_name="SFe * E3T"                                unit="mmol/m2"  > SFe * e3t </field > 
     937       <field id="SFe_e3t"      long_name="SFe * e3t"                                unit="mmol/m2"  > SFe * e3t </field > 
    821938       <field id="DFe"          long_name="Diatoms iron  Concentration"              unit="mmol/m3" /> 
    822        <field id="DFe_E3T"      long_name="DFe * E3T"                                unit="mmol/m2"  > DFe * e3t </field > 
     939       <field id="DFe_e3t"      long_name="DFe * e3t"                                unit="mmol/m2"  > DFe * e3t </field > 
    823940       <field id="GSi"          long_name="Sinking biogenic Silicate Concentration"  unit="mmol/m3" /> 
    824        <field id="GSi_E3T"      long_name="GSi * E3T"                                unit="mmol/m2"  > GSi * e3t </field > 
     941       <field id="GSi_e3t"      long_name="GSi * e3t"                                unit="mmol/m2"  > GSi * e3t </field > 
    825942       <field id="NFe"          long_name="Nano iron Concentration"                  unit="mmol/m3" /> 
    826        <field id="NFe_E3T"      long_name="NFe * E3T"                                unit="mmol/m2"  > NFe * e3t </field > 
     943       <field id="NFe_e3t"      long_name="NFe * e3t"                                unit="mmol/m2"  > NFe * e3t </field > 
    827944       <field id="NCHL"         long_name="Nano chlorophyl Concentration"            unit="mg/m3"   /> 
    828        <field id="NCHL_E3T"     long_name="NCHL * E3T"                               unit="mmol/m2"  > NCHL * e3t </field > 
     945       <field id="NCHL_e3t"     long_name="NCHL * e3t"                               unit="mmol/m2"  > NCHL * e3t </field > 
    829946       <field id="DCHL"         long_name="Diatoms chlorophyl Concentration"         unit="mg/m3"   /> 
    830        <field id="DCHL_E3T"     long_name="DCHL * E3T"                               unit="mmol/m2"  > DCHL * e3t </field > 
     947       <field id="DCHL_e3t"     long_name="DCHL * e3t"                               unit="mmol/m2"  > DCHL * e3t </field > 
    831948       <field id="NO3"          long_name="Nitrate Concentration"                    unit="mmol/m3" /> 
    832        <field id="NO3_E3T"      long_name="NO3 * E3T"                                unit="mmol/m2"  > NO3 * e3t </field > 
     949       <field id="NO3_e3t"      long_name="NO3 * e3t"                                unit="mmol/m2"  > NO3 * e3t </field > 
    833950       <field id="NH4"          long_name="Ammonium Concentration"                   unit="mmol/m3" /> 
    834        <field id="NH4_E3T"      long_name="NH4 * E3T"                                unit="mmol/m2"  > NH4 * e3t </field > 
    835  
    836        <!-- PISCES with Kriest parametisation : variables available with key_kriest --> 
    837        <field id="Num"         long_name="Number of organic particles"              unit="1" /> 
    838        <field id="Num_E3T"     long_name="Num * E3T"                                unit="m"  > Num * e3t </field > 
    839  
    840        <!-- PISCES light : variables available with key_pisces_reduced --> 
     951       <field id="NH4_e3t"      long_name="NH4 * e3t"                                unit="mmol/m2"  > NH4 * e3t </field > 
     952 
     953       <!-- PISCES quota : variables available with ln_p5z  --> 
     954 
     955       <field id="DON"          long_name="Dissolved organic N Concentration"        unit="mmol/m3" /> 
     956       <field id="DON_e3t"      long_name="DON * e3t"                                unit="mmol/m2"  > DON * e3t </field > 
     957       <field id="DOP"          long_name="Dissolved organic P Concentration"        unit="mmol/m3" /> 
     958       <field id="DOP_e3t"      long_name="DOP * e3t"                                unit="mmol/m2"  > DOP * e3t </field > 
     959       <field id="PON"          long_name="Small PON Concentration"                  unit="mmol/m3" /> 
     960       <field id="PON_e3t"      long_name="PON * e3t"                                unit="mmol/m2"  > PON * e3t </field > 
     961       <field id="POP"          long_name="Small POP Concentration"                  unit="mmol/m3" /> 
     962       <field id="POP_e3t"      long_name="POP * e3t"                                unit="mmol/m2"  > POP * e3t </field > 
     963       <field id="GON"          long_name="Big PON Concentration"                    unit="mmol/m3" /> 
     964       <field id="GON_e3t"      long_name="GON * e3t"                                unit="mmol/m2"  > GON * e3t </field > 
     965       <field id="GOP"          long_name="Big POP Concentration"                    unit="mmol/m3" /> 
     966       <field id="GOP_e3t"      long_name="GOP * e3t"                                unit="mmol/m2"  > GOP * e3t </field > 
     967       <field id="PHYN"         long_name="Nanophytoplankton N biomass"              unit="mmol/m3" /> 
     968       <field id="PHYN_e3t"     long_name="PHYN * e3t"                               unit="mmol/m2"  > PHYN * e3t </field > 
     969       <field id="PHYP"         long_name="Nanophytoplankton P biomass"              unit="mmol/m3" /> 
     970       <field id="PHYP_e3t"     long_name="PHYP * e3t"                               unit="mmol/m2"  > PHYP * e3t </field > 
     971       <field id="DIAN"         long_name="Diatoms N biomass"                        unit="mmol/m3" /> 
     972       <field id="DIAN_e3t"     long_name="DIAN * e3t"                               unit="mmol/m2"  > DIAN * e3t </field > 
     973       <field id="DIAP"         long_name="Diatoms P biomass"                        unit="mmol/m3" /> 
     974       <field id="DIAP_e3t"     long_name="DIAP * e3t"                               unit="mmol/m2"  > DIAP * e3t </field > 
     975       <field id="PIC"          long_name="Picophytoplankton C biomass"              unit="mmol/m3" /> 
     976       <field id="PIC_e3t"      long_name="PIC * e3t"                                unit="mmol/m2"  > PIC * e3t </field > 
     977       <field id="PICN"         long_name="Picophytoplankton N biomass"              unit="mmol/m3" /> 
     978       <field id="PICN_e3t"     long_name="PICN * e3t"                               unit="mmol/m2"  > PICN * e3t </field > 
     979       <field id="PICP"         long_name="Picophytoplankton P biomass"              unit="mmol/m3" /> 
     980       <field id="PICP_e3t"     long_name="PICP * e3t"                               unit="mmol/m2"  > PICP * e3t </field > 
     981       <field id="PFe"          long_name="Picophytoplankton Fe biomass"             unit="mmol/m3" /> 
     982       <field id="PFe_e3t"      long_name="PFe * e3t"                                unit="mmol/m2"  > PFe * e3t </field > 
     983       <field id="PCHL"         long_name="Picophytoplankton Chl biomass"            unit="mg/m3" /> 
     984       <field id="PCHL_e3t"     long_name="PCHL * e3t"                               unit="mmol/m2"  > PCHL * e3t </field > 
     985 
     986      <!-- PISCES with ligand parametisation : variables available namelist paramter ln_ligand --> 
     987       <field id="LGW"         long_name="Weak ligands concentration"                unit="mmol/m3" /> 
     988       <field id="LGW_e3t"     long_name="LGW * e3t"                                 unit="mmol/m2"  > LGW * e3t </field > 
     989       <field id="LFe"         long_name="Lithogenic iron concentration"             unit="mmol/m3" /> 
     990       <field id="LFe_e3t"     long_name="LFe * e3t"                                 unit="mmol/m2"  > LFe * e3t </field > 
     991 
     992       <!-- PISCES light : variables available with ln_p2z  --> 
    841993       <field id="DET"         long_name="Detritus"                                 unit="mmol-N/m3" /> 
    842        <field id="DET_E3T"     long_name="DET * E3T"                                unit="mmol-N/m2"  > DET * e3t </field > 
     994       <field id="DET_e3t"     long_name="DET * e3t"                                unit="mmol-N/m2"  > DET * e3t </field > 
    843995       <field id="DOM"         long_name="Dissolved Organic Matter"                 unit="mmol-N/m3" /> 
    844        <field id="DOM_E3T"     long_name="DOM * E3T"                                unit="mmol-N/m2"  > DOM * e3t </field > 
    845  
    846        <!-- CFC11 : variables available with key_cfc --> 
    847        <field id="CFC11"       long_name="CFC-11 Concentration"                     unit="umol/m3" /> 
    848        <field id="CFC11_E3T"   long_name="CFC11 * E3T"                              unit="umol/m2"  > CFC11 * e3t </field > 
    849        <!-- Bomb C14 : variables available with key_c14b --> 
    850        <field id="C14B"     long_name="Bomb C14 Concentration"                      unit="1"         /> 
    851        <field id="C14B_E3T"    long_name="C14B * E3T"                               unit="m"  > C14B * e3t </field > 
     996       <field id="DOM_e3t"     long_name="DOM * e3t"                                unit="mmol-N/m2"  > DOM * e3t </field > 
     997 
     998       <!-- CFC11 : variables available with ln_cfc11 --> 
     999       <field id="CFC11"       long_name="Chlorofluoro carbon11 Concentration"      unit="umol/m3" /> 
     1000       <field id="CFC11_e3t"   long_name="CFC11 * e3t"                              unit="umol/m2"  > CFC11 * e3t </field > 
     1001 
     1002       <!-- CFC12 : variables available with ln_cfc12 --> 
     1003       <field id="CFC12"       long_name="Chlorofluoro carbon12 Concentration"      unit="umol/m3" /> 
     1004       <field id="CFC12_e3t"   long_name="CFC12 * e3t"                              unit="umol/m2"  > CFC12 * e3t </field > 
     1005 
     1006       <!-- SF6 : variables available with ln_sf6 --> 
     1007       <field id="SF6"       long_name="Sulfur hexafluoride Concentration"      unit="umol/m3" /> 
     1008       <field id="SF6_e3t"   long_name="SF6 * e3t"                              unit="umol/m2"  > SF6 * e3t </field > 
     1009 
     1010       <!-- C14 : variables available with ln_c14 --> 
     1011       <field id="RC14"        long_name="Radiocarbon ratio"                        unit="-"         /> 
     1012       <field id="RC14_e3t"    long_name="RC14 * e3t"                               unit="m"  > RC14 * e3t </field > 
     1013 
     1014       <!-- AGE : variables available with ln_age --> 
     1015       <field id="Age"        long_name="Sea water age since surface contact"       unit="yr"         /> 
     1016       <field id="Age_e3t"    long_name="Age * e3t"                                 unit="yr * m"  > Age * e3t </field > 
     1017 
    8521018     </field_group> 
    8531019 
     
    8601026       <field id="PAR"         long_name="Photosynthetically Available Radiation"  unit="W/m2"       grid_ref="grid_T_3D" /> 
    8611027       <field id="PARDM"       long_name="Daily mean PAR"                          unit="W/m2"       grid_ref="grid_T_3D" /> 
    862        <field id="PPPHY"       long_name="Primary production of nanophyto"         unit="mol/m3/s"   grid_ref="grid_T_3D" /> 
    863        <field id="PPPHY2"      long_name="Primary production of diatoms"           unit="mol/m3/s"   grid_ref="grid_T_3D" /> 
    864        <field id="PPNEWN"      long_name="New Primary production of nanophyto"     unit="mol/m3/s"   grid_ref="grid_T_3D" /> 
    865        <field id="PPNEWD"      long_name="New Primary production of diatoms"       unit="mol/m3/s"   grid_ref="grid_T_3D" /> 
    866        <field id="PBSi"        long_name="Primary production of Si diatoms"        unit="mol/m3/s"   grid_ref="grid_T_3D" /> 
    867        <field id="PFeN"        long_name="Primary production of nano iron"         unit="mol/m3/s"   grid_ref="grid_T_3D" /> 
     1028       <field id="PPPHYN"      long_name="Primary production of nanophyto"         unit="molC/m3/s"  grid_ref="grid_T_3D" /> 
     1029       <field id="PPPHYP"      long_name="Primary production of picophyto"         unit="molC/m3/s"  grid_ref="grid_T_3D" /> 
     1030       <field id="PPPHYD"      long_name="Primary production of diatoms"           unit="molC/m3/s"  grid_ref="grid_T_3D" /> 
     1031       <field id="PPNEWN"      long_name="New Primary production of nanophyto"     unit="molC/m3/s"  grid_ref="grid_T_3D" /> 
     1032       <field id="PPNEWP"      long_name="New Primary production of picophyto"     unit="molC/m3/s"  grid_ref="grid_T_3D" /> 
     1033       <field id="PPNEWD"      long_name="New Primary production of diatoms"       unit="molC/m3/s"  grid_ref="grid_T_3D" /> 
     1034       <field id="PBSi"        long_name="Primary production of Si diatoms"        unit="molC/m3/s"  grid_ref="grid_T_3D" /> 
     1035       <field id="PFeN"        long_name="Primary production of nano iron"         unit="molC/m3/s"  grid_ref="grid_T_3D" /> 
     1036       <field id="PFeP"        long_name="Primary production of pico iron"         unit="molC/m3/s"  grid_ref="grid_T_3D" /> 
    8681037       <field id="PFeD"        long_name="Primary production of diatoms iron"      unit="mol/m3/s"   grid_ref="grid_T_3D" /> 
    8691038       <field id="xfracal"     long_name="Calcifying fraction"                     unit="1"          grid_ref="grid_T_3D" /> 
     
    8741043       <field id="REMIN"       long_name="Oxic remineralization of OM"             unit="mol/m3/s"   grid_ref="grid_T_3D" /> 
    8751044       <field id="DENIT"       long_name="Anoxic remineralization of OM"           unit="mol/m3/s"   grid_ref="grid_T_3D" /> 
     1045       <field id="REMINP"       long_name="Oxic remineralization rate of POC"      unit="d-1"        grid_ref="grid_T_3D" /> 
     1046       <field id="REMING"       long_name="Oxic remineralization rate of GOC"      unit="d-1"        grid_ref="grid_T_3D" /> 
    8761047       <field id="Nfix"        long_name="Nitrogen fixation"                       unit="mol/m3/s"   grid_ref="grid_T_3D" /> 
    8771048       <field id="Mumax"       long_name="Maximum growth rate"                     unit="s-1"        grid_ref="grid_T_3D" /> 
    8781049       <field id="MuN"         long_name="Realized growth rate for nanophyto"      unit="s-1"        grid_ref="grid_T_3D" /> 
     1050       <field id="MuP"         long_name="Realized growth rate for picophyto"      unit="s-1"        grid_ref="grid_T_3D" /> 
    8791051       <field id="MuD"         long_name="Realized growth rate for diatomes"       unit="s-1"        grid_ref="grid_T_3D" /> 
     1052       <field id="MunetN"      long_name="Net growth rate for nanophyto"           unit="s-1"        grid_ref="grid_T_3D" /> 
     1053       <field id="MunetP"      long_name="Net growth rate for picophyto"           unit="s-1"        grid_ref="grid_T_3D" /> 
     1054       <field id="MunetD"      long_name="Net growth rate for diatomes"            unit="s-1"        grid_ref="grid_T_3D" /> 
    8801055       <field id="LNnut"       long_name="Nutrient limitation term in Nanophyto"   unit=""           grid_ref="grid_T_3D" /> 
     1056       <field id="LPnut"       long_name="Nutrient limitation term in Picophyto"   unit="-"          grid_ref="grid_T_3D" /> 
    8811057       <field id="LDnut"       long_name="Nutrient limitation term in Diatoms"     unit=""           grid_ref="grid_T_3D" /> 
    8821058       <field id="LNFe"        long_name="Iron limitation term in Nanophyto"       unit=""           grid_ref="grid_T_3D" /> 
     1059       <field id="LPFe"        long_name="Iron limitation term in Picophyto"       unit="-"          grid_ref="grid_T_3D" /> 
    8831060       <field id="LDFe"        long_name="Iron limitation term in Diatoms"         unit=""           grid_ref="grid_T_3D" /> 
    8841061       <field id="LNlight"     long_name="Light limitation term in Nanophyto"      unit=""           grid_ref="grid_T_3D" /> 
     1062       <field id="LPlight"     long_name="Light limitation term in Picophyto"      unit="-"          grid_ref="grid_T_3D" /> 
    8851063       <field id="LDlight"     long_name="Light limitation term in Diatoms"        unit=""           grid_ref="grid_T_3D" /> 
     1064       <field id="SIZEN"       long_name="Mean relative size of nanophyto."        unit="-"          grid_ref="grid_T_3D" /> 
     1065       <field id="SIZEP"       long_name="Mean relative size of picophyto."        unit="-"          grid_ref="grid_T_3D" /> 
     1066       <field id="SIZED"       long_name="Mean relative size of diatoms"           unit="-"          grid_ref="grid_T_3D" /> 
    8861067       <field id="Fe2"         long_name="Iron II concentration"                   unit="nmol/m3"    grid_ref="grid_T_3D" /> 
    8871068       <field id="Fe3"         long_name="Iron III concentration"                  unit="nmol/m3"    grid_ref="grid_T_3D" /> 
     
    8961077       <field id="Sdenit"      long_name="Nitrate reduction in the sediments"      unit="mol/m2/s"                        /> 
    8971078       <field id="Ironice"     long_name="Iron input/uptake due to sea ice"        unit="mol/m2/s"                        /> 
     1079       <field id="SedCal"      long_name="Calcite burial in the sediments"         unit="molC/m2/s"                       /> 
     1080       <field id="SedSi"       long_name="Silicon burial in the sediments"         unit="molSi/m2/s"                      /> 
     1081       <field id="SedC"        long_name="Organic C burial in the sediments"       unit="molC/m2/s"                       /> 
    8981082       <field id="HYDR"        long_name="Iron input from hydrothemal vents"       unit="mol/m2/s"   grid_ref="grid_T_3D" /> 
    8991083       <field id="EPC100"      long_name="Export of carbon particles at 100 m"     unit="mol/m2/s"                        /> 
     
    9141098       <field id="Ironsed"     long_name="Iron deposition from sediment"           unit="mol/m2/s"   grid_ref="grid_T_3D" /> 
    9151099 
    916  
    917        <!-- PISCES with Kriest parametisation : variables available with key_kriest --> 
    918        <field id="EPN100"      long_name="Particulate number flux at 100 m"        unit="mol/m2/s"                        /> 
    919        <field id="EXPN"        long_name="Particulate number flux"                 unit="mol/m2/s"   grid_ref="grid_T_3D" /> 
    920        <field id="XNUM"        long_name="Number of particles in aggregats"        unit="1"          grid_ref="grid_T_3D" /> 
    921        <field id="WSC"         long_name="sinking speed of mass flux"              unit="m2/s"       grid_ref="grid_T_3D" /> 
    922        <field id="WSN"         long_name="sinking speed of number flux"            unit="m2/s"       grid_ref="grid_T_3D" /> 
    923  
    9241100       <!-- dbio_T on T grid : variables available with key_diaar5 --> 
    9251101       <field id="TPP"         long_name="Total Primary production of phyto"                   unit="mol/m3/s"  grid_ref="grid_T_3D" /> 
     
    9301106       <field id="ZO2MIN"      long_name="Depth of oxygen minimum concentration"               unit="m"                              /> 
    9311107       <field id="INTNFIX"     long_name="Nitrogen fixation rate : vert. integrated"           unit="mol/m2/s"                       /> 
    932        <field id="INTPPPHY"    long_name="Vertically integrated primary production by nanophy" unit="mol/m2/s"                       /> 
    933        <field id="INTPPPHY2"   long_name="Vertically integrated primary production by diatom"  unit="mol/m2/s"                       /> 
     1108       <field id="INTPPPHYN"    long_name="Vertically integrated primary production by nanophy" unit="mol/m2/s"                       /> 
     1109       <field id="INTPPPHYD"   long_name="Vertically integrated primary production by diatom"  unit="mol/m2/s"                       /> 
    9341110       <field id="INTPP"       long_name="Vertically integrated primary production by phyto"   unit="mol/m2/s"                       /> 
    9351111       <field id="INTPNEW"     long_name="Vertically integrated new primary production"        unit="mol/m2/s"                       /> 
     
    9611137       <field id="TDETSED"     long_name="TDETSED"                                 unit=""  />  
    9621138 
    963        <!-- CFC11 : variables available with key_cfc --> 
    964        <field id="qtrCFC11"    long_name="Air-sea flux of CFC-11"                  unit="mol/m2/s"   /> 
    965        <field id="qintCFC11"   long_name="Cumulative air-sea flux of CFC-11"       unit="mol/m2"     /> 
    966  
    967        <!-- Bomb C14 : variables available with key_c14b --> 
    968        <field id="qtrC14b"     long_name="Air-sea flux of Bomb C14"                unit="mol/m2/s"   /> 
    969        <field id="qintC14b"    long_name="Cumulative air-sea flux of Bomb C14"     unit="mol/m2"     /> 
    970        <field id="fdecay"      long_name="Radiactive decay of Bomb C14"            unit="mol/m3"  grid_ref="grid_T_3D"  /> 
     1139       <!-- CFC11 : variables available with ln_cfc11 --> 
     1140       <field id="qtr_CFC11"    long_name="Air-sea flux of CFC-11"                  unit="mol/m2/s"   /> 
     1141       <field id="qint_CFC11"   long_name="Cumulative air-sea flux of CFC-11"       unit="mol/m2"     /> 
     1142 
     1143       <!-- CFC12 : variables available with ln_cfc12 --> 
     1144       <field id="qtr_CFC12"    long_name="Air-sea flux of CFC12"                  unit="mol/m2/s"   /> 
     1145       <field id="qint_CFC12"   long_name="Cumulative air-sea flux of CFC12"       unit="mol/m2"     /> 
     1146 
     1147       <!-- SF6 : variables available with ln_sf6 --> 
     1148       <field id="qtr_SF6"      long_name="Air-sea flux of SF6"                    unit="mol/m2/s"   /> 
     1149       <field id="qint_SF6"     long_name="Cumulative air-sea flux of SF6"         unit="mol/m2"     /> 
     1150 
     1151       <!--  C14 : variables available with ln_c14 --> 
     1152       <field id="DeltaC14"     long_name="Delta C14"                              unit="permil" grid_ref="grid_T_3D"   /> 
     1153       <field id="C14Age"       long_name="Radiocarbon age"                        unit="yr"     grid_ref="grid_T_3D"   /> 
     1154       <field id="RAge"         long_name="Reservoir Age"                          unit="yr"     /> 
     1155       <field id="qtr_C14"      long_name="Air-sea flux of C14"                    unit="1/m2/s"   /> 
     1156       <field id="qint_C14"     long_name="Cumulative air-sea flux of C14"         unit="1/m2"     /> 
    9711157     </field_group> 
    9721158 
    973      <field_group id="PISCES_scalar"  domain_ref="1point" > 
    974        <field id="pno3tot"         long_name="global mean nitrate concentration"                  unit="mol/m3"   /> 
     1159     <field_group id="tracer_scalar"  domain_ref="1point" > 
     1160     <!-- PISCES scalar  --> 
     1161       <field id="pno3tot"         long_name="Global mean nitrate concentration"                  unit="mol/m3"   /> 
    9751162       <field id="ppo4tot"         long_name="global mean phosphorus concentration"               unit="mol/m3"   /> 
    976        <field id="psiltot"         long_name="global mean silicate concentration"                 unit="mol/m3"   /> 
    977        <field id="palktot"         long_name="global mean alkalinity concentration"               unit="mol/m3"   /> 
    978        <field id="pfertot"         long_name="global mean iron concentration"                     unit="mol/m3"   /> 
    979        <field id="tcflx"           long_name="total Flux of Carbon out of the ocean"              unit="mol/s"   /> 
    980        <field id="tcflxcum"        long_name="cumulative total Flux of Carbon out of the ocean"   unit="mol/s"   /> 
    981        <field id="tcexp"           long_name="total Carbon export at 100m"                        unit="mol/s"   /> 
    982        <field id="tintpp"          long_name="global total integrated primary production"         unit="mol/s"   /> 
    983        <field id="tnfix"           long_name="global total nitrogen fixation"                     unit="mol/s"   /> 
     1163       <field id="psiltot"         long_name="Global mean silicate concentration"                 unit="mol/m3"   /> 
     1164       <field id="palktot"         long_name="Global mean alkalinity concentration"               unit="mol/m3"   /> 
     1165       <field id="pfertot"         long_name="Global mean iron concentration"                     unit="mol/m3"   /> 
     1166       <field id="tcflx"           long_name="Total Flux of Carbon out of the ocean"              unit="mol/s"   /> 
     1167       <field id="tcflxcum"        long_name="Cumulative total Flux of Carbon out of the ocean"   unit="mol/s"   /> 
     1168       <field id="tcexp"           long_name="Total Carbon export at 100m"                        unit="mol/s"   /> 
     1169       <field id="tintpp"          long_name="Global total integrated primary production"         unit="mol/s"   /> 
     1170       <field id="tnfix"           long_name="Global total nitrogen fixation"                     unit="mol/s"   /> 
    9841171       <field id="tdenit"          long_name="Total denitrification"                              unit="mol/s"   /> 
     1172     <!-- C14 scalar  --> 
     1173       <field id="AtmCO2"          long_name="Global atmospheric CO2"                             unit="ppm"   /> 
     1174       <field id="AtmC14"          long_name="Global atmospheric DeltaC14"                        unit="permil"   /> 
     1175       <field id="K_C14"           long_name="Global 14C/C exchange velocity"                     unit="m/yr"   /> 
     1176       <field id="K_CO2"           long_name="Global CO2 piston velocity"                         unit="cm/h"   /> 
     1177       <field id="C14Inv"          long_name="global Radiocarbon ocean inventory"                 unit="10^26 atoms"   /> 
    9851178     </field_group> 
    9861179 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/SHARED/namelist_pisces_ref

    r6945 r7403  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! PISCES  (key_pisces) reference namelist (see below for key_pisces_reduced) 
     2!! PISCES reference namelist  
    33!!              1  - air-sea exchange                         (nampisext) 
    44!!              2  - biological parameters                    (nampisbio) 
     
    99!!              7  - parameters for calcite chemistry         (nampiscal) 
    1010!!              8  - parameters for inputs deposition         (nampissed) 
    11 !!              9  - parameters for Kriest parameterization   (nampiskrp, nampiskrs) 
    12 !!              10 - additional 2D/3D  diagnostics            (nampisdia) 
    1311!!              11 - Damping                                  (nampisdmp) 
    14 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    15 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     12!----------------------------------------------------------------------- 
     13&nampismod     !  Model used  
     14!----------------------------------------------------------------------- 
     15  ln_p2z    = .false.        !  LOBSTER model used 
     16  ln_p4z    = .true.         !  PISCES model used 
     17  ln_p5z    = .false.        !  PISCES QUOTA model used 
     18  ln_ligand = .false.        !  Enable  organic ligands 
     19/ 
     20!----------------------------------------------------------------------- 
    1621&nampisext     !   air-sea exchange 
    17 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     22!----------------------------------------------------------------------- 
    1823   ln_co2int  =  .false. ! read atm pco2 from a file (T) or constant (F) 
    1924   atcco2     =  280.    ! Constant value atmospheric pCO2 - ln_co2int = F 
     
    2328!                        ! then the first atmospheric CO2 record read is at years(1) 
    2429/ 
    25 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     30!----------------------------------------------------------------------- 
    2631&nampisatm     !  Atmospheric prrssure  
    27 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     32!----------------------------------------------------------------------- 
    2833!              !  file name   ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    2934!              !              !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    3035   sn_patm     = 'presatm'    ,     -1            , 'patm'    ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
     36   sn_atmco2   = 'presatmco2' ,     -1            , 'xco2'    ,  .true.      , .true. , 'yearly'  , ''       , ''       , '' 
    3137   cn_dir      = './'     !  root directory for the location of the dynamical files 
    3238! 
    33    ln_presatm  = .false.   ! constant atmopsheric pressure (F) or from a file (T) 
    34 / 
    35 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     39   ln_presatm    = .false.   ! constant atmopsheric pressure (F) or from a file (T) 
     40   ln_presatmco2 = .false.   ! Read spatialized atm co2 files [ppm] if TRUE 
     41/ 
     42!----------------------------------------------------------------------- 
    3643&nampisbio     !   biological parameters 
    37 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     44!----------------------------------------------------------------------- 
    3845   nrdttrc    =  1        ! time step frequency for biology 
    3946   wsbio      =  2.       ! POC sinking speed 
    4047   xkmort     =  2.E-7    ! half saturation constant for mortality 
    4148   ferat3     =  10.E-6   ! Fe/C in zooplankton  
    42    wsbio2     =  30.      ! Big particles sinking speed 
     49   wsbio2     =  50.      ! Big particles sinking speed 
     50   wsbio2max  =  50.      ! Big particles maximum sinking speed 
     51   wsbio2scale=  5000.    ! Big particles length scale of sinking 
    4352   niter1max  =  1        ! Maximum number of iterations for POC 
    44    niter2max  =  1        ! Maximum number of iterations for GOC 
    45 / 
    46 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    47 &nampislim     !   parameters for nutrient limitations 
    48 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     53   niter2max  =  2        ! Maximum number of iterations for GOC 
     54!                         !  ln_ligand enabled 
     55   wfep       =  0.2      ! FeP sinking speed  
     56   ldocp      =  1.E-5    ! Phyto ligand production per unit doc  
     57   ldocz      =  1.E-5    ! Zoo ligand production per unit doc  
     58   lthet      =  0.5      ! Proportional loss of ligands due to Fe uptake  
     59!                         !  ln_p5z enabled 
     60   no3rat3    =  0.182    ! N/C ratio in zooplankton 
     61   po4rat3    =  0.0094   ! P/C ratio in zooplankton 
     62/ 
     63!----------------------------------------------------------------------- 
     64&namp4zlim     !   parameters for nutrient limitations for PISCES std - ln_p4z 
     65!----------------------------------------------------------------------- 
    4966   concnno3   =  1.e-6    ! Nitrate half saturation of nanophytoplankton 
    5067   concdno3   =  3.E-6    ! Nitrate half saturation for diatoms 
     
    6683   qdfelim    =  7.E-6    ! Optimal quota of diatoms 
    6784   caco3r     =  0.3      ! mean rain ratio 
    68    oxymin    =  1.E-6     ! Half-saturation constant for anoxia 
    69 / 
    70 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     85   oxymin     =  1.E-6    ! Half-saturation constant for anoxia 
     86/ 
     87!----------------------------------------------------------------------- 
     88&namp5zlim     !   parameters for nutrient limitations PISCES QUOTA - ln_p5z 
     89!----------------------------------------------------------------------- 
     90   concnno3   =  3e-6     ! Nitrate half saturation of nanophytoplankton 
     91   concpno3   =  1e-6 
     92   concdno3   =  4E-6     ! Phosphate half saturation for diatoms 
     93   concnnh4   =  1.5E-6   ! NH4 half saturation for phyto 
     94   concpnh4   =  4E-7 
     95   concdnh4   =  2E-6     ! NH4 half saturation for diatoms 
     96   concnpo4   =  3E-6     ! PO4 half saturation for phyto 
     97   concppo4   =  1.5E-6 
     98   concdpo4   =  4E-6     ! PO4 half saturation for diatoms 
     99   concnfer   =  3E-9   ! Iron half saturation for phyto 
     100   concpfer   =  1.5E-9 
     101   concdfer   =  4E-9   ! Iron half saturation for diatoms 
     102   concbfe    =  1.E-11   ! Half-saturation for Fe limitation of Bacteria 
     103   concbnh4   =  1.E-7    ! NH4 half saturation for phyto 
     104   concbno3   =  5.E-7    ! Phosphate half saturation for diatoms 
     105   concbpo4   =  1E-7     ! Phosphate half saturation for bacteria 
     106   xsizedia   =  1.E-6    ! Minimum size criteria for diatoms 
     107   xsizephy   =  1.E-6    ! Minimum size criteria for phyto 
     108   xsizepic   =  1.E-6 
     109   xsizern    =  1.0      ! Size ratio for nanophytoplankton 
     110   xsizerp    =  1.0 
     111   xsizerd    =  4.0      ! Size ratio for diatoms 
     112   xksi1      =  2.E-6    ! half saturation constant for Si uptake 
     113   xksi2      =  20E-6  ! half saturation constant for Si/C 
     114   xkdoc      =  417.E-6  ! half-saturation constant of DOC remineralization 
     115   caco3r     =  0.35     ! mean rain ratio 
     116   oxymin     =  1.E-6    ! Half-saturation constant for anoxia 
     117/ 
     118!----------------------------------------------------------------------- 
     119&namp5zquota    !   parameters for nutrient limitations PISCES quota - ln_p5z 
     120!----------------------------------------------------------------------- 
     121   qfnopt     =  7.E-6    ! Optimal Fe quota of nanophyto 
     122   qfpopt     =  7.E-6    ! Optimal Fe quota of picophyto 
     123   qfdopt     =  7.E-6    ! Optimal quota of diatoms 
     124   qnnmin     =  0.29     ! Minimal N quota for nano 
     125   qnnmax     =  1.39     ! Maximal N quota for nano 
     126   qpnmin     =  0.28     ! Minimal P quota for nano 
     127   qpnmax     =  1.06     ! Maximal P quota for nano 
     128   qnpmin     =  0.42     ! Minimal N quota for pico 
     129   qnpmax     =  1.39     ! Maximal N quota for pico 
     130   qppmin     =  0.25     ! Minimal P quota for pico 
     131   qppmax     =  0.7      ! Maximal P quota for pico 
     132   qndmin     =  0.25     ! Minimal N quota for diatoms 
     133   qndmax     =  1.39     ! Maximal N quota for diatoms 
     134   qpdmin     =  0.29     ! Minimal P quota for diatoms 
     135   qpdmax     =  1.32     ! Maximal P quota for diatoms 
     136   qfnmax     =  40E-6    ! Maximal Fe quota for nano 
     137   qfpmax     =  40E-6    ! Maximal Fe quota for pico 
     138   qfdmax     =  40E-6    ! Maximal Fe quota for diatoms 
     139/ 
     140!----------------------------------------------------------------------- 
    71141&nampisopt     !   parameters for optics 
    72 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     142!----------------------------------------------------------------------- 
    73143!              !  file name       ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    74144!              !                  !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     
    78148   parlux      =  0.43      ! Fraction of shortwave as PAR 
    79149/  
    80 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    81 &nampisprod     !   parameters for phytoplankton growth 
    82 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    83    pislope    =  2.       ! P-I slope 
    84    pislope2   =  2.       ! P-I slope  for diatoms 
     150!----------------------------------------------------------------------- 
     151&namp4zprod     !   parameters for phytoplankton growth for PISCES std - ln_p4z 
     152!----------------------------------------------------------------------- 
     153   pislopen   =  2.       ! P-I slope 
     154   pisloped   =  2.       ! P-I slope  for diatoms 
    85155   xadap      =  0.       ! Adaptation factor to low light 
    86    excret     =  0.05     ! excretion ratio of phytoplankton 
    87    excret2    =  0.05     ! excretion ratio of diatoms 
     156   excretn    =  0.05     ! excretion ratio of phytoplankton 
     157   excretd    =  0.05     ! excretion ratio of diatoms 
    88158   ln_newprod =  .true.   ! Enable new parame. of production (T/F)  
    89159   bresp      =  0.033    ! Basal respiration rate 
     
    95165   grosip     =  0.159    ! mean Si/C ratio 
    96166/ 
    97 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    98 &nampismort     !   parameters for phytoplankton sinks 
    99 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    100    wchl       =  0.01    ! quadratic mortality of phytoplankton 
     167!----------------------------------------------------------------------- 
     168&namp5zprod     !   parameters for phytoplankton growth for PISCES quota - ln_p5z 
     169!----------------------------------------------------------------------- 
     170   pislopen   =  3.       ! P-I slope 
     171   pislopep   =  3.       ! P-I slope for picophytoplankton 
     172   pisloped   =  3.       ! P-I slope  for diatoms 
     173   excretn    =  0.05     ! excretion ratio of phytoplankton 
     174   excretp    =  0.05     ! excretion ratio of picophytoplankton 
     175   excretd    =  0.05     ! excretion ratio of diatoms 
     176   xadap      =  0.       ! Adaptation factor to low light 
     177   bresp      =  0.02     ! Basal respiration rate 
     178   thetannm   =  0.25     ! Maximum Chl/N in nanophytoplankton 
     179   thetanpm   =  0.25     ! Maximum Chl/N in picophytoplankton 
     180   thetandm   =  0.3      ! Maximum Chl/N in diatoms 
     181   chlcmin    =  0.004    ! Minimum Chl/c in phytoplankton 
     182   grosip     =  0.131    ! mean Si/C ratio 
     183/ 
     184!----------------------------------------------------------------------- 
     185&namp4zmort     !   parameters for phytoplankton sinks for PISCES std - ln_p4z 
     186!----------------------------------------------------------------------- 
     187   wchl       =  0.01     ! quadratic mortality of phytoplankton 
    101188   wchld      =  0.01     ! maximum quadratic mortality of diatoms 
    102189   wchldm     =  0.03     ! maximum quadratic mortality of diatoms 
     
    104191   mprat2     =  0.01     ! Diatoms mortality rate 
    105192/ 
    106 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    107 &nampismes     !   parameters for mesozooplankton 
    108 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     193!----------------------------------------------------------------------- 
     194&namp5zmort     !   parameters for phytoplankton sinks for PISCES quota - ln_p5z 
     195!----------------------------------------------------------------------- 
     196   wchln      =  0.01     ! quadratic mortality of nanophytoplankton 
     197   wchlp      =  0.01     ! quadratic mortality of picophytoplankton 
     198   wchld      =  0.01     ! maximum quadratic mortality of diatoms 
     199   wchldm     =  0.02     ! maximum quadratic mortality of diatoms 
     200   mpratn     =  0.01     ! nanophytoplankton mortality rate 
     201   mpratp     =  0.01     ! picophytoplankton mortality rate 
     202   mprat2     =  0.01     ! Diatoms mortality rate 
     203/ 
     204!----------------------------------------------------------------------- 
     205&namp4zmes     !   parameters for mesozooplankton for PISCES std - ln_p4z 
     206!----------------------------------------------------------------------- 
    109207   part2      =  0.75     ! part of calcite not dissolved in mesozoo guts 
    110208   grazrat2   =  0.75     ! maximal mesozoo grazing rate 
     
    126224   grazflux   =  2.e3     ! flux-feeding rate 
    127225/ 
    128 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    129 &nampiszoo     !   parameters for microzooplankton 
    130 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    131    part       =  0.5      ! part of calcite not dissolved in microzoo gutsa 
     226!----------------------------------------------------------------------- 
     227&namp5zmes     !   parameters for mesozooplankton 
     228!----------------------------------------------------------------------- 
     229   part2      =  0.75     ! part of calcite not dissolved in mesozoo guts 
     230   grazrat2   =  0.85     ! maximal mesozoo grazing rate 
     231   bmetexc2   =  .true.   ! Metabolic use of excess carbon  
     232   resrat2    =  0.005    ! exsudation rate of mesozooplankton 
     233   mzrat2     =  0.02     ! mesozooplankton mortality rate 
     234   xpref2d    =  1.       ! zoo preference for phyto 
     235   xpref2p    =  1.       ! zoo preference for POC 
     236   xpref2z    =  1.       ! zoo preference for zoo 
     237   xpref2m    =  0.2      ! meso preference for zoo 
     238   xpref2c    =  0.3      ! zoo preference for poc 
     239   xthresh2zoo = 1E-8     ! zoo feeding threshold for mesozooplankton 
     240   xthresh2dia = 1E-8     ! diatoms feeding threshold for mesozooplankton 
     241   xthresh2phy = 1E-8     ! nanophyto feeding threshold for mesozooplankton 
     242   xthresh2mes = 1E-8     ! meso feeding threshold for mesozooplankton 
     243   xthresh2poc = 1E-8     ! poc feeding threshold for mesozooplankton 
     244   xthresh2    =  3E-7     ! Food threshold for grazing 
     245   xkgraz2     =  20.E-6   ! half sturation constant for meso grazing 
     246   epsher2     =  0.5      ! Efficicency of Mesozoo growth 
     247   ssigma2     =  0.5     ! Fraction excreted as semi-labile DOM 
     248   srespir2    =  0.2     ! Active respiration 
     249   unass2c     =  0.3     ! non assimilated fraction of P by mesozoo 
     250   unass2n     =  0.3     ! non assimilated fraction of N by mesozoo 
     251   unass2p     =  0.3     ! non assimilated fraction of P by mesozoo 
     252   grazflux   =  3.e3     ! flux-feeding rate 
     253/ 
     254!----------------------------------------------------------------------- 
     255&namp4zzoo     !   parameters for microzooplankton for PISCES std - ln_p4z 
     256!----------------------------------------------------------------------- 
     257   part       =  0.5      ! part of calcite not dissolved in microzoo guts 
    132258   grazrat    =  3.0      ! maximal zoo grazing rate 
    133259   resrat     =  0.03     ! exsudation rate of zooplankton 
     
    145271   unass      =  0.3      ! non assimilated fraction of phyto by zoo 
    146272/ 
    147 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     273!----------------------------------------------------------------------- 
     274&namp5zzoo     !   parameters for microzooplankton 
     275!----------------------------------------------------------------------- 
     276   part       =  0.5      ! part of calcite not dissolved in microzoo gutsa 
     277   grazrat    =  2.75     ! maximal zoo grazing rate 
     278   bmetexc    =  .true.   ! Metabolic use of excess carbon 
     279   resrat     =  0.03     ! exsudation rate of zooplankton 
     280   mzrat      =  0.005    ! zooplankton mortality rate 
     281   xprefc     =  0.1      ! Microzoo preference for POM 
     282   xprefn     =  1.       ! Microzoo preference for Nanophyto 
     283   xprefp     =  1.6      ! Microzoo preference for picophyto 
     284   xprefd     =  1.0      ! Microzoo preference for Diatoms 
     285   xprefz     =  0.3      ! Microzoo preference for microzooplankton 
     286   xthreshdia =  1.E-8    ! Diatoms feeding threshold for microzooplankton 
     287   xthreshphy =  1.E-8    ! Nanophyto feeding threshold for microzooplankton 
     288   xthreshpic =  1.E-8 
     289   xthreshzoo =  1.E-8    ! Nanophyto feeding threshold for microzooplankton 
     290   xthreshpoc =  1.E-8    ! POC feeding threshold for microzooplankton 
     291   xthresh    =  3.E-7    ! Food threshold for feeding 
     292   xkgraz     =  20.E-6   ! half sturation constant for grazing 
     293   epsher     =  0.5      ! Efficiency of microzoo growth 
     294   ssigma     =  0.5      ! Fraction excreted as semi-labile DOM 
     295   srespir    =  0.2      ! Active respiration 
     296   unassc     =  0.3      ! non assimilated fraction of C by zoo 
     297   unassn     =  0.3      ! non assimilated fraction of C by zoo 
     298   unassp     =  0.3      ! non assimilated fraction of C by zoo 
     299/ 
     300!----------------------------------------------------------------------- 
    148301&nampisfer     !   parameters for iron chemistry 
    149 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     302!----------------------------------------------------------------------- 
    150303   ln_fechem =  .false.   ! complex iron chemistry ( T/F ) 
    151304   ln_ligvar =  .false.   ! variable ligand concentration 
    152    xlam1     =  0.005     ! scavenging rate of Iron 
    153    xlamdust  =  150.0     ! Scavenging rate of dust 
    154    ligand    =  0.6E-9    ! Ligands concentration  
    155  
    156 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     305   ln_fecolloid = .false. ! variable colloidal fraction 
     306   xlam1        =  0.005  ! scavenging rate of Iron 
     307   xlamdust     =  150.0  ! Scavenging rate of dust 
     308   ligand       =  0.6E-9 ! Ligands concentration  
     309   kfep         =  0.     ! Nanoparticle formation rate constant 
     310/ 
     311!-----------------------------------------------------------------------   
    157312&nampisrem     !   parameters for remineralization 
    158 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     313!----------------------------------------------------------------------- 
    159314   xremik    =  0.3       ! remineralization rate of DOC 
    160    xremip    =  0.025     ! remineralisation rate of POC 
    161315   nitrif    =  0.05      ! NH4 nitrification rate 
    162316   xsirem    =  0.003     ! remineralization rate of Si 
    163317   xsiremlab =  0.03      ! fast remineralization rate of Si 
    164318   xsilab    =  0.5       ! Fraction of labile biogenic silica 
    165 / 
    166 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     319   feratb    =  10.E-6    ! Fe/C quota in bacteria 
     320   xkferb    =  2.5E-10   ! Half-saturation constant for bacteria Fe/C 
     321!                         ! ln_p5z 
     322   xremikc   =  0.25      ! remineralization rate of DOC 
     323   xremikn   =  0.35      ! remineralization rate of DON 
     324   xremikp   =  0.4       ! remineralization rate of DOP 
     325!   feratb    =  20E-6     ! Bacterial Fe/C ratio 
     326!   xkferb    =  3E-10     ! Half-saturation constant for bact. Fe/C 
     327/ 
     328!----------------------------------------------------------------------- 
     329&nampispoc     !   parameters for organic particles 
     330!----------------------------------------------------------------------- 
     331   xremip    =  0.035     ! remineralisation rate of PON 
     332   jcpoc     =  15        ! Number of lability classes 
     333   rshape    =  1.0       ! Shape of the gamma function 
     334!                         ! ln_p5z 
     335   xremipc   =  0.02      ! remineralisation rate of POC 
     336   xremipn   =  0.025     ! remineralisation rate of PON 
     337   xremipp   =  0.03      ! remineralisation rate of POP 
     338/ 
     339!----------------------------------------------------------------------- 
    167340&nampiscal     !   parameters for Calcite chemistry 
    168 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     341!----------------------------------------------------------------------- 
    169342   kdca       =  6.       ! calcite dissolution rate constant (1/time) 
    170343   nca        =  1.       ! order of dissolution reaction (dimensionless) 
    171344/ 
    172 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     345!----------------------------------------------------------------------- 
    173346&nampissbc     !   parameters for inputs deposition 
    174 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     347!----------------------------------------------------------------------- 
    175348!              !  file name        ! frequency (hours) ! variable      ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    176349!              !                   !  (if <0  months)  !   name        !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     
    205378   concfediaz  =  1.e-10   ! Diazotrophs half-saturation Cste for Iron 
    206379   hratio      =  1.e+7    ! Fe to 3He ratio assumed for vent iron supply  
    207 / 
    208 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     380!                          ! ln_ligand 
     381   fep_rats    =  1.       ! Fep/Fer ratio from sed sources  
     382   fep_rath    =  1.       ! Fep/Fer ratio from sed hydro sources  
     383   lgw_rath    =  0.5      ! Weak ligand ratio from sed hydro sources  
     384/ 
     385!----------------------------------------------------------------------- 
     386&nampislig      !  Namelist parameters for ligands, nampislig 
     387!----------------------------------------------------------------------- 
     388   rfep        =  0.001    ! Dissolution rate of FeP 
     389   rlgw        =  1.       ! Lifetime (years) of weak ligands 
     390   rlig        =  1.E-4    ! Remin ligand production per unit C 
     391   prlgw       =  1.E-4    ! Photolysis of weak ligand 
     392   rlgs        =  1000.    ! Lifetime (years) of strong ligands 
     393/ 
     394!----------------------------------------------------------------------- 
    209395&nampisice      !  Prescribed sea ice tracers 
    210 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     396!----------------------------------------------------------------------- 
    211397! constant ocean tracer concentrations are defined in trcice_pisces.F90 (Global, Arctic, Antarctic and Baltic) 
    212398! trc_ice_ratio     * betw 0 and 1: prescribed ice/ocean tracer concentration ratio 
     
    219405! cn_trc_o          * 'GL' use global ocean values making the Baltic distinction only 
    220406!                     'AA' use specific Arctic/Antarctic/Baltic values 
    221 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     407!----------------------------------------------------------------------- 
    222408!    sn_tri_ ! trc_ice_ratio ! trc_ice_prescr !     cn_trc_o 
    223409   sn_tri_dic =           -1.,           -99.,          'AA' 
     
    247433   sn_tri_nh4 =            1.,           -99.,          'AA' 
    248434/ 
    249 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    250 &nampiskrp     !   Kriest parameterization : parameters     "key_kriest" 
    251 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    252    xkr_eta      = 1.17    ! Sinking  exponent 
    253    xkr_zeta     = 2.28    ! N content exponent 
    254    xkr_ncontent = 5.7E-6  ! N content factor 
    255    xkr_mass_min = 0.0002  ! Minimum mass for Aggregates 
    256    xkr_mass_max = 1.      ! Maximum mass for Aggregates 
    257 / 
    258 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    259 &nampiskrs     !   Kriest parameterization : size classes  "key_kriest" 
    260 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    261    xkr_sfact    = 942.    ! Sinking factor 
    262    xkr_stick    = 0.5     ! Stickiness 
    263    xkr_nnano    = 2.337   ! Nbr of cell in nano size class 
    264    xkr_ndiat    = 3.718   ! Nbr of cell in diatoms size class 
    265    xkr_nmeso    = 7.147   ! Nbr of cell in mesozoo size class 
    266    xkr_naggr    = 9.877   ! Nbr of cell in aggregates  size class 
    267 / 
    268 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    269 &nampisdia     !   additional 2D/3D tracers diagnostics  
    270 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    271 !              !    name   !           title of the field          !     units      ! 
    272 !              !           !                                       !                !   
    273    pisdia2d(1)  = 'Cflx     ' , 'DIC flux                          ',  'molC/m2/s    ' 
    274    pisdia2d(2)  = 'Oflx     ' , 'Oxygen flux                       ',  'molC/m2/s    ' 
    275    pisdia2d(3)  = 'Kg       ' , 'Gas transfer                      ',  'mol/m2/s/uatm' 
    276    pisdia2d(4)  = 'Delc     ' , 'Delta CO2                         ',  'uatm         ' 
    277    pisdia2d(5)  = 'PMO      ' , 'POC export                        ',  'molC/m2/s    ' 
    278    pisdia2d(6)  = 'PMO2     ' , 'GOC export                        ',  'molC/m2/s    ' 
    279    pisdia2d(7)  = 'ExpFe1   ' , 'Nano iron export                  ',  'molFe/m2/s   ' 
    280    pisdia2d(8)  = 'ExpFe2   ' , 'Diatoms iron export               ',  'molFe/m2/s   ' 
    281    pisdia2d(9)  = 'ExpSi    ' , 'Silicate export                   ',  'molSi/m2/s   ' 
    282    pisdia2d(10) = 'ExpCaCO3 ' , 'Calcite export                    ',  'molC/m2/s    ' 
    283    pisdia2d(11) = 'heup     ' , 'euphotic layer depth              ',  'm            ' 
    284    pisdia2d(12) = 'Fedep    ' , 'Iron dep                          ',  'molFe/m2/s   ' 
    285    pisdia2d(13) = 'Nfix     ' , 'Nitrogen Fixation                 ',  'molN/m2/s    ' 
    286    pisdia3d(1)  = 'PH       ' , 'PH                                ',  '-            ' 
    287    pisdia3d(2)  = 'CO3      ' , 'Bicarbonates                      ',  'mol/l        ' 
    288    pisdia3d(3)  = 'CO3sat   ' , 'CO3 saturation                    ',  'mol/l        ' 
    289    pisdia3d(4)  = 'PAR      ' , 'light penetration                 ',  'W/m2         ' 
    290    pisdia3d(5)  = 'PPPHY    ' , 'Primary production of nanophyto   ',  'molC/m3/s    ' 
    291    pisdia3d(6)  = 'PPPHY2   ' , 'Primary production of diatoms     ',  'molC/m3/s    ' 
    292    pisdia3d(7)  = 'PPNEWN   ' , 'New Primary production of nano    ',  'molC/m3/s    ' 
    293    pisdia3d(8)  = 'PPNEWD   ' , 'New Primary production of diat    ',  'molC/m3/s    ' 
    294    pisdia3d(9)  = 'PBSi     ' , 'Primary production of Si diatoms  ',  'molSi/m3/s   ' 
    295    pisdia3d(10) = 'PFeN     ' , 'Primary production of nano iron   ',  'molFe/m3/s   ' 
    296    pisdia3d(11) = 'PFeD     ' , 'Primary production of diatoms iron',  'molFe/m3/s   ' 
    297 / 
    298 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     435!----------------------------------------------------------------------- 
    299436&nampisdmp     !  Damping  
    300 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     437!----------------------------------------------------------------------- 
    301438   ln_pisdmp    =  .true.     !  Relaxation fo some tracers to a mean value 
    302439   nn_pisdmp    =  5475       !  Frequency of Relaxation  
    303440/ 
    304 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     441!----------------------------------------------------------------------- 
    305442&nampismass     !  Mass conservation 
    306 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     443!----------------------------------------------------------------------- 
    307444   ln_check_mass =  .false.    !  Check mass conservation 
    308445/ 
     
    317454!!              7  - general coefficients                       (namlobrat) 
    318455!!              8  - optical parameters                         (namlobopt) 
    319  
    320 !!              10 - biological diagnostics trends              (namlobdbi)  
    321456!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    322 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     457!----------------------------------------------------------------------- 
    323458&namlobphy     !   biological parameters for phytoplankton 
    324 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     459!----------------------------------------------------------------------- 
    325460   tmumax  =  1.21e-5   ! maximal phytoplankton growth rate            [s-1]  
    326461   rgamma  =  0.05      ! phytoplankton exudation fraction             [%] 
     
    329464   aki     =  33.       ! light photosynthesis half saturation constant[W/m2] 
    330465/ 
    331 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     466!----------------------------------------------------------------------- 
    332467&namlobnut     !   biological parameters for nutrients 
    333 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     468!----------------------------------------------------------------------- 
    334469   akno3   =  0.7       ! nitrate limitation half-saturation value     [mmol/m3] 
    335470   aknh4   =  0.001     ! ammonium limitation half-saturation value    [mmol/m3] 
     
    337472   psinut  =  3.        ! inhibition of nitrate uptake by ammonium 
    338473/ 
    339 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     474!----------------------------------------------------------------------- 
    340475&namlobzoo     !   biological parameters for zooplankton 
    341 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     476!----------------------------------------------------------------------- 
    342477   rppz    = 0.8        ! zooplankton nominal preference for phytoplancton food  [%] 
    343478   taus    = 9.26E-6    ! specific zooplankton maximal grazing rate              [s-1]  
     
    351486   tmminz  = 2.31e-6    ! minimal zooplankton mortality rate                     [(mmolN/m3)-1 d-1] 
    352487/ 
    353 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     488!----------------------------------------------------------------------- 
    354489&namlobdet     !   biological parameters for detritus 
    355 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     490!----------------------------------------------------------------------- 
    356491   taudn   = 5.80e-7    ! detritus breakdown rate                        [0.1/86400 s-1=10 days] 
    357492   fdetlab = 0.         ! NH4 fraction of detritus dissolution            
    358493/ 
    359 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     494!----------------------------------------------------------------------- 
    360495&namlobdom     !   biological parameters for DOM 
    361 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     496!----------------------------------------------------------------------- 
    362497   taudomn = 6.43e-8    ! DOM breakdown rate                             [s-1]  
    363498!                       ! slow remineralization rate of semi-labile dom to nh4 (1 month) 
    364499/ 
    365 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     500!----------------------------------------------------------------------- 
    366501&namlobsed     !   parameters from aphotic layers to sediment 
    367 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     502!----------------------------------------------------------------------- 
    368503   sedlam     = 3.86e-7    ! time coefficient of POC remineralization in sediments [s-1] 
    369504   sedlostpoc = 0.         ! mass of POC lost in sediments 
     
    371506   xhr        = -0.858     ! coeff for martin''s remineralisation profile 
    372507/ 
    373 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     508!----------------------------------------------------------------------- 
    374509&namlobrat     !   general coefficients 
    375 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     510!----------------------------------------------------------------------- 
    376511   rcchl   = 60.       ! Carbone/Chlorophyl ratio                     [mgC.mgChla-1] 
    377512   redf    = 6.56       ! redfield ratio (C:N) for phyto 
    378513   reddom  = 6.56       ! redfield ratio (C:N) for DOM 
    379514/ 
    380 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
     515!----------------------------------------------------------------------- 
    381516&namlobopt     !   optical parameters 
    382 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
     517!----------------------------------------------------------------------- 
    383518   xkg0   = 0.0232     ! green absorption coefficient of water 
    384519   xkr0   = 0.225      ! red absorption coefficent of water 
     
    389524   rpig   = 0.7        ! chla/chla+pheo ratio 
    390525/ 
    391 !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' 
    392 &nampisdbi     !   biological diagnostics trends      
    393 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, 
    394 !                !  2D bio diagnostics   units : mmole/m2/s   ("key_trdmld_trc") 
    395 !                !  name    !       title of the field      !     units      ! 
    396    pisdiabio(1)  = 'NO3PHY' , 'Flux from NO3 to PHY          ',  'mmole/m3/s' 
    397    pisdiabio(2)  = 'NH4PHY' , 'Flux from NH4 to PHY          ',  'mmole/m3/s' 
    398    pisdiabio(3)  = 'PHYNH4' , 'Flux from PHY to NH4          ',  'mmole/m3/s' 
    399    pisdiabio(4)  = 'PHYDOM' , 'Flux from PHY to DOM          ',  'mmole/m3/s' 
    400    pisdiabio(5)  = 'PHYZOO' , 'Flux from PHY to ZOO          ',  'mmole/m3/s' 
    401    pisdiabio(6)  = 'PHYDET' , 'Flux from PHY to DET          ',  'mmole/m3/s' 
    402    pisdiabio(7)  = 'DETZOO' , 'Flux from DET to ZOO          ',  'mmole/m3/s' 
    403    pisdiabio(8)  = 'DETSED' , 'Flux from DET to SED          ',  'mmole/m3/s' 
    404    pisdiabio(9)  = 'ZOODET' , 'Flux from ZOO to DET          ',  'mmole/m3/s' 
    405    pisdiabio(10)  = 'ZOOBOD' , 'Zooplankton closure          ',  'mmole/m3/s' 
    406    pisdiabio(11)  = 'ZOONH4' , 'Flux from ZOO to NH4         ',  'mmole/m3/s' 
    407    pisdiabio(12)  = 'ZOODOM' , 'Flux from ZOO to DOM         ',  'mmole/m3/s' 
    408    pisdiabio(13)  = 'NH4NO3' , 'Flux from NH4 to NO3         ',  'mmole/m3/s' 
    409    pisdiabio(14)  = 'DOMNH4' , 'Flux from DOM to NH4         ',  'mmole/m3/s' 
    410    pisdiabio(15)  = 'DETNH4' , 'Flux from DET to NH4         ',  'mmole/m3/s' 
    411    pisdiabio(16)  = 'DETDOM' , 'Flux from DET to DOM         ',  'mmole/m3/s' 
    412    pisdiabio(17)  = 'SEDNO3' , 'NO3 remineralization from SED',  'mmole/m3/s' 
    413 / 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/SHARED/namelist_ref

    r6497 r7403  
    288288   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    289289   ln_isf      = .false.   !  ice shelf                                 (T   => fill namsbc_isf) 
    290    ln_wave     = .false.   !  coupling with surface wave                (T => fill namsbc_wave) 
     290   ln_wave     = .false.   !  Activate coupling with wave  (T => fill namsbc_wave) 
     291   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 
     292   ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)  
     293   ln_tauoc    = .false.   !  Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave) 
     294   ln_stcor    = .false.   !  Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave) 
    291295   nn_lsm      = 0         !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 
    292296                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) 
     
    380384   sn_snd_crt    =   'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
    381385   sn_snd_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     386   sn_snd_crtw   =       'none'                 ,    'no'    ,     ''      ,         ''           , 'U,V' 
     387   sn_snd_ifrac  =       'none'                 ,    'no'    ,     ''      ,         ''           ,   '' 
     388   sn_snd_wlev   =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
    382389! receive 
    383390   sn_rcv_w10m   =   'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     
    391398   sn_rcv_cal    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    392399   sn_rcv_co2    =   'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     400   sn_rcv_hsig   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     401   sn_rcv_iceflx =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     402   sn_rcv_mslp   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     403   sn_rcv_phioc  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     404   sn_rcv_sdrfx  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     405   sn_rcv_sdrfy  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     406   sn_rcv_wper   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     407   sn_rcv_wnum   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     408   sn_rcv_wstrf  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     409   sn_rcv_wdrag  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    393410! 
    394411   nn_cplmodel   =     1   !  Maximum number of models to/from which NEMO is potentialy sending/receiving data 
     
    535552/ 
    536553!----------------------------------------------------------------------- 
    537 &namsbc_wave   ! External fields from wave model                        (ln_wave=T) 
    538 !----------------------------------------------------------------------- 
    539 !              ! file name ! frequency (hours) ! variable    ! time interp.!  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
    540 !              !           !  (if <0  months)  !   name      !  (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
    541    sn_cdg      = 'cdg_wave',        1          , 'drag_coeff',   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
    542    sn_usd      = 'sdw_wave',        1          , 'u_sd2d'    ,   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
    543    sn_vsd      = 'sdw_wave',        1          , 'v_sd2d'    ,   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
    544    sn_wn       = 'sdw_wave',        1          , 'wave_num'  ,   .true.    , .false., 'daily'   ,    ''      , ''       , '' 
    545 ! 
    546    cn_dir_cdg  = './'      !  root directory for the location of drag coefficient files 
    547    ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model 
    548    ln_sdw      = .false.   !  Computation of 3D stokes drift                
     554&namsbc_wave   ! External fields from wave model 
     555!----------------------------------------------------------------------- 
     556!              !  file name  ! frequency (hours) ! variable     ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     557!              !             !  (if <0  months)  !   name       !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      ! 
     558   sn_cdg      =  'sdw_wave' ,        1          , 'drag_coeff' ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     559   sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     560   sn_vsd      =  'sdw_wave' ,        1          , 'v_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     561   sn_swh      =  'sdw_wave' ,        1          , 'hs'         ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     562   sn_wmp      =  'sdw_wave' ,        1          , 'wmp'        ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     563   sn_wnum     =  'sdw_wave' ,        1          , 'wave_num'   ,     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     564   sn_tauoc    =  'sdw_wave' ,        1          , 'wave_stress',     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
     565! 
     566   cn_dir  = './'  !  root directory for the location of drag coefficient files 
    549567/ 
    550568!----------------------------------------------------------------------- 
     
    973991   ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping 
    974992      nn_zdfexp   =    3        ! number of sub-timestep for ln_zdfexp=T 
     993   ln_zdfqiao  = .false.   !  Enhanced wave vertical mixing Qiao (2010) (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave) 
    975994/ 
    976995!----------------------------------------------------------------------- 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/SHARED/namelist_top_ref

    r6403 r7403  
    88!!               - tracer newtonian damping              (namtrc_dmp) 
    99!!               - dynamical tracer trends               (namtrc_trd) 
    10 !!               - tracer output diagonstics             (namtrc_dia) 
    1110!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    1211!----------------------------------------------------------------------- 
     
    1413!----------------------------------------------------------------------- 
    1514   nn_dttrc      =  1        !  time step frequency for passive sn_tracers 
    16    nn_writetrc   =  5475     !  time step frequency for sn_tracer outputs 
    1715   ln_top_euler  = .false.   !  use Euler time-stepping for TOP 
    1816   ln_rsttr      = .false.   !  start from a restart file (T) or not (F) 
     
    2826&namtrc          !   tracers definition 
    2927!----------------------------------------------------------------------- 
    30    ln_trcdta     =   .true.  !  Initialisation from data input file (T) or not (F) 
     28   jp_bgc        =  0           !  Number of passive tracers of the BGC model 
     29! 
     30   ln_pisces     =  .false.     !  Run PISCES BGC model  
     31   ln_my_trc     =  .false.     !  Run MY_TRC BGC model 
     32   ln_age        =  .false.     !  Run the sea water age tracer 
     33   ln_cfc11      =  .false.     !  Run the CFC11 passive tracer 
     34   ln_cfc12      =  .false.     !  Run the CFC12 passive tracer 
     35   ln_sf6        =  .false.     !  Run the SF6 passive tracer 
     36   ln_c14        =  .false.     !  Run the Radiocarbon passive tracer 
     37! 
     38   ln_trcdta     =  .false.  !  Initialisation from data input file (T) or not (F) 
    3139   ln_trcdmp     =  .false.  !  add a damping termn (T) or not (F) 
    3240   ln_trcdmp_clo =  .false.  !  damping term (T) or not (F) on closed seas 
     41! 
     42   jp_dia3d      = 0         ! Number of 3D diagnostic variables 
     43   jp_dia2d      = 0         ! Number of 2D diagnostic variables 
     44!                !           !                                         !            !                               ! 
     45!                !    name   !           title of the field            !   units    ! initial data from file or not !  
     46!  sn_tracer(1)  = 'tracer  ' , 'Tracer  Concentration                 ',   ' - '    ,           .false. 
     47/ 
     48!----------------------------------------------------------------------- 
     49&namage         !   AGE  
     50!----------------------------------------------------------------------- 
     51   rn_age_depth      = 10            ! depth over which age tracer reset to zero 
     52   rn_age_kill_rate  = -0.000138888  !  = -1/7200 recip of relaxation timescale (s) for  age tracer shallower than age_depth 
    3353/ 
    3454!----------------------------------------------------------------------- 
     
    3656!----------------------------------------------------------------------- 
    3757   cn_dir        =  './'     !  root directory for the location of the data files 
     58!          !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     59!          !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      ! 
     60   sn_trcdta(1)  = 'data_TRC_nomask'        ,        -12        ,  'TRC'     ,    .false.   , .true. , 'yearly'  , ''       , ''   , '' 
    3861/ 
    3962!----------------------------------------------------------------------- 
     
    111134   ln_trdtrc(23) = .true. 
    112135/ 
    113 !----------------------------------------------------------------------- 
    114 &namtrc_dia      !   parameters for passive tracer additional diagnostics 
    115 !---------------------------------------------------------------------- 
    116    ln_diatrc     =  .true.   !  save additional diag. (T) or not (F) 
    117    ln_diabio     =  .true.   !  output biological trends 
    118    nn_writedia   =  5475     !  time step frequency for diagnostics 
    119    nn_writebio   =    10     !  frequency of biological outputs 
    120 / 
    121136!---------------------------------------------------------------------- 
    122137&namtrc_bc       !   data for boundary conditions 
     
    125140   cn_dir_cbc    =  './'     !  root directory for the location of COASTAL data files 
    126141   cn_dir_obc    =  './'     !  root directory for the location of OPEN data files 
     142   ln_rnf_ctl    = .false.   !  Remove runoff dilution on tracers with absent river load 
     143   rn_bc_time    =  86400.   !  Time scaling factor for SBC and CBC data (seconds in a day) 
    127144/ 
    128145!---------------------------------------------------------------------- 
  • branches/2016/dev_merge_2016/NEMOGCM/CONFIG/cfg.txt

    r6403 r7403  
    1111GYRE OPA_SRC 
    1212ORCA2_LIM_PISCES OPA_SRC LIM_SRC_2 NST_SRC TOP_SRC 
     13ORCA2_LIM3_TRC OPA_SRC LIM_SRC_3 NST_SRC TOP_SRC 
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/LIM_SRC_3/limistate.F90

    r6695 r7403  
    313313 
    314314                  IF(lwp) THEN  
    315                      WRITE(numout,*) ' ztests : ', ztests 
    316315                     IF( ztests .NE. 4 )THEN 
    317316                        WRITE(numout,*) 
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/DIA/diaar5.F90

    r6665 r7403  
    66   !! History :  3.2  !  2009-11  (S. Masson)  Original code 
    77   !!            3.3  !  2010-10  (C. Ethe, G. Madec) reorganisation of initialisation phase + merge TRC-TRA 
    8    !!---------------------------------------------------------------------- 
    9 #if defined key_diaar5 
    10    !!---------------------------------------------------------------------- 
    11    !!   'key_diaar5'  :                           activate ar5 diagnotics 
    128   !!---------------------------------------------------------------------- 
    139   !!   dia_ar5       : AR5 diagnostics 
     
    2420   USE phycst         ! physical constant 
    2521   USE in_out_manager  ! I/O manager 
     22   USE zdfddm 
     23   USE zdf_oce 
    2624 
    2725   IMPLICIT NONE 
     
    2927 
    3028   PUBLIC   dia_ar5        ! routine called in step.F90 module 
    31    PUBLIC   dia_ar5_init   ! routine called in opa.F90 module 
    3229   PUBLIC   dia_ar5_alloc  ! routine called in nemogcm.F90 module 
    33  
    34    LOGICAL, PUBLIC, PARAMETER :: lk_diaar5 = .TRUE.   ! coupled flag 
     30   PUBLIC   dia_ar5_hst    ! heat/salt transport 
    3531 
    3632   REAL(wp)                         ::   vol0         ! ocean volume (interior domain) 
     
    3935   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:  ) ::   thick0       ! ocean thickness (interior domain) 
    4036   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   sn0          ! initial salinity 
     37 
     38   LOGICAL  :: l_ar5 
    4139       
     40   !! * Substitutions 
     41#  include "zdfddm_substitute.h90" 
     42#  include "vectopt_loop_substitute.h90" 
    4243   !!---------------------------------------------------------------------- 
    4344   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
     
    7374      INTEGER  ::   ji, jj, jk                      ! dummy loop arguments 
    7475      REAL(wp) ::   zvolssh, zvol, zssh_steric, zztmp, zarho, ztemp, zsal, zmass 
     76      REAL(wp) ::   zaw, zbw, zrw 
    7577      ! 
    7678      REAL(wp), POINTER, DIMENSION(:,:)     :: zarea_ssh , zbotpres       ! 2D workspace  
     79      REAL(wp), POINTER, DIMENSION(:,:)     :: zpe                         ! 2D workspace  
    7780      REAL(wp), POINTER, DIMENSION(:,:,:)   :: zrhd , zrhop               ! 3D workspace 
    7881      REAL(wp), POINTER, DIMENSION(:,:,:,:) :: ztsn                       ! 4D workspace 
     
    8083      IF( nn_timing == 1 )   CALL timing_start('dia_ar5') 
    8184  
    82       CALL wrk_alloc( jpi , jpj              , zarea_ssh , zbotpres ) 
    83       CALL wrk_alloc( jpi , jpj , jpk        , zrhd      , zrhop    ) 
    84       CALL wrk_alloc( jpi , jpj , jpk , jpts , ztsn                 ) 
    85  
    86       zarea_ssh(:,:) = area(:,:) * sshn(:,:) 
    87  
    88       !                                         ! total volume of liquid seawater 
    89       zvolssh = SUM( zarea_ssh(:,:) )  
    90       IF( lk_mpp )   CALL mpp_sum( zvolssh ) 
    91       zvol = vol0 + zvolssh 
     85      IF( kt == nit000 )     CALL dia_ar5_init 
     86 
     87      IF( l_ar5 ) THEN 
     88         CALL wrk_alloc( jpi , jpj              , zarea_ssh , zbotpres ) 
     89         CALL wrk_alloc( jpi , jpj , jpk        , zrhd      , zrhop    ) 
     90         CALL wrk_alloc( jpi , jpj , jpk , jpts , ztsn                 ) 
     91         zarea_ssh(:,:) = area(:,:) * sshn(:,:) 
     92      ENDIF 
     93      ! 
     94      IF( iom_use( 'voltot' ) .OR. iom_use( 'sshtot' )  .OR. iom_use( 'sshdyn' )  ) THEN     
     95         !                                         ! total volume of liquid seawater 
     96         zvolssh = SUM( zarea_ssh(:,:) )  
     97         IF( lk_mpp )   CALL mpp_sum( zvolssh ) 
     98         zvol = vol0 + zvolssh 
    9299       
    93       CALL iom_put( 'voltot', zvol               ) 
    94       CALL iom_put( 'sshtot', zvolssh / area_tot ) 
    95  
    96       !                      
    97       ztsn(:,:,:,jp_tem) = tsn(:,:,:,jp_tem)                    ! thermosteric ssh 
    98       ztsn(:,:,:,jp_sal) = sn0(:,:,:) 
    99       CALL eos( ztsn, zrhd, gdept_n(:,:,:) )                       ! now in situ density using initial salinity 
    100       ! 
    101       zbotpres(:,:) = 0._wp                        ! no atmospheric surface pressure, levitating sea-ice 
    102       DO jk = 1, jpkm1 
    103          zbotpres(:,:) = zbotpres(:,:) + e3t_n(:,:,jk) * zrhd(:,:,jk) 
    104       END DO 
    105       IF( ln_linssh ) THEN 
    106          IF( ln_isfcav ) THEN 
    107             DO ji=1,jpi 
    108                DO jj=1,jpj 
    109                   zbotpres(ji,jj) = zbotpres(ji,jj) + sshn(ji,jj) * zrhd(ji,jj,mikt(ji,jj)) + riceload(ji,jj) 
    110                END DO 
    111             END DO 
    112          ELSE 
    113             zbotpres(:,:) = zbotpres(:,:) + sshn(:,:) * zrhd(:,:,1) 
    114          END IF 
     100         CALL iom_put( 'voltot', zvol               ) 
     101         CALL iom_put( 'sshtot', zvolssh / area_tot ) 
     102         CALL iom_put( 'sshdyn', sshn(:,:) - (zvolssh / area_tot) ) 
     103         ! 
     104      ENDIF 
     105 
     106      IF( iom_use( 'botpres' ) .OR. iom_use( 'sshthster' )  .OR. iom_use( 'sshsteric' )  ) THEN     
     107         !                      
     108         ztsn(:,:,:,jp_tem) = tsn(:,:,:,jp_tem)                    ! thermosteric ssh 
     109         ztsn(:,:,:,jp_sal) = sn0(:,:,:) 
     110         CALL eos( ztsn, zrhd, gdept_n(:,:,:) )                       ! now in situ density using initial salinity 
     111         ! 
     112         zbotpres(:,:) = 0._wp                        ! no atmospheric surface pressure, levitating sea-ice 
     113         DO jk = 1, jpkm1 
     114            zbotpres(:,:) = zbotpres(:,:) + e3t_n(:,:,jk) * zrhd(:,:,jk) 
     115         END DO 
     116         IF( ln_linssh ) THEN 
     117            IF( ln_isfcav ) THEN 
     118               DO ji = 1, jpi 
     119                  DO jj = 1, jpj 
     120                     zbotpres(ji,jj) = zbotpres(ji,jj) + sshn(ji,jj) * zrhd(ji,jj,mikt(ji,jj)) + riceload(ji,jj) 
     121                  END DO 
     122               END DO 
     123            ELSE 
     124               zbotpres(:,:) = zbotpres(:,:) + sshn(:,:) * zrhd(:,:,1) 
     125            END IF 
    115126!!gm 
    116127!!gm   riceload should be added in both ln_linssh=T or F, no? 
    117128!!gm 
    118       END IF 
    119       !                                          
    120       zarho = SUM( area(:,:) * zbotpres(:,:) )  
    121       IF( lk_mpp )   CALL mpp_sum( zarho ) 
    122       zssh_steric = - zarho / area_tot 
    123       CALL iom_put( 'sshthster', zssh_steric ) 
     129         END IF 
     130         !                                          
     131         zarho = SUM( area(:,:) * zbotpres(:,:) )  
     132         IF( lk_mpp )   CALL mpp_sum( zarho ) 
     133         zssh_steric = - zarho / area_tot 
     134         CALL iom_put( 'sshthster', zssh_steric ) 
    124135       
    125       !                                         ! steric sea surface height 
    126       CALL eos( tsn, zrhd, zrhop, gdept_n(:,:,:) )                 ! now in situ and potential density 
    127       zrhop(:,:,jpk) = 0._wp 
    128       CALL iom_put( 'rhop', zrhop ) 
    129       ! 
    130       zbotpres(:,:) = 0._wp                        ! no atmospheric surface pressure, levitating sea-ice 
     136         !                                         ! steric sea surface height 
     137         CALL eos( tsn, zrhd, zrhop, gdept_n(:,:,:) )                 ! now in situ and potential density 
     138         zrhop(:,:,jpk) = 0._wp 
     139         CALL iom_put( 'rhop', zrhop ) 
     140         ! 
     141         zbotpres(:,:) = 0._wp                        ! no atmospheric surface pressure, levitating sea-ice 
     142         DO jk = 1, jpkm1 
     143            zbotpres(:,:) = zbotpres(:,:) + e3t_n(:,:,jk) * zrhd(:,:,jk) 
     144         END DO 
     145         IF( ln_linssh ) THEN 
     146            IF ( ln_isfcav ) THEN 
     147               DO ji = 1,jpi 
     148                  DO jj = 1,jpj 
     149                     zbotpres(ji,jj) = zbotpres(ji,jj) + sshn(ji,jj) * zrhd(ji,jj,mikt(ji,jj)) + riceload(ji,jj) 
     150                  END DO 
     151               END DO 
     152            ELSE 
     153               zbotpres(:,:) = zbotpres(:,:) + sshn(:,:) * zrhd(:,:,1) 
     154            END IF 
     155         END IF 
     156         !     
     157         zarho = SUM( area(:,:) * zbotpres(:,:) )  
     158         IF( lk_mpp )   CALL mpp_sum( zarho ) 
     159         zssh_steric = - zarho / area_tot 
     160         CALL iom_put( 'sshsteric', zssh_steric ) 
     161       
     162         !                                         ! ocean bottom pressure 
     163         zztmp = rau0 * grav * 1.e-4_wp               ! recover pressure from pressure anomaly and cover to dbar = 1.e4 Pa 
     164         zbotpres(:,:) = zztmp * ( zbotpres(:,:) + sshn(:,:) + thick0(:,:) ) 
     165         CALL iom_put( 'botpres', zbotpres ) 
     166         ! 
     167      ENDIF 
     168 
     169      IF( iom_use( 'masstot' ) .OR. iom_use( 'temptot' )  .OR. iom_use( 'saltot' )  ) THEN     
     170         !                                         ! Mean density anomalie, temperature and salinity 
     171         ztemp = 0._wp 
     172         zsal  = 0._wp 
     173         DO jk = 1, jpkm1 
     174            DO jj = 1, jpj 
     175               DO ji = 1, jpi 
     176                  zztmp = area(ji,jj) * e3t_n(ji,jj,jk) 
     177                  ztemp = ztemp + zztmp * tsn(ji,jj,jk,jp_tem) 
     178                  zsal  = zsal  + zztmp * tsn(ji,jj,jk,jp_sal) 
     179               END DO 
     180            END DO 
     181         END DO 
     182         IF( ln_linssh ) THEN 
     183            IF( ln_isfcav ) THEN 
     184               DO ji = 1, jpi 
     185                  DO jj = 1, jpj 
     186                     ztemp = ztemp + zarea_ssh(ji,jj) * tsn(ji,jj,mikt(ji,jj),jp_tem)  
     187                     zsal  = zsal  + zarea_ssh(ji,jj) * tsn(ji,jj,mikt(ji,jj),jp_sal)  
     188                  END DO 
     189               END DO 
     190            ELSE 
     191               ztemp = ztemp + SUM( zarea_ssh(:,:) * tsn(:,:,1,jp_tem) ) 
     192               zsal  = zsal  + SUM( zarea_ssh(:,:) * tsn(:,:,1,jp_sal) ) 
     193            END IF 
     194         ENDIF 
     195         IF( lk_mpp ) THEN   
     196            CALL mpp_sum( ztemp ) 
     197            CALL mpp_sum( zsal  ) 
     198         END IF 
     199         ! 
     200         zmass = rau0 * ( zarho + zvol )                 ! total mass of liquid seawater 
     201         ztemp = ztemp / zvol                            ! potential temperature in liquid seawater 
     202         zsal  = zsal  / zvol                            ! Salinity of liquid seawater 
     203         ! 
     204         CALL iom_put( 'masstot', zmass ) 
     205         CALL iom_put( 'temptot', ztemp ) 
     206         CALL iom_put( 'saltot' , zsal  ) 
     207         ! 
     208      ENDIF 
     209 
     210      IF( iom_use( 'tnpeo' )) THEN     
     211      ! Work done against stratification by vertical mixing 
     212      ! Exclude points where rn2 is negative as convection kicks in here and 
     213      ! work is not being done against stratification 
     214          CALL wrk_alloc( jpi, jpj, zpe ) 
     215          zpe(:,:) = 0._wp 
     216          IF( lk_zdfddm ) THEN 
     217             DO ji=1,jpi 
     218                DO jj=1,jpj 
     219                   DO jk=1,jpk 
     220                      zrw =   ( gdepw_n(ji,jj,jk  ) - gdept_n(ji,jj,jk) )   & 
     221                         &  / ( gdept_n(ji,jj,jk-1) - gdept_n(ji,jj,jk) ) 
     222                      ! 
     223                      zaw = rab_n(ji,jj,jk,jp_tem) * (1. - zrw) + rab_n(ji,jj,jk-1,jp_tem)* zrw 
     224                      zbw = rab_n(ji,jj,jk,jp_sal) * (1. - zrw) + rab_n(ji,jj,jk-1,jp_sal)* zrw 
     225                      ! 
     226                      zpe(ji, jj) = zpe(ji, jj) - MIN(0._wp, rn2(ji,jj,jk)) * & 
     227                           &       grav * (avt(ji,jj,jk) * zaw * (tsn(ji,jj,jk-1,jp_tem) - tsn(ji,jj,jk,jp_tem) )  & 
     228                           &       - fsavs(ji,jj,jk) * zbw * (tsn(ji,jj,jk-1,jp_sal) - tsn(ji,jj,jk,jp_sal) ) ) 
     229 
     230                   ENDDO 
     231                ENDDO 
     232             ENDDO 
     233          ELSE 
     234             DO ji = 1, jpi 
     235                DO jj = 1, jpj 
     236                   DO jk = 1, jpk 
     237                       zpe(ji,jj) = zpe(ji,jj) + avt(ji, jj, jk) * MIN(0._wp,rn2(ji, jj, jk)) * rau0 * e3w_n(ji, jj, jk) 
     238                   ENDDO 
     239                ENDDO 
     240             ENDDO 
     241          ENDIF 
     242          CALL lbc_lnk( zpe, 'T', 1._wp)          
     243          CALL iom_put( 'tnpeo', zpe ) 
     244          CALL wrk_dealloc( jpi, jpj, zpe ) 
     245      ENDIF 
     246      ! 
     247      IF( l_ar5 ) THEN 
     248        CALL wrk_dealloc( jpi , jpj              , zarea_ssh , zbotpres ) 
     249        CALL wrk_dealloc( jpi , jpj , jpk        , zrhd      , zrhop    ) 
     250        CALL wrk_dealloc( jpi , jpj , jpk , jpts , ztsn                 ) 
     251      ENDIF 
     252      ! 
     253      IF( nn_timing == 1 )   CALL timing_stop('dia_ar5') 
     254      ! 
     255   END SUBROUTINE dia_ar5 
     256 
     257   SUBROUTINE dia_ar5_hst( ktra, cptr, pua, pva )  
     258      !!---------------------------------------------------------------------- 
     259      !!                    ***  ROUTINE dia_ar5_htr *** 
     260      !!---------------------------------------------------------------------- 
     261      !! Wrapper for heat transport calculations 
     262      !! Called from all advection and/or diffusion routines 
     263      !!---------------------------------------------------------------------- 
     264      INTEGER                         , INTENT(in )  :: ktra  ! tracer index 
     265      CHARACTER(len=3)                , INTENT(in)   :: cptr  ! transport type  'adv'/'ldf' 
     266      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in)   :: pua   ! 3D input array of advection/diffusion 
     267      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in)   :: pva   ! 3D input array of advection/diffusion 
     268      ! 
     269      INTEGER    ::  ji, jj, jk 
     270      REAL(wp), POINTER, DIMENSION(:,:)  :: z2d 
     271 
     272     
     273 
     274      CALL wrk_alloc( jpi, jpj, z2d ) 
     275      z2d(:,:) = pua(:,:,1)  
    131276      DO jk = 1, jpkm1 
    132          zbotpres(:,:) = zbotpres(:,:) + e3t_n(:,:,jk) * zrhd(:,:,jk) 
    133       END DO 
    134       IF( ln_linssh ) THEN 
    135          IF ( ln_isfcav ) THEN 
    136             DO ji=1,jpi 
    137                DO jj=1,jpj 
    138                   zbotpres(ji,jj) = zbotpres(ji,jj) + sshn(ji,jj) * zrhd(ji,jj,mikt(ji,jj)) + riceload(ji,jj) 
    139                END DO 
     277         DO jj = 2, jpjm1 
     278            DO ji = fs_2, fs_jpim1   ! vector opt. 
     279               z2d(ji,jj) = z2d(ji,jj) + pua(ji,jj,jk)  
    140280            END DO 
    141          ELSE 
    142             zbotpres(:,:) = zbotpres(:,:) + sshn(:,:) * zrhd(:,:,1) 
    143          END IF 
    144       END IF 
    145       !     
    146       zarho = SUM( area(:,:) * zbotpres(:,:) )  
    147       IF( lk_mpp )   CALL mpp_sum( zarho ) 
    148       zssh_steric = - zarho / area_tot 
    149       CALL iom_put( 'sshsteric', zssh_steric ) 
    150        
    151       !                                         ! ocean bottom pressure 
    152       zztmp = rau0 * grav * 1.e-4_wp               ! recover pressure from pressure anomaly and cover to dbar = 1.e4 Pa 
    153       zbotpres(:,:) = zztmp * ( zbotpres(:,:) + sshn(:,:) + thick0(:,:) ) 
    154       CALL iom_put( 'botpres', zbotpres ) 
    155  
    156       !                                         ! Mean density anomalie, temperature and salinity 
    157       ztemp = 0._wp 
    158       zsal  = 0._wp 
    159       DO jk = 1, jpkm1 
    160          DO jj = 1, jpj 
    161             DO ji = 1, jpi 
    162                zztmp = area(ji,jj) * e3t_n(ji,jj,jk) 
    163                ztemp = ztemp + zztmp * tsn(ji,jj,jk,jp_tem) 
    164                zsal  = zsal  + zztmp * tsn(ji,jj,jk,jp_sal) 
    165             END DO 
    166          END DO 
    167       END DO 
    168       IF( ln_linssh ) THEN 
    169          IF( ln_isfcav ) THEN 
    170             DO ji=1,jpi 
    171                DO jj=1,jpj 
    172                   ztemp = ztemp + zarea_ssh(ji,jj) * tsn(ji,jj,mikt(ji,jj),jp_tem)  
    173                   zsal  = zsal  + zarea_ssh(ji,jj) * tsn(ji,jj,mikt(ji,jj),jp_sal)  
    174                END DO 
    175             END DO 
    176          ELSE 
    177             ztemp = ztemp + SUM( zarea_ssh(:,:) * tsn(:,:,1,jp_tem) ) 
    178             zsal  = zsal  + SUM( zarea_ssh(:,:) * tsn(:,:,1,jp_sal) ) 
    179          END IF 
    180       ENDIF 
    181       IF( lk_mpp ) THEN   
    182          CALL mpp_sum( ztemp ) 
    183          CALL mpp_sum( zsal  ) 
    184       END IF 
    185       ! 
    186       zmass = rau0 * ( zarho + zvol )                 ! total mass of liquid seawater 
    187       ztemp = ztemp / zvol                            ! potential temperature in liquid seawater 
    188       zsal  = zsal  / zvol                            ! Salinity of liquid seawater 
    189       ! 
    190       CALL iom_put( 'masstot', zmass ) 
    191       CALL iom_put( 'temptot', ztemp ) 
    192       CALL iom_put( 'saltot' , zsal  ) 
    193       ! 
    194       CALL wrk_dealloc( jpi , jpj              , zarea_ssh , zbotpres ) 
    195       CALL wrk_dealloc( jpi , jpj , jpk        , zrhd      , zrhop    ) 
    196       CALL wrk_dealloc( jpi , jpj , jpk , jpts , ztsn                 ) 
    197       ! 
    198       IF( nn_timing == 1 )   CALL timing_stop('dia_ar5') 
    199       ! 
    200    END SUBROUTINE dia_ar5 
     281         END DO 
     282       END DO 
     283       CALL lbc_lnk( z2d, 'U', -1. ) 
     284       IF( cptr == 'adv' ) THEN 
     285          IF( ktra == jp_tem ) CALL iom_put( "uadv_heattr" , rau0_rcp * z2d )  ! advective heat transport in i-direction 
     286          IF( ktra == jp_sal ) CALL iom_put( "uadv_salttr" , rau0     * z2d )  ! advective salt transport in i-direction 
     287       ENDIF 
     288       IF( cptr == 'ldf' ) THEN 
     289          IF( ktra == jp_tem ) CALL iom_put( "udiff_heattr" , rau0_rcp * z2d ) ! diffusive heat transport in i-direction 
     290          IF( ktra == jp_sal ) CALL iom_put( "udiff_salttr" , rau0     * z2d ) ! diffusive salt transport in i-direction 
     291       ENDIF 
     292       ! 
     293       z2d(:,:) = pva(:,:,1)  
     294       DO jk = 1, jpkm1 
     295          DO jj = 2, jpjm1 
     296             DO ji = fs_2, fs_jpim1   ! vector opt. 
     297                z2d(ji,jj) = z2d(ji,jj) + pva(ji,jj,jk)  
     298             END DO 
     299          END DO 
     300       END DO 
     301       CALL lbc_lnk( z2d, 'V', -1. ) 
     302       IF( cptr == 'adv' ) THEN 
     303          IF( ktra == jp_tem ) CALL iom_put( "vadv_heattr" , rau0_rcp * z2d )  ! advective heat transport in j-direction 
     304          IF( ktra == jp_sal ) CALL iom_put( "vadv_salttr" , rau0     * z2d )  ! advective salt transport in j-direction 
     305       ENDIF 
     306       IF( cptr == 'ldf' ) THEN 
     307          IF( ktra == jp_tem ) CALL iom_put( "vdiff_heattr" , rau0_rcp * z2d ) ! diffusive heat transport in j-direction 
     308          IF( ktra == jp_sal ) CALL iom_put( "vdiff_salttr" , rau0     * z2d ) ! diffusive salt transport in j-direction 
     309       ENDIF 
     310           
     311       CALL wrk_dealloc( jpi, jpj, z2d ) 
     312 
     313   END SUBROUTINE dia_ar5_hst 
    201314 
    202315 
     
    217330      IF( nn_timing == 1 )   CALL timing_start('dia_ar5_init') 
    218331      ! 
    219       CALL wrk_alloc( jpi , jpj , jpk, jpts, zsaldta ) 
    220       !                                      ! allocate dia_ar5 arrays 
    221       IF( dia_ar5_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'dia_ar5_init : unable to allocate arrays' ) 
    222  
    223       area(:,:) = e1e2t(:,:) * tmask_i(:,:) 
    224  
    225       area_tot = SUM( area(:,:) )   ;   IF( lk_mpp )   CALL mpp_sum( area_tot ) 
    226  
    227       vol0        = 0._wp 
    228       thick0(:,:) = 0._wp 
    229       DO jk = 1, jpkm1 
    230          vol0        = vol0        + SUM( area (:,:) * tmask(:,:,jk) * e3t_0(:,:,jk) ) 
    231          thick0(:,:) = thick0(:,:) +    tmask_i(:,:) * tmask(:,:,jk) * e3t_0(:,:,jk) 
    232       END DO 
    233       IF( lk_mpp )   CALL mpp_sum( vol0 ) 
    234  
    235  
    236       CALL iom_open ( 'sali_ref_clim_monthly', inum ) 
    237       CALL iom_get  ( inum, jpdom_data, 'vosaline' , zsaldta(:,:,:,1), 1  ) 
    238       CALL iom_get  ( inum, jpdom_data, 'vosaline' , zsaldta(:,:,:,2), 12 ) 
    239       CALL iom_close( inum ) 
    240  
    241       sn0(:,:,:) = 0.5_wp * ( zsaldta(:,:,:,1) + zsaldta(:,:,:,2) )         
    242       sn0(:,:,:) = sn0(:,:,:) * tmask(:,:,:) 
    243       IF( ln_zps ) THEN               ! z-coord. partial steps 
    244          DO jj = 1, jpj               ! interpolation of salinity at the last ocean level (i.e. the partial step) 
    245             DO ji = 1, jpi 
    246                ik = mbkt(ji,jj) 
    247                IF( ik > 1 ) THEN 
    248                   zztmp = ( gdept_1d(ik) - gdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 
    249                   sn0(ji,jj,ik) = ( 1._wp - zztmp ) * sn0(ji,jj,ik) + zztmp * sn0(ji,jj,ik-1) 
    250                ENDIF 
     332      l_ar5 = .FALSE. 
     333      IF(   iom_use( 'voltot'  ) .OR. iom_use( 'sshtot'    )  .OR. iom_use( 'sshdyn' )  .OR.  &  
     334         &  iom_use( 'masstot' ) .OR. iom_use( 'temptot'   )  .OR. iom_use( 'saltot' ) .OR.  &     
     335         &  iom_use( 'botpres' ) .OR. iom_use( 'sshthster' )  .OR. iom_use( 'sshsteric' )  ) L_ar5 = .TRUE. 
     336   
     337      IF( l_ar5 ) THEN 
     338         ! 
     339         CALL wrk_alloc( jpi , jpj , jpk, jpts, zsaldta ) 
     340         !                                      ! allocate dia_ar5 arrays 
     341         IF( dia_ar5_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'dia_ar5_init : unable to allocate arrays' ) 
     342 
     343         area(:,:) = e1e2t(:,:) * tmask_i(:,:) 
     344 
     345         area_tot = SUM( area(:,:) )   ;   IF( lk_mpp )   CALL mpp_sum( area_tot ) 
     346 
     347         vol0        = 0._wp 
     348         thick0(:,:) = 0._wp 
     349         DO jk = 1, jpkm1 
     350            vol0        = vol0        + SUM( area (:,:) * tmask(:,:,jk) * e3t_0(:,:,jk) ) 
     351            thick0(:,:) = thick0(:,:) +    tmask_i(:,:) * tmask(:,:,jk) * e3t_0(:,:,jk) 
     352         END DO 
     353         IF( lk_mpp )   CALL mpp_sum( vol0 ) 
     354 
     355 
     356         CALL iom_open ( 'sali_ref_clim_monthly', inum ) 
     357         CALL iom_get  ( inum, jpdom_data, 'vosaline' , zsaldta(:,:,:,1), 1  ) 
     358         CALL iom_get  ( inum, jpdom_data, 'vosaline' , zsaldta(:,:,:,2), 12 ) 
     359         CALL iom_close( inum ) 
     360 
     361         sn0(:,:,:) = 0.5_wp * ( zsaldta(:,:,:,1) + zsaldta(:,:,:,2) )         
     362         sn0(:,:,:) = sn0(:,:,:) * tmask(:,:,:) 
     363         IF( ln_zps ) THEN               ! z-coord. partial steps 
     364            DO jj = 1, jpj               ! interpolation of salinity at the last ocean level (i.e. the partial step) 
     365               DO ji = 1, jpi 
     366                  ik = mbkt(ji,jj) 
     367                  IF( ik > 1 ) THEN 
     368                     zztmp = ( gdept_1d(ik) - gdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 
     369                     sn0(ji,jj,ik) = ( 1._wp - zztmp ) * sn0(ji,jj,ik) + zztmp * sn0(ji,jj,ik-1) 
     370                  ENDIF 
     371               END DO 
    251372            END DO 
    252          END DO 
    253       ENDIF 
    254       ! 
    255       CALL wrk_dealloc( jpi , jpj , jpk, jpts, zsaldta ) 
     373         ENDIF 
     374         ! 
     375         CALL wrk_dealloc( jpi , jpj , jpk, jpts, zsaldta ) 
     376         ! 
     377      ENDIF 
    256378      ! 
    257379      IF( nn_timing == 1 )   CALL timing_stop('dia_ar5_init') 
    258380      ! 
    259381   END SUBROUTINE dia_ar5_init 
    260  
    261 #else 
    262    !!---------------------------------------------------------------------- 
    263    !!   Default option :                                         NO diaar5 
    264    !!---------------------------------------------------------------------- 
    265    LOGICAL, PUBLIC, PARAMETER :: lk_diaar5 = .FALSE.   ! coupled flag 
    266 CONTAINS 
    267    SUBROUTINE dia_ar5_init    ! Dummy routine 
    268    END SUBROUTINE dia_ar5_init 
    269    SUBROUTINE dia_ar5( kt )   ! Empty routine 
    270       INTEGER ::   kt 
    271       WRITE(*,*) 'dia_ar5: You should not have seen this print! error?', kt 
    272    END SUBROUTINE dia_ar5 
    273 #endif 
    274382 
    275383   !!====================================================================== 
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/DIA/diaptr.F90

    r6140 r7403  
    99   !!            3.3  ! 2010-10  (G. Madec)  dynamical allocation 
    1010   !!            3.6  ! 2014-12  (C. Ethe) use of IOM 
     11   !!            3.6  ! 2016-06  (T. Graham) Addition of diagnostics for CMIP6 
    1112   !!---------------------------------------------------------------------- 
    1213 
     
    3839   PUBLIC   dia_ptr_init   ! call in step module 
    3940   PUBLIC   dia_ptr        ! call in step module 
     41   PUBLIC   dia_ptr_hst    ! called from tra_ldf/tra_adv routines 
    4042 
    4143   !                                  !!** namelist  namptr  ** 
    42    REAL(wp), ALLOCATABLE, SAVE, PUBLIC, DIMENSION(:) ::   htr_adv, htr_ldf   !: Heat TRansports (adv, diff, overturn.) 
    43    REAL(wp), ALLOCATABLE, SAVE, PUBLIC, DIMENSION(:) ::   str_adv, str_ldf   !: Salt TRansports (adv, diff, overturn.) 
    44     
     44   REAL(wp), ALLOCATABLE, SAVE, PUBLIC, DIMENSION(:,:) ::   htr_adv, htr_ldf, htr_eiv   !: Heat TRansports (adv, diff, Bolus.) 
     45   REAL(wp), ALLOCATABLE, SAVE, PUBLIC, DIMENSION(:,:) ::   str_adv, str_ldf, str_eiv   !: Salt TRansports (adv, diff, Bolus.) 
     46   REAL(wp), ALLOCATABLE, SAVE, PUBLIC, DIMENSION(:,:) ::   htr_ove, str_ove   !: heat Salt TRansports ( overturn.) 
     47   REAL(wp), ALLOCATABLE, SAVE, PUBLIC, DIMENSION(:,:) ::   htr_btr, str_btr   !: heat Salt TRansports ( barotropic ) 
    4548 
    4649   LOGICAL, PUBLIC ::   ln_diaptr   !  Poleward transport flag (T) or not (F) 
    4750   LOGICAL, PUBLIC ::   ln_subbas   !  Atlantic/Pacific/Indian basins calculation 
    48    INTEGER        ::   nptr        ! = 1 (l_subbas=F) or = 5 (glo, atl, pac, ind, ipc) (l_subbas=T)  
     51   INTEGER, PUBLIC ::   nptr        ! = 1 (l_subbas=F) or = 5 (glo, atl, pac, ind, ipc) (l_subbas=T)  
    4952 
    5053   REAL(wp) ::   rc_sv    = 1.e-6_wp   ! conversion from m3/s to Sverdrup 
     
    7578      ! 
    7679      INTEGER  ::   ji, jj, jk, jn   ! dummy loop indices 
    77       REAL(wp) ::   zv, zsfc               ! local scalar 
     80      REAL(wp) ::   zsfc,zvfc               ! local scalar 
    7881      REAL(wp), DIMENSION(jpi,jpj)     ::  z2d   ! 2D workspace 
    7982      REAL(wp), DIMENSION(jpi,jpj,jpk) ::  z3d   ! 3D workspace 
    8083      REAL(wp), DIMENSION(jpi,jpj,jpk) ::  zmask   ! 3D workspace 
    8184      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts) ::  zts   ! 3D workspace 
    82       CHARACTER( len = 10 )  :: cl1 
     85      REAL(wp), DIMENSION(jpj)     ::  vsum   ! 1D workspace 
     86      REAL(wp), DIMENSION(jpj,jpts)     ::  tssum   ! 1D workspace 
     87  
     88      ! 
     89      !overturning calculation 
     90      REAL(wp), DIMENSION(jpj,jpk,nptr) ::   sjk  , r1_sjk ! i-mean i-k-surface and its inverse 
     91      REAL(wp), DIMENSION(jpj,jpk,nptr) ::   v_msf, sn_jk  , tn_jk ! i-mean T and S, j-Stream-Function 
     92      REAL(wp), DIMENSION(jpi,jpj,jpk) ::  zvn   ! 3D workspace 
     93 
     94 
     95      CHARACTER( len = 12 )  :: cl1 
    8396      !!---------------------------------------------------------------------- 
    8497      ! 
     
    109122            END DO 
    110123         ENDIF 
     124         IF( iom_use("sopstove") .OR. iom_use("sophtove") .OR. iom_use("sopstbtr") .OR. iom_use("sophtbtr") ) THEN 
     125            ! define fields multiplied by scalar 
     126            zmask(:,:,:) = 0._wp 
     127            zts(:,:,:,:) = 0._wp 
     128            zvn(:,:,:) = 0._wp 
     129            DO jk = 1, jpkm1 
     130               DO jj = 1, jpjm1 
     131                  DO ji = 1, jpi 
     132                     zvfc = e1v(ji,jj) * e3v_n(ji,jj,jk) 
     133                     zmask(ji,jj,jk)      = vmask(ji,jj,jk)      * zvfc 
     134                     zts(ji,jj,jk,jp_tem) = (tsn(ji,jj,jk,jp_tem)+tsn(ji,jj+1,jk,jp_tem)) * 0.5 * zvfc  !Tracers averaged onto V grid 
     135                     zts(ji,jj,jk,jp_sal) = (tsn(ji,jj,jk,jp_sal)+tsn(ji,jj+1,jk,jp_sal)) * 0.5 * zvfc 
     136                     zvn(ji,jj,jk)        = vn(ji,jj,jk)         * zvfc 
     137                  ENDDO 
     138               ENDDO 
     139             ENDDO 
     140         ENDIF 
     141         IF( iom_use("sopstove") .OR. iom_use("sophtove") ) THEN 
     142             sjk(:,:,1) = ptr_sjk( zmask(:,:,:), btmsk(:,:,1) ) 
     143             r1_sjk(:,:,1) = 0._wp 
     144             WHERE( sjk(:,:,1) /= 0._wp )   r1_sjk(:,:,1) = 1._wp / sjk(:,:,1) 
     145 
     146             ! i-mean T and S, j-Stream-Function, global 
     147             tn_jk(:,:,1) = ptr_sjk( zts(:,:,:,jp_tem) ) * r1_sjk(:,:,1) 
     148             sn_jk(:,:,1) = ptr_sjk( zts(:,:,:,jp_sal) ) * r1_sjk(:,:,1) 
     149             v_msf(:,:,1) = ptr_sjk( zvn(:,:,:) ) 
     150 
     151             htr_ove(:,1) = SUM( v_msf(:,:,1)*tn_jk(:,:,1) ,2 ) 
     152             str_ove(:,1) = SUM( v_msf(:,:,1)*sn_jk(:,:,1) ,2 ) 
     153 
     154             z2d(1,:) = htr_ove(:,1) * rc_pwatt        !  (conversion in PW) 
     155             DO ji = 1, jpi 
     156               z2d(ji,:) = z2d(1,:) 
     157             ENDDO 
     158             cl1 = 'sophtove' 
     159             CALL iom_put( TRIM(cl1), z2d ) 
     160             z2d(1,:) = str_ove(:,1) * rc_ggram        !  (conversion in Gg) 
     161             DO ji = 1, jpi 
     162               z2d(ji,:) = z2d(1,:) 
     163             ENDDO 
     164             cl1 = 'sopstove' 
     165             CALL iom_put( TRIM(cl1), z2d ) 
     166             IF( ln_subbas ) THEN 
     167                DO jn = 2, nptr 
     168                    sjk(:,:,jn) = ptr_sjk( zmask(:,:,:), btmsk(:,:,jn) ) 
     169                    r1_sjk(:,:,jn) = 0._wp 
     170                    WHERE( sjk(:,:,jn) /= 0._wp )   r1_sjk(:,:,jn) = 1._wp / sjk(:,:,jn) 
     171 
     172                    ! i-mean T and S, j-Stream-Function, basin 
     173                    tn_jk(:,:,jn) = ptr_sjk( zts(:,:,:,jp_tem), btmsk(:,:,jn) ) * r1_sjk(:,:,jn) 
     174                    sn_jk(:,:,jn) = ptr_sjk( zts(:,:,:,jp_sal), btmsk(:,:,jn) ) * r1_sjk(:,:,jn) 
     175                    v_msf(:,:,jn) = ptr_sjk( zvn(:,:,:), btmsk(:,:,jn) )  
     176                    htr_ove(:,jn) = SUM( v_msf(:,:,jn)*tn_jk(:,:,jn) ,2 ) 
     177                    str_ove(:,jn) = SUM( v_msf(:,:,jn)*sn_jk(:,:,jn) ,2 ) 
     178 
     179                    z2d(1,:) = htr_ove(:,jn) * rc_pwatt !  (conversion in PW) 
     180                    DO ji = 1, jpi 
     181                        z2d(ji,:) = z2d(1,:) 
     182                    ENDDO 
     183                    cl1 = TRIM('sophtove_'//clsubb(jn)) 
     184                    CALL iom_put( cl1, z2d ) 
     185                    z2d(1,:) = str_ove(:,jn) * rc_ggram        ! (conversion in Gg) 
     186                    DO ji = 1, jpi 
     187                        z2d(ji,:) = z2d(1,:) 
     188                    ENDDO 
     189                    cl1 = TRIM('sopstove_'//clsubb(jn)) 
     190                    CALL iom_put( cl1, z2d ) 
     191                END DO 
     192             ENDIF 
     193         ENDIF 
     194         IF( iom_use("sopstbtr") .OR. iom_use("sophtbtr") ) THEN 
     195         ! Calculate barotropic heat and salt transport here  
     196             sjk(:,1,1) = ptr_sj( zmask(:,:,:), btmsk(:,:,1) ) 
     197             r1_sjk(:,1,1) = 0._wp 
     198             WHERE( sjk(:,1,1) /= 0._wp )   r1_sjk(:,1,1) = 1._wp / sjk(:,1,1) 
     199             
     200            vsum = ptr_sj( zvn(:,:,:), btmsk(:,:,1)) 
     201            tssum(:,jp_tem) = ptr_sj( zts(:,:,:,jp_tem), btmsk(:,:,1) ) 
     202            tssum(:,jp_sal) = ptr_sj( zts(:,:,:,jp_sal), btmsk(:,:,1) ) 
     203            htr_btr(:,1) = vsum * tssum(:,jp_tem) * r1_sjk(:,1,1) 
     204            str_btr(:,1) = vsum * tssum(:,jp_sal) * r1_sjk(:,1,1) 
     205            z2d(1,:) = htr_btr(:,1) * rc_pwatt        !  (conversion in PW) 
     206            DO ji = 2, jpi 
     207               z2d(ji,:) = z2d(1,:) 
     208            ENDDO 
     209            cl1 = 'sophtbtr' 
     210            CALL iom_put( TRIM(cl1), z2d ) 
     211            z2d(1,:) = str_btr(:,1) * rc_ggram        !  (conversion in Gg) 
     212            DO ji = 2, jpi 
     213              z2d(ji,:) = z2d(1,:) 
     214            ENDDO 
     215            cl1 = 'sopstbtr' 
     216            CALL iom_put( TRIM(cl1), z2d ) 
     217            IF( ln_subbas ) THEN 
     218                DO jn = 2, nptr 
     219                    sjk(:,1,jn) = ptr_sj( zmask(:,:,:), btmsk(:,:,jn) ) 
     220                    r1_sjk(:,1,jn) = 0._wp 
     221                    WHERE( sjk(:,1,jn) /= 0._wp )   r1_sjk(:,1,jn) = 1._wp / sjk(:,1,jn) 
     222                    vsum = ptr_sj( zvn(:,:,:), btmsk(:,:,jn)) 
     223                    tssum(:,jp_tem) = ptr_sj( zts(:,:,:,jp_tem), btmsk(:,:,jn) ) 
     224                    tssum(:,jp_sal) = ptr_sj( zts(:,:,:,jp_sal), btmsk(:,:,jn) ) 
     225                    htr_btr(:,jn) = vsum * tssum(:,jp_tem) * r1_sjk(:,1,jn) 
     226                    str_btr(:,jn) = vsum * tssum(:,jp_sal) * r1_sjk(:,1,jn) 
     227                    z2d(1,:) = htr_btr(:,jn) * rc_pwatt !  (conversion in PW) 
     228                    DO ji = 1, jpi 
     229                        z2d(ji,:) = z2d(1,:) 
     230                    ENDDO 
     231                    cl1 = TRIM('sophtbtr_'//clsubb(jn)) 
     232                    CALL iom_put( cl1, z2d ) 
     233                    z2d(1,:) = str_btr(:,jn) * rc_ggram        ! (conversion in Gg) 
     234                    DO ji = 1, jpi 
     235                        z2d(ji,:) = z2d(1,:) 
     236                    ENDDO 
     237                    cl1 = TRIM('sopstbtr_'//clsubb(jn)) 
     238                    CALL iom_put( cl1, z2d ) 
     239               ENDDO 
     240            ENDIF !ln_subbas 
     241         ENDIF !iom_use("sopstbtr....) 
    111242         ! 
    112243      ELSE 
     
    148279         !                                ! Advective and diffusive heat and salt transport 
    149280         IF( iom_use("sophtadv") .OR. iom_use("sopstadv") ) THEN    
    150             z2d(1,:) = htr_adv(:) * rc_pwatt        !  (conversion in PW) 
     281            z2d(1,:) = htr_adv(:,1) * rc_pwatt        !  (conversion in PW) 
    151282            DO ji = 1, jpi 
    152283               z2d(ji,:) = z2d(1,:) 
     
    154285            cl1 = 'sophtadv'                  
    155286            CALL iom_put( TRIM(cl1), z2d ) 
    156             z2d(1,:) = str_adv(:) * rc_ggram        ! (conversion in Gg) 
     287            z2d(1,:) = str_adv(:,1) * rc_ggram        ! (conversion in Gg) 
    157288            DO ji = 1, jpi 
    158289               z2d(ji,:) = z2d(1,:) 
     
    160291            cl1 = 'sopstadv' 
    161292            CALL iom_put( TRIM(cl1), z2d ) 
     293            IF( ln_subbas ) THEN 
     294              DO jn=2,nptr 
     295               z2d(1,:) = htr_adv(:,jn) * rc_pwatt        !  (conversion in PW) 
     296               DO ji = 1, jpi 
     297                 z2d(ji,:) = z2d(1,:) 
     298               ENDDO 
     299               cl1 = TRIM('sophtadv_'//clsubb(jn))                  
     300               CALL iom_put( cl1, z2d ) 
     301               z2d(1,:) = str_adv(:,jn) * rc_ggram        ! (conversion in Gg) 
     302               DO ji = 1, jpi 
     303                  z2d(ji,:) = z2d(1,:) 
     304               ENDDO 
     305               cl1 = TRIM('sopstadv_'//clsubb(jn))                  
     306               CALL iom_put( cl1, z2d )               
     307              ENDDO 
     308            ENDIF 
    162309         ENDIF 
    163310         ! 
    164311         IF( iom_use("sophtldf") .OR. iom_use("sopstldf") ) THEN    
    165             z2d(1,:) = htr_ldf(:) * rc_pwatt        !  (conversion in PW)  
     312            z2d(1,:) = htr_ldf(:,1) * rc_pwatt        !  (conversion in PW)  
    166313            DO ji = 1, jpi 
    167314               z2d(ji,:) = z2d(1,:) 
     
    169316            cl1 = 'sophtldf' 
    170317            CALL iom_put( TRIM(cl1), z2d ) 
    171             z2d(1,:) = str_ldf(:) * rc_ggram        !  (conversion in Gg) 
     318            z2d(1,:) = str_ldf(:,1) * rc_ggram        !  (conversion in Gg) 
    172319            DO ji = 1, jpi 
    173320               z2d(ji,:) = z2d(1,:) 
     
    175322            cl1 = 'sopstldf' 
    176323            CALL iom_put( TRIM(cl1), z2d ) 
     324            IF( ln_subbas ) THEN 
     325              DO jn=2,nptr 
     326               z2d(1,:) = htr_ldf(:,jn) * rc_pwatt        !  (conversion in PW) 
     327               DO ji = 1, jpi 
     328                 z2d(ji,:) = z2d(1,:) 
     329               ENDDO 
     330               cl1 = TRIM('sophtldf_'//clsubb(jn))                  
     331               CALL iom_put( cl1, z2d ) 
     332               z2d(1,:) = str_ldf(:,jn) * rc_ggram        ! (conversion in Gg) 
     333               DO ji = 1, jpi 
     334                  z2d(ji,:) = z2d(1,:) 
     335               ENDDO 
     336               cl1 = TRIM('sopstldf_'//clsubb(jn))                  
     337               CALL iom_put( cl1, z2d )               
     338              ENDDO 
     339            ENDIF 
     340         ENDIF 
     341 
     342         IF( iom_use("sophteiv") .OR. iom_use("sopsteiv") ) THEN  
     343            z2d(1,:) = htr_eiv(:,1) * rc_pwatt        !  (conversion in PW)  
     344            DO ji = 1, jpi 
     345               z2d(ji,:) = z2d(1,:) 
     346            ENDDO 
     347            cl1 = 'sophteiv' 
     348            CALL iom_put( TRIM(cl1), z2d ) 
     349            z2d(1,:) = str_eiv(:,1) * rc_ggram        !  (conversion in Gg) 
     350            DO ji = 1, jpi 
     351               z2d(ji,:) = z2d(1,:) 
     352            ENDDO 
     353            cl1 = 'sopsteiv' 
     354            CALL iom_put( TRIM(cl1), z2d ) 
     355            IF( ln_subbas ) THEN 
     356               DO jn=2,nptr 
     357                  z2d(1,:) = htr_eiv(:,jn) * rc_pwatt        !  (conversion in PW) 
     358                  DO ji = 1, jpi 
     359                     z2d(ji,:) = z2d(1,:) 
     360                  ENDDO 
     361                  cl1 = TRIM('sophteiv_'//clsubb(jn))                  
     362                  CALL iom_put( cl1, z2d ) 
     363                  z2d(1,:) = str_eiv(:,jn) * rc_ggram        ! (conversion in Gg) 
     364                  DO ji = 1, jpi 
     365                     z2d(ji,:) = z2d(1,:) 
     366                  ENDDO 
     367                  cl1 = TRIM('sopsteiv_'//clsubb(jn))  
     368                  CALL iom_put( cl1, z2d )               
     369               ENDDO 
     370            ENDIF 
    177371         ENDIF 
    178372         ! 
     
    254448         ! Initialise arrays to zero because diatpr is called before they are first calculated 
    255449         ! Note that this means diagnostics will not be exactly correct when model run is restarted. 
    256          htr_adv(:) = 0._wp  ;  str_adv(:) =  0._wp   
    257          htr_ldf(:) = 0._wp  ;  str_ldf(:) =  0._wp  
     450         htr_adv(:,:) = 0._wp  ;  str_adv(:,:) =  0._wp  
     451         htr_ldf(:,:) = 0._wp  ;  str_ldf(:,:) =  0._wp  
     452         htr_eiv(:,:) = 0._wp  ;  str_eiv(:,:) =  0._wp  
     453         htr_ove(:,:) = 0._wp  ;   str_ove(:,:) =  0._wp 
     454         htr_btr(:,:) = 0._wp  ;   str_btr(:,:) =  0._wp 
    258455         ! 
    259456      ENDIF  
     
    261458   END SUBROUTINE dia_ptr_init 
    262459 
     460   SUBROUTINE dia_ptr_hst( ktra, cptr, pva )  
     461      !!---------------------------------------------------------------------- 
     462      !!                    ***  ROUTINE dia_ptr_hst *** 
     463      !!---------------------------------------------------------------------- 
     464      !! Wrapper for heat and salt transport calculations to calculate them for each basin 
     465      !! Called from all advection and/or diffusion routines 
     466      !!---------------------------------------------------------------------- 
     467      INTEGER                         , INTENT(in )  :: ktra  ! tracer index 
     468      CHARACTER(len=3)                , INTENT(in)   :: cptr  ! transport type  'adv'/'ldf'/'eiv' 
     469      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in)   :: pva   ! 3D input array of advection/diffusion 
     470      INTEGER                                        :: jn    ! 
     471 
     472      IF( cptr == 'adv' ) THEN 
     473         IF( ktra == jp_tem )  htr_adv(:,1) = ptr_sj( pva(:,:,:) ) 
     474         IF( ktra == jp_sal )  str_adv(:,1) = ptr_sj( pva(:,:,:) ) 
     475      ENDIF 
     476      IF( cptr == 'ldf' ) THEN 
     477         IF( ktra == jp_tem )  htr_ldf(:,1) = ptr_sj( pva(:,:,:) ) 
     478         IF( ktra == jp_sal )  str_ldf(:,1) = ptr_sj( pva(:,:,:) ) 
     479      ENDIF 
     480      IF( cptr == 'eiv' ) THEN 
     481         IF( ktra == jp_tem )  htr_eiv(:,1) = ptr_sj( pva(:,:,:) ) 
     482         IF( ktra == jp_sal )  str_eiv(:,1) = ptr_sj( pva(:,:,:) ) 
     483      ENDIF 
     484      ! 
     485      IF( ln_subbas ) THEN 
     486         ! 
     487         IF( cptr == 'adv' ) THEN 
     488             IF( ktra == jp_tem ) THEN  
     489                DO jn = 2, nptr 
     490                   htr_adv(:,jn) = ptr_sj( pva(:,:,:), btmsk(:,:,jn) ) 
     491                END DO 
     492             ENDIF 
     493             IF( ktra == jp_sal ) THEN  
     494                DO jn = 2, nptr 
     495                   str_adv(:,jn) = ptr_sj( pva(:,:,:), btmsk(:,:,jn) ) 
     496                END DO 
     497             ENDIF 
     498         ENDIF 
     499         IF( cptr == 'ldf' ) THEN 
     500             IF( ktra == jp_tem ) THEN  
     501                DO jn = 2, nptr 
     502                    htr_ldf(:,jn) = ptr_sj( pva(:,:,:), btmsk(:,:,jn) ) 
     503                 END DO 
     504             ENDIF 
     505             IF( ktra == jp_sal ) THEN  
     506                DO jn = 2, nptr 
     507                   str_ldf(:,jn) = ptr_sj( pva(:,:,:), btmsk(:,:,jn) ) 
     508                END DO 
     509             ENDIF 
     510         ENDIF 
     511         IF( cptr == 'eiv' ) THEN 
     512             IF( ktra == jp_tem ) THEN  
     513                DO jn = 2, nptr 
     514                    htr_eiv(:,jn) = ptr_sj( pva(:,:,:), btmsk(:,:,jn) ) 
     515                 END DO 
     516             ENDIF 
     517             IF( ktra == jp_sal ) THEN  
     518                DO jn = 2, nptr 
     519                   str_eiv(:,jn) = ptr_sj( pva(:,:,:), btmsk(:,:,jn) ) 
     520                END DO 
     521             ENDIF 
     522         ENDIF 
     523         ! 
     524      ENDIF 
     525   END SUBROUTINE dia_ptr_hst 
     526 
    263527 
    264528   FUNCTION dia_ptr_alloc() 
     
    271535      ierr(:) = 0 
    272536      ! 
    273       ALLOCATE( btmsk(jpi,jpj,nptr) ,           & 
    274          &      htr_adv(jpj) , str_adv(jpj) ,   & 
    275          &      htr_ldf(jpj) , str_ldf(jpj) , STAT=ierr(1)  ) 
     537      ALLOCATE( btmsk(jpi,jpj,nptr) ,              & 
     538         &      htr_adv(jpj,nptr) , str_adv(jpj,nptr) ,   & 
     539         &      htr_eiv(jpj,nptr) , str_eiv(jpj,nptr) ,   & 
     540         &      htr_ove(jpj,nptr) , str_ove(jpj,nptr) ,   & 
     541         &      htr_btr(jpj,nptr) , str_btr(jpj,nptr) ,   & 
     542         &      htr_ldf(jpj,nptr) , str_ldf(jpj,nptr) , STAT=ierr(1)  ) 
    276543         ! 
    277544      ALLOCATE( p_fval1d(jpj), p_fval2d(jpj,jpk), Stat=ierr(2)) 
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r6387 r7403  
    302302      CALL iom_put( "hdiv", hdivn )                  ! Horizontal divergence 
    303303      ! 
    304       IF( iom_use("u_masstr") .OR. iom_use("u_heattr") .OR. iom_use("u_salttr") ) THEN 
     304      IF( iom_use("u_masstr") .OR. iom_use("u_masstr_vint") .OR. iom_use("u_heattr") .OR. iom_use("u_salttr") ) THEN 
    305305         z3d(:,:,jpk) = 0.e0 
     306         z2d(:,:) = 0.e0 
    306307         DO jk = 1, jpkm1 
    307308            z3d(:,:,jk) = rau0 * un(:,:,jk) * e2u(:,:) * e3u_n(:,:,jk) * umask(:,:,jk) 
     309            z2d(:,:) = z2d(:,:) + z3d(:,:,jk) 
    308310         END DO 
    309311         CALL iom_put( "u_masstr", z3d )                  ! mass transport in i-direction 
     312         CALL iom_put( "u_masstr_vint", z2d )             ! mass transport in i-direction vertical sum 
    310313      ENDIF 
    311314       
     
    370373         CALL iom_put( "v_salttr", 0.5 * z2d )            !  heat transport in j-direction 
    371374      ENDIF 
     375 
     376      ! Vertical integral of temperature 
     377      IF( iom_use("tosmint") ) THEN 
     378         z2d(:,:)=0._wp 
     379         DO jk = 1, jpkm1 
     380            DO jj = 2, jpjm1 
     381               DO ji = fs_2, fs_jpim1   ! vector opt. 
     382                  z2d(ji,jj) = z2d(ji,jj) + rau0 * e3t_n(ji,jj,jk) *  tsn(ji,jj,jk,jp_tem) 
     383               END DO 
     384            END DO 
     385         END DO 
     386         CALL lbc_lnk( z2d, 'T', -1. ) 
     387         CALL iom_put( "tosmint", z2d )  
     388      ENDIF 
     389 
     390      ! Vertical integral of salinity 
     391      IF( iom_use("somint") ) THEN 
     392         z2d(:,:)=0._wp 
     393         DO jk = 1, jpkm1 
     394            DO jj = 2, jpjm1 
     395               DO ji = fs_2, fs_jpim1   ! vector opt. 
     396                  z2d(ji,jj) = z2d(ji,jj) + rau0 * e3t_n(ji,jj,jk) * tsn(ji,jj,jk,jp_sal) 
     397               END DO 
     398            END DO 
     399         END DO 
     400         CALL lbc_lnk( z2d, 'T', -1. ) 
     401         CALL iom_put( "somint", z2d )  
     402      ENDIF 
     403 
     404      CALL iom_put( "bn2", rn2 )  !Brunt-Vaisala buoyancy frequency (N^2) 
    372405      ! 
    373406      CALL wrk_dealloc( jpi , jpj      , z2d ) 
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90

    r6140 r7403  
    2424   USE ldfslp          ! lateral diffusion: slope of iso-neutral surfaces 
    2525   USE ldfc1d_c2d      ! lateral diffusion: 1D & 2D cases  
    26    USE diaar5, ONLY:   lk_diaar5 
     26   USE diaptr 
    2727   ! 
    2828   USE trc_oce, ONLY: lk_offline ! offline flag 
     
    730730      CALL iom_put( "woce_eiv", zw3d ) 
    731731      ! 
     732      !       
     733      ! 
     734      CALL wrk_alloc( jpi,jpj,   zw2d ) 
     735      ! 
     736      zztmp = 0.5_wp * rau0 * rcp  
     737      IF( iom_use('ueiv_heattr') .OR. iom_use('ueiv_heattr3d') ) THEN 
     738        zw2d(:,:)   = 0._wp  
     739        zw3d(:,:,:) = 0._wp  
     740        DO jk = 1, jpkm1 
     741           DO jj = 2, jpjm1 
     742              DO ji = fs_2, fs_jpim1   ! vector opt. 
     743                 zw3d(ji,jj,jk) = zw3d(ji,jj,jk) + ( psi_uw(ji,jj,jk+1)      - psi_uw(ji,jj,jk)          )   & 
     744                    &                            * ( tsn   (ji,jj,jk,jp_tem) + tsn   (ji+1,jj,jk,jp_tem) )  
     745                 zw2d(ji,jj) = zw2d(ji,jj) + zw3d(ji,jj,jk) 
     746              END DO 
     747           END DO 
     748        END DO 
     749        CALL lbc_lnk( zw2d, 'U', -1. ) 
     750        CALL lbc_lnk( zw3d, 'U', -1. ) 
     751        CALL iom_put( "ueiv_heattr"  , zztmp * zw2d )                  ! heat transport in i-direction 
     752        CALL iom_put( "ueiv_heattr3d", zztmp * zw3d )                  ! heat transport in i-direction 
     753      ENDIF 
     754      zw2d(:,:)   = 0._wp  
     755      zw3d(:,:,:) = 0._wp  
     756      DO jk = 1, jpkm1 
     757         DO jj = 2, jpjm1 
     758            DO ji = fs_2, fs_jpim1   ! vector opt. 
     759               zw3d(ji,jj,jk) = zw3d(ji,jj,jk) + ( psi_vw(ji,jj,jk+1)      - psi_vw(ji,jj,jk)          )   & 
     760                  &                            * ( tsn   (ji,jj,jk,jp_tem) + tsn   (ji,jj+1,jk,jp_tem) )  
     761               zw2d(ji,jj) = zw2d(ji,jj) + zw3d(ji,jj,jk) 
     762            END DO 
     763         END DO 
     764      END DO 
     765      CALL lbc_lnk( zw2d, 'V', -1. ) 
     766      CALL iom_put( "veiv_heattr", zztmp * zw2d )                  !  heat transport in j-direction 
     767      CALL iom_put( "veiv_heattr", zztmp * zw3d )                  !  heat transport in j-direction 
     768      ! 
     769      IF( ln_diaptr )  CALL dia_ptr_hst( jp_tem, 'eiv', 0.5 * zw3d ) 
     770      ! 
     771      zztmp = 0.5_wp * 0.5 
     772      IF( iom_use('ueiv_salttr') .OR. iom_use('ueiv_salttr3d')) THEN 
     773        zw2d(:,:) = 0._wp  
     774        zw3d(:,:,:) = 0._wp  
     775        DO jk = 1, jpkm1 
     776           DO jj = 2, jpjm1 
     777              DO ji = fs_2, fs_jpim1   ! vector opt. 
     778                 zw3d(ji,jj,jk) = zw3d(ji,jj,jk) * ( psi_uw(ji,jj,jk+1)      - psi_uw(ji,jj,jk)          )   & 
     779                    &                            * ( tsn   (ji,jj,jk,jp_sal) + tsn   (ji+1,jj,jk,jp_sal) )  
     780                 zw2d(ji,jj) = zw2d(ji,jj) + zw3d(ji,jj,jk) 
     781              END DO 
     782           END DO 
     783        END DO 
     784        CALL lbc_lnk( zw2d, 'U', -1. ) 
     785        CALL lbc_lnk( zw3d, 'U', -1. ) 
     786        CALL iom_put( "ueiv_salttr", zztmp * zw2d )                  ! salt transport in i-direction 
     787        CALL iom_put( "ueiv_salttr3d", zztmp * zw3d )                  ! salt transport in i-direction 
     788      ENDIF 
     789      zw2d(:,:) = 0._wp  
     790      zw3d(:,:,:) = 0._wp  
     791      DO jk = 1, jpkm1 
     792         DO jj = 2, jpjm1 
     793            DO ji = fs_2, fs_jpim1   ! vector opt. 
     794               zw3d(ji,jj,jk) = zw3d(ji,jj,jk) + ( psi_vw(ji,jj,jk+1)      - psi_vw(ji,jj,jk)          )   & 
     795                  &                            * ( tsn   (ji,jj,jk,jp_sal) + tsn   (ji,jj+1,jk,jp_sal) )  
     796               zw2d(ji,jj) = zw2d(ji,jj) + zw3d(ji,jj,jk) 
     797            END DO 
     798         END DO 
     799      END DO 
     800      CALL lbc_lnk( zw2d, 'V', -1. ) 
     801      CALL iom_put( "veiv_salttr", zztmp * zw2d )                  !  salt transport in j-direction 
     802      CALL iom_put( "veiv_salttr", zztmp * zw3d )                  !  salt transport in j-direction 
     803      ! 
     804      IF( ln_diaptr ) CALL dia_ptr_hst( jp_sal, 'eiv', 0.5 * zw3d ) 
     805      ! 
     806      CALL wrk_dealloc( jpi,jpj,   zw2d ) 
    732807      CALL wrk_dealloc( jpi,jpj,jpk,   zw3d ) 
    733       !       
    734       ! 
    735       IF( lk_diaar5 ) THEN                               !==  eiv heat transport: calculate and output  ==! 
    736          CALL wrk_alloc( jpi,jpj,   zw2d ) 
    737          ! 
    738          zztmp = 0.5_wp * rau0 * rcp  
    739          zw2d(:,:) = 0._wp  
    740          DO jk = 1, jpkm1 
    741             DO jj = 2, jpjm1 
    742                DO ji = fs_2, fs_jpim1   ! vector opt. 
    743                   zw2d(ji,jj) = zw2d(ji,jj) + zztmp * ( psi_uw(ji,jj,jk+1)      - psi_uw(ji,jj,jk)          )   & 
    744                      &                              * ( tsn   (ji,jj,jk,jp_tem) + tsn   (ji+1,jj,jk,jp_tem) )  
    745                END DO 
    746             END DO 
    747          END DO 
    748          CALL lbc_lnk( zw2d, 'U', -1. ) 
    749          CALL iom_put( "ueiv_heattr", zw2d )                  ! heat transport in i-direction 
    750          zw2d(:,:) = 0._wp  
    751          DO jk = 1, jpkm1 
    752             DO jj = 2, jpjm1 
    753                DO ji = fs_2, fs_jpim1   ! vector opt. 
    754                   zw2d(ji,jj) = zw2d(ji,jj) + zztmp * ( psi_vw(ji,jj,jk+1)      - psi_vw(ji,jj,jk)          )   & 
    755                      &                              * ( tsn   (ji,jj,jk,jp_tem) + tsn   (ji,jj+1,jk,jp_tem) )  
    756                END DO 
    757             END DO 
    758          END DO 
    759          CALL lbc_lnk( zw2d, 'V', -1. ) 
    760          CALL iom_put( "veiv_heattr", zw2d )                  !  heat transport in i-direction 
    761          ! 
    762          CALL wrk_dealloc( jpi,jpj,   zw2d ) 
    763       ENDIF 
    764808      ! 
    765809      IF( nn_timing == 1 )  CALL timing_stop( 'ldf_eiv_dia')       
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_oasis3.F90

    r5836 r7403  
    6666   INTEGER                    ::   nsnd         ! total number of fields sent  
    6767   INTEGER                    ::   ncplmodel    ! Maximum number of models to/from which NEMO is potentialy sending/receiving data 
    68    INTEGER, PUBLIC, PARAMETER ::   nmaxfld=50   ! Maximum number of coupling fields 
     68   INTEGER, PUBLIC, PARAMETER ::   nmaxfld=55   ! Maximum number of coupling fields 
    6969   INTEGER, PUBLIC, PARAMETER ::   nmaxcat=5    ! Maximum number of coupling fields 
    7070   INTEGER, PUBLIC, PARAMETER ::   nmaxcpl=5    ! Maximum number of coupling fields 
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90

    r6140 r7403  
    6565   LOGICAL , PUBLIC ::   ln_cdgw        !: true if neutral drag coefficient from wave model 
    6666   LOGICAL , PUBLIC ::   ln_sdw         !: true if 3d stokes drift from wave model 
     67   LOGICAL , PUBLIC ::   ln_tauoc       !: true if normalized stress from wave is used 
     68   LOGICAL , PUBLIC ::   ln_stcor       !: true if Stokes-Coriolis term is used 
    6769   ! 
    6870   LOGICAL , PUBLIC ::   ln_icebergs    !: Icebergs 
     
    120122   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   sprecip           !: solid precipitation                          [Kg/m2/s] 
    121123   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   fr_i              !: ice fraction = 1 - lead fraction      (between 0 to 1) 
    122 #if defined key_cpl_carbon_cycle 
    123124   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   atm_co2           !: atmospheric pCO2                             [ppm] 
    124 #endif 
    125125   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: xcplmask          !: coupling mask for ln_mixcpl (warning: allocated in sbccpl) 
    126126 
     
    166166         ! 
    167167      ALLOCATE( tprecip(jpi,jpj) , sprecip(jpi,jpj) , fr_i(jpi,jpj) ,     & 
    168 #if defined key_cpl_carbon_cycle 
    169168         &      atm_co2(jpi,jpj) ,                                        & 
    170 #endif 
    171169         &      ssu_m  (jpi,jpj) , sst_m(jpi,jpj) , frq_m(jpi,jpj) ,      & 
    172170         &      ssv_m  (jpi,jpj) , sss_m(jpi,jpj) , ssh_m(jpi,jpj) , STAT=ierr(4) ) 
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

    r6813 r7403  
    745745 
    746746      !! Neutral coefficients at 10m: 
    747       IF( ln_cdgw ) THEN      ! wave drag case 
     747      IF( ln_wave .AND. ln_cdgw ) THEN      ! wave drag case 
    748748         cdn_wave(:,:) = cdn_wave(:,:) + rsmall * ( 1._wp - tmask(:,:,1) ) 
    749749         ztmp0   (:,:) = cdn_wave(:,:) 
     
    791791         END IF 
    792792        
    793          IF( ln_cdgw ) THEN      ! surface wave case 
     793         IF( ln_wave .AND. ln_cdgw ) THEN      ! surface wave case 
    794794            sqrt_Cd = vkarmn / ( vkarmn / sqrt_Cd_n10 - zpsi_m_u )  
    795795            Cd      = sqrt_Cd * sqrt_Cd 
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90

    r6140 r7403  
    1717   USE fldread        ! read input fields 
    1818   USE sbc_oce        ! Surface boundary condition: ocean fields 
    19    USE sbcwave  ,ONLY :   cdn_wave !wave module 
    2019   ! 
    2120   USE iom            ! I/O manager library 
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

    r6722 r7403  
    1818   !!   sbc_cpl_snd     : send     fields to the atmosphere 
    1919   !!---------------------------------------------------------------------- 
    20    USE dom_oce        ! ocean space and time domain 
    21    USE sbc_oce        ! Surface boundary condition: ocean fields 
    22    USE sbc_ice        ! Surface boundary condition: ice fields 
    23    USE sbcapr         ! Stochastic param. : ??? 
    24    USE sbcdcy         ! surface boundary condition: diurnal cycle 
    25    USE phycst         ! physical constants 
     20   USE dom_oce         ! ocean space and time domain 
     21   USE sbc_oce         ! Surface boundary condition: ocean fields 
     22   USE trc_oce         ! share SMS/Ocean variables 
     23   USE sbc_ice         ! Surface boundary condition: ice fields 
     24   USE sbcapr          ! Stochastic param. : ??? 
     25   USE sbcdcy          ! surface boundary condition: diurnal cycle 
     26   USE sbcwave         ! surface boundary condition: waves 
     27   USE phycst          ! physical constants 
    2628#if defined key_lim3 
    2729   USE ice            ! ice variables 
     
    3638   USE albedo         !  
    3739   USE eosbn2         !  
    38    USE sbcrnf  , ONLY : l_rnfcpl 
    39 #if defined key_cpl_carbon_cycle 
    40    USE p4zflx, ONLY : oce_co2 
    41 #endif 
     40   USE sbcrnf, ONLY : l_rnfcpl 
    4241#if defined key_cice 
    4342   USE ice_domain_size, only: ncat 
     
    106105   INTEGER, PARAMETER ::   jpr_e3t1st = 41   ! first T level thickness  
    107106   INTEGER, PARAMETER ::   jpr_fraqsr = 42   ! fraction of solar net radiation absorbed in the first ocean level 
    108    INTEGER, PARAMETER ::   jprcv      = 42   ! total number of fields received 
     107   INTEGER, PARAMETER ::   jpr_mslp   = 43   ! mean sea level pressure  
     108   INTEGER, PARAMETER ::   jpr_hsig   = 44   ! Hsig  
     109   INTEGER, PARAMETER ::   jpr_phioc  = 45   ! Wave=>ocean energy flux  
     110   INTEGER, PARAMETER ::   jpr_sdrftx = 46   ! Stokes drift on grid 1  
     111   INTEGER, PARAMETER ::   jpr_sdrfty = 47   ! Stokes drift on grid 2  
     112   INTEGER, PARAMETER ::   jpr_wper   = 48   ! Mean wave period 
     113   INTEGER, PARAMETER ::   jpr_wnum   = 49   ! Mean wavenumber 
     114   INTEGER, PARAMETER ::   jpr_wstrf  = 50   ! Stress fraction adsorbed by waves 
     115   INTEGER, PARAMETER ::   jpr_wdrag  = 51   ! Neutral surface drag coefficient 
     116   INTEGER, PARAMETER ::   jprcv      = 51   ! total number of fields received   
    109117 
    110118   INTEGER, PARAMETER ::   jps_fice   =  1   ! ice fraction sent to the atmosphere 
     
    136144   INTEGER, PARAMETER ::   jps_e3t1st = 27   ! first level depth (vvl) 
    137145   INTEGER, PARAMETER ::   jps_fraqsr = 28   ! fraction of solar net radiation absorbed in the first ocean level 
    138    INTEGER, PARAMETER ::   jpsnd      = 28   ! total number of fields sended 
     146   INTEGER, PARAMETER ::   jps_ficet  = 29   ! total ice fraction   
     147   INTEGER, PARAMETER ::   jps_ocxw   = 30   ! currents on grid 1   
     148   INTEGER, PARAMETER ::   jps_ocyw   = 31   ! currents on grid 2 
     149   INTEGER, PARAMETER ::   jps_wlev   = 32   ! water level  
     150   INTEGER, PARAMETER ::   jpsnd      = 32   ! total number of fields sent  
    139151 
    140152   !                                  !!** namelist namsbc_cpl ** 
     
    150162   !                                   ! Received from the atmosphere 
    151163   TYPE(FLD_C) ::   sn_rcv_w10m, sn_rcv_taumod, sn_rcv_tau, sn_rcv_dqnsdt, sn_rcv_qsr, sn_rcv_qns, sn_rcv_emp, sn_rcv_rnf 
    152    TYPE(FLD_C) ::   sn_rcv_cal, sn_rcv_iceflx, sn_rcv_co2                         
     164   TYPE(FLD_C) ::   sn_rcv_cal, sn_rcv_iceflx, sn_rcv_co2, sn_rcv_mslp                            
     165   ! Send to waves  
     166   TYPE(FLD_C) ::   sn_snd_ifrac, sn_snd_crtw, sn_snd_wlev  
     167   ! Received from waves  
     168   TYPE(FLD_C) ::   sn_rcv_hsig,sn_rcv_phioc,sn_rcv_sdrfx,sn_rcv_sdrfy,sn_rcv_wper,sn_rcv_wnum,sn_rcv_wstrf,sn_rcv_wdrag 
    153169   !                                   ! Other namelist parameters 
    154170   INTEGER     ::   nn_cplmodel           ! Maximum number of models to/from which NEMO is potentialy sending/receiving data 
     
    163179   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   albedo_oce_mix    ! ocean albedo sent to atmosphere (mix clear/overcast sky) 
    164180 
    165    INTEGER , ALLOCATABLE, SAVE, DIMENSION(    :) ::   nrcvinfo          ! OASIS info argument 
     181   REAL(wp) ::   rpref = 101000._wp   ! reference atmospheric pressure[N/m2]  
     182   REAL(wp) ::   r1_grau              ! = 1.e0 / (grav * rau0)  
     183 
     184   INTEGER , ALLOCATABLE, SAVE, DIMENSION(    :) ::   nrcvinfo           ! OASIS info argument 
    166185 
    167186   !! Substitution 
     
    178197      !!             ***  FUNCTION sbc_cpl_alloc  *** 
    179198      !!---------------------------------------------------------------------- 
    180       INTEGER :: ierr(3) 
     199      INTEGER :: ierr(4) 
    181200      !!---------------------------------------------------------------------- 
    182201      ierr(:) = 0 
     
    189208      ALLOCATE( xcplmask(jpi,jpj,0:nn_cplmodel) , STAT=ierr(3) ) 
    190209      ! 
     210      IF( .NOT. ln_apr_dyn ) ALLOCATE( ssh_ib(jpi,jpj), ssh_ibb(jpi,jpj), apr(jpi, jpj), STAT=ierr(4) )  
     211 
    191212      sbc_cpl_alloc = MAXVAL( ierr ) 
    192213      IF( lk_mpp            )   CALL mpp_sum ( sbc_cpl_alloc ) 
     
    214235      REAL(wp), POINTER, DIMENSION(:,:) ::   zacs, zaos 
    215236      !! 
    216       NAMELIST/namsbc_cpl/  sn_snd_temp, sn_snd_alb   , sn_snd_thick, sn_snd_crt   , sn_snd_co2,      & 
    217          &                  sn_rcv_w10m, sn_rcv_taumod, sn_rcv_tau  , sn_rcv_dqnsdt, sn_rcv_qsr,      & 
    218          &                  sn_rcv_qns , sn_rcv_emp   , sn_rcv_rnf  , sn_rcv_cal   , sn_rcv_iceflx,   & 
    219          &                  sn_rcv_co2 , nn_cplmodel  , ln_usecplmask 
     237      NAMELIST/namsbc_cpl/  sn_snd_temp , sn_snd_alb  , sn_snd_thick , sn_snd_crt   , sn_snd_co2,      &  
     238         &                  sn_rcv_w10m, sn_rcv_taumod, sn_rcv_tau   , sn_rcv_dqnsdt, sn_rcv_qsr,      &  
     239         &                  sn_snd_ifrac, sn_snd_crtw , sn_snd_wlev  , sn_rcv_hsig  , sn_rcv_phioc ,   &  
     240         &                  sn_rcv_sdrfx, sn_rcv_sdrfy, sn_rcv_wper  , sn_rcv_wnum  , sn_rcv_wstrf ,   & 
     241         &                  sn_rcv_wdrag, sn_rcv_qns  , sn_rcv_emp   , sn_rcv_rnf   , sn_rcv_cal   ,   & 
     242         &                  sn_rcv_iceflx,sn_rcv_co2  , nn_cplmodel  , ln_usecplmask, sn_rcv_mslp  
    220243      !!--------------------------------------------------------------------- 
    221244      ! 
     
    258281         WRITE(numout,*)'      sea ice heat fluxes             = ', TRIM(sn_rcv_iceflx%cldes), ' (', TRIM(sn_rcv_iceflx%clcat), ')' 
    259282         WRITE(numout,*)'      atm co2                         = ', TRIM(sn_rcv_co2%cldes   ), ' (', TRIM(sn_rcv_co2%clcat   ), ')' 
     283         WRITE(numout,*)'      significant wave heigth         = ', TRIM(sn_rcv_hsig%cldes  ), ' (', TRIM(sn_rcv_hsig%clcat  ), ')'  
     284         WRITE(numout,*)'      wave to oce energy flux         = ', TRIM(sn_rcv_phioc%cldes ), ' (', TRIM(sn_rcv_phioc%clcat ), ')'  
     285         WRITE(numout,*)'      Surface Stokes drift grid u     = ', TRIM(sn_rcv_sdrfx%cldes ), ' (', TRIM(sn_rcv_sdrfx%clcat ), ')'  
     286         WRITE(numout,*)'      Surface Stokes drift grid v     = ', TRIM(sn_rcv_sdrfy%cldes ), ' (', TRIM(sn_rcv_sdrfy%clcat ), ')'  
     287         WRITE(numout,*)'      Mean wave period                = ', TRIM(sn_rcv_wper%cldes  ), ' (', TRIM(sn_rcv_wper%clcat  ), ')'  
     288         WRITE(numout,*)'      Mean wave number                = ', TRIM(sn_rcv_wnum%cldes  ), ' (', TRIM(sn_rcv_wnum%clcat  ), ')'  
     289         WRITE(numout,*)'      Stress frac adsorbed by waves   = ', TRIM(sn_rcv_wstrf%cldes ), ' (', TRIM(sn_rcv_wstrf%clcat ), ')'  
     290         WRITE(numout,*)'      Neutral surf drag coefficient   = ', TRIM(sn_rcv_wdrag%cldes ), ' (', TRIM(sn_rcv_wdrag%clcat ), ')'  
    260291         WRITE(numout,*)'  sent fields (multiple ice categories)' 
    261292         WRITE(numout,*)'      surface temperature             = ', TRIM(sn_snd_temp%cldes  ), ' (', TRIM(sn_snd_temp%clcat  ), ')' 
    262293         WRITE(numout,*)'      albedo                          = ', TRIM(sn_snd_alb%cldes   ), ' (', TRIM(sn_snd_alb%clcat   ), ')' 
    263294         WRITE(numout,*)'      ice/snow thickness              = ', TRIM(sn_snd_thick%cldes ), ' (', TRIM(sn_snd_thick%clcat ), ')' 
     295         WRITE(numout,*)'      total ice fraction              = ', TRIM(sn_snd_ifrac%cldes ), ' (', TRIM(sn_snd_ifrac%clcat ), ')'  
    264296         WRITE(numout,*)'      surface current                 = ', TRIM(sn_snd_crt%cldes   ), ' (', TRIM(sn_snd_crt%clcat   ), ')' 
    265297         WRITE(numout,*)'                      - referential   = ', sn_snd_crt%clvref  
     
    267299         WRITE(numout,*)'                      - mesh          = ', sn_snd_crt%clvgrd 
    268300         WRITE(numout,*)'      oce co2 flux                    = ', TRIM(sn_snd_co2%cldes   ), ' (', TRIM(sn_snd_co2%clcat   ), ')' 
     301         WRITE(numout,*)'      water level                     = ', TRIM(sn_snd_wlev%cldes  ), ' (', TRIM(sn_snd_wlev%clcat  ), ')'  
     302         WRITE(numout,*)'      mean sea level pressure         = ', TRIM(sn_rcv_mslp%cldes  ), ' (', TRIM(sn_rcv_mslp%clcat  ), ')'  
     303         WRITE(numout,*)'      surface current to waves        = ', TRIM(sn_snd_crtw%cldes  ), ' (', TRIM(sn_snd_crtw%clcat  ), ')'  
     304         WRITE(numout,*)'                      - referential   = ', sn_snd_crtw%clvref  
     305         WRITE(numout,*)'                      - orientation   = ', sn_snd_crtw%clvor  
     306         WRITE(numout,*)'                      - mesh          = ', sn_snd_crtw%clvgrd  
    269307         WRITE(numout,*)'  nn_cplmodel                         = ', nn_cplmodel 
    270308         WRITE(numout,*)'  ln_usecplmask                       = ', ln_usecplmask 
     
    305343      !  
    306344      ! Vectors: change of sign at north fold ONLY if on the local grid 
     345      IF( TRIM( sn_rcv_tau%cldes ) == 'oce only' .OR. TRIM(sn_rcv_tau%cldes ) == 'oce and ice') THEN ! avoid working with the atmospheric fields if they are not coupled 
    307346      IF( TRIM( sn_rcv_tau%clvor ) == 'local grid' )   srcv(jpr_otx1:jpr_itz2)%nsgn = -1. 
    308347       
     
    372411         srcv(jpr_ity1)%clgrid = 'V'                  ! i.e. it is always at U- & V-points for i- & j-comp. resp. 
    373412      ENDIF 
    374       ! 
     413      ENDIF 
     414 
    375415      !                                                      ! ------------------------- ! 
    376416      !                                                      !    freshwater budget      !   E-P 
     
    467507      !                                                      !      Atmospheric CO2      ! 
    468508      !                                                      ! ------------------------- ! 
    469       srcv(jpr_co2 )%clname = 'O_AtmCO2'   ;   IF( TRIM(sn_rcv_co2%cldes   ) == 'coupled' )    srcv(jpr_co2 )%laction = .TRUE. 
     509      srcv(jpr_co2 )%clname = 'O_AtmCO2'    
     510      IF( TRIM(sn_rcv_co2%cldes   ) == 'coupled' )  THEN 
     511         srcv(jpr_co2 )%laction = .TRUE. 
     512         l_co2cpl = .TRUE. 
     513         IF(lwp) WRITE(numout,*) 
     514         IF(lwp) WRITE(numout,*) '   Atmospheric pco2 received from oasis ' 
     515         IF(lwp) WRITE(numout,*) 
     516      ENDIF 
     517 
     518      !                                                      ! ------------------------- !  
     519      !                                                      ! Mean Sea Level Pressure   !  
     520      !                                                      ! ------------------------- !  
     521      srcv(jpr_mslp)%clname = 'O_MSLP'     ;   IF( TRIM(sn_rcv_mslp%cldes  ) == 'coupled' )    srcv(jpr_mslp)%laction = .TRUE.  
     522 
    470523      !                                                      ! ------------------------- ! 
    471524      !                                                      !   topmelt and botmelt     !    
     
    481534         srcv(jpr_topm:jpr_botm)%laction = .TRUE. 
    482535      ENDIF 
     536      !                                                      ! ------------------------- ! 
     537      !                                                      !      Wave breaking        !     
     538      !                                                      ! ------------------------- !  
     539      srcv(jpr_hsig)%clname  = 'O_Hsigwa'    ! significant wave height 
     540      IF( TRIM(sn_rcv_hsig%cldes  ) == 'coupled' )  THEN 
     541         srcv(jpr_hsig)%laction = .TRUE. 
     542         cpl_hsig = .TRUE. 
     543      ENDIF 
     544      srcv(jpr_phioc)%clname = 'O_PhiOce'    ! wave to ocean energy 
     545      IF( TRIM(sn_rcv_phioc%cldes ) == 'coupled' )  THEN 
     546         srcv(jpr_phioc)%laction = .TRUE. 
     547         cpl_phioc = .TRUE. 
     548      ENDIF 
     549      srcv(jpr_sdrftx)%clname = 'O_Sdrfx'    ! Stokes drift in the u direction 
     550      IF( TRIM(sn_rcv_sdrfx%cldes ) == 'coupled' )  THEN 
     551         srcv(jpr_sdrftx)%laction = .TRUE. 
     552         cpl_sdrftx = .TRUE. 
     553      ENDIF 
     554      srcv(jpr_sdrfty)%clname = 'O_Sdrfy'    ! Stokes drift in the v direction 
     555      IF( TRIM(sn_rcv_sdrfy%cldes ) == 'coupled' )  THEN 
     556         srcv(jpr_sdrfty)%laction = .TRUE. 
     557         cpl_sdrfty = .TRUE. 
     558      ENDIF 
     559      srcv(jpr_wper)%clname = 'O_WPer'       ! mean wave period 
     560      IF( TRIM(sn_rcv_wper%cldes  ) == 'coupled' )  THEN 
     561         srcv(jpr_wper)%laction = .TRUE. 
     562         cpl_wper = .TRUE. 
     563      ENDIF 
     564      srcv(jpr_wnum)%clname = 'O_WNum'       ! mean wave number 
     565      IF( TRIM(sn_rcv_wnum%cldes ) == 'coupled' )  THEN 
     566         srcv(jpr_wnum)%laction = .TRUE. 
     567         cpl_wnum = .TRUE. 
     568      ENDIF 
     569      srcv(jpr_wstrf)%clname = 'O_WStrf'     ! stress fraction adsorbed by the wave 
     570      IF( TRIM(sn_rcv_wstrf%cldes ) == 'coupled' )  THEN 
     571         srcv(jpr_wstrf)%laction = .TRUE. 
     572         cpl_wstrf = .TRUE. 
     573      ENDIF 
     574      srcv(jpr_wdrag)%clname = 'O_WDrag'     ! neutral surface drag coefficient 
     575      IF( TRIM(sn_rcv_wdrag%cldes ) == 'coupled' )  THEN 
     576         srcv(jpr_wdrag)%laction = .TRUE. 
     577         cpl_wdrag = .TRUE. 
     578      ENDIF 
     579      !  
    483580      !                                                      ! ------------------------------- ! 
    484581      !                                                      !   OPA-SAS coupling - rcv by opa !    
     
    635732      !                                                      ! ------------------------- ! 
    636733      ssnd(jps_fice)%clname = 'OIceFrc' 
     734      ssnd(jps_ficet)%clname = 'OIceFrcT'  
    637735      ssnd(jps_hice)%clname = 'OIceTck' 
    638736      ssnd(jps_hsnw)%clname = 'OSnwTck' 
     
    643741      ENDIF 
    644742       
     743      IF (TRIM( sn_snd_ifrac%cldes )  == 'coupled') ssnd(jps_ficet)%laction = .TRUE.  
     744 
    645745      SELECT CASE ( TRIM( sn_snd_thick%cldes ) ) 
    646746      CASE( 'none'         )       ! nothing to do 
     
    663763      ssnd(jps_ocy1)%clname = 'O_OCury1'   ;   ssnd(jps_ivy1)%clname = 'O_IVely1' 
    664764      ssnd(jps_ocz1)%clname = 'O_OCurz1'   ;   ssnd(jps_ivz1)%clname = 'O_IVelz1' 
     765      ssnd(jps_ocxw)%clname = 'O_OCurxw'  
     766      ssnd(jps_ocyw)%clname = 'O_OCuryw'  
    665767      ! 
    666768      ssnd(jps_ocx1:jps_ivz1)%nsgn = -1.   ! vectors: change of the sign at the north fold 
     
    683785      END SELECT 
    684786 
     787      ssnd(jps_ocxw:jps_ocyw)%nsgn = -1.   ! vectors: change of the sign at the north fold  
     788         
     789      IF( sn_snd_crtw%clvgrd == 'U,V' ) THEN  
     790         ssnd(jps_ocxw)%clgrid = 'U' ; ssnd(jps_ocyw)%clgrid = 'V'  
     791      ELSE IF( sn_snd_crtw%clvgrd /= 'T' ) THEN  
     792         CALL ctl_stop( 'sn_snd_crtw%clvgrd must be equal to T' )  
     793      ENDIF  
     794      IF( TRIM( sn_snd_crtw%clvor ) == 'eastward-northward' ) ssnd(jps_ocxw:jps_ocyw)%nsgn = 1.  
     795      SELECT CASE( TRIM( sn_snd_crtw%cldes ) )  
     796         CASE( 'none'                 )   ; ssnd(jps_ocxw:jps_ocyw)%laction = .FALSE.  
     797         CASE( 'oce only'             )   ; ssnd(jps_ocxw:jps_ocyw)%laction = .TRUE.  
     798         CASE( 'weighted oce and ice' )   !   nothing to do  
     799         CASE( 'mixed oce-ice'        )   ; ssnd(jps_ivx1:jps_ivz1)%laction = .FALSE.  
     800         CASE default   ;   CALL ctl_stop( 'sbc_cpl_init: wrong definition of sn_snd_crtw%cldes' )  
     801      END SELECT  
     802 
    685803      !                                                      ! ------------------------- ! 
    686804      !                                                      !          CO2 flux         ! 
    687805      !                                                      ! ------------------------- ! 
    688806      ssnd(jps_co2)%clname = 'O_CO2FLX' ;  IF( TRIM(sn_snd_co2%cldes) == 'coupled' )    ssnd(jps_co2 )%laction = .TRUE. 
     807 
     808      !                                                      ! ------------------------- !  
     809      !                                                      !     Sea surface height    !  
     810      !                                                      ! ------------------------- !  
     811      ssnd(jps_wlev)%clname = 'O_Wlevel' ;  IF( TRIM(sn_snd_wlev%cldes) == 'coupled' )   ssnd(jps_wlev)%laction = .TRUE.  
    689812 
    690813      !                                                      ! ------------------------------- ! 
     
    781904      IF( ln_dm2dc .AND. ln_cpl .AND. ncpl_qsr_freq /= 86400 )   & 
    782905         &   CALL ctl_stop( 'sbc_cpl_init: diurnal cycle reconstruction (ln_dm2dc) needs daily couping for solar radiation' ) 
    783       ncpl_qsr_freq = 86400 / ncpl_qsr_freq 
     906      IF( ln_dm2dc .AND. ln_cpl ) ncpl_qsr_freq = 86400 / ncpl_qsr_freq 
    784907 
    785908      CALL wrk_dealloc( jpi,jpj,   zacs, zaos ) 
     
    835958      !!                        emp          upward mass flux [evap. - precip. (- runoffs) (- calving)] (ocean only case) 
    836959      !!---------------------------------------------------------------------- 
    837       INTEGER, INTENT(in) ::   kt       ! ocean model time step index 
    838       INTEGER, INTENT(in) ::   k_fsbc   ! frequency of sbc (-> ice model) computation  
    839       INTEGER, INTENT(in) ::   k_ice    ! ice management in the sbc (=0/1/2/3) 
    840  
     960      USE zdf_oce,  ONLY : ln_zdfqiao 
     961 
     962      IMPLICIT NONE 
     963 
     964      INTEGER, INTENT(in)           ::   kt          ! ocean model time step index 
     965      INTEGER, INTENT(in)           ::   k_fsbc      ! frequency of sbc (-> ice model) computation  
     966      INTEGER, INTENT(in)           ::   k_ice       ! ice management in the sbc (=0/1/2/3) 
    841967      !! 
    842968      LOGICAL  ::   llnewtx, llnewtau      ! update wind stress components and module?? 
     
    9841110      ENDIF 
    9851111 
    986 #if defined key_cpl_carbon_cycle 
    9871112      !                                                      ! ================== ! 
    9881113      !                                                      ! atmosph. CO2 (ppm) ! 
    9891114      !                                                      ! ================== ! 
    9901115      IF( srcv(jpr_co2)%laction )   atm_co2(:,:) = frcv(jpr_co2)%z3(:,:,1) 
    991 #endif 
     1116      !  
     1117      !                                                      ! ========================= !  
     1118      !                                                      ! Mean Sea Level Pressure   !   (taum)  
     1119      !                                                      ! ========================= !  
     1120      !  
     1121      IF( srcv(jpr_mslp)%laction ) THEN                    ! UKMO SHELF effect of atmospheric pressure on SSH  
     1122          IF( kt /= nit000 )   ssh_ibb(:,:) = ssh_ib(:,:)    !* Swap of ssh_ib fields  
     1123 
     1124          r1_grau = 1.e0 / (grav * rau0)               !* constant for optimization  
     1125          ssh_ib(:,:) = - ( frcv(jpr_mslp)%z3(:,:,1) - rpref ) * r1_grau    ! equivalent ssh (inverse barometer)  
     1126          apr   (:,:) =     frcv(jpr_mslp)%z3(:,:,1)                         !atmospheric pressure  
     1127     
     1128          IF( kt == nit000 ) ssh_ibb(:,:) = ssh_ib(:,:)  ! correct this later (read from restart if possible)  
     1129      END IF  
     1130      ! 
     1131      IF( ln_sdw ) THEN  ! Stokes Drift correction activated 
     1132      !                                                      ! ========================= !  
     1133      !                                                      !       Stokes drift u      ! 
     1134      !                                                      ! ========================= !  
     1135         IF( srcv(jpr_sdrftx)%laction ) zusd2dt(:,:) = frcv(jpr_sdrftx)%z3(:,:,1) 
     1136      ! 
     1137      !                                                      ! ========================= !  
     1138      !                                                      !       Stokes drift v      ! 
     1139      !                                                      ! ========================= !  
     1140         IF( srcv(jpr_sdrfty)%laction ) zvsd2dt(:,:) = frcv(jpr_sdrfty)%z3(:,:,1) 
     1141      ! 
     1142      !                                                      ! ========================= !  
     1143      !                                                      !      Wave mean period     ! 
     1144      !                                                      ! ========================= !  
     1145         IF( srcv(jpr_wper)%laction ) wmp(:,:) = frcv(jpr_wper)%z3(:,:,1) 
     1146      ! 
     1147      !                                                      ! ========================= !  
     1148      !                                                      !  Significant wave height  ! 
     1149      !                                                      ! ========================= !  
     1150         IF( srcv(jpr_hsig)%laction ) swh(:,:) = frcv(jpr_hsig)%z3(:,:,1) 
     1151      ! 
     1152      !                                                      ! ========================= !  
     1153      !                                                      !    Vertical mixing Qiao   ! 
     1154      !                                                      ! ========================= !  
     1155         IF( srcv(jpr_wnum)%laction .AND. ln_zdfqiao ) wnum(:,:) = frcv(jpr_wnum)%z3(:,:,1) 
     1156 
     1157         ! Calculate the 3D Stokes drift both in coupled and not fully uncoupled mode 
     1158         IF( srcv(jpr_sdrftx)%laction .OR. srcv(jpr_sdrfty)%laction .OR. srcv(jpr_wper)%laction & 
     1159                                                                    .OR. srcv(jpr_hsig)%laction ) THEN 
     1160            CALL sbc_stokes() 
     1161            IF( ln_zdfqiao .AND. .NOT. srcv(jpr_wnum)%laction ) CALL sbc_qiao() 
     1162         ENDIF 
     1163         IF( ln_zdfqiao .AND. srcv(jpr_wnum)%laction ) CALL sbc_qiao() 
     1164      ENDIF 
     1165      !                                                      ! ========================= !  
     1166      !                                                      ! Stress adsorbed by waves  ! 
     1167      !                                                      ! ========================= !  
     1168      IF( srcv(jpr_wstrf)%laction .AND. ln_tauoc ) tauoc_wave(:,:) = frcv(jpr_wstrf)%z3(:,:,1) 
     1169 
     1170      !                                                      ! ========================= !  
     1171      !                                                      !   Wave drag coefficient   ! 
     1172      !                                                      ! ========================= !  
     1173      IF( srcv(jpr_wdrag)%laction .AND. ln_cdgw ) cdn_wave(:,:) = frcv(jpr_wdrag)%z3(:,:,1) 
    9921174 
    9931175      !  Fields received by SAS when OASIS coupling 
     
    19192101         IF( ssnd(jps_hsnw)%laction )   CALL cpl_snd( jps_hsnw, isec, ztmp4, info ) 
    19202102      ENDIF 
    1921       ! 
    1922 #if defined key_cpl_carbon_cycle 
    19232103      !                                                      ! ------------------------- ! 
    19242104      !                                                      !  CO2 flux from PISCES     !  
    19252105      !                                                      ! ------------------------- ! 
    1926       IF( ssnd(jps_co2)%laction )   CALL cpl_snd( jps_co2, isec, RESHAPE ( oce_co2, (/jpi,jpj,1/) ) , info ) 
    1927       ! 
    1928 #endif 
     2106      IF( ssnd(jps_co2)%laction .AND. l_co2cpl )   CALL cpl_snd( jps_co2, isec, RESHAPE ( oce_co2, (/jpi,jpj,1/) ) , info ) 
     2107      ! 
    19292108      !                                                      ! ------------------------- ! 
    19302109      IF( ssnd(jps_ocx1)%laction ) THEN                      !      Surface current      ! 
     
    20632242      ENDIF 
    20642243      ! 
     2244      !                                                      ! ------------------------- !  
     2245      !                                                      !  Surface current to waves !  
     2246      !                                                      ! ------------------------- !  
     2247      IF( ssnd(jps_ocxw)%laction .OR. ssnd(jps_ocyw)%laction ) THEN  
     2248          !      
     2249          !                                                  j+1  j     -----V---F  
     2250          ! surface velocity always sent from T point                    !       |  
     2251          !                                                       j      |   T   U  
     2252          !                                                              |       |  
     2253          !                                                   j   j-1   -I-------|  
     2254          !                                               (for I)        |       |  
     2255          !                                                             i-1  i   i  
     2256          !                                                              i      i+1 (for I)  
     2257          SELECT CASE( TRIM( sn_snd_crtw%cldes ) )  
     2258          CASE( 'oce only'             )      ! C-grid ==> T  
     2259             DO jj = 2, jpjm1  
     2260                DO ji = fs_2, fs_jpim1   ! vector opt.  
     2261                   zotx1(ji,jj) = 0.5 * ( un(ji,jj,1) + un(ji-1,jj  ,1) )  
     2262                   zoty1(ji,jj) = 0.5 * ( vn(ji,jj,1) + vn(ji , jj-1,1) )   
     2263                END DO  
     2264             END DO  
     2265          CASE( 'weighted oce and ice' )     
     2266             SELECT CASE ( cp_ice_msh )  
     2267             CASE( 'C' )                      ! Ocean and Ice on C-grid ==> T  
     2268                DO jj = 2, jpjm1  
     2269                   DO ji = fs_2, fs_jpim1   ! vector opt.  
     2270                      zotx1(ji,jj) = 0.5 * ( un   (ji,jj,1) + un   (ji-1,jj  ,1) ) * zfr_l(ji,jj)    
     2271                      zoty1(ji,jj) = 0.5 * ( vn   (ji,jj,1) + vn   (ji  ,jj-1,1) ) * zfr_l(ji,jj)  
     2272                      zitx1(ji,jj) = 0.5 * ( u_ice(ji,jj  ) + u_ice(ji-1,jj    ) ) *  fr_i(ji,jj)  
     2273                      zity1(ji,jj) = 0.5 * ( v_ice(ji,jj  ) + v_ice(ji  ,jj-1  ) ) *  fr_i(ji,jj)  
     2274                   END DO  
     2275                END DO  
     2276             CASE( 'I' )                      ! Ocean on C grid, Ice on I-point (B-grid) ==> T  
     2277                DO jj = 2, jpjm1  
     2278                   DO ji = 2, jpim1   ! NO vector opt.  
     2279                      zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj  ,1) ) * zfr_l(ji,jj)    
     2280                      zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)      + vn(ji  ,jj-1,1) ) * zfr_l(ji,jj)    
     2281                      zitx1(ji,jj) = 0.25 * ( u_ice(ji+1,jj+1) + u_ice(ji,jj+1)                     &  
     2282                         &                  + u_ice(ji+1,jj  ) + u_ice(ji,jj  )  ) *  fr_i(ji,jj)  
     2283                      zity1(ji,jj) = 0.25 * ( v_ice(ji+1,jj+1) + v_ice(ji,jj+1)                     &  
     2284                         &                  + v_ice(ji+1,jj  ) + v_ice(ji,jj  )  ) *  fr_i(ji,jj)  
     2285                   END DO  
     2286                END DO  
     2287             CASE( 'F' )                      ! Ocean on C grid, Ice on F-point (B-grid) ==> T  
     2288                DO jj = 2, jpjm1  
     2289                   DO ji = 2, jpim1   ! NO vector opt.  
     2290                      zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj  ,1) ) * zfr_l(ji,jj)    
     2291                      zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)      + vn(ji  ,jj-1,1) ) * zfr_l(ji,jj)    
     2292                      zitx1(ji,jj) = 0.25 * ( u_ice(ji-1,jj-1) + u_ice(ji,jj-1)                     &  
     2293                         &                  + u_ice(ji-1,jj  ) + u_ice(ji,jj  )  ) *  fr_i(ji,jj)  
     2294                      zity1(ji,jj) = 0.25 * ( v_ice(ji-1,jj-1) + v_ice(ji,jj-1)                     &  
     2295                         &                  + v_ice(ji-1,jj  ) + v_ice(ji,jj  )  ) *  fr_i(ji,jj)  
     2296                   END DO  
     2297                END DO  
     2298             END SELECT  
     2299             CALL lbc_lnk( zitx1, 'T', -1. )   ;   CALL lbc_lnk( zity1, 'T', -1. )  
     2300          CASE( 'mixed oce-ice'        )  
     2301             SELECT CASE ( cp_ice_msh )  
     2302             CASE( 'C' )                      ! Ocean and Ice on C-grid ==> T  
     2303                DO jj = 2, jpjm1  
     2304                   DO ji = fs_2, fs_jpim1   ! vector opt.  
     2305                      zotx1(ji,jj) = 0.5 * ( un   (ji,jj,1) + un   (ji-1,jj  ,1) ) * zfr_l(ji,jj)   &  
     2306                         &         + 0.5 * ( u_ice(ji,jj  ) + u_ice(ji-1,jj    ) ) *  fr_i(ji,jj)  
     2307                      zoty1(ji,jj) = 0.5 * ( vn   (ji,jj,1) + vn   (ji  ,jj-1,1) ) * zfr_l(ji,jj)   &  
     2308                         &         + 0.5 * ( v_ice(ji,jj  ) + v_ice(ji  ,jj-1  ) ) *  fr_i(ji,jj)  
     2309                   END DO  
     2310                END DO  
     2311             CASE( 'I' )                      ! Ocean on C grid, Ice on I-point (B-grid) ==> T  
     2312                DO jj = 2, jpjm1  
     2313                   DO ji = 2, jpim1   ! NO vector opt.  
     2314                      zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj  ,1) ) * zfr_l(ji,jj)   &     
     2315                         &         + 0.25 * ( u_ice(ji+1,jj+1) + u_ice(ji,jj+1)                     &  
     2316                         &                  + u_ice(ji+1,jj  ) + u_ice(ji,jj  )  ) *  fr_i(ji,jj)  
     2317                      zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)      + vn(ji  ,jj-1,1) ) * zfr_l(ji,jj)   &   
     2318                         &         + 0.25 * ( v_ice(ji+1,jj+1) + v_ice(ji,jj+1)                     &  
     2319                         &                  + v_ice(ji+1,jj  ) + v_ice(ji,jj  )  ) *  fr_i(ji,jj)  
     2320                   END DO  
     2321                END DO  
     2322             CASE( 'F' )                      ! Ocean on C grid, Ice on F-point (B-grid) ==> T  
     2323                DO jj = 2, jpjm1  
     2324                   DO ji = 2, jpim1   ! NO vector opt.  
     2325                      zotx1(ji,jj) = 0.5  * ( un(ji,jj,1)      + un(ji-1,jj  ,1) ) * zfr_l(ji,jj)   &     
     2326                         &         + 0.25 * ( u_ice(ji-1,jj-1) + u_ice(ji,jj-1)                     &  
     2327                         &                  + u_ice(ji-1,jj  ) + u_ice(ji,jj  )  ) *  fr_i(ji,jj)  
     2328                      zoty1(ji,jj) = 0.5  * ( vn(ji,jj,1)      + vn(ji  ,jj-1,1) ) * zfr_l(ji,jj)   &   
     2329                         &         + 0.25 * ( v_ice(ji-1,jj-1) + v_ice(ji,jj-1)                     &  
     2330                         &                  + v_ice(ji-1,jj  ) + v_ice(ji,jj  )  ) *  fr_i(ji,jj)  
     2331                   END DO  
     2332                END DO  
     2333             END SELECT  
     2334          END SELECT  
     2335         CALL lbc_lnk( zotx1, ssnd(jps_ocxw)%clgrid, -1. )   ;   CALL lbc_lnk( zoty1, ssnd(jps_ocyw)%clgrid, -1. )  
     2336         !  
     2337         !  
     2338         IF( TRIM( sn_snd_crtw%clvor ) == 'eastward-northward' ) THEN             ! Rotation of the components  
     2339         !                                                                        ! Ocean component  
     2340            CALL rot_rep( zotx1, zoty1, ssnd(jps_ocxw)%clgrid, 'ij->e', ztmp1 )       ! 1st component   
     2341            CALL rot_rep( zotx1, zoty1, ssnd(jps_ocxw)%clgrid, 'ij->n', ztmp2 )       ! 2nd component   
     2342            zotx1(:,:) = ztmp1(:,:)                                                   ! overwrite the components   
     2343            zoty1(:,:) = ztmp2(:,:)   
     2344            IF( ssnd(jps_ivx1)%laction ) THEN                                     ! Ice component  
     2345               CALL rot_rep( zitx1, zity1, ssnd(jps_ivx1)%clgrid, 'ij->e', ztmp1 )    ! 1st component   
     2346               CALL rot_rep( zitx1, zity1, ssnd(jps_ivx1)%clgrid, 'ij->n', ztmp2 )    ! 2nd component   
     2347               zitx1(:,:) = ztmp1(:,:)                                                ! overwrite the components   
     2348               zity1(:,:) = ztmp2(:,:)  
     2349            ENDIF  
     2350         ENDIF  
     2351         !  
     2352!         ! spherical coordinates to cartesian -> 2 components to 3 components  
     2353!         IF( TRIM( sn_snd_crtw%clvref ) == 'cartesian' ) THEN  
     2354!            ztmp1(:,:) = zotx1(:,:)                     ! ocean currents  
     2355!            ztmp2(:,:) = zoty1(:,:)  
     2356!            CALL oce2geo ( ztmp1, ztmp2, 'T', zotx1, zoty1, zotz1 )  
     2357!            !  
     2358!            IF( ssnd(jps_ivx1)%laction ) THEN           ! ice velocities  
     2359!               ztmp1(:,:) = zitx1(:,:)  
     2360!               ztmp1(:,:) = zity1(:,:)  
     2361!               CALL oce2geo ( ztmp1, ztmp2, 'T', zitx1, zity1, zitz1 )  
     2362!            ENDIF  
     2363!         ENDIF  
     2364         !  
     2365         IF( ssnd(jps_ocxw)%laction )   CALL cpl_snd( jps_ocxw, isec, RESHAPE ( zotx1, (/jpi,jpj,1/) ), info )   ! ocean x current 1st grid  
     2366         IF( ssnd(jps_ocyw)%laction )   CALL cpl_snd( jps_ocyw, isec, RESHAPE ( zoty1, (/jpi,jpj,1/) ), info )   ! ocean y current 1st grid  
     2367         !   
     2368      ENDIF  
     2369      !  
     2370      IF( ssnd(jps_ficet)%laction ) THEN  
     2371         CALL cpl_snd( jps_ficet, isec, RESHAPE ( fr_i, (/jpi,jpj,1/) ), info )  
     2372      END IF  
     2373      !                                                      ! ------------------------- !  
     2374      !                                                      !   Water levels to waves   !  
     2375      !                                                      ! ------------------------- !  
     2376      IF( ssnd(jps_wlev)%laction ) THEN  
     2377         IF( ln_apr_dyn ) THEN   
     2378            IF( kt /= nit000 ) THEN   
     2379               ztmp1(:,:) = sshb(:,:) - 0.5 * ( ssh_ib(:,:) + ssh_ibb(:,:) )   
     2380            ELSE   
     2381               ztmp1(:,:) = sshb(:,:)   
     2382            ENDIF   
     2383         ELSE   
     2384            ztmp1(:,:) = sshn(:,:)   
     2385         ENDIF   
     2386         CALL cpl_snd( jps_wlev  , isec, RESHAPE ( ztmp1, (/jpi,jpj,1/) ), info )  
     2387      END IF  
    20652388      ! 
    20662389      !  Fields sent by OPA to SAS when doing OPA<->SAS coupling 
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/SBC/sbcisf.F90

    r6140 r7403  
    9090      INTEGER, INTENT( in ) :: kt                   ! ocean time step 
    9191      ! 
    92       INTEGER               :: ji, jj               ! loop index 
     92      INTEGER               :: ji, jj, jk           ! loop index 
     93      INTEGER               :: ikt, ikb             ! loop index 
    9394      REAL(wp), DIMENSION (:,:), POINTER :: zt_frz, zdep ! freezing temperature (zt_frz) at depth (zdep)  
     95      REAL(wp), DIMENSION(:,:,:), POINTER :: zfwfisf3d, zqhcisf3d, zqlatisf3d 
     96      REAL(wp), DIMENSION(:,:  ), POINTER :: zqhcisf2d 
    9497      !!--------------------------------------------------------------------- 
    9598      ! 
     
    161164         CALL lbc_lnk(risf_tsc(:,:,jp_tem),'T',1.) 
    162165         CALL lbc_lnk(risf_tsc(:,:,jp_sal),'T',1.) 
    163          CALL lbc_lnk(fwfisf(:,:)         ,'T',1.) 
    164          CALL lbc_lnk(qisf(:,:)           ,'T',1.) 
     166         CALL lbc_lnk(fwfisf(:,:)   ,'T',1.) 
     167         CALL lbc_lnk(qisf(:,:)     ,'T',1.) 
     168 
     169!============================================================================================================================================= 
     170         IF ( iom_use('fwfisf3d') .OR. iom_use('qlatisf3d') .OR. iom_use('qhcisf3d') .OR. iom_use('qhcisf')) THEN 
     171            CALL wrk_alloc( jpi,jpj,jpk, zfwfisf3d, zqhcisf3d, zqlatisf3d ) 
     172            CALL wrk_alloc( jpi,jpj,     zqhcisf2d                        ) 
     173 
     174            zfwfisf3d(:,:,:) = 0.0_wp                         ! 3d ice shelf melting (kg/m2/s) 
     175            zqhcisf3d(:,:,:) = 0.0_wp                         ! 3d heat content flux (W/m2) 
     176            zqlatisf3d(:,:,:)= 0.0_wp                         ! 3d ice shelf melting latent heat flux (W/m2) 
     177            zqhcisf2d(:,:)   = fwfisf(:,:) * zt_frz * rcp     ! 2d heat content flux (W/m2) 
     178 
     179            DO jj = 1,jpj 
     180               DO ji = 1,jpi 
     181                  ikt = misfkt(ji,jj) 
     182                  ikb = misfkb(ji,jj) 
     183                  DO jk = ikt, ikb - 1 
     184                     zfwfisf3d (ji,jj,jk) = zfwfisf3d (ji,jj,jk) + fwfisf   (ji,jj) * r1_hisf_tbl(ji,jj) * e3t_n(ji,jj,jk) 
     185                     zqhcisf3d (ji,jj,jk) = zqhcisf3d (ji,jj,jk) + zqhcisf2d(ji,jj) * r1_hisf_tbl(ji,jj) * e3t_n(ji,jj,jk) 
     186                     zqlatisf3d(ji,jj,jk) = zqlatisf3d(ji,jj,jk) + qisf     (ji,jj) * r1_hisf_tbl(ji,jj) * e3t_n(ji,jj,jk) 
     187                  END DO 
     188                  zfwfisf3d (ji,jj,jk) = zfwfisf3d (ji,jj,jk) + fwfisf   (ji,jj) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) * e3t_n(ji,jj,jk) 
     189                  zqhcisf3d (ji,jj,jk) = zqhcisf3d (ji,jj,jk) + zqhcisf2d(ji,jj) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) * e3t_n(ji,jj,jk) 
     190                  zqlatisf3d(ji,jj,jk) = zqlatisf3d(ji,jj,jk) + qisf     (ji,jj) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) * e3t_n(ji,jj,jk) 
     191               END DO 
     192            END DO 
     193 
     194            CALL iom_put('fwfisf3d' , zfwfisf3d (:,:,:)) 
     195            CALL iom_put('qlatisf3d', zqlatisf3d(:,:,:)) 
     196            CALL iom_put('qhcisf3d' , zqhcisf3d (:,:,:)) 
     197            CALL iom_put('qhcisf'   , zqhcisf2d (:,:  )) 
     198 
     199            CALL wrk_dealloc( jpi,jpj,jpk, zfwfisf3d, zqhcisf3d, zqlatisf3d ) 
     200            CALL wrk_dealloc( jpi,jpj,     zqhcisf2d                        ) 
     201         END IF 
     202 
     203         ! output 
     204         CALL iom_put('qlatisf'  , qisf) 
     205         CALL iom_put('fwfisf', fwfisf) 
     206!============================================================================================================================================= 
    165207 
    166208         IF( kt == nit000 ) THEN                         !   set the forcing field at nit000 - 1    ! 
     
    177219         END IF 
    178220         !  
    179          ! output 
    180          CALL iom_put('qisf'  , qisf) 
    181          CALL iom_put('fwfisf', fwfisf) 
    182  
    183221         ! deallocation 
    184222         CALL wrk_dealloc( jpi,jpj, zt_frz, zdep  ) 
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

    r6460 r7403  
    8989      NAMELIST/namsbc/ nn_fsbc  , ln_ana   , ln_flx, ln_blk_clio, ln_blk_core, ln_blk_mfs,   & 
    9090         &             ln_cpl   , ln_mixcpl, nn_components      , nn_limflx  ,               & 
    91          &             ln_traqsr, ln_dm2dc ,                                                 &   
     91         &             ln_traqsr, ln_dm2dc ,                                                 & 
    9292         &             nn_ice   , nn_ice_embd,                                               & 
    9393         &             ln_rnf   , ln_ssr   , ln_isf   , nn_fwb    , ln_apr_dyn,              & 
    94          &             ln_wave  ,                                                            & 
    95          &             nn_lsm    
     94         &             ln_wave  , ln_cdgw  , ln_sdw   , ln_tauoc  , ln_stcor  ,              & 
     95         &             nn_lsm 
    9696      INTEGER  ::   ios 
    9797      INTEGER  ::   ierr, ierr0, ierr1, ierr2, ierr3, jpm 
     
    153153         WRITE(numout,*) '              closed sea (=0/1) (set in namdom)          nn_closea     = ', nn_closea 
    154154         WRITE(numout,*) '              nb of iterations if land-sea-mask applied  nn_lsm        = ', nn_lsm 
    155          WRITE(numout,*) '              surface wave                               ln_wave       = ', ln_wave   
     155         WRITE(numout,*) '              surface wave                               ln_wave       = ', ln_wave 
     156         WRITE(numout,*) '                 Stokes drift corr. to vert. velocity    ln_sdw        = ', ln_sdw 
     157         WRITE(numout,*) '                 wave modified ocean stress              ln_tauoc      = ', ln_tauoc 
     158         WRITE(numout,*) '                 Stokes coriolis term                    ln_stcor      = ', ln_stcor 
     159         WRITE(numout,*) '                 neutral drag coefficient (CORE, MFS)    ln_cdgw       = ', ln_cdgw 
    156160      ENDIF 
    157161      ! 
     
    220224         &   CALL ctl_stop( 'diurnal cycle into qsr field from daily values requires a flux or core-bulk formulation' ) 
    221225       
     226      IF ( ln_wave ) THEN 
     227      !Activated wave module but neither drag nor stokes drift activated 
     228         IF ( .NOT.(ln_cdgw .OR. ln_sdw .OR. ln_tauoc .OR. ln_stcor ) )   THEN 
     229            CALL ctl_warn( 'Ask for wave coupling but ln_cdgw=F, ln_sdw=F, ln_tauoc=F, ln_stcor=F') 
     230      !drag coefficient read from wave model definable only with mfs bulk formulae and core  
     231         ELSEIF (ln_cdgw .AND. .NOT.(ln_blk_mfs .OR. ln_blk_core) )       THEN        
     232             CALL ctl_stop( 'drag coefficient read from wave model definable only with mfs bulk formulae and core') 
     233         ELSEIF (ln_stcor .AND. .NOT. ln_sdw)                             THEN 
     234             CALL ctl_stop( 'Stokes-Coriolis term calculated only if activated Stokes Drift ln_sdw=T') 
     235         ENDIF 
     236      ELSE 
     237      IF ( ln_cdgw .OR. ln_sdw .OR. ln_tauoc .OR. ln_stcor )                &  
     238         &   CALL ctl_stop( 'Not Activated Wave Module (ln_wave=F) but asked coupling ',    & 
     239         &                  'with drag coefficient (ln_cdgw =T) '  ,                        & 
     240         &                  'or Stokes Drift (ln_sdw=T) ' ,                                 & 
     241         &                  'or ocean stress modification due to waves (ln_tauoc=T) ',      &   
     242         &                  'or Stokes-Coriolis term (ln_stcori=T)'  ) 
     243      ENDIF  
    222244      !                          ! Choice of the Surface Boudary Condition (set nsbc) 
    223245      ll_purecpl = ln_cpl .AND. .NOT. ln_mixcpl 
     
    357379            &                CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice )   ! OPA-SAS coupling: OPA receiving fields from SAS 
    358380      END SELECT 
    359  
     381      IF ( ln_wave .AND. ln_tauoc) THEN                                 ! Wave stress subctracted 
     382            utau(:,:) = utau(:,:)*tauoc_wave(:,:) 
     383            vtau(:,:) = vtau(:,:)*tauoc_wave(:,:) 
     384            taum(:,:) = taum(:,:)*tauoc_wave(:,:) 
     385      ! 
     386            SELECT CASE( nsbc ) 
     387            CASE(  0,1,2,3,5,-1 )  ; 
     388                IF(lwp .AND. kt == nit000 ) WRITE(numout,*) 'WARNING: You are subtracting the wave stress to the ocean. & 
     389                        & If not requested select ln_tauoc=.false' 
     390            END SELECT 
     391      ! 
     392      END IF 
    360393      IF( ln_mixcpl )        CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice )   ! forced-coupled mixed formulation after forcing 
    361394 
  • branches/2016/dev_merge_2016/NEMOGCM/NEMO/OPA_SRC/SBC/sbcwave.F90

    r6140 r7403  
    44   !! Wave module  
    55   !!====================================================================== 
    6    !! History :  3.3  !   2011-09  (Adani M)  Original code: Drag Coefficient  
    7    !!         :  3.4  !   2012-10  (Adani M)                 Stokes Drift  
    8    !!---------------------------------------------------------------------- 
    9  
    10    !!---------------------------------------------------------------------- 
    11    !!   sbc_wave      : read drag coefficient from wave model in netcdf files  
     6   !! History :  3.3  !   2011-09  (M. Adani)  Original code: Drag Coefficient  
     7   !!         :  3.4  !   2012-10  (M. Adani)  Stokes Drift  
     8   !!            3.6  !   2014-09  (E. Clementi,P. Oddo) New Stokes Drift Computation 
     9   !!---------------------------------------------------------------------- 
     10 
     11   !!---------------------------------------------------------------------- 
     12   !!   sbc_wave      : wave data from wave model in netcdf files  
    1213   !!---------------------------------------------------------------------- 
    1314   USE oce            !  
    14    USE sbc_oce        ! Surface boundary condition: ocean fields 
     15   USE sbc_oce       ! Surface boundary condition: ocean fields 
    1516   USE bdy_oce        ! 
    1617   USE domvvl         ! 
    17    ! 
    1818   USE iom            ! I/O manager library 
    1919   USE in_out_manager ! I/O manager 
    2020   USE lib_mpp        ! distribued memory computing library 
    21    USE fldread        ! read input fields 
     21   USE fldread       ! read input fields 
    2222   USE wrk_nemo       ! 
     23   USE phycst         ! physical constants  
    2324 
    2425   IMPLICIT NONE 
    2526   PRIVATE 
    2627 
    27    PUBLIC   sbc_wave    ! routine called in sbc_blk_core or sbc_blk_mfs 
     28   PUBLIC   sbc_stokes, sbc_qiao  ! routines called in sbccpl 
     29   PUBLIC   sbc_wave    ! routine called in sbcmod 
    2830    
    29    INTEGER , PARAMETER ::   jpfld  = 3   ! maximum number of files to read for srokes drift 
    30    INTEGER , PARAMETER ::   jp_usd = 1   ! index of stokes drift  (i-component) (m/s)    at T-point 
    31    INTEGER , PARAMETER ::   jp_vsd = 2   ! index of stokes drift  (j-component) (m/s)    at T-point 
    32    INTEGER , PARAMETER ::   jp_wn  = 3   ! index of wave number                 (1/m)    at T-point 
    33  
    34    TYPE(FLD), ALLOCATABLE, DIMENSION(:)  :: sf_cd    ! structure of input fields (file informations, fields read) Drag Coefficient 
    35    TYPE(FLD), ALLOCATABLE, DIMENSION(:)  :: sf_sd    ! structure of input fields (file informations, fields read) Stokes Drift 
    36  
    37    REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION (:,:)   :: cdn_wave  
    38    REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION (:,:,:) :: usd3d, vsd3d, wsd3d  
    39    REAL(wp),         ALLOCATABLE, DIMENSION (:,:)   :: usd2d, vsd2d, uwavenum, vwavenum  
     31   ! Variables checking if the wave parameters are coupled (if not, they are read from file) 
     32   LOGICAL, PUBLIC     ::   cpl_hsig=.FALSE. 
     33   LOGICAL, PUBLIC     ::   cpl_phioc=.FALSE. 
     34   LOGICAL, PUBLIC     ::   cpl_sdrftx=.FALSE. 
     35   LOGICAL, PUBLIC     ::   cpl_sdrfty=.FALSE. 
     36   LOGICAL, PUBLIC     ::   cpl_wper=.FALSE. 
     37   LOGICAL, PUBLIC     ::   cpl_wnum=.FALSE. 
     38   LOGICAL, PUBLIC     ::   cpl_wstrf=.FALSE. 
     39   LOGICAL, PUBLIC     ::   cpl_wdrag=.FALSE. 
     40 
     41   INTEGER ::   jpfld                ! number of files to read for stokes drift 
     42   INTEGER ::   jp_usd               ! index of stokes drift  (i-component) (m/s)    at T-point 
     43   INTEGER ::   jp_vsd               ! index of stokes drift  (j-component) (m/s)    at T-point 
     44   INTEGER ::   jp_swh               ! index of significant wave hight      (m)      at T-point 
     45   INTEGER ::   jp_wmp               ! index of mean wave period            (s)      at T-point 
     46 
     47   TYPE(FLD), ALLOCATABLE, DIMENSION(:)  :: sf_cd    ! structure of input fields (file informations, fields read) Drag Coefficient 
     48   TYPE(FLD), ALLOCATABLE, DIMENSION(:)  :: sf_sd    ! structure of input fields (file informations, fields read) Stokes Drift 
     49   TYPE(FLD), ALLOCATABLE, DIMENSION(:)  :: sf_wn    ! structure of input fields (file informations, fields read) wave number for Qiao 
     50   TYPE(FLD), ALLOCATABLE, DIMENSION(:)  :: sf_tauoc ! structure of input fields (file informations, fields read) normalized wave stress into the ocean 
     51   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:)       :: cdn_wave  
     52   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:)       :: swh,wmp, wnum 
     53   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:)       :: tauoc_wave 
     54   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:)       :: tsd2d 
     55   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:)       :: zusd2dt, zvsd2dt 
     56   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:)     :: usd3d, vsd3d, wsd3d  
     57   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:)     :: usd3dt, vsd3dt 
    4058 
    4159   !! * Substitutions 
     
    4866CONTAINS 
    4967 
     68   SUBROUTINE sbc_stokes( ) 
     69      !!-------------------------------