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 5034 for branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG – NEMO

Ignore:
Timestamp:
2015-01-15T14:48:42+01:00 (9 years ago)
Author:
andrewryan
Message:

merge with trunk

Location:
branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG
Files:
40 edited
5 copied

Legend:

Unmodified
Added
Removed
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/AMM12/EXP00/iodef.xml

    r4343 r5034  
    8585   </file> 
    8686    
    87    <file id="file8" name_suffix="_icemod" description="ice variables" > 
    88      <field field_ref="ice_pres"                     /> 
    89      <field field_ref="snowthic_cea" name="snd"     long_name="surface_snow_thickness"   /> 
    90      <field field_ref="icethic_cea"  name="sit"     long_name="sea_ice_thickness"        /> 
    91      <field field_ref="iceprod_cea"  name="sip"     long_name="sea_ice_thickness"        /> 
    92      <field field_ref="ist_ipa"      /> 
    93      <field field_ref="ioceflxb"     /> 
    94      <field field_ref="uice_ipa"     /> 
    95      <field field_ref="vice_ipa"     /> 
    96      <field field_ref="utau_ice"     /> 
    97      <field field_ref="vtau_ice"     /> 
    98      <field field_ref="qsr_io_cea"   /> 
    99      <field field_ref="qns_io_cea"   /> 
    100      <field field_ref="snowpre"      /> 
    101    </file> 
    102     
    10387      </file_group> 
    10488 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/AMM12/EXP00/namelist_cfg

    r4493 r5034  
    1414   cn_ocerst_in   = "amm12_restart_oce"   !  suffix of ocean restart name (input) 
    1515   cn_ocerst_out  = "restart_oce_out"   !  suffix of ocean restart name (input) 
    16    nn_istate   =       1   !  output the initial state (1) or not (0) 
    1716   nn_stock    =    1296   !  frequency of creation of a restart file (modulo referenced to 1) 
    1817   nn_write    =     144   !  frequency of write in the output file   (modulo referenced to nit000)  
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/C1D_PAPA/EXP00/iodef.xml

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

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

    r4247 r5034  
    1  bld::tool::fppkeys key_c1d key_dynspg_flt key_diahth key_zdfgls key_iomput key_mpp_mpi 
     1 bld::tool::fppkeys key_c1d key_zdfgls 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE/EXP00/namelist_cfg

    r4370 r5034  
    1111   nn_stock    =    4320   !  frequency of creation of a restart file (modulo referenced to 1) 
    1212   nn_write    =      60   !  frequency of write in the output file   (modulo referenced to nn_it000) 
     13 
     14   ln_clobber  = .true.    !  clobber (overwrite) an existing file 
     15 
    1316/ 
    1417!----------------------------------------------------------------------- 
     
    8285!----------------------------------------------------------------------- 
    8386   nn_fsbc     = 1         !  frequency of surface boundary condition computation 
    84                            !     (also = the frequency of sea-ice model call) 
     87   !                       !     (also = the frequency of sea-ice model call) 
    8588   ln_ana      = .true.    !  analytical formulation                    (T => fill namsbc_ana ) 
    8689   ln_blk_core = .false.   !  CORE bulk formulation                     (T => fill namsbc_core) 
     
    198201&nameos        !   ocean physical parameters 
    199202!----------------------------------------------------------------------- 
    200    nn_eos      =   2       !  type of equation of state and Brunt-Vaisala frequency 
     203   nn_eos      =  0       !  type of equation of state and Brunt-Vaisala frequency 
     204                                 !  =-1, TEOS-10  
     205                                 !  = 0, EOS-80  
     206                                 !  = 1, S-EOS   (simplified eos) 
     207   ln_useCT    = .false.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 
     208   !                             ! 
     209   !                      ! S-EOS coefficients : 
     210   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
     211   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1) 
     212   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1) 
     213   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos) 
     214   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos) 
     215   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos) 
     216   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos) 
     217   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos) 
     218!!org GYRE   rn_alpha    =   2.0e-4  !  thermal expension coefficient (nn_eos= 1 or 2) 
     219!!org GYRE   rn_beta     =   7.7e-4  !  saline  expension coefficient (nn_eos= 2) 
     220!!org  caution  now a0 = alpha / rau0   with rau0 = 1026 
    201221/ 
    202222!----------------------------------------------------------------------- 
    203223&namtra_adv    !   advection scheme for tracer 
    204224!----------------------------------------------------------------------- 
    205    ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl  
     225   ln_traadv_cen2   =  .false.   !  2nd order centered scheme 
     226   ln_traadv_tvd    =  .true.    !  TVD scheme 
     227   ln_traadv_muscl  =  .false.   !  MUSCL scheme 
     228   ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries 
     229   ln_traadv_ubs    =  .false.   !  UBS scheme 
     230   ln_traadv_qck    =  .false.   !  QUICKEST scheme 
     231   ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl 
    206232/ 
    207233!----------------------------------------------------------------------- 
     
    242268!namdyn_spg    !   surface pressure gradient   (CPP key only) 
    243269!----------------------------------------------------------------------- 
     270 
    244271!----------------------------------------------------------------------- 
    245272&namdyn_ldf    !   lateral diffusion on momentum 
     
    301328!              !       or mixed-layer trends or barotropic vorticity    ("key_trdmld" or     "key_trdvor") 
    302329!----------------------------------------------------------------------- 
    303 / 
     330   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE 
     331   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output 
     332   ln_dyn_mxl  = .FALSE.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 
     333   ln_vor_trd  = .FALSE.   ! (T) 2D barotropic vorticity trends (not coded yet) 
     334   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends 
     335   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends 
     336   ln_tra_trd  = .false.    ! (T) 3D tracer trend output 
     337   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 
     338   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step) 
     339/ 
     340!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk) 
     341!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
     342!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input) 
     343!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output) 
     344!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics 
     345!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
     346!!gm 
    304347!----------------------------------------------------------------------- 
    305348&namflo       !   float parameters                                      ("key_float") 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE/cpp_GYRE.fcm

    r4230 r5034  
    1  bld::tool::fppkeys key_dynspg_flt key_ldfslp key_zdftke key_iomput key_mpp_mpi 
     1 bld::tool::fppkeys key_dynspg_flt key_ldfslp key_zdftke key_iomput key_mpp_mpi key_nosignedzero 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_cfg

    r4370 r5034  
    206206&nameos        !   ocean physical parameters 
    207207!----------------------------------------------------------------------- 
    208    nn_eos      =   2       !  type of equation of state and Brunt-Vaisala frequency 
     208   nn_eos      =  0       !  type of equation of state and Brunt-Vaisala frequency 
     209                                 !  =-1, TEOS-10 
     210                                 !  = 0, EOS-80 
     211                                 !  = 1, S-EOS   (simplified eos) 
     212   ln_useCT    = .false.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 
     213   !                             ! 
     214   !                      ! S-EOS coefficients : 
     215   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
     216   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1) 
     217   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1) 
     218   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos) 
     219   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos) 
     220   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos) 
     221   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos) 
     222   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos) 
     223!!org GYRE   rn_alpha    =   2.0e-4  !  thermal expension coefficient (nn_eos= 1 or 2) 
     224!!org GYRE   rn_beta     =   7.7e-4  !  saline  expension coefficient (nn_eos= 2) 
     225!!org  caution  now a0 = alpha / rau0   with rau0 = 1026 
    209226/ 
    210227!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_cfg

    r4370 r5034  
    121121&nameos        !   ocean physical parameters 
    122122!----------------------------------------------------------------------- 
    123    nn_eos      =   2       !  type of equation of state and Brunt-Vaisala frequency 
     123   nn_eos      =  0       !  type of equation of state and Brunt-Vaisala frequency 
     124                                 !  =-1, TEOS-10 
     125                                 !  = 0, EOS-80 
     126                                 !  = 1, S-EOS   (simplified eos) 
     127   ln_useCT    = .false.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 
     128   !                             ! 
     129   !                      ! S-EOS coefficients : 
     130   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
     131   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1) 
     132   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1) 
     133   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos) 
     134   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos) 
     135   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos) 
     136   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos) 
     137   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos) 
     138!!org GYRE   rn_alpha    =   2.0e-4  !  thermal expension coefficient (nn_eos= 1 or 2) 
     139!!org GYRE   rn_beta     =   7.7e-4  !  saline  expension coefficient (nn_eos= 2) 
     140!!org  caution  now a0 = alpha / rau0   with rau0 = 1026 
    124141/ 
    125142!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/GYRE_PISCES/cpp_GYRE_PISCES.fcm

    r4230 r5034  
    1 bld::tool::fppkeys  key_dynspg_flt key_ldfslp key_zdftke key_top key_pisces_reduced key_iomput key_mpp_mpi  
     1bld::tool::fppkeys  key_dynspg_flt key_ldfslp key_zdftke key_top key_pisces_reduced key_iomput key_mpp_mpi 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/1_namelist_cfg

    r4329 r5034  
    5555/ 
    5656!----------------------------------------------------------------------- 
     57&namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
     58!----------------------------------------------------------------------- 
     59/ 
     60!----------------------------------------------------------------------- 
     61&namcrs        !   Grid coarsening for dynamics output and/or 
     62               !   passive tracer coarsened online simulations 
     63!----------------------------------------------------------------------- 
     64/ 
     65!----------------------------------------------------------------------- 
    5766&namtsd    !   data : Temperature  & Salinity 
    5867!----------------------------------------------------------------------- 
     
    7281/ 
    7382!----------------------------------------------------------------------- 
    74 &namsbc_core   !   namsbc_core  CORE bulk formulea 
    75 !----------------------------------------------------------------------- 
     83&namsbc_core   !   namsbc_core  CORE bulk formulae 
     84!----------------------------------------------------------------------- 
     85!              !  file name                    ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                               ! rotation ! land/sea mask ! 
     86!              !                               !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                              ! pairing  ! filename      ! 
     87   sn_wndi     = 'u_10.15JUNE2009_fill'        ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bicubic.nc'      , 'Uwnd'   , '' 
     88   sn_wndj     = 'v_10.15JUNE2009_fill'        ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bicubic.nc'      , 'Vwnd'   , '' 
     89   sn_qsr      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     90   sn_qlw      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     91   sn_tair     = 't_10.15JUNE2009_fill'        ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     92   sn_humi     = 'q_10.15JUNE2009_fill'        ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     93   sn_prec     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     94   sn_snow     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     95   sn_tdif     = 'taudif_core'                 ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     96   cn_dir      = './'      !  root directory for the location of the bulk files 
     97   ln_2m       = .false.   !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
     98   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
     99   ln_bulk2z   = .false.   !  Air temperature/humidity and wind vectors are referenced at heights rn_zqt and rn_zu 
     100   rn_zqt      = 3.        !  Air temperature and humidity reference height (m) (ln_bulk2z) 
     101   rn_zu       = 4.        !  Wind vector reference height (m)                  (ln_bulk2z) 
     102   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
     103   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
     104   rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity  
     105                           !  in the calculation of the wind stress (0.=absolute winds or 1.=relative winds) 
    76106/ 
    77107!----------------------------------------------------------------------- 
    78108&namtra_qsr    !   penetrative solar radiation 
    79109!----------------------------------------------------------------------- 
     110!              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
     111!              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
     112   sn_chl      ='chlorophyll',        -1         , 'CHLA'    ,   .true.     , .true. , 'yearly'  , 'weights_bilin.nc'       , '' 
    80113/ 
    81114!----------------------------------------------------------------------- 
     
    106139&nameos        !   ocean physical parameters 
    107140!----------------------------------------------------------------------- 
     141/ 
     142!----------------------------------------------------------------------- 
     143&namtra_dmp    !   tracer: T & S newtonian damping 
     144!----------------------------------------------------------------------- 
     145   ln_tradmp   =  .false.   !  add a damping termn (T) or not (F) 
    108146/ 
    109147!----------------------------------------------------------------------- 
     
    183221!----------------------------------------------------------------------- 
    184222/ 
     223!----------------------------------------------------------------------- 
     224&namobs       !  observation usage                                      ('key_diaobs') 
     225!----------------------------------------------------------------------- 
     226/ 
     227!----------------------------------------------------------------------- 
     228&nam_asminc   !   assimilation increments                               ('key_asminc') 
     229!----------------------------------------------------------------------- 
     230/ 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef.xml

    r4272 r5034  
    8282     <field field_ref="iceprod_cea"  name="sip"     long_name="sea_ice_thickness"        /> 
    8383     <field field_ref="ist_ipa"      /> 
    84      <field field_ref="ioceflxb"     /> 
    8584     <field field_ref="uice_ipa"     /> 
    8685     <field field_ref="vice_ipa"     /> 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist_cfg

    r4370 r5034  
    188188!----------------------------------------------------------------------- 
    189189/ 
     190!----------------------------------------------------------------------- 
     191&namobs       !  observation usage                                      ('key_diaobs') 
     192!----------------------------------------------------------------------- 
     193/ 
     194!----------------------------------------------------------------------- 
     195&nam_asminc   !   assimilation increments                               ('key_asminc') 
     196!----------------------------------------------------------------------- 
     197/ 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/README

    r4394 r5034  
    2424 
    2525 
    26  
    27  
     26#---------------------------------------------------------------------- 
     27# HOW TO USE  
     28#---------------------------------------------------------------------- 
     29./std_main.sh -plot -pdf 
     30or 
     31./std_main.sh -ts -pdf 
    2832 
    2933 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceFrac.pro

    r4559 r5034  
    2121   
    2222  IF keyword_set(arc) THEN BEGIN  
    23     domdef, 20, 380, 50, 90 
     23    domdef, 20, 380, 40, 90 
    2424    map = [90, 0, 0] 
    2525  ENDIF 
    2626  IF keyword_set(ant) THEN BEGIN  
    27     domdef, 20, 380, -90, -50 
     27    domdef, 20, 380, -90, -40 
    2828    map = [-90, 0, 0] 
    2929  ENDIF 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro

    r4559 r5034  
    172172  Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 
    173173  jan = where(mm EQ 1, cnt) 
     174  Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    174175  sept = where(mm EQ 9, cnt) 
    175   Ithi_jan_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(jan)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    176176  Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit}  
    177177  undefine, Ithi_1 
     
    288288     Isal_sept_2 = {arr:-1} 
    289289   ENDELSE 
    290 ; 
     290 
    291291  IvelU_1 = read_ncdf(getenv('VAR1_IvelU'), allrecords = allrec, filename = std_file1_I )  
    292292  IvelV_1 = read_ncdf(getenv('VAR1_IvelV'), allrecords = allrec, filename = std_file1_I )  
     
    335335    Ivelo_sept_2 = {arr:-1} 
    336336  ENDELSE 
    337 ; 
     337 
    338338  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )  
    339339  caldat, time, mm 
     
    396396  cnt = cnt+1   &   blabla = 'Erp salinity damping term' 
    397397  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex 
    398 ; 
     398;; 
    399399  cnt = cnt+1   &   blabla = 'Evaporation - Precipitation - Runoff term' 
    400400  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex 
     
    439439  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD2, MLD, _extra = ex 
    440440; 
     441  cnt = cnt+1   &   blabla = 'Mixed layer depth ortho plan ARCTIC' 
     442  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_MLD_ortho, MLD1, MLD2, /ARC, _extra = ex 
     443;sf 
     444  cnt = cnt+1   &   blabla = 'Mixed layer depth ortho plan ANTARTIC' 
     445  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_MLD_ortho, MLD1, MLD2, /ANT, _extra = ex 
     446;sf 
    441447  cnt = cnt+1   &   blabla = 'Zonal mean Mixed layer depth' 
    442448  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ZonMld, MLD1, MLD2, MLD, _extra = ex 
     
    564570  cnt = cnt+1   &   blabla = 'Antartic Ice Velocity: SEPT' 
    565571  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceVel, IvelU_sept_1, IvelU_sept_2, IvelV_sept_1, IvelV_sept_2, Ivelo_sept_1, Ivelo_sept_2, /ANT,  /SEPT, _extra = ex 
    566 ; 
     572 
    567573  cnt = cnt+1   &   blabla = 'Arctic Ice Age: MARCH' 
    568574  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh

    r4562 r5034  
    2121#idl_command=idl71 
    2222# 
    23 PS_DIR=$( pwd )/OR2LPv36_ps_plot 
    24 PDF_DIR=$( pwd )/OR2LPv36_pdf_plot 
     23PS_DIR=$( pwd )/p4H25a50-testht_ps_plot 
     24PDF_DIR=$( pwd )/p4H25a50-testht_pdf_plot 
    2525HTML_DIR=$( pwd )/html_plot 
    2626SAXO_DIR=/Users/sflod/SAXO_DIR 
    2727# 
    28 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/500y_pisces     # path of data in NetCDF format 
     28DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format 
    2929DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
    3030DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
     
    5454# 
    5555#===================== EXP1 ===================== 
    56 # 23610101_23701231 
    57 FILE1_T=OR2LPv36_23610101_23701231_1Y_grid_T.nc   # exp1 grid T input file 
    58 FILE1_U=OR2LPv36_23610101_23701231_1Y_grid_U.nc   # exp1 grid U input file 
    59 FILE1_V=OR2LPv36_23610101_23701231_1Y_grid_V.nc   # exp1 grid V input file 
    60 FILE1_I=OR2LPv36_23610101_23701231_1M_icemod.nc   # exp1 ice    input file 
     56# 
     57FILE1_T=p4H25a50_20410101_20501231_1Y_grid_T.nc   # exp1 grid T input file 
     58FILE1_U=p4H25a50_20410101_20501231_1Y_grid_U.nc   # exp1 grid U input file 
     59FILE1_V=p4H25a50_20410101_20501231_1Y_grid_V.nc   # exp1 grid V input file 
     60FILE1_I=p4H25a50_20410101_20501231_1M_icemod.nc   # exp1 ice    input file 
    6161VAR1_T=thetao 
    6262VAR1_S=so         
    63 VAR1_QNET=tohfls 
    64 VAR1_ERP=sowafldp       
     63VAR1_QNET=qt 
     64VAR1_ERP=wfcorr 
    6565VAR1_EMP=wfo 
    6666VAR1_MLD=mldr10_1 
    6767VAR1_U=uocetr_eff 
    6868VAR1_V=vocetr_eff 
    69 VAR1_Ithick=iicethic 
    70 VAR1_Ifrac=soicecov 
    71 VAR1_Isnow=isnowthi 
    72 #VAR1_Isal=sisali 
    73 #VAR1_Iage=siages 
    74 VAR1_IvelU=iicevelu 
    75 VAR1_IvelV=iicevelv 
    76 #VAR1_Ivelo=sivelo 
     69VAR1_Ithick=sithic 
     70VAR1_Ifrac=siconc 
     71VAR1_Isnow=snvolu 
     72VAR1_Isal=sisali 
     73VAR1_Iage=siages 
     74VAR1_IvelU=sivelu 
     75VAR1_IvelV=sivelv 
     76VAR1_Ivelo=sivelo 
    7777# 
    7878#===================== EXP2 ===================== 
    7979# 
    80 FILE2_T=OR2LPv36_23610101_23701231_1Y_grid_T.nc   # exp1 grid T input file 
    81 FILE2_U=OR2LPv36_23610101_23701231_1Y_grid_U.nc   # exp1 grid U input file 
    82 FILE2_V=OR2LPv36_23610101_23701231_1Y_grid_V.nc   # exp1 grid V input file 
    83 FILE2_I=OR2LPv36_23610101_23701231_1M_icemod.nc   # exp1 ice    input file 
     80FILE2_T=testht_20410101_20501231_1Y_grid_T.nc   # exp1 grid T input file 
     81FILE2_U=testht_20410101_20501231_1Y_grid_U.nc   # exp1 grid U input file 
     82FILE2_V=testht_20410101_20501231_1Y_grid_V.nc   # exp1 grid V input file 
     83FILE2_I=testht_20410101_20501231_1M_icemod.nc   # exp1 ice    input file 
    8484VAR2_T=thetao 
    8585VAR2_S=so 
    86 VAR2_QNET=tohfls 
    87 VAR2_ERP=sowafldp 
     86VAR2_QNET=qt 
     87VAR2_ERP=wfcorr 
    8888VAR2_EMP=wfo 
    8989VAR2_MLD=mldr10_1 
    9090VAR2_U=uocetr_eff 
    9191VAR2_V=vocetr_eff 
    92 VAR2_Ithick=iicethic 
    93 VAR2_Ifrac=soicecov 
    94 VAR2_Isnow=isnowthi 
    95 #VAR2_Isal=sisali 
    96 #VAR2_Iage=siages 
    97 VAR2_IvelU=iicevelu 
    98 VAR2_IvelV=iicevelv 
    99 #VAR2_Ivelo=sivelo 
     92VAR2_Ithick=sithic 
     93VAR2_Ifrac=siconc 
     94VAR2_Isnow=snvolu 
     95VAR2_Isal=sisali 
     96VAR2_Iage=siages 
     97VAR2_IvelU=sivelu 
     98VAR2_IvelV=sivelv 
     99VAR2_Ivelo=sivelo 
    100100# 
    101101######################### Export Variables ############################### 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.old_names

    r4559 r5034  
    7474VAR1_IvelU=iicevelu 
    7575VAR1_IvelV=iicevelv 
    76 VAR1_Ivelo=iicevelo 
     76#VAR1_Ivelo=iicevelo 
    7777# 
    7878#===================== EXP2 ===================== 
     
    9999VAR2_IvelU=iicevelu 
    100100VAR2_IvelV=iicevelv 
    101 VAR2_Ivelo=iicevelo 
     101#VAR2_Ivelo=iicevelo 
    102102# 
    103103######################### Export Variables ############################### 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE.pro

    r4497 r5034  
    8888  ICE_S_15 = {arr:ICE_S_15 * 1.e-12, unit : '10^12 m^2'} 
    8989  ; 
    90   title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
     90  ;;title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
     91  title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (DASHED) '+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 
    9192  jpt=12 
    9293  time=julday(1,15,1900)+30*lindgen(12) 
    93   pltt, ICE_N, 't', MIN = 4., MAX = 16., /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
     94  pltt, ICE_N, 't', MIN = 4., MAX = 16., /REMPLI, /PORTRAIT, LINESTYLE=2, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
    9495       , COLOR = 000 , small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    95   pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $       ;;; dashed lines is LINESTYLE=2  $ 
     96  pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT $       ;;; dashed lines is LINESTYLE=2  $ 
    9697        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    97   pltt, vICE_area_NH, 't',  /REMPLI, /PORTRAIT  $  
     98  pltt, vICE_area_NH, 't',  /REMPLI, /PORTRAIT, LINESTYLE=2  $  
    9899        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    99   pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $   ;;; dashed lines is LINESTYLE=2  $ 
     100  pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT $   ;;; dashed lines is LINESTYLE=2  $ 
    100101        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    101102; 
    102   title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
    103   pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, /NOERASE , XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
     103  title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) '+d1_d2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (DASHED) '+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 
     104  pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, LINESTYLE=2, /NOERASE , XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
    104105        ,COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    105   pltt, ICE_S_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     106  pltt, ICE_S_15, 't', /REMPLI, /PORTRAIT $  
    106107        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    107   pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT  $  
     108  pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2  $  
    108109         , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    109   pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     110  pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT $  
    110111        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    111112; 
     
    159160 
    160161 
    161   title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
     162  ;;title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (CONTINUOUS) '+'!C'+ 'and Extend minus 15% (DASHED)' 
     163  title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+' Global Annual Mean Ice Area (DASHED) '+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 
    162164  jpt=12 
    163165  time=julday(1,15,1900)+30*lindgen(12) 
    164   pltt, ICE_N, 't', MIN = 4, MAX = 16,  /REMPLI, /PORTRAIT, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $ 
     166  pltt, ICE_N, 't', MIN = 4, MAX = 16,  /REMPLI, /PORTRAIT, LINESTYLE=2, XGRIDSTYLE = 1, window = 2, DATE_FORMAT = '%M' $ 
    165167        , COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex ; BLACK 
    166   pltt, ICE_N2, 't', /REMPLI, /PORTRAIT $ 
     168  pltt, ICE_N2, 't', /REMPLI, /PORTRAIT , LINESTYLE=2 $ 
    167169        , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; RED 
    168   pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     170  pltt, ICE_N_15, 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    169171        , /ov1d, COLOR = 000, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  
    170   pltt, ICE_N2_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     172  pltt, ICE_N2_15, 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    171173        , /ov1d, COLOR = 250, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    172   pltt, vICE_area_NH, 't', /REMPLI, /PORTRAIT  $  
     174  pltt, vICE_area_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2  $  
    173175        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; light blue 
    174   pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     176  pltt, vICE_ext_NH, 't', /REMPLI, /PORTRAIT $  
    175177        , /ov1d, COLOR = 100, small = [1, 2, 1], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex  ; blu scuro 
    176178; 
    177   title ='Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+'Global Annual Mean Ice Area (CONTINUS)'+'!C'+ 'and Extend minus 15% (DASHED)' 
     179  title ='Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'OBSERVATION (light blue) '+'!C'+'Global Annual Mean Ice Area (DASHED)'+'!C'+ 'and Extend minus 15% (CONTINUOUS)' 
    178180;  title ='Southern Hemisphere'+'!C' 
    179   pltt, ICE_S, 't', MIN = 0., MAX = 20.,  /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
     181  pltt, ICE_S, 't', MIN = 0., MAX = 20., /REMPLI, LINESTYLE=2, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
    180182         , COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    181   pltt, ICE_S2, 't', /REMPLI, /NOERASE $ 
     183  pltt, ICE_S2, 't', /REMPLI, /NOERASE, LINESTYLE=2 $ 
    182184        , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    183   pltt, ICE_S_15 , 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     185  pltt, ICE_S_15 , 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    184186        , /ov1d, COLOR = 000, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    185   pltt, ICE_S2_15, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $ ; linee tratteggiate LINESTYLE=2  $ 
     187  pltt, ICE_S2_15, 't', /REMPLI, /PORTRAIT $ ; linee tratteggiate LINESTYLE=2  $ 
    186188        , /ov1d, COLOR = 250, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    187   pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT $  
     189  pltt,  vICE_area_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
    188190          , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    189   pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT, LINESTYLE=2 $  
     191  pltt,  vICE_ext_SH, 't', /REMPLI, /PORTRAIT $  
    190192        , /ov1d, COLOR = 100, small = [1, 2, 2], YTITLE = '10^12 m^2 ', TITLE = title, _extra = ex 
    191193; 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro

    r4514 r5034  
    55arr2d = arr2d[t1:t2] 
    66arr2d = reform(arr2d,12,nyear) ; put in 2D array 
    7 arr2d = total(arr2d,2)/nyear ; total over 2th dimension (i.e.years) 
     7;arr2d = total(arr2d,2)/nyear ; total over 2th dimension (i.e.years) 
     8arr2d = arr2d[*, nyear-1] ;  select last year 
    89 
    910return, arr2d 
     
    3435; 
    3536  d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')' 
     37  d2 = '('+strtrim(date2, 1)+')' 
    3638; 
    3739  iodir = std_iodir_data 
     
    7476  ;; AREA 
    7577  ICE_area_export = reform(ICE_area_export, 12, nyr) 
    76   ICE_area_export = total(ICE_area_export,2)/nyr 
     78  ;ICE_area_export = total(ICE_area_export,2)/nyr  ; old version monthly mean average  over all years 
     79  ICE_area_export = ICE_area_export[*,nyr-1]  ; dim= 12, index 0-11; last year choosen 
    7780  ; ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'}  ; annual mean 
    7881  ICE_area_export = {arr:ICE_area_export * 1.e-12 * 86400 * 30 , unit : '10^6 Km^2/month'}    ; monthly mean 
     
    8083  ; 
    8184  ICE_vol_export = reform(ICE_vol_export, 12, nyr) 
    82   ICE_vol_export = total(ICE_vol_export,2)/nyr 
     85  ;ICE_vol_export = total(ICE_vol_export,2)/nyr 
     86  ICE_vol_export = ICE_vol_export[*,nyr-1] 
    8387  ; ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 365 , unit : '10^3 Km^3/year'}      ; annual mean 
    8488  ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'}      ; monthly mean 
    8589 
    8690  ; 
    87   title = 'Fram Strait Areal Export'+'!C'+prefix+' '+d1_d2 
     91  ;title = 'Fram Strait Areal Export: LAST YEAR'+'!C'+prefix+' '+d1_d2 
     92  title = 'Fram Strait Areal Export: LAST YEAR'+'!C'+prefix+' '+d2 
    8893  jpt=12 
    8994  time=julday(1,15,1900)+30*lindgen(12) 
     
    99104  xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
    100105  ; 
    101   title = 'Fram Strait Volume Export'+'!C'+prefix+' '+d1_d2 
     106  ;title = 'Fram Strait Volume Export LAST YEAR'+'!C'+prefix+' '+d1_d2 
     107  title = 'Fram Strait Volume Export LAST YEAR'+'!C'+prefix+' '+d2 
    102108  pltt, ICE_vol_export, 't',  /REMPLI, MIN = 60., MAX = 500. , /NOERASE, XGRIDSTYLE = 1 $ 
    103109       , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 
     
    118124 
    119125    d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')' 
     126    d2_2 = '('+strtrim(date2_2, 1)+')' 
    120127    tsave = time 
    121128    domdef, indx1, indx2, indy1, indy2, /xindex, /yindex,/memeindices 
     
    140147    ;; AREA 
    141148    ICE_area_export_2 = reform(ICE_area_export_2, 12, nyr) 
    142     ICE_area_export_2 = total(ICE_area_export_2,2)/nyr 
     149    ;ICE_area_export_2 = total(ICE_area_export_2,2)/nyr 
     150    ICE_area_export_2 = ICE_area_export_2[*,nyr-1] 
    143151    ; ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 365 , unit : '10^6 Km^2/year'}   ; annual mean 
    144152    ICE_area_export_2 = {arr:ICE_area_export_2 * 1.e-12 * 86400 * 30 , unit : '10^6 Km^2/month'}     ; monthly mean 
    145153    ; 
    146154    ICE_vol_export_2 = reform(ICE_vol_export_2, 12, nyr) 
    147     ICE_vol_export_2 = total(ICE_vol_export_2,2)/nyr 
     155    ;ICE_vol_export_2 = total(ICE_vol_export_2,2)/nyr 
     156    ICE_vol_export_2 = ICE_vol_export_2[*,nyr-1] 
    148157    ; ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'}     ; annual mean 
    149158    ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'}       ; monthly mean 
     
    152161    if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1 
    153162 
    154     title = 'Fram Strait Areal Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
     163    ;title = 'Fram Strait Areal Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
     164    title = 'Fram Strait Areal Export LAST YEAR'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d2_2 
    155165    jpt=12 
    156166    time=julday(1,15,1900)+30*lindgen(12) 
     
    170180    xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
    171181    ; 
    172     title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
     182    ;title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
     183    title = 'Fram Strait Volume Export LAST YEAR'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d2_2 
    173184    pltt, ICE_vol_export, 't', /REMPLI,  MIN = 60., MAX = 500., /NOERASE, XGRIDSTYLE = 1 $ 
    174185        , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vol.pro

    r4504 r5034  
    7070; 
    7171  title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (BLACK)' 
    72   pltt, ICE_vol_S, 't', MIN = 0., MAX = 11000., /REMPLI, /NOERASE, XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
     72  pltt, ICE_vol_S, 't', MIN = 0., MAX = 12000., /REMPLI, /NOERASE, XGRIDSTYLE = 1 , DATE_FORMAT = '%M' $ 
    7373        , small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex 
    7474; 
     
    119119    ; 
    120120    title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2+'!C'+'Global Annual Mean Ice Volume' 
    121     pltt, ICE_vol_S, 't', MIN = 0., MAX = 11000., /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
     121    pltt, ICE_vol_S, 't', MIN = 0., MAX = 12000., /REMPLI, /NOERASE, XGRIDSTYLE = 1, DATE_FORMAT = '%M' $ 
    122122          , small = [1, 2, 2], YTITLE = '10^9 Km^3 ', TITLE = title, _extra = ex    
    123123    pltt, ICE_vol_S2, 't',  /REMPLI, /NOERASE $ 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro

    r4559 r5034  
    6060  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_SSH, masknp, POSTSCRIPT = postscript, _extra = ex 
    6161; 
    62   cnt = cnt+1   &   blabla = 'Global Mean Q net' 
    63   IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_Q, masknp, POSTSCRIPT = postscript, _extra = ex 
     62;;;;  cnt = cnt+1   &   blabla = 'Global Mean Q net' 
     63;;;;  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_Q, masknp, POSTSCRIPT = postscript, _extra = ex 
    6464; 
    6565  cnt = cnt+1   &   blabla = 'Global Mean EMP' 
     
    7777  cnt = cnt+1   &   blabla = 'Sea-Ice Volume' 
    7878  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 
    79 ; 
    80   cnt = cnt+1   &   blabla = 'ICE Velocity ' 
    81   IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex 
    82 ; 
     79;SF; and probably useless 
     80;SF  cnt = cnt+1   &   blabla = 'ICE Velocity ' 
     81;SF  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex 
     82;SF; 
    8383  cnt = cnt+1   &   blabla = 'Sea-Ice Volume Export at Fram Strait' 
    8484  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_FRAM, masknp, POSTSCRIPT = postscript, _extra = ex 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_vardef.sh

    r4559 r5034  
    1818idl_command=/Applications/itt/idl64/bin/idl 
    1919# 
    20 PS_DIR=$( pwd )/ts_heatbug_ps 
    21 PDF_DIR=$( pwd )/ts_heatbug_pdf 
     20PS_DIR=$( pwd )/ts_p4H25a50-testht_ps 
     21PDF_DIR=$( pwd )/ts_p4H25a50-testht_pdf 
    2222HTML_DIR=$( pwd )/html 
    2323SAXO_DIR=/Users/sflod/SAXO_DIR 
    2424# 
    25 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada      # path of data in NetCDF format 
     25DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format 
    2626DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
    2727DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
     
    5252READ_ONLY_FIRST_RECORD=0   # if 0 then read all records in files else read only the first reacord in each file 
    5353# 
    54 #DATE1=20010101   ;   DATE2=25001231 
    55 DATE1=20310101   ;   DATE2=20401231 
    56 #DATE1=00010101      ;   DATE2=00101231 
     54DATE1=20010101         ;   DATE2=20501231 
    5755# 
    58 VAR1_T=thetao        ;   V1T_PREF=heatbug       ;   V1T_SUFF=_1Y_grid_T.nc 
    59 VAR1_S=so            ;   V1S_PREF=heatbug      ;   V1S_SUFF=_1Y_grid_T.nc 
    60 VAR1_SSH=zos         ;   V1SSH_PREF=heatbug    ;   V1SSH_SUFF=_1Y_grid_T.nc 
    61 VAR1_Q=tohfls        ;   V1Q_PREF=heatbug      ;   V1Q_SUFF=_1Y_grid_T.nc 
    62 VAR1_EMP=wfo         ;   V1EMP_PREF=heatbug    ;   V1EMP_SUFF=_1Y_grid_T.nc 
    63 VAR1_U=uocetr_eff    ;   V1U_PREF=heatbug      ;   V1U_SUFF=_1Y_grid_U.nc 
    64 VAR1_V=vocetr_eff    ;   V1V_PREF=heatbug      ;   V1V_SUFF=_1Y_grid_V.nc 
    65 VAR1_ICE=siconc      ;   V1ICE_PREF=heatbug    ;   V1ICE_SUFF=_1M_icemod.nc 
    66 VAR1_Ithick=sithic   ;   V1It_PREF=heatbug     ;   V1It_SUFF=_1M_icemod.nc 
    67 VAR1_SNOW=sndept     ;   V1SNOW_PREF=heatbug   ;   V1SNOW_SUFF=_1M_icemod.nc 
    68 VAR1_IvelV=sivelv    ;   V1IvV_PREF=heatbug    ;   V1IvV_SUFF=_1M_icemod.nc 
    69 VAR1_Ivel=sivelo     ;   V1Iv_PREF=heatbug     ;   V1Iv_SUFF=_1M_icemod.nc 
     56VAR1_T=thetao          ;   V1T_PREF=p4H25a50      ;   V1T_SUFF=_1Y_grid_T.nc 
     57VAR1_S=so              ;   V1S_PREF=p4H25a50      ;   V1S_SUFF=_1Y_grid_T.nc 
     58VAR1_SSH=zos         ;   V1SSH_PREF=p4H25a50    ;   V1SSH_SUFF=_1Y_grid_T.nc 
     59##VAR1_Q=qt              ;   V1Q_PREF=p4H25a50      ;   V1Q_SUFF=_1Y_grid_T.nc 
     60VAR1_EMP=wfo         ;   V1EMP_PREF=p4H25a50    ;   V1EMP_SUFF=_1Y_grid_T.nc 
     61VAR1_U=uocetr_eff      ;   V1U_PREF=p4H25a50      ;   V1U_SUFF=_1Y_grid_U.nc 
     62VAR1_V=vocetr_eff      ;   V1V_PREF=p4H25a50      ;   V1V_SUFF=_1Y_grid_V.nc 
     63VAR1_ICE=siconc      ;   V1ICE_PREF=p4H25a50    ;   V1ICE_SUFF=_1M_icemod.nc 
     64VAR1_Ithick=sithic    ;   V1It_PREF=p4H25a50     ;   V1It_SUFF=_1M_icemod.nc 
     65VAR1_SNOW=snthic    ;   V1SNOW_PREF=p4H25a50   ;   V1SNOW_SUFF=_1M_icemod.nc 
     66VAR1_IvelV=sivelv    ;   V1IvV_PREF=p4H25a50    ;   V1IvV_SUFF=_1M_icemod.nc  
     67VAR1_Ivel=sivelo      ;   V1Iv_PREF=p4H25a50     ;   V1Iv_SUFF=_1M_icemod.nc 
    7068# 
    7169#===================== EXP2 ===================== 
    7270# 
    73 #DATE1_2=20010101   ;   DATE2_2=20101231 
    74 DATE1_2=20310101   ;   DATE2_2=20401231 
    75 #DATE1_2=20010101   ;   DATE2_2=25001231 
     71DATE1_2=20010101       ;   DATE2_2=20501231 
    7672# 
    77 VAR2_T=thetao      ;   V2T_PREF=heatbug      ;   V2T_SUFF=_1Y_grid_T.nc 
    78 VAR2_S=so          ;   V2S_PREF=heatbug      ;   V2S_SUFF=_1Y_grid_T.nc 
    79 VAR2_SSH=zos       ;   V2SSH_PREF=heatbug    ;   V2SSH_SUFF=_1Y_grid_T.nc 
    80 VAR2_Q=tohfls      ;   V2Q_PREF=heatbug      ;   V2Q_SUFF=_1Y_grid_T.nc 
    81 VAR2_EMP=wfo       ;   V2EMP_PREF=heatbug    ;   V2EMP_SUFF=_1Y_grid_T.nc 
    82 VAR2_U=uocetr_eff  ;   V2U_PREF=heatbug      ;   V2U_SUFF=_1Y_grid_U.nc 
    83 VAR2_V=vocetr_eff  ;   V2V_PREF=heatbug      ;   V2V_SUFF=_1Y_grid_V.nc 
    84 VAR2_ICE=siconc    ;   V2ICE_PREF=heatbug    ;   V2ICE_SUFF=_1M_icemod.nc 
    85 VAR2_Ithick=sithic ;   V2It_PREF=heatbug     ;   V2It_SUFF=_1M_icemod.nc 
    86 VAR2_SNOW=sndept   ;   V2SNOW_PREF=heatbug   ;   V2SNOW_SUFF=_1M_icemod.nc 
    87 VAR2_IvelV=sivelv  ;   V2IvV_PREF=heatbug    ;   V2IvV_SUFF=_1M_icemod.nc 
    88 VAR2_Ivel=sivelo   ;   V2Iv_PREF=heatbug     ;   V2Iv_SUFF=_1M_icemod.nc 
     73VAR2_T=thetao          ;   V2T_PREF=testht     ;   V2T_SUFF=_1Y_grid_T.nc 
     74VAR2_S=so           ;   V2S_PREF=testht     ;   V2S_SUFF=_1Y_grid_T.nc 
     75VAR2_SSH=zos         ;   V2SSH_PREF=testht   ;   V2SSH_SUFF=_1Y_grid_T.nc 
     76##VAR2_Q=qt            ;   V2Q_PREF=testht     ;   V2Q_SUFF=_1Y_grid_T.nc 
     77VAR2_EMP=wfo            ;   V2EMP_PREF=testht   ;   V2EMP_SUFF=_1Y_grid_T.nc 
     78VAR2_U=uocetr_eff      ;   V2U_PREF=testht     ;   V2U_SUFF=_1Y_grid_U.nc 
     79VAR2_V=vocetr_eff      ;   V2V_PREF=testht     ;   V2V_SUFF=_1Y_grid_V.nc 
     80VAR2_ICE=siconc         ;   V2ICE_PREF=testht   ;   V2ICE_SUFF=_1M_icemod.nc 
     81VAR2_Ithick=sithic       ;   V2It_PREF=testht    ;   V2It_SUFF=_1M_icemod.nc 
     82VAR2_SNOW=snthic       ;   V2SNOW_PREF=testht  ;   V2SNOW_SUFF=_1M_icemod.nc 
     83VAR2_IvelV=sivelv       ;   V2IvV_PREF=testht   ;   V2IvV_SUFF=_1M_icemod.nc  
     84VAR2_Ivel=sivelo        ;   V2IvV_PREF=testht   ;   V2IvV_SUFF=_1M_icemod.nc  
    8985# 
    9086######################### Export Variables ############################### 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM/cpp_ORCA2_LIM.fcm

    r4523 r5034  
    1  bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi key_diaobs key_asminc 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/1_namelist_cfg

    r4161 r5034  
    5555/ 
    5656!----------------------------------------------------------------------- 
     57&namsplit      !   time splitting parameters                            ("key_dynspg_ts") 
     58!----------------------------------------------------------------------- 
     59/ 
     60!----------------------------------------------------------------------- 
     61&namcrs        !   Grid coarsening for dynamics output and/or 
     62               !   passive tracer coarsened online simulations 
     63!----------------------------------------------------------------------- 
     64/ 
     65!----------------------------------------------------------------------- 
    5766&namtsd    !   data : Temperature  & Salinity 
    5867!----------------------------------------------------------------------- 
     
    7281/ 
    7382!----------------------------------------------------------------------- 
    74 &namsbc_core   !   namsbc_core  CORE bulk formulea 
    75 !----------------------------------------------------------------------- 
    76 !              !  file name  ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! 
    77 !              !             !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! 
    78    sn_wndi = 'u_10.15JUNE2009_orca2'       ,  6  , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_bicub.nc'       , '' 
    79    sn_wndj = 'v_10.15JUNE2009_orca2'       ,  6  , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_bicub.nc'       , '' 
    80    sn_qsr  = 'ncar_rad.15JUNE2009_orca2'   , 24  , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_bilin.nc'       , '' 
    81    sn_qlw  = 'ncar_rad.15JUNE2009_orca2'   , 24  , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_bilin.nc'       , '' 
    82    sn_tair = 't_10.15JUNE2009_orca2'       ,  6  , 'T_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_bilin.nc'       , '' 
    83    sn_humi = 'q_10.15JUNE2009_orca2'       ,  6  , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_bilin.nc'       , '' 
    84    sn_prec = 'ncar_precip.15JUNE2009_orca2', -1  , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_bilin.nc'       , '' 
    85    sn_snow = 'ncar_precip.15JUNE2009_orca2', -1  , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_bilin.nc'       , '' 
    86    sn_tdif = 'taudif_core'                 , 24  , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_bilin.nc'       , '' 
     83&namsbc_core   !   namsbc_core  CORE bulk formulae 
     84!----------------------------------------------------------------------- 
     85!              !  file name                    ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights                               ! rotation ! land/sea mask ! 
     86!              !                               !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename                              ! pairing  ! filename      ! 
     87   sn_wndi     = 'u_10.15JUNE2009_fill'        ,         6         , 'U_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bicubic.nc'      , 'Uwnd'   , '' 
     88   sn_wndj     = 'v_10.15JUNE2009_fill'        ,         6         , 'V_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bicubic.nc'      , 'Vwnd'   , '' 
     89   sn_qsr      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'SWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     90   sn_qlw      = 'ncar_rad.15JUNE2009_fill'    ,        24         , 'LWDN_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     91   sn_tair     = 't_10.15JUNE2009_fill'        ,         6         , 'T_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     92   sn_humi     = 'q_10.15JUNE2009_fill'        ,         6         , 'Q_10_MOD',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     93   sn_prec     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'PRC_MOD1',   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     94   sn_snow     = 'ncar_precip.15JUNE2009_fill' ,        -1         , 'SNOW'    ,   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     95   sn_tdif     = 'taudif_core'                 ,        24         , 'taudif'  ,   .false.    , .true. , 'yearly'  , 'weights_core2_agrif_bilinear.nc'     , ''       , '' 
     96   cn_dir      = './'      !  root directory for the location of the bulk files 
     97   ln_2m       = .false.   !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
     98   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
     99   ln_bulk2z   = .false.   !  Air temperature/humidity and wind vectors are referenced at heights rn_zqt and rn_zu 
     100   rn_zqt      = 3.        !  Air temperature and humidity reference height (m) (ln_bulk2z) 
     101   rn_zu       = 4.        !  Wind vector reference height (m)                  (ln_bulk2z) 
     102   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
     103   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
     104   rn_vfac     = 0.        !  multiplicative factor for ocean/ice velocity  
     105                           !  in the calculation of the wind stress (0.=absolute winds or 1.=relative winds) 
    87106/ 
    88107!----------------------------------------------------------------------- 
     
    120139&nameos        !   ocean physical parameters 
    121140!----------------------------------------------------------------------- 
     141/ 
     142!----------------------------------------------------------------------- 
     143&namtra_dmp    !   tracer: T & S newtonian damping 
     144!----------------------------------------------------------------------- 
     145   ln_tradmp   =  .false.   !  add a damping termn (T) or not (F) 
    122146/ 
    123147!----------------------------------------------------------------------- 
     
    197221!----------------------------------------------------------------------- 
    198222/ 
     223!----------------------------------------------------------------------- 
     224&namobs       !  observation usage                                      ('key_diaobs') 
     225!----------------------------------------------------------------------- 
     226/ 
     227!----------------------------------------------------------------------- 
     228&nam_asminc   !   assimilation increments                               ('key_asminc') 
     229!----------------------------------------------------------------------- 
     230/ 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/iodef.xml

    r4100 r5034  
    4747     <field field_ref="mldkz5"       /> 
    4848     <field field_ref="mldr10_1"     /> 
     49      <field field_ref="heatc"        name="heatc"    long_name="Heat content vertically integrated"            /> 
     50      <field field_ref="saltc"        name="saltc"    long_name="Salt content vertically integrated"            /> 
    4951   </file> 
    5052 
    5153   <file id="file2" name_suffix="_SBC" description="surface fluxes variables" > <!-- time step automaticaly defined based on nn_fsbc --> 
    5254     <field field_ref="empmr"        name="wfo"      long_name="water_flux_into_sea_water"                     /> 
    53      <field field_ref="qsr"          name="rsntds"   long_name="surface_net_downward_shortwave_flux"           /> 
    54      <field field_ref="qt"           name="tohfls"   long_name="surface_net_downward_total_heat_flux"          /> 
    55      <field field_ref="saltflx"      name="sosflxdo"  /> 
     55     <field field_ref="qsr_oce"      name="qsr_oce"  long_name="downward shortwave flux at ocean surface"           /> 
     56     <field field_ref="qns_oce"      name="qns_oce"  long_name="downward non solar flux at ocean surface"           /> 
     57     <field field_ref="qt_oce"       name="qt_oce"   long_name="downward total flux at ocean surface"           /> 
     58     <field field_ref="qsr_ice"      name="qsr_ice"  long_name="downward shortwave flux at ice surface"           /> 
     59     <field field_ref="qns_ice"      name="qns_ice"  long_name="downward non solar flux at ice surface"           /> 
     60     <field field_ref="qtr_ice"      name="qtr_ice"  long_name="shortwave flux transmitted thru the ice"           /> 
     61     <field field_ref="qt_ice"       name="qt_ice"   long_name="downward total flux at ice surface"           /> 
     62     <field field_ref="saltflx"      name="sfx"  /> 
    5663     <field field_ref="taum"         name="taum" /> 
    57      <field field_ref="wspd"         name="sowindsp"  /> 
    58           <field field_ref="precip"       name="soprecip" /> 
     64     <field field_ref="wspd"         name="windsp"  /> 
     65      <field field_ref="precip"       name="precip" /> 
    5966     <!-- ice and snow --> 
    6067     <field field_ref="snowpre"      /> 
    61           <field field_ref="qsr_io"           name="iicesflx" /> 
    62           <field field_ref="qns_io"           name="iicenflx" /> 
    63           <field field_ref="utau_ice"         name="iicestru" /> 
    64           <field field_ref="vtau_ice"         name="iicestrv" /> 
     68      <field field_ref="utau_ice"         name="utau_ice" /> 
     69      <field field_ref="vtau_ice"         name="vtau_ice" /> 
    6570 
    6671   </file> 
    6772 
    6873   <file id="file3" name_suffix="_grid_U" description="ocean U grid variables" > 
    69      <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
     74     <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    7075     <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
    7176     <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
     
    7378     <field field_ref="u_masstr"     name="vozomatr"  /> 
    7479     <field field_ref="u_heattr"     name="sozohetr"  /> 
     80      <field field_ref="u_salttr"     name="sozosatr"  /> 
    7581   </file> 
    7682    
    7783   <file id="file4" name_suffix="_grid_V" description="ocean V grid variables" > 
    78      <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
     84     <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    7985     <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
    8086     <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
     
    8288     <field field_ref="v_masstr"     name="vomematr"  /> 
    8389     <field field_ref="v_heattr"     name="somehetr"  /> 
     90      <field field_ref="v_salttr"     name="somesatr"  /> 
    8491   </file> 
    8592    
     
    9097   </file> 
    9198 
    92    <file id="file6" name_suffix="_icemod" description="ice variables" > 
    93      <field field_ref="snowthic_cea" name="snd"     long_name="surface_snow_thickness"   /> 
    94      <field field_ref="icethic_cea"  name="sit"     long_name="sea_ice_thickness"        /> 
    95           <field field_ref="icevolu"      name="iicevolu" /> 
    96           <field field_ref="snowvol"      name="isnowvol" /> 
    97           <field field_ref="iceconc"      name="iiceconc"  /> 
    98  
    99           <field field_ref="icebopr"          name="iicebopr" /> 
    100           <field field_ref="icedypr"          name="iicedypr" /> 
    101           <field field_ref="ioceflxb"         name="ioceflxb" /> 
    102           <field field_ref="uice_ipa"         name="iicevelu" /> 
    103           <field field_ref="vice_ipa"         name="iicevelv" /> 
    104           <field field_ref="isst"             name="isstempe" /> 
    105           <field field_ref="isss"             name="isssalin" /> 
    106           <field field_ref="micesalt"         name="iicesali" /> 
    107           <field field_ref="miceage"          name="iiceages" /> 
    108           <field field_ref="icelapr"          name="iicelapr" /> 
    109           <field field_ref="icesipr"          name="iicesipr" /> 
    110           <field field_ref="micet"            name="iicetemp" /> 
    111           <field field_ref="icehc"            name="iiceheco" /> 
    112           <field field_ref="isnowhc"          name="isnoheco" /> 
    113           <field field_ref="icest"            name="iicesurt" /> 
    114           <field field_ref="sfxbri"            name="iicefsbr" /> 
    115           <field field_ref="sfxthd"            name="iicefseq" /> 
    116           <field field_ref="ibrinv"           name="ibrinvol" /> 
    117           <field field_ref="icecolf"          name="iicecolf" /> 
    118           <field field_ref="icestr"           name="iicestre" /> 
    119           <field field_ref="icevel"           name="iicevelo" /> 
    120           <field field_ref="isume"            name="iicesume" /> 
    121           <field field_ref="ibome"            name="iicebome" /> 
    122           <field field_ref="idive"            name="iicedive" /> 
    123           <field field_ref="ishear"           name="iiceshea" /> 
    124           <field field_ref="icerepr"          name="iicerepr" /> 
    125           <field field_ref="sfxmec"            name="iicefsrp" /> 
    126           <field field_ref="sfxres"            name="iicefsre" /> 
    127           <field field_ref="icetrp"           name="iicevtrp" /> 
    128    </file> 
    129  
    130  
    131    <file id="file7" name_suffix="_scalar" description="scalar variables" > 
     99   <file id="file6" name_suffix="_icemod" description="ice variables" enabled=".true." > 
     100     <field field_ref="snowthic_cea"     name="snthic"     long_name="surface_snow_thickness"   /> 
     101     <field field_ref="icethic_cea"      name="sithic"     long_name="sea_ice_thickness"        /> 
     102          <field field_ref="icevolu"          name="sivolu" /> 
     103          <field field_ref="snowvol"          name="snvolu" /> 
     104          <field field_ref="iceconc"          name="siconc"  /> 
     105 
     106          <field field_ref="vfxbog"          name="vfxbog" /> 
     107          <field field_ref="vfxdyn"          name="vfxdyn" /> 
     108          <field field_ref="vfxopw"          name="vfxopw" /> 
     109          <field field_ref="vfxsni"          name="vfxsni" /> 
     110          <field field_ref="vfxsum"          name="vfxsum" /> 
     111          <field field_ref="vfxbom"          name="vfxbom" /> 
     112          <field field_ref="vfxres"          name="vfxres" /> 
     113          <field field_ref="vfxice"          name="vfxice" /> 
     114          <field field_ref="vfxsnw"          name="vfxsnw" /> 
     115          <field field_ref="vfxsub"          name="vfxsub" /> 
     116          <field field_ref="vfxspr"          name="vfxspr" /> 
     117 
     118          <field field_ref="icetrp"          name="sivtrp" /> 
     119          <field field_ref="snwtrp"          name="snvtrp" /> 
     120          <field field_ref="deitrp"          name="deitrp" /> 
     121          <field field_ref="destrp"          name="destrp" /> 
     122 
     123          <field field_ref="sfxbri"           name="sfxbri" /> 
     124          <field field_ref="sfxdyn"           name="sfxdyn" /> 
     125          <field field_ref="sfxres"           name="sfxres" /> 
     126          <field field_ref="sfxbog"           name="sfxbog" /> 
     127          <field field_ref="sfxbom"           name="sfxbom" /> 
     128          <field field_ref="sfxsum"           name="sfxsum" /> 
     129          <field field_ref="sfxsni"           name="sfxsni" /> 
     130          <field field_ref="sfxopw"           name="sfxopw" /> 
     131          <field field_ref="sfx"              name="sfx" /> 
     132 
     133          <field field_ref="hfxsum"          name="hfxsum"    /> 
     134          <field field_ref="hfxbom"          name="hfxbom"    /> 
     135          <field field_ref="hfxbog"          name="hfxbog"    /> 
     136          <field field_ref="hfxdif"          name="hfxdif"    /> 
     137          <field field_ref="hfxopw"          name="hfxopw"    /> 
     138          <field field_ref="hfxout"          name="hfxout"    /> 
     139          <field field_ref="hfxin"           name="hfxin"    /> 
     140          <field field_ref="hfxsnw"          name="hfxsnw"    /> 
     141          <field field_ref="hfxerr"          name="hfxerr"    /> 
     142          <field field_ref="hfxerr_rem"      name="hfxerr_rem"    /> 
     143 
     144     <!-- ice-ocean heat flux from mass exchange --> 
     145          <field field_ref="hfxdyn"          name="hfxdyn"    /> 
     146          <field field_ref="hfxres"          name="hfxres"    /> 
     147          <field field_ref="hfxthd"          name="hfxthd"    /> 
     148     <!-- ice-atm. heat flux from mass exchange --> 
     149          <field field_ref="hfxsub"          name="hfxsub"    /> 
     150          <field field_ref="hfxspr"          name="hfxspr"    /> 
     151 
     152     <!-- diags --> 
     153          <field field_ref="hfxdhc"          name="hfxdhc"    /> 
     154          <field field_ref="hfxtur"          name="hfxtur"    /> 
     155 
     156          <field field_ref="isst"             name="sst" /> 
     157          <field field_ref="isss"             name="sss" /> 
     158          <field field_ref="micesalt"         name="sisali" /> 
     159          <field field_ref="micet"            name="sitemp" /> 
     160          <field field_ref="icest"            name="sistem" /> 
     161          <field field_ref="icehc"            name="siheco" /> 
     162          <field field_ref="isnowhc"          name="snheco" /> 
     163          <field field_ref="miceage"          name="siages" /> 
     164 
     165          <field field_ref="uice_ipa"         name="sivelu" /> 
     166          <field field_ref="vice_ipa"         name="sivelv" /> 
     167          <field field_ref="icevel"           name="sivelo" /> 
     168          <field field_ref="idive"            name="sidive" /> 
     169          <field field_ref="ishear"           name="sishea" /> 
     170          <field field_ref="icestr"           name="sistre" /> 
     171 
     172          <field field_ref="ibrinv"           name="sibrin" /> 
     173          <field field_ref="icecolf"          name="sicolf" /> 
     174 
     175          <field field_ref="iceage_cat"       name="siagecat"/> 
     176          <field field_ref="iceconc_cat"      name="siconcat"/> 
     177          <field field_ref="icethic_cat"      name="sithicat"/> 
     178          <field field_ref="snowthic_cat"     name="snthicat"/> 
     179          <field field_ref="salinity_cat"     name="salincat"/> 
     180          <field field_ref="brinevol_cat"     name="sibricat"/> 
     181 
     182   </file> 
     183 
     184        <file id="file7" name_suffix="_scalar" description="scalar variables" enabled=".true." > 
    132185          <field field_ref="voltot"       name="scvoltot"  /> 
    133186          <field field_ref="sshtot"       name="scsshtot"  /> 
     
    137190          <field field_ref="temptot"      name="sctemtot"  /> 
    138191          <field field_ref="saltot"       name="scsaltot"  /> 
    139     
    140           <!-- available with ln_diahsb --> 
    141      <field field_ref="bgtemper"     name="bgtemper"    /> 
    142      <field field_ref="bgsaline"     name="bgsaline"    /> 
    143      <field field_ref="bgheatco"     name="bgheatco"    /> 
    144      <field field_ref="bgsaltco"     name="bgsaltco"    /> 
    145      <field field_ref="bgvolssh"     name="bgvolssh"    /> 
    146      <field field_ref="bgvoltot"     name="bgvoltot"    /> 
    147      <field field_ref="bgsshtot"     name="bgsshtot"    /> 
    148      <field field_ref="bgfrcvol"     name="bgfrcvol"    /> 
    149      <field field_ref="bgfrctem"     name="bgfrctem"    /> 
    150      <field field_ref="bgfrcsal"     name="bgfrcsal"    /> 
    151      <!-- available with ln_limdiahsb --> 
    152      <field field_ref="ibgvoltot"    name="ibgvoltot"   /> 
    153      <field field_ref="sbgvoltot"    name="sbgvoltot"   /> 
    154      <field field_ref="ibgarea"      name="ibgarea"     /> 
    155      <field field_ref="ibgsaline"    name="ibgsaline"   /> 
    156      <field field_ref="ibgtemper"    name="ibgtemper"   /> 
    157      <field field_ref="ibgheatco"    name="ibgheatco"   /> 
    158      <field field_ref="ibgsaltco"    name="ibgsaltco"   /> 
    159      <field field_ref="sbgheatco"    name="sbgheatco"   /> 
    160      <field field_ref="ibgfrcsfx"    name="ibgfrcsfx"  /> 
    161      <field field_ref="ibgfrcemp"    name="ibgfrcemp"   /> 
    162      <field field_ref="ibgsfx"       name="ibgsfx"     /> 
    163      <field field_ref="ibgemp"       name="ibgemp"      /> 
    164      <field field_ref="ibgsfxbri"    name="ibgsfxbri"    /> 
    165      <field field_ref="ibgsfxthd"    name="ibgsfxthd"    /> 
    166      <field field_ref="ibgsfxres"    name="ibgsfxres" /> 
    167      <field field_ref="ibgsfxmec"    name="ibgsfxmec" /> 
    168      <field field_ref="ibggrpme"     name="ibggrpme"    /> 
    169  
    170    </file> 
     192 
     193          <field field_ref="bgtemper"     name="bgtemper"    /> 
     194          <field field_ref="bgsaline"     name="bgsaline"    /> 
     195          <field field_ref="bgheatco"     name="bgheatco"    /> 
     196          <field field_ref="bgsaltco"     name="bgsaltco"    /> 
     197          <field field_ref="bgvolssh"     name="bgvolssh"    />  
     198          <field field_ref="bgvole3t"     name="bgvole3t"    /> 
     199          <field field_ref="bgfrcvol"     name="bgfrcvol"    /> 
     200          <field field_ref="bgfrctem"     name="bgfrctem"    /> 
     201          <field field_ref="bgfrcsal"     name="bgfrcsal"    /> 
     202 
     203          <field field_ref="ibgvoltot"    name="ibgvoltot"   /> 
     204          <field field_ref="sbgvoltot"    name="sbgvoltot"   /> 
     205          <field field_ref="ibgarea"      name="ibgarea"     /> 
     206          <field field_ref="ibgsaline"    name="ibgsaline"   /> 
     207          <field field_ref="ibgtemper"    name="ibgtemper"   /> 
     208          <field field_ref="ibgheatco"    name="ibgheatco"   /> 
     209          <field field_ref="sbgheatco"    name="sbgheatco"   /> 
     210          <field field_ref="ibgsaltco"    name="ibgsaltco"   /> 
     211 
     212          <field field_ref="ibgvfx"       name="ibgvfx"      /> 
     213          <field field_ref="ibgvfxbog"    name="ibgvfxbog"      /> 
     214          <field field_ref="ibgvfxopw"    name="ibgvfxopw"      /> 
     215          <field field_ref="ibgvfxsni"    name="ibgvfxsni"      /> 
     216          <field field_ref="ibgvfxdyn"    name="ibgvfxdyn"      /> 
     217          <field field_ref="ibgvfxbom"    name="ibgvfxbom"      /> 
     218          <field field_ref="ibgvfxsum"    name="ibgvfxsum"      /> 
     219          <field field_ref="ibgvfxres"    name="ibgvfxres"      /> 
     220          <field field_ref="ibgvfxspr"    name="ibgvfxspr"      /> 
     221          <field field_ref="ibgvfxsnw"    name="ibgvfxsnw"      /> 
     222          <field field_ref="ibgvfxsub"    name="ibgvfxsub"      /> 
     223 
     224          <field field_ref="ibgsfx"       name="ibgsfx"     /> 
     225          <field field_ref="ibgsfxbri"    name="ibgsfxbri"    /> 
     226          <field field_ref="ibgsfxdyn"    name="ibgsfxdyn" /> 
     227          <field field_ref="ibgsfxres"    name="ibgsfxres" /> 
     228          <field field_ref="ibgsfxbog"    name="ibgsfxbog"      /> 
     229          <field field_ref="ibgsfxopw"    name="ibgsfxopw"      /> 
     230          <field field_ref="ibgsfxsni"    name="ibgsfxsni"      /> 
     231          <field field_ref="ibgsfxbom"    name="ibgsfxbom"      /> 
     232          <field field_ref="ibgsfxsum"    name="ibgsfxsum"      /> 
     233 
     234          <field field_ref="ibghfxdhc"    name="ibghfxdhc"    /> 
     235          <field field_ref="ibghfxspr"    name="ibghfxspr"    /> 
     236 
     237          <field field_ref="ibghfxres"    name="ibghfxres"    /> 
     238          <field field_ref="ibghfxsub"    name="ibghfxsub"    /> 
     239          <field field_ref="ibghfxdyn"    name="ibghfxdyn"    /> 
     240          <field field_ref="ibghfxthd"    name="ibghfxthd"    /> 
     241          <field field_ref="ibghfxsum"    name="ibghfxsum"    /> 
     242          <field field_ref="ibghfxbom"    name="ibghfxbom"    /> 
     243          <field field_ref="ibghfxbog"    name="ibghfxbog"    /> 
     244          <field field_ref="ibghfxdif"    name="ibghfxdif"    /> 
     245          <field field_ref="ibghfxopw"    name="ibghfxopw"    /> 
     246          <field field_ref="ibghfxout"    name="ibghfxout"    /> 
     247          <field field_ref="ibghfxin"    name="ibghfxin"    /> 
     248          <field field_ref="ibghfxsnw"    name="ibghfxsnw"    /> 
     249 
     250          <field field_ref="ibgfrcvol"    name="ibgfrcvol"   /> 
     251          <field field_ref="ibgfrcsfx"    name="ibgfrcsfx"  /> 
     252          <field field_ref="ibgvolgrm"    name="ibgvolgrm"    /> 
     253 
     254        </file> 
    171255 
    172256   <!-- 
     
    214298      <axis id="nfloat" long_name="Float number"      unit="-"  /> 
    215299      <axis id="icbcla" long_name="Iceberg class"     unit="-"  /> 
     300      <axis id="ncatice" long_name="Ice categories"   unit="-"  /> 
    216301   </axis_definition>  
    217302     
     
    241326     <variable id="buffer_server_factor_size" type="integer">2</variable> 
    242327     <variable id="info_level"                type="integer">0</variable> 
    243      <variable id="using_server"              type="boolean">true</variable> 
     328     <variable id="using_server"              type="boolean">false</variable> 
    244329     <variable id="using_oasis"               type="boolean">false</variable> 
    245330     <variable id="oasis_codes_id"            type="string" >oceanx</variable> 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/iodef_default.xml

    r4100 r5034  
    4747     <field field_ref="mldkz5"       /> 
    4848     <field field_ref="mldr10_1"     /> 
     49      <field field_ref="heatc"        name="heatc"    long_name="Heat content vertically integrated"            /> 
     50      <field field_ref="saltc"        name="saltc"    long_name="Salt content vertically integrated"            /> 
    4951   </file> 
    5052 
    5153   <file id="file2" name_suffix="_SBC" description="surface fluxes variables" > <!-- time step automaticaly defined based on nn_fsbc --> 
    5254     <field field_ref="empmr"        name="wfo"      long_name="water_flux_into_sea_water"                     /> 
    53      <field field_ref="qsr"          name="rsntds"   long_name="surface_net_downward_shortwave_flux"           /> 
    54      <field field_ref="qt"           name="tohfls"   long_name="surface_net_downward_total_heat_flux"          /> 
    55      <field field_ref="saltflx"      name="sosflxdo"  /> 
     55     <field field_ref="qsr_oce"      name="qsr_oce"  long_name="downward shortwave flux at ocean surface"           /> 
     56     <field field_ref="qns_oce"      name="qns_oce"  long_name="downward non solar flux at ocean surface"           /> 
     57     <field field_ref="qt_oce"       name="qt_oce"   long_name="downward total flux at ocean surface"           /> 
     58     <field field_ref="qsr_ice"      name="qsr_ice"  long_name="downward shortwave flux at ice surface"           /> 
     59     <field field_ref="qns_ice"      name="qns_ice"  long_name="downward non solar flux at ice surface"           /> 
     60     <field field_ref="qtr_ice"      name="qtr_ice"  long_name="shortwave flux transmitted thru the ice"           /> 
     61     <field field_ref="qt_ice"       name="qt_ice"   long_name="downward total flux at ice surface"           /> 
     62     <field field_ref="saltflx"      name="sfx"  /> 
    5663     <field field_ref="taum"         name="taum" /> 
    57      <field field_ref="wspd"         name="sowindsp"  /> 
    58           <field field_ref="precip"       name="soprecip" /> 
     64     <field field_ref="wspd"         name="windsp"  /> 
     65      <field field_ref="precip"       name="precip" /> 
    5966     <!-- ice and snow --> 
    6067     <field field_ref="snowpre"      /> 
    61           <field field_ref="qsr_io"           name="iicesflx" /> 
    62           <field field_ref="qns_io"           name="iicenflx" /> 
    63           <field field_ref="utau_ice"         name="iicestru" /> 
    64           <field field_ref="vtau_ice"         name="iicestrv" /> 
     68      <field field_ref="utau_ice"         name="utau_ice" /> 
     69      <field field_ref="vtau_ice"         name="vtau_ice" /> 
    6570 
    6671   </file> 
     
    7378     <field field_ref="u_masstr"     name="vozomatr"  /> 
    7479     <field field_ref="u_heattr"     name="sozohetr"  /> 
     80      <field field_ref="u_salttr"     name="sozosatr"  /> 
    7581   </file> 
    7682    
     
    8288     <field field_ref="v_masstr"     name="vomematr"  /> 
    8389     <field field_ref="v_heattr"     name="somehetr"  /> 
     90      <field field_ref="v_salttr"     name="somesatr"  /> 
    8491   </file> 
    8592    
     
    9097   </file> 
    9198 
    92    <file id="file6" name_suffix="_icemod" description="ice variables" > 
    93      <field field_ref="snowthic_cea" name="snd"     long_name="surface_snow_thickness"   /> 
    94      <field field_ref="icethic_cea"  name="sit"     long_name="sea_ice_thickness"        /> 
    95           <field field_ref="icevolu"      name="iicevolu" /> 
    96           <field field_ref="snowvol"      name="isnowvol" /> 
    97           <field field_ref="iceconc"      name="iiceconc"  /> 
    98  
    99           <field field_ref="icebopr"          name="iicebopr" /> 
    100           <field field_ref="icedypr"          name="iicedypr" /> 
    101           <field field_ref="ioceflxb"         name="ioceflxb" /> 
    102           <field field_ref="uice_ipa"         name="iicevelu" /> 
    103           <field field_ref="vice_ipa"         name="iicevelv" /> 
    104           <field field_ref="isst"             name="isstempe" /> 
    105           <field field_ref="isss"             name="isssalin" /> 
    106           <field field_ref="micesalt"         name="iicesali" /> 
    107           <field field_ref="miceage"          name="iiceages" /> 
    108           <field field_ref="icelapr"          name="iicelapr" /> 
    109           <field field_ref="icesipr"          name="iicesipr" /> 
    110           <field field_ref="micet"            name="iicetemp" /> 
    111           <field field_ref="icehc"            name="iiceheco" /> 
    112           <field field_ref="isnowhc"          name="isnoheco" /> 
    113           <field field_ref="icest"            name="iicesurt" /> 
    114           <field field_ref="sfxbri"            name="iicefsbr" /> 
    115           <field field_ref="sfxthd"            name="iicefseq" /> 
    116           <field field_ref="ibrinv"           name="ibrinvol" /> 
    117           <field field_ref="icecolf"          name="iicecolf" /> 
    118           <field field_ref="icestr"           name="iicestre" /> 
    119           <field field_ref="icevel"           name="iicevelo" /> 
    120           <field field_ref="isume"            name="iicesume" /> 
    121           <field field_ref="ibome"            name="iicebome" /> 
    122           <field field_ref="idive"            name="iicedive" /> 
    123           <field field_ref="ishear"           name="iiceshea" /> 
    124           <field field_ref="icerepr"          name="iicerepr" /> 
    125           <field field_ref="sfxmec"            name="iicefsrp" /> 
    126           <field field_ref="sfxres"            name="iicefsre" /> 
    127           <field field_ref="icetrp"           name="iicevtrp" /> 
    128    </file> 
    129  
    130  
    131    <file id="file7" name_suffix="_scalar" description="scalar variables" > 
     99   <file id="file6" name_suffix="_icemod" description="ice variables" enabled=".true." > 
     100     <field field_ref="snowthic_cea"     name="snthic"     long_name="surface_snow_thickness"   /> 
     101     <field field_ref="icethic_cea"      name="sithic"     long_name="sea_ice_thickness"        /> 
     102          <field field_ref="icevolu"          name="sivolu" /> 
     103          <field field_ref="snowvol"          name="snvolu" /> 
     104          <field field_ref="iceconc"          name="siconc"  /> 
     105 
     106          <field field_ref="vfxbog"          name="vfxbog" /> 
     107          <field field_ref="vfxdyn"          name="vfxdyn" /> 
     108          <field field_ref="vfxopw"          name="vfxopw" /> 
     109          <field field_ref="vfxsni"          name="vfxsni" /> 
     110          <field field_ref="vfxsum"          name="vfxsum" /> 
     111          <field field_ref="vfxbom"          name="vfxbom" /> 
     112          <field field_ref="vfxres"          name="vfxres" /> 
     113          <field field_ref="vfxice"          name="vfxice" /> 
     114          <field field_ref="vfxsnw"          name="vfxsnw" /> 
     115          <field field_ref="vfxsub"          name="vfxsub" /> 
     116          <field field_ref="vfxspr"          name="vfxspr" /> 
     117 
     118          <field field_ref="icetrp"          name="sivtrp" /> 
     119          <field field_ref="snwtrp"          name="snvtrp" /> 
     120          <field field_ref="deitrp"          name="deitrp" /> 
     121          <field field_ref="destrp"          name="destrp" /> 
     122 
     123          <field field_ref="sfxbri"           name="sfxbri" /> 
     124          <field field_ref="sfxdyn"           name="sfxdyn" /> 
     125          <field field_ref="sfxres"           name="sfxres" /> 
     126          <field field_ref="sfxbog"           name="sfxbog" /> 
     127          <field field_ref="sfxbom"           name="sfxbom" /> 
     128          <field field_ref="sfxsum"           name="sfxsum" /> 
     129          <field field_ref="sfxsni"           name="sfxsni" /> 
     130          <field field_ref="sfxopw"           name="sfxopw" /> 
     131          <field field_ref="sfx"              name="sfx" /> 
     132 
     133          <field field_ref="hfxsum"          name="hfxsum"    /> 
     134          <field field_ref="hfxbom"          name="hfxbom"    /> 
     135          <field field_ref="hfxbog"          name="hfxbog"    /> 
     136          <field field_ref="hfxdif"          name="hfxdif"    /> 
     137          <field field_ref="hfxopw"          name="hfxopw"    /> 
     138          <field field_ref="hfxout"          name="hfxout"    /> 
     139          <field field_ref="hfxin"           name="hfxin"    /> 
     140          <field field_ref="hfxsnw"          name="hfxsnw"    /> 
     141          <field field_ref="hfxerr"          name="hfxerr"    /> 
     142          <field field_ref="hfxerr_rem"      name="hfxerr_rem"    /> 
     143 
     144     <!-- ice-ocean heat flux from mass exchange --> 
     145          <field field_ref="hfxdyn"          name="hfxdyn"    /> 
     146          <field field_ref="hfxres"          name="hfxres"    /> 
     147          <field field_ref="hfxthd"          name="hfxthd"    /> 
     148     <!-- ice-atm. heat flux from mass exchange --> 
     149          <field field_ref="hfxsub"          name="hfxsub"    /> 
     150          <field field_ref="hfxspr"          name="hfxspr"    /> 
     151 
     152     <!-- diags --> 
     153          <field field_ref="hfxdhc"          name="hfxdhc"    /> 
     154          <field field_ref="hfxtur"          name="hfxtur"    /> 
     155 
     156          <field field_ref="isst"             name="sst" /> 
     157          <field field_ref="isss"             name="sss" /> 
     158          <field field_ref="micesalt"         name="sisali" /> 
     159          <field field_ref="micet"            name="sitemp" /> 
     160          <field field_ref="icest"            name="sistem" /> 
     161          <field field_ref="icehc"            name="siheco" /> 
     162          <field field_ref="isnowhc"          name="snheco" /> 
     163          <field field_ref="miceage"          name="siages" /> 
     164 
     165          <field field_ref="uice_ipa"         name="sivelu" /> 
     166          <field field_ref="vice_ipa"         name="sivelv" /> 
     167          <field field_ref="icevel"           name="sivelo" /> 
     168          <field field_ref="idive"            name="sidive" /> 
     169          <field field_ref="ishear"           name="sishea" /> 
     170          <field field_ref="icestr"           name="sistre" /> 
     171 
     172          <field field_ref="ibrinv"           name="sibrin" /> 
     173          <field field_ref="icecolf"          name="sicolf" /> 
     174 
     175          <field field_ref="iceage_cat"       name="siagecat"/> 
     176          <field field_ref="iceconc_cat"      name="siconcat"/> 
     177          <field field_ref="icethic_cat"      name="sithicat"/> 
     178          <field field_ref="snowthic_cat"     name="snthicat"/> 
     179          <field field_ref="salinity_cat"     name="salincat"/> 
     180          <field field_ref="brinevol_cat"     name="sibricat"/> 
     181 
     182   </file> 
     183 
     184        <file id="file7" name_suffix="_scalar" description="scalar variables" enabled=".true." > 
    132185          <field field_ref="voltot"       name="scvoltot"  /> 
    133186          <field field_ref="sshtot"       name="scsshtot"  /> 
     
    137190          <field field_ref="temptot"      name="sctemtot"  /> 
    138191          <field field_ref="saltot"       name="scsaltot"  /> 
    139     
    140           <!-- available with ln_diahsb --> 
    141      <field field_ref="bgtemper"     name="bgtemper"    /> 
    142      <field field_ref="bgsaline"     name="bgsaline"    /> 
    143      <field field_ref="bgheatco"     name="bgheatco"    /> 
    144      <field field_ref="bgsaltco"     name="bgsaltco"    /> 
    145      <field field_ref="bgvolssh"     name="bgvolssh"    /> 
    146      <field field_ref="bgvoltot"     name="bgvoltot"    /> 
    147      <field field_ref="bgsshtot"     name="bgsshtot"    /> 
    148      <field field_ref="bgfrcvol"     name="bgfrcvol"    /> 
    149      <field field_ref="bgfrctem"     name="bgfrctem"    /> 
    150      <field field_ref="bgfrcsal"     name="bgfrcsal"    /> 
    151      <!-- available with ln_limdiahsb --> 
    152      <field field_ref="ibgvoltot"    name="ibgvoltot"   /> 
    153      <field field_ref="sbgvoltot"    name="sbgvoltot"   /> 
    154      <field field_ref="ibgarea"      name="ibgarea"     /> 
    155      <field field_ref="ibgsaline"    name="ibgsaline"   /> 
    156      <field field_ref="ibgtemper"    name="ibgtemper"   /> 
    157      <field field_ref="ibgheatco"    name="ibgheatco"   /> 
    158      <field field_ref="ibgsaltco"    name="ibgsaltco"   /> 
    159      <field field_ref="sbgheatco"    name="sbgheatco"   /> 
    160      <field field_ref="ibgfrcsfx"    name="ibgfrcsfx"  /> 
    161      <field field_ref="ibgfrcemp"    name="ibgfrcemp"   /> 
    162      <field field_ref="ibgsfx"       name="ibgsfx"     /> 
    163      <field field_ref="ibgemp"       name="ibgemp"      /> 
    164      <field field_ref="ibgsfxbri"    name="ibgsfxbri"    /> 
    165      <field field_ref="ibgsfxthd"    name="ibgsfxthd"    /> 
    166      <field field_ref="ibgsfxres"    name="ibgsfxres" /> 
    167      <field field_ref="ibgsfxmec"    name="ibgsfxmec" /> 
    168      <field field_ref="ibggrpme"     name="ibggrpme"    /> 
    169  
    170    </file> 
     192 
     193          <field field_ref="bgtemper"     name="bgtemper"    /> 
     194          <field field_ref="bgsaline"     name="bgsaline"    /> 
     195          <field field_ref="bgheatco"     name="bgheatco"    /> 
     196          <field field_ref="bgsaltco"     name="bgsaltco"    /> 
     197          <field field_ref="bgvolssh"     name="bgvolssh"    />  
     198          <field field_ref="bgvole3t"     name="bgvole3t"    /> 
     199          <field field_ref="bgfrcvol"     name="bgfrcvol"    /> 
     200          <field field_ref="bgfrctem"     name="bgfrctem"    /> 
     201          <field field_ref="bgfrcsal"     name="bgfrcsal"    /> 
     202 
     203          <field field_ref="ibgvoltot"    name="ibgvoltot"   /> 
     204          <field field_ref="sbgvoltot"    name="sbgvoltot"   /> 
     205          <field field_ref="ibgarea"      name="ibgarea"     /> 
     206          <field field_ref="ibgsaline"    name="ibgsaline"   /> 
     207          <field field_ref="ibgtemper"    name="ibgtemper"   /> 
     208          <field field_ref="ibgheatco"    name="ibgheatco"   /> 
     209          <field field_ref="sbgheatco"    name="sbgheatco"   /> 
     210          <field field_ref="ibgsaltco"    name="ibgsaltco"   /> 
     211 
     212          <field field_ref="ibgvfx"       name="ibgvfx"      /> 
     213          <field field_ref="ibgvfxbog"    name="ibgvfxbog"      /> 
     214          <field field_ref="ibgvfxopw"    name="ibgvfxopw"      /> 
     215          <field field_ref="ibgvfxsni"    name="ibgvfxsni"      /> 
     216          <field field_ref="ibgvfxdyn"    name="ibgvfxdyn"      /> 
     217          <field field_ref="ibgvfxbom"    name="ibgvfxbom"      /> 
     218          <field field_ref="ibgvfxsum"    name="ibgvfxsum"      /> 
     219          <field field_ref="ibgvfxres"    name="ibgvfxres"      /> 
     220          <field field_ref="ibgvfxspr"    name="ibgvfxspr"      /> 
     221          <field field_ref="ibgvfxsnw"    name="ibgvfxsnw"      /> 
     222          <field field_ref="ibgvfxsub"    name="ibgvfxsub"      /> 
     223 
     224          <field field_ref="ibgsfx"       name="ibgsfx"     /> 
     225          <field field_ref="ibgsfxbri"    name="ibgsfxbri"    /> 
     226          <field field_ref="ibgsfxdyn"    name="ibgsfxdyn" /> 
     227          <field field_ref="ibgsfxres"    name="ibgsfxres" /> 
     228          <field field_ref="ibgsfxbog"    name="ibgsfxbog"      /> 
     229          <field field_ref="ibgsfxopw"    name="ibgsfxopw"      /> 
     230          <field field_ref="ibgsfxsni"    name="ibgsfxsni"      /> 
     231          <field field_ref="ibgsfxbom"    name="ibgsfxbom"      /> 
     232          <field field_ref="ibgsfxsum"    name="ibgsfxsum"      /> 
     233 
     234          <field field_ref="ibghfxdhc"    name="ibghfxdhc"    /> 
     235          <field field_ref="ibghfxspr"    name="ibghfxspr"    /> 
     236 
     237          <field field_ref="ibghfxres"    name="ibghfxres"    /> 
     238          <field field_ref="ibghfxsub"    name="ibghfxsub"    /> 
     239          <field field_ref="ibghfxdyn"    name="ibghfxdyn"    /> 
     240          <field field_ref="ibghfxthd"    name="ibghfxthd"    /> 
     241          <field field_ref="ibghfxsum"    name="ibghfxsum"    /> 
     242          <field field_ref="ibghfxbom"    name="ibghfxbom"    /> 
     243          <field field_ref="ibghfxbog"    name="ibghfxbog"    /> 
     244          <field field_ref="ibghfxdif"    name="ibghfxdif"    /> 
     245          <field field_ref="ibghfxopw"    name="ibghfxopw"    /> 
     246          <field field_ref="ibghfxout"    name="ibghfxout"    /> 
     247          <field field_ref="ibghfxin"    name="ibghfxin"    /> 
     248          <field field_ref="ibghfxsnw"    name="ibghfxsnw"    /> 
     249 
     250          <field field_ref="ibgfrcvol"    name="ibgfrcvol"   /> 
     251          <field field_ref="ibgfrcsfx"    name="ibgfrcsfx"  /> 
     252          <field field_ref="ibgvolgrm"    name="ibgvolgrm"    /> 
     253 
     254        </file> 
    171255 
    172256   <!-- 
     
    214298      <axis id="nfloat" long_name="Float number"      unit="-"  /> 
    215299      <axis id="icbcla" long_name="Iceberg class"     unit="-"  /> 
     300      <axis id="ncatice" long_name="Ice categories"   unit="-"  /> 
    216301   </axis_definition>  
    217302     
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_cfg

    r4370 r5034  
    139139&namdyn_hpg    !   Hydrostatic pressure gradient option 
    140140!----------------------------------------------------------------------- 
     141   ln_hpg_zco  = .false.   !  z-coordinate - full steps 
     142   ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation) 
     143   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation) 
     144   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
     145   ln_hpg_prj  = .true.    !  s-coordinate (Pressure Jacobian scheme) 
     146   ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T) 
     147                                 !           centered      time scheme  (F) 
    141148/ 
    142149!----------------------------------------------------------------------- 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_ice_cfg

    r4161 r5034  
    1616/ 
    1717!----------------------------------------------------------------------- 
    18 &namicetrp     !   ice transport 
     18&namicethd     !   ice thermodynamic 
    1919!----------------------------------------------------------------------- 
    2020/ 
    2121!----------------------------------------------------------------------- 
    22 &namicethd     !   ice thermodynamic 
     22&namicesal     !   ice salinity 
    2323!----------------------------------------------------------------------- 
    2424/ 
     25!----------------------------------------------------------------------- 
     26&namiceitdme   !   parameters for mechanical redistribution of ice  
     27!----------------------------------------------------------------------- 
     28/ 
     29!----------------------------------------------------------------------- 
     30&namicedia     !   ice diagnostics 
     31!----------------------------------------------------------------------- 
     32/ 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM3/cpp_ORCA2_LIM3.fcm

    r4161 r5034  
    1  bld::tool::fppkeys key_trabbl key_lim3 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi 
     1 bld::tool::fppkeys key_trabbl key_lim3 key_vvl key_dynspg_ts key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi key_diaobs key_asminc 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/iodef.xml

    r4245 r5034  
    4646     <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
    4747   </file> 
    48     
     48 
     49        <file id="file4" name_suffix="_bioscalar" description="pisces sms variables" > 
     50          <field field_ref="tdenit"   name="tdenit"   unit="TgN/yr" operation="instant" > tdenit * 14. * 86400. * 365. / 1e12 </field> 
     51          <field field_ref="tnfix"    name="tnfix"    unit="TgN/yr" operation="instant" > tnfix * 14. * 86400. * 365. / 1e12 </field> 
     52          <field field_ref="tcflx"    name="tcflx"    unit="PgC/yr" operation="instant" > tcflx * -1. * 12. * 86400. * 365. / 1e15 </field> 
     53          <field field_ref="tcflxcum" name="tcflxcum" unit="PgC"    operation="instant" > tcflxcum * -1. * 12. / 1e15 </field> 
     54          <field field_ref="tcexp"    name="tcexp"    unit="PgC/yr" operation="instant" > tcexp * 12. * 86400. * 365. / 1e15 </field> 
     55          <field field_ref="tintpp"   name="tintpp"   unit="PgC/yr" operation="instant" > tintpp * 12. * 86400. * 365. / 1e15 </field> 
     56          <field field_ref="pno3tot"  name="pno3tot"  unit="umolN"  > pno3tot * 16. / 122. * 1e6 </field> 
     57          <field field_ref="psiltot"  name="psiltot"  unit="umolC"  > psiltot * 1e6  </field> 
     58          <field field_ref="palktot"  name="palktot"  unit="umolC"  > palktot * 1e6  </field> 
     59        <field field_ref="pfertot"  name="pfertot"  unit="nmolFe" > pfertot * 1e9  </field> 
     60                                                                                                                             </file> 
    4961      </file_group> 
    5062 
     
    5365      <file_group id="5d" output_freq="5d"  output_level="10" enabled=".TRUE.">  <!-- 5d files -->    
    5466 
    55    <file id="file4" name_suffix="_grid_T" description="ocean T grid variables" > 
     67   <file id="file5" name_suffix="_grid_T" description="ocean T grid variables" > 
    5668     <field field_ref="toce"         name="thetao"   long_name="sea_water_potential_temperature"               /> 
    5769     <field field_ref="soce"         name="so"       long_name="sea_water_salinity"                            /> 
     
    7284   </file> 
    7385    
    74    <file id="file5" name_suffix="_grid_U" description="ocean U grid variables" > 
     86   <file id="file6" name_suffix="_grid_U" description="ocean U grid variables" > 
    7587     <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
    7688     <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
     
    8193   </file> 
    8294    
    83    <file id="file6" name_suffix="_grid_V" description="ocean V grid variables" > 
     95   <file id="file7" name_suffix="_grid_V" description="ocean V grid variables" > 
    8496     <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
    8597     <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
     
    90102   </file> 
    91103    
    92    <file id="file7" name_suffix="_grid_W" description="ocean W grid variables" > 
     104   <file id="file8" name_suffix="_grid_W" description="ocean W grid variables" > 
    93105     <field field_ref="woce"         name="wo"      long_name="ocean vertical velocity"         /> 
    94106     <field field_ref="avt"          name="difvho"  long_name="ocean_vertical_heat_diffusivity" /> 
    95107   </file> 
    96108    
    97    <file id="file8" name_suffix="_icemod" description="ice variables" > 
     109   <file id="file9" name_suffix="_icemod" description="ice variables" > 
    98110     <field field_ref="ice_pres"                     /> 
    99111     <field field_ref="snowthic_cea" name="snd"     long_name="surface_snow_thickness"   /> 
     
    101113     <field field_ref="iceprod_cea"  name="sip"     long_name="sea_ice_thickness"        /> 
    102114     <field field_ref="ist_ipa"      /> 
    103      <field field_ref="ioceflxb"     /> 
    104115     <field field_ref="uice_ipa"     /> 
    105116     <field field_ref="vice_ipa"     /> 
     
    115126      <file_group id="1m" output_freq="1mo" output_level="10" enabled=".TRUE."> <!-- real monthly files --> 
    116127 
    117    <file id="file9" name_suffix="_ptrc_T" description="pisces sms variables" > 
     128   <file id="file10" name_suffix="_ptrc_T" description="pisces sms variables" > 
    118129          <field field_ref="DIC"      /> 
    119130          <field field_ref="Alkalini" /> 
     
    127138   </file> 
    128139    
    129    <file id="file10" name_suffix="_diad_T" description="additional pisces diagnostics" > 
     140   <file id="file11" name_suffix="_diad_T" description="additional pisces diagnostics" > 
    130141          <field field_ref="Cflx"     /> 
    131142          <field field_ref="Dpco2"    /> 
     
    141152      <file_group id="1y"  output_freq="1y" output_level="10" enabled=".TRUE."> <!-- real yearly files --> 
    142153 
    143    <file id="file11" name_suffix="_ptrc_T" description="pisces sms variables" > 
     154   <file id="file12" name_suffix="_ptrc_T" description="pisces sms variables" > 
    144155          <field field_ref="DIC"      /> 
    145156          <field field_ref="Alkalini" /> 
     
    168179   </file> 
    169180 
    170    <file id="file12" name_suffix="_diad_T" description="additional pisces diagnostics" > 
     181   <file id="file13" name_suffix="_diad_T" description="additional pisces diagnostics" > 
    171182          <field field_ref="PH"       /> 
    172183          <field field_ref="CO3"      /> 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/iodef.xml

    r4153 r5034  
    3030      <file_group id="6h" output_freq="6h"  output_level="10" enabled=".TRUE."/> <!-- 6h files --> 
    3131      
    32       <file_group id="1d" output_freq="1d"  output_level="10" enabled=".TRUE."/> <!-- 1d files --> 
     32      <file_group id="1d" output_freq="1d"  output_level="10" enabled=".TRUE."> <!-- 1d files --> 
     33        <file id="file1" name_suffix="_bioscalar" description="pisces sms variables" > 
     34           <field field_ref="tdenit"   name="tdenit"   unit="TgN/yr" operation="instant" > tdenit * 14. * 86400. * 365. / 1e12 </field> 
     35           <field field_ref="tnfix"    name="tnfix"    unit="TgN/yr" operation="instant" > tnfix * 14. * 86400. * 365. / 1e12 </field> 
     36           <field field_ref="tcflx"    name="tcflx"    unit="PgC/yr" operation="instant" > tcflx * -1. * 12. * 86400. * 365. / 1e15 </field> 
     37           <field field_ref="tcflxcum" name="tcflxcum" unit="PgC"    operation="instant" > tcflxcum * -1. * 12. / 1e15 </field> 
     38           <field field_ref="tcexp"    name="tcexp"    unit="PgC/yr" operation="instant" > tcexp * 12. * 86400. * 365. / 1e15 </field> 
     39           <field field_ref="tintpp"   name="tintpp"   unit="PgC/yr" operation="instant" > tintpp * 12. * 86400. * 365. / 1e15 </field> 
     40           <field field_ref="pno3tot"  name="pno3tot"  unit="umolN"  > pno3tot * 16. / 122. * 1e6 </field> 
     41           <field field_ref="psiltot"  name="psiltot"  unit="umolC"  > psiltot * 1e6  </field> 
     42           <field field_ref="palktot"  name="palktot"  unit="umolC"  > palktot * 1e6  </field> 
     43           <field field_ref="pfertot"  name="pfertot"  unit="nmolFe" > pfertot * 1e9  </field> 
     44        </file> 
     45      </file_group> 
     46 
    3347      <file_group id="3d" output_freq="3d"  output_level="10" enabled=".TRUE."/> <!-- 3d files -->     
    3448      <file_group id="5d" output_freq="5d"  output_level="10" enabled=".TRUE."/>  <!-- 5d files -->    
     
    3650      <file_group id="1m" output_freq="1mo" output_level="10" enabled=".TRUE."> <!-- real monthly files --> 
    3751 
    38    <file id="file1" name_suffix="_ptrc_T" description="pisces sms variables" > 
     52   <file id="file2" name_suffix="_ptrc_T" description="pisces sms variables" > 
    3953          <field field_ref="DIC"      /> 
    4054          <field field_ref="Alkalini" /> 
     
    4862   </file> 
    4963    
    50    <file id="file2" name_suffix="_diad_T" description="additional pisces diagnostics" > 
     64   <file id="file3" name_suffix="_diad_T" description="additional pisces diagnostics" > 
    5165          <field field_ref="Cflx"     /> 
    5266          <field field_ref="Dpco2"    /> 
     
    6175      <file_group id="1y"  output_freq="1y" output_level="10" enabled=".TRUE."> <!-- real yearly files --> 
    6276 
    63    <file id="file3" name_suffix="_ptrc_T" description="pisces sms variables" > 
     77   <file id="file4" name_suffix="_ptrc_T" description="pisces sms variables" > 
    6478          <field field_ref="DIC"      /> 
    6579          <field field_ref="Alkalini" /> 
     
    88102   </file> 
    89103 
    90    <file id="file4" name_suffix="_diad_T" description="additional pisces diagnostics" > 
     104   <file id="file5" name_suffix="_diad_T" description="additional pisces diagnostics" > 
    91105          <field field_ref="PH"       /> 
    92106          <field field_ref="CO3"      /> 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/ORCA2_SAS_LIM/EXP00/iodef.xml

    r4186 r5034  
    8484     <field field_ref="iceprod_cea"  name="sip"     long_name="sea_ice_thickness"        /> 
    8585     <field field_ref="ist_ipa"      /> 
    86      <field field_ref="ioceflxb"     /> 
    8786     <field field_ref="uice_ipa"     /> 
    8887     <field field_ref="vice_ipa"     /> 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/README_configs_namcfg_namdom

    r4147 r5034  
    331331!----------------------------------------------------------------------- 
    332332   cp_cfg      =  "orca"               !  name of the configuration 
    333    jp_cfg      =    1442               !  resolution of the configuration 
    334    jpidta      =    1021               !  1st lateral dimension ( >= jpi ) 
    335    jpjdta      =     511               !  2nd    "         "    ( >= jpj ) 
     333   jp_cfg      =     025               !  resolution of the configuration 
     334   jpidta      =    1442               !  1st lateral dimension ( >= jpi ) 
     335   jpjdta      =    1021               !  2nd    "         "    ( >= jpj ) 
    336336   jpkdta      =      75               !  number of levels      ( >= jpk ) 
    337337   jpiglo      =    1442               !  1st dimension of global domain --> i =jpidta 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/domain_def.xml

    r3771 r5034  
    55    <!--   My zoom: example of hand defined zoom   --> 
    66         <domain id="myzoom" zoom_ibegin="10" zoom_jbegin="10" zoom_ni="5" zoom_nj="5" /> 
    7     <!--   Eq section   --> 
     7         <domain id="1point" zoom_ibegin="10" zoom_jbegin="10" zoom_ni="1" zoom_nj="1" /> 
     8         <!--   Eq section   --> 
    89         <domain id="EqT" zoom_ibegin="1" zoom_jbegin="0000" zoom_ni="0000" zoom_nj="1" /> 
    910         <!--   TAO   --> 
     
    170171      <domain_group id="scalarpoint"> 
    171172         <domain id="scalarpoint" long_name="scalar"/> 
     173         <!-- 
    172174         <domain id="1point" zoom_ibegin="1" zoom_jbegin="1" zoom_ni="1" zoom_nj="1" /> 
     175         --> 
    173176      </domain_group> 
    174177 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/field_def.xml

    r4565 r5034  
    2828         <field id="mldkz5"       long_name="mixing layer depth (Turbocline)"           unit="m"                        /> 
    2929         <field id="mldr10_1"     long_name="Mixed Layer Depth 0.01 ref.10m"            unit="m"                        /> 
    30          <field id="rhop"         long_name="potential density (sigma0)"                unit="kg/m3" grid_ref="grid_T_3D"/> 
    31     <field id="eken"         long_name="kinetic energy"                            unit="m2/s2" grid_ref="grid_T_3D"/> 
    32     <field id="hdiv"         long_name="horizontal divergence"                     unit="s-1"   grid_ref="grid_T_3D"/> 
     30         <field id="heatc"        long_name="Heat content vertically integrated"        unit="J/m2"                      /> 
     31         <field id="saltc"        long_name="Salt content vertically integrated"        unit="PSU*kg/m2"                /> 
     32         <!-- EOS --> 
     33         <field id="alpha"        long_name="thermal expansion"                         unit="1/degC" grid_ref="grid_T_3D"/> 
     34         <field id="beta"         long_name="haline contraction"                        unit="1/psu"  grid_ref="grid_T_3D"/> 
     35         <field id="bn2"          long_name="squared Brunt-Vaisala frequency"           unit="1/s"    grid_ref="grid_T_3D"/> 
     36         <field id="rhop"         long_name="potential density (sigma0)"                unit="kg/m3"  grid_ref="grid_T_3D"/> 
     37         <!-- Energy - horizontal divergence --> 
     38         <field id="eken"         long_name="kinetic energy"                            unit="m2/s2"  grid_ref="grid_T_3D"/> 
     39         <field id="hdiv"         long_name="horizontal divergence"                     unit="s-1"    grid_ref="grid_T_3D"/> 
    3340         <!-- variables available with MLE --> 
    3441         <field id="Lf_NHpf"      long_name="MLE: Lf = N H / f"                         unit="m"                        /> 
     
    156163    <!-- * variable relative to atmospheric pressure forcing : available with ln_apr_dyn --> 
    157164    <field id="ssh_ib"       long_name="Inverse barometer sea surface height"                         unit="m"        /> 
    158     
     165 
     166         <!-- * variable related to ice shelf forcing * --> 
     167         <field id="fwfisf"       long_name="Ice shelf melting"                                            unit="Kg/m2/s"  /> 
     168         <field id="qisf"         long_name="Ice Shelf Heat Flux"                                          unit="W/m2"     /> 
     169         <field id="isfgammat"    long_name="transfert coefficient for isf (temperature) "                 unit="m/s"      /> 
     170         <field id="isfgammas"    long_name="transfert coefficient for isf (salinity)    "                 unit="m/s"      /> 
     171        <field id="stbl"         long_name="salinity in the Losh tbl                    "                 unit="PSU"      /> 
     172        <field id="ttbl"         long_name="temperature in the Losh tbl                 "                 unit="C"        /> 
     173 
    159174         <!-- *_oce variables available with ln_blk_clio or ln_blk_core --> 
    160          <field id="qns_oce"      long_name="Non solar Downward Heat Flux over open ocean"                 unit="W/m2"     /> 
    161175         <field id="qlw_oce"      long_name="Longwave Downward Heat Flux over open ocean"                  unit="W/m2"     /> 
    162176         <field id="qsb_oce"      long_name="Sensible Downward Heat Flux over open ocean"                  unit="W/m2"     /> 
     
    194208         <field id="ice_cover"    long_name="Ice fraction"                                                 unit="1"        /> 
    195209 
    196          <field id="ioceflxb"     long_name="Oceanic flux at the ice base"                                 unit="W/m2"     /> 
    197210         <field id="qsr_ai_cea"   long_name="Air-Ice downward solar heat flux (cell average)"              unit="W/m2"     /> 
    198211         <field id="qns_ai_cea"   long_name="Air-Ice downward non-solar heat flux (cell average)"          unit="W/m2"     /> 
     
    205218         <field id="icethic_cea"  long_name="Ice thickness (cell average)"                                 unit="m"        /> 
    206219         <field id="iceprod_cea"  long_name="Ice production (cell average)"                                unit="m/s"      /> 
     220         <field id="iiceconc"     long_name="Ice concentration"                                            unit=""         /> 
    207221          
    208222         <field id="ice_pres"     long_name="Ice presence"                                                 unit="-"        /> 
     
    219233         <field id="emp_x_sst"    long_name="Concentration/Dilution term on SST"                           unit="kgC/m2/s" /> 
    220234         <field id="emp_x_sss"    long_name="Concentration/Dilution term on SSS"                         unit="kgPSU/m2/s" />         
    221           
    222           
     235        
    223236         <field id="iceconc"      long_name="ice concentration"                                            unit="%"        /> 
    224          <field id="icebopr"      long_name="daily bottom thermo ice prod."                                unit="km3/day"   /> 
    225          <field id="icedypr"      long_name="daily  dynamic ice prod."                                     unit="km3/day"   /> 
    226     <field id="ioceflxb"     long_name="Oceanic flux at the ice base"                                 unit="W/m2"     /> 
    227237         <field id="uice_ipa"     long_name="Ice velocity along i-axis at I-point (ice presence average)"  unit="m/s"      /> 
    228238         <field id="vice_ipa"     long_name="Ice velocity along j-axis at I-point (ice presence average)"  unit="m/s"      /> 
     
    232242         <field id="qsr_oce"      long_name="solar heat flux at ocean surface"                             unit="W/m2"     /> 
    233243         <field id="qns_oce"      long_name="non-solar heat flux at ocean surface"                         unit="W/m2"     /> 
    234          <field id="hfbri"        long_name="heat flux due to brine release"                               unit="W/m2"     /> 
     244         <field id="qt_ice"       long_name="total flux at ice surface"                                    unit="W/m2"     /> 
     245         <field id="qsr_ice"      long_name="solar heat flux at ice surface"                               unit="W/m2"     /> 
     246         <field id="qns_ice"      long_name="non-solar heat flux at ice surface"                           unit="W/m2"     /> 
     247         <field id="qtr_ice"      long_name="solar heat flux transmitted thru the ice"                     unit="W/m2"     /> 
    235248         <field id="utau_ice"     long_name="Wind stress along i-axis over the ice at i-point"             unit="N/m2"     /> 
    236249         <field id="vtau_ice"     long_name="Wind stress along j-axis over the ice at i-point"             unit="N/m2"     /> 
    237     <field id="qsr_io"       long_name="Ice-Oce downward solar heat flux"                             unit="W/m2"     /> 
    238     <field id="qns_io"       long_name="Ice-Oce downward non-solar heat flux"                         unit="W/m2"     /> 
    239250         <field id="micesalt"     long_name="Mean ice salinity"                                            unit="psu"      /> 
    240251         <field id="miceage"      long_name="Mean ice age"                                                 unit="years"    /> 
    241          <field id="icelapr"      long_name="daily lateral thermo ice prod."                               unit="km3/day"   /> 
    242          <field id="icesipr"      long_name="daily snowice ice prod."                                      unit="km3/day"   /> 
     252 
     253         <field id="iceage_cat"   long_name="Ice age for categories"                                       unit="days"   axis_ref="ncatice" /> 
     254         <field id="iceconc_cat"  long_name="Ice concentration for categories"                             unit="%"      axis_ref="ncatice" /> 
     255         <field id="icethic_cat"  long_name="Ice thickness for categories"                                 unit="m"      axis_ref="ncatice" /> 
     256         <field id="snowthic_cat" long_name="Snow thicknessi for categories"                               unit="m"      axis_ref="ncatice" /> 
     257         <field id="salinity_cat" long_name="Sea-Ice Bulk salinity for categories"                         unit="ppt"    axis_ref="ncatice" /> 
     258         <field id="brinevol_cat" long_name="Brine volume for categories"                                  unit="%"      axis_ref="ncatice" /> 
     259 
    243260         <field id="micet"        long_name="Mean ice temperature"                                         unit="degC"     /> 
    244          <field id="icehc"        long_name="ice total heat content"                                       unit="10^9 J"   />  
     261         <field id="icehc"        long_name="ice total heat content"                                       unit="10^9J"   />  
    245262         <field id="isnowhc"      long_name="snow total heat content"                                      unit="10^9J"    /> 
    246263         <field id="icest"        long_name="ice surface temperature"                                      unit="degC"     /> 
    247          <field id="sfxbri"       long_name="brine salt flux"                                              unit="psu*kg/m2/day" /> 
    248          <field id="sfxthd"       long_name="equivalent FW salt flux"                                      unit="psu*kg/m2/day" /> 
    249264         <field id="ibrinv"       long_name="brine volume"                                                 unit="%"        /> 
    250265         <field id="icecolf"      long_name="frazil ice collection thickness"                              unit="m"        /> 
    251266         <field id="icestr"       long_name="ice strength"                                                 unit="N/m"      /> 
    252267         <field id="icevel"       long_name="ice velocity"                                                 unit="m/s"      /> 
    253          <field id="isume"        long_name="surface melt"                                                 unit="km3/day"   /> 
    254          <field id="ibome"        long_name="bottom melt"                                                  unit="km3/day"   /> 
    255268         <field id="idive"        long_name="divergence"                                                   unit="10-8s-1"  /> 
    256269         <field id="ishear"       long_name="shear"                                                        unit="10-8s-1"  /> 
    257          <field id="icerepr"      long_name="daily resultant ice prod./melting from limupdate"             unit="km3/day"   /> 
    258          <field id="icevolu"      long_name="ice volume"                                                   unit="km3"      /> 
    259          <field id="snowvol"      long_name="snow volume"                                                  unit="km3"      /> 
    260          <field id="sfxmec"       long_name="salt flux from ridging rafting"                               unit="psu*kg/m2/day" /> 
     270         <field id="icevolu"      long_name="ice volume"                                                   unit="m"      /> 
     271         <field id="snowvol"      long_name="snow volume"                                                  unit="m"      /> 
     272 
     273         <field id="icetrp"       long_name="ice volume transport"                                         unit="m/day"   /> 
     274         <field id="snwtrp"       long_name="snw volume transport"                                         unit="m/day"   /> 
     275         <field id="deitrp"       long_name="advected ice enhalpy"                                         unit="W/m2"   /> 
     276         <field id="destrp"       long_name="advected snw enhalpy"                                         unit="W/m2"   /> 
     277 
     278         <field id="sfxbri"       long_name="brine salt flux"                                              unit="psu*kg/m2/day" /> 
     279         <field id="sfxdyn"       long_name="salt flux from ridging rafting"                               unit="psu*kg/m2/day" /> 
    261280         <field id="sfxres"       long_name="salt flux from lipupdate (resultant)"                         unit="psu*kg/m2/day" /> 
    262          <field id="icetrp"       long_name="ice volume transport"                                         unit="km3/day"   /> 
    263  
     281         <field id="sfxbog"       long_name="salt flux from bot growth"                                    unit="psu*kg/m2/day" /> 
     282         <field id="sfxbom"       long_name="salt flux from bot melt"                                      unit="psu*kg/m2/day" /> 
     283         <field id="sfxsum"       long_name="salt flux from surf melt"                                     unit="psu*kg/m2/day" /> 
     284         <field id="sfxsni"       long_name="salt flux from snow-ice formation"                            unit="psu*kg/m2/day" /> 
     285         <field id="sfxopw"       long_name="salt flux from open water ice formation"                      unit="psu*kg/m2/day" /> 
     286         <field id="sfx"          long_name="salt flux total"                                              unit="psu*kg/m2/day" /> 
     287 
     288         <field id="vfxbog"       long_name="daily bottom thermo ice prod."                                unit="m/day"   /> 
     289         <field id="vfxdyn"       long_name="daily  dynamic ice prod."                                     unit="m/day"   /> 
     290         <field id="vfxopw"       long_name="daily lateral thermo ice prod."                               unit="m/day"   /> 
     291         <field id="vfxsni"       long_name="daily snowice ice prod."                                      unit="m/day"   /> 
     292         <field id="vfxsum"       long_name="surface melt"                                                 unit="m/day"   /> 
     293         <field id="vfxbom"       long_name="bottom melt"                                                  unit="m/day"   /> 
     294         <field id="vfxres"       long_name="daily resultant ice prod./melting from limupdate"             unit="m/day"   /> 
     295         <field id="vfxice"       long_name="ice melt/growth"                                              unit="m/day"   /> 
     296         <field id="vfxsnw"       long_name="snw melt/growth"                                              unit="m/day"   /> 
     297         <field id="vfxsub"       long_name="snw sublimation"                                              unit="m/day"   /> 
     298         <field id="vfxspr"       long_name="snw precipitation on ice"                                     unit="m/day"   /> 
     299 
     300         <field id="hfxsum"   long_name="heat fluxes causing surface ice melt"            unit="W/m2"  /> 
     301         <field id="hfxbom"   long_name="heat fluxes causing bottom ice melt"             unit="W/m2"  /> 
     302         <field id="hfxbog"   long_name="heat fluxes causing bottom ice growth"           unit="W/m2"  /> 
     303         <field id="hfxdif"   long_name="heat fluxes causing ice temperature change"      unit="W/m2"  /> 
     304         <field id="hfxopw"   long_name="heat fluxes causing open water ice formation"    unit="W/m2"  /> 
     305         <field id="hfxsnw"   long_name="heat fluxes causing snow melt"                   unit="W/m2"  /> 
     306         <field id="hfxerr"   long_name="heat fluxes error after heat diffusion"          unit="W/m2"  /> 
     307         <field id="hfxerr_rem" long_name="heat fluxes error after remapping"             unit="W/m2"  /> 
     308         <field id="hfxout"   long_name="total heat fluxes received by the ocean"         unit="W/m2"  /> 
     309         <field id="hfxin"    long_name="total heat fluxes at the ice/ocean surface"      unit="W/m2"  /> 
     310 
     311    <!-- heat flux associated with mass exchange --> 
     312         <field id="hfxthd"   long_name="heat fluxes from ice-ocean mass exchange during thermo"              unit="W/m2"  /> 
     313         <field id="hfxdyn"   long_name="heat fluxes from ice-ocean mass exchange during dynamic"             unit="W/m2"  /> 
     314         <field id="hfxres"   long_name="heat fluxes from ice-ocean mass exchange during resultant"           unit="W/m2"  /> 
     315         <field id="hfxsub"   long_name="heat fluxes from ice-atm. mass exchange during sublimation"          unit="W/m2"  /> 
     316         <field id="hfxspr"   long_name="heat fluxes from ice-atm. mass exchange during snow precip"          unit="W/m2" /> 
     317 
     318    <!-- diags --> 
     319         <field id="hfxdhc"    long_name="Heat content variation in snow and ice"   unit="W/m2" /> 
     320         <field id="hfxtur"    long_name="turbulent heat flux at the ice base"      unit="W/m2"  /> 
    264321 
    265322      </field_group> 
     
    269326      <field_group id="grid_U"   grid_ref="grid_U_2D"> 
    270327         <field id="utau"         long_name="Wind Stress along i-axis"                    unit="N/m2" /> 
    271          <field id="suoce"        long_name="ocean surface current along i-axis"          unit="m/s"  /> 
    272328         <field id="uoce"         long_name="ocean current along i-axis"                  unit="m/s"  grid_ref="grid_U_3D" /> 
     329         <field id="ssu"          long_name="ocean surface current along i-axis"          unit="m/s"     /> 
    273330         <field id="uocetr_eff"   long_name="Effective ocean transport along i-axis"      unit="m3/s" grid_ref="grid_U_3D" /> 
    274331         <field id="uocet"        long_name="ocean transport along i-axis times temperature" unit="degC.m/s" grid_ref="grid_U_3D" /> 
     
    280337         <!-- uoce_eiv: available with key_trabbl --> 
    281338         <field id="uoce_bbl"     long_name="BBL ocean current along i-axis"              unit="m/s"  grid_ref="grid_U_3D" /> 
    282     <field id="ahu_bbl"      long_name="BBL diffusive flux along i-axis"             unit="m3/s" /> 
     339        <field id="ahu_bbl"      long_name="BBL diffusive flux along i-axis"             unit="m3/s" /> 
     340         <!-- variable for ice shelves --> 
     341        <field id="utbl"         long_name="zonal current in the Losh tbl"               unit="m/s"  axis_ref="none" /> 
    283342         <!-- variables available with key_diaar5 --> 
    284343         <field id="u_masstr"     long_name="ocean eulerian mass transport along i-axis"  unit="kg/s" grid_ref="grid_U_3D" /> 
    285344         <field id="u_heattr"     long_name="ocean eulerian heat transport along i-axis"  unit="W"    /> 
     345         <field id="u_salttr"     long_name="ocean eulerian salt transport along i-axis"  unit="PSU*kg/s"    /> 
    286346         <field id="ueiv_heattr"  long_name="ocean bolus heat transport along i-axis"     unit="W"    /> 
    287347         <field id="udiff_heattr" long_name="ocean diffusion heat transport along i-axis" unit="W"    /> 
     
    292352      <field_group id="grid_V"   grid_ref="grid_V_2D"> 
    293353         <field id="vtau"         long_name="Wind Stress along j-axis"                    unit="N/m2" /> 
    294          <field id="svoce"        long_name="ocean surface current along j-axis"          unit="m/s"  /> 
    295354         <field id="voce"         long_name="ocean current along j-axis"                  unit="m/s"  grid_ref="grid_V_3D" /> 
     355         <field id="ssv"          long_name="ocean surface current along j-axis"          unit="m/s"  /> 
    296356         <field id="vocetr_eff"   long_name="Effective ocean transport along j-axis"      unit="m3/s" grid_ref="grid_V_3D" /> 
    297357         <field id="vocet"        long_name="ocean transport along j-axis times temperature" unit="degC.m/s" grid_ref="grid_V_3D" /> 
     
    303363         <!-- voce_eiv: available with key_trabbl --> 
    304364         <field id="voce_bbl"     long_name="BBL ocean current along j-axis"              unit="m/s"  grid_ref="grid_V_3D" /> 
    305     <field id="ahv_bbl"      long_name="BBL diffusive flux along j-axis"             unit="m3/s"                   /> 
     365        <field id="ahv_bbl"      long_name="BBL diffusive flux along j-axis"             unit="m3/s"                   /> 
     366         <!-- variable for ice shelves --> 
     367        <field id="vtbl"         long_name="meridional current in the Losh tbl"        unit="m/s"  axis_ref="none" /> 
    306368         <!-- variables available with key_diaar5 --> 
    307369         <field id="v_masstr"     long_name="ocean eulerian mass transport along j-axis"  unit="kg/s" grid_ref="grid_V_3D" /> 
    308370         <field id="v_heattr"     long_name="ocean eulerian heat transport along j-axis"  unit="W"     /> 
     371         <field id="v_salttr"     long_name="ocean eulerian salt transport along i-axis"  unit="PSU*kg/s"    /> 
    309372         <field id="veiv_heattr"  long_name="ocean bolus heat transport along j-axis"     unit="W"     /> 
    310373         <field id="vdiff_heattr" long_name="ocean diffusion heat transport along j-axis" unit="W"     /> 
     
    350413 
    351414       <!-- available with ln_diahsb --> 
    352     <field id="bgtemper"     long_name="global mean temperature"                  unit="degC"   /> 
    353     <field id="bgsaline"     long_name="global mean salinity"                     unit="psu"    /> 
    354     <field id="bgheatco"     long_name="global mean heat content"                 unit="10^9J"  /> 
    355     <field id="bgsaltco"     long_name="global mean salt content"                 unit="psu*m3" /> 
    356     <field id="bgvolssh"     long_name="global mean ssh volume"                   unit="km3"     /> 
    357     <field id="bgvoltot"     long_name="global mean volume"                       unit="km3"     /> 
    358     <field id="bgsshtot"     long_name="global mean ssh"                          unit="m"      /> 
    359     <field id="bgfrcvol"     long_name="global mean volume from forcing"          unit="km3"     /> 
    360     <field id="bgfrctem"     long_name="global mean heat content from forcing"    unit="10^9J"  /> 
    361     <field id="bgfrcsal"     long_name="global mean salt content from forcing"    unit="psu*km3" /> 
    362     <field id="bgmistem"     long_name="global mean temperature error due to free surface" unit="degC" /> 
    363     <field id="bgmissal"     long_name="global mean salinity error due to free surface"    unit="psu"  /> 
     415       <field id="bgtemper"     long_name="global mean temperature"                  unit="degC"   /> 
     416       <field id="bgsaline"     long_name="global mean salinity"                     unit="psu"    /> 
     417       <field id="bgheatco"     long_name="global mean heat content"                 unit="10^9J"  /> 
     418       <field id="bgsaltco"     long_name="global mean salt content"                 unit="psu*m3" /> 
     419       <field id="bgvolssh"     long_name="global mean ssh volume"                   unit="km3"     /> 
     420        <field id="bgvole3t"     long_name="global mean volume variation (e3t)"           unit="km3"/> 
     421       <field id="bgvoltot"     long_name="global mean volume"                       unit="km3"     /> 
     422       <field id="bgsshtot"     long_name="global mean ssh"                          unit="m"      /> 
     423       <field id="bgfrcvol"     long_name="global mean volume from forcing"          unit="km3"     /> 
     424       <field id="bgfrctem"     long_name="global mean heat content from forcing"    unit="10^9J"  /> 
     425       <field id="bgfrcsal"     long_name="global mean salt content from forcing"    unit="psu*km3" /> 
     426       <field id="bgmistem"     long_name="global mean temperature error due to free surface" unit="degC" /> 
     427       <field id="bgmissal"     long_name="global mean salinity error due to free surface"    unit="psu"  /> 
    364428      </field_group> 
    365429 
    366430      <field_group id="SBC_scalar"  domain_ref="1point" > 
    367          <!-- available with ln_limdiahsb --> 
     431         <!-- available with ln_limdiaout --> 
    368432    <field id="ibgvoltot"    long_name="global mean ice volume"                   unit="km3"   /> 
    369433    <field id="sbgvoltot"    long_name="global mean snow volume"                  unit="km3"   /> 
     
    371435    <field id="ibgsaline"    long_name="global mean ice salinity"                 unit="psu"   /> 
    372436    <field id="ibgtemper"    long_name="global mean ice temperature"              unit="degC"   /> 
    373     <field id="ibgheatco"    long_name="global mean ice heat content"             unit="10^9J"   /> 
     437    <field id="ibgheatco"    long_name="global mean ice heat content"             unit="10^20J"   /> 
     438    <field id="sbgheatco"    long_name="global mean snow heat content"            unit="10^20J"   /> 
    374439    <field id="ibgsaltco"    long_name="global mean ice salt content"             unit="psu*km3"   /> 
    375     <field id="sbgheatco"    long_name="global mean snow heat content"            unit="10^9J"   /> 
    376     <field id="ibgfrcsfx"    long_name="global mean salt content from sfx"        unit="psu*km3"   /> 
    377     <field id="ibgfrcemp"    long_name="global mean volume from emp"              unit="km3"      /> 
    378     <field id="ibgsfx"       long_name="global mean emps"                         unit="psu*kg/m2/day"   /> 
    379     <field id="ibgemp"       long_name="global mean emp"                          unit="kg/m2/day"   /> 
    380     <field id="ibgsfxbri"    long_name="global mean ice sfx_bri"                  unit="psu*kg/m2/day"   /> 
    381     <field id="ibgsfxthd"    long_name="global mean ice sfx_thd"                  unit="psu*kg/m2/day"   /> 
    382     <field id="ibgsfxres"    long_name="global mean ice sfx_res"                  unit="psu*kg/m2/day"   /> 
    383     <field id="ibgsfxmec"    long_name="global mean ice fsalt_rpo"                unit="psu*kg/m2/day"   /> 
    384     <field id="ibggrpme"     long_name="global mean ice growth+melt volume"       unit="km3"      /> 
     440 
     441    <field id="ibgvfx"       long_name="global mean volume flux (emp)"               unit="m/day"   /> 
     442    <field id="ibgvfxbog"    long_name="global mean volume flux (bottom growth)"     unit="m/day"   /> 
     443    <field id="ibgvfxopw"    long_name="global mean volume flux (open water growth)" unit="m/day"   /> 
     444    <field id="ibgvfxsni"    long_name="global mean volume flux (snow-ice growth)"   unit="m/day"   /> 
     445    <field id="ibgvfxdyn"    long_name="global mean volume flux (dynamic growth)"    unit="m/day"   /> 
     446    <field id="ibgvfxbom"    long_name="global mean volume flux (bottom melt)"       unit="m/day"   /> 
     447    <field id="ibgvfxsum"    long_name="global mean volume flux (surface melt)"      unit="m/day"   /> 
     448    <field id="ibgvfxres"    long_name="global mean volume flux (resultant)"         unit="m/day"   /> 
     449    <field id="ibgvfxspr"    long_name="global mean volume flux (snow precip)"       unit="m/day"   /> 
     450    <field id="ibgvfxsnw"    long_name="global mean volume flux (snow melt)"         unit="m/day"   /> 
     451    <field id="ibgvfxsub"    long_name="global mean volume flux (snow sublimation)"  unit="m/day"   /> 
     452 
     453    <field id="ibgsfx"       long_name="global mean salt flux (total)"            unit="psu*m/day"   /> 
     454    <field id="ibgsfxbri"    long_name="global mean salt flux (brines)"           unit="psu*m/day"   /> 
     455    <field id="ibgsfxdyn"    long_name="global mean salt flux (dynamic)"          unit="psu*m/day"   /> 
     456    <field id="ibgsfxres"    long_name="global mean salt flux (resultant)"        unit="psu*m/day"   /> 
     457    <field id="ibgsfxbog"    long_name="global mean salt flux (thermo)"           unit="psu*m/day"   /> 
     458    <field id="ibgsfxopw"    long_name="global mean salt flux (thermo)"           unit="psu*m/day"   /> 
     459    <field id="ibgsfxsni"    long_name="global mean salt flux (thermo)"           unit="psu*m/day"   /> 
     460    <field id="ibgsfxbom"    long_name="global mean salt flux (thermo)"           unit="psu*m/day"   /> 
     461    <field id="ibgsfxsum"    long_name="global mean salt flux (thermo)"           unit="psu*m/day"   /> 
     462 
     463 
     464        <field id="ibghfxdhc"    long_name="Heat content variation in snow and ice"   unit="W" /> 
     465        <field id="ibghfxspr"    long_name="Heat content of snow precip"              unit="W" /> 
     466 
     467        <field id="ibghfxthd"   long_name="heat fluxes from ice-ocean exchange during thermo"              unit="W"  /> 
     468        <field id="ibghfxsum"   long_name="heat fluxes causing surface ice melt"              unit="W"  /> 
     469        <field id="ibghfxbom"   long_name="heat fluxes causing bottom ice melt"              unit="W"  /> 
     470        <field id="ibghfxbog"   long_name="heat fluxes causing bottom ice growth"              unit="W"  /> 
     471        <field id="ibghfxdif"   long_name="heat fluxes causing ice temperature change"              unit="W"  /> 
     472        <field id="ibghfxopw"   long_name="heat fluxes causing open water ice formation"              unit="W"  /> 
     473        <field id="ibghfxdyn"   long_name="heat fluxes from ice-ocean exchange during dynamic"             unit="W"  /> 
     474        <field id="ibghfxres"   long_name="heat fluxes from ice-ocean exchange during resultant"           unit="W"  /> 
     475        <field id="ibghfxsub"   long_name="heat fluxes from sublimation"                                   unit="W"  /> 
     476        <field id="ibghfxsnw"   long_name="heat fluxes from snow-ocean exchange"                           unit="W"  /> 
     477        <field id="ibghfxout"   long_name="non solar heat fluxes received by the ocean"                    unit="W"  /> 
     478        <field id="ibghfxin"    long_name="total heat fluxes at the ice surface"                           unit="W"  /> 
     479 
     480    <field id="ibgfrcvol"    long_name="global mean forcing volume (emp)"         unit="km3"      /> 
     481    <field id="ibgfrcsfx"    long_name="global mean forcing salt   (sfx)"         unit="psu*km3"   /> 
     482    <field id="ibgvolgrm"    long_name="global mean ice growth+melt volume"       unit="km3"      /> 
    385483      </field_group> 
    386484   
     
    501599       <field id="EPSI100"     long_name="Export of Silicate at 100 m"             unit="mol/m2/s"                     /> 
    502600       <field id="EPCAL100"    long_name="Export of Calcite at 100 m"              unit="mol/m2/s"                     /> 
     601       <field id="EXPC"        long_name="Export of carbon"                        unit="molC/m2/s"  grid_ref="grid_T_3D" /> 
     602       <field id="EXPFE"       long_name="Export of biogenic iron"                 unit="molFe/m2/s" grid_ref="grid_T_3D" /> 
     603       <field id="EXPSI"       long_name="Export of Silicate"                      unit="molSi/m2/s" grid_ref="grid_T_3D" /> 
     604       <field id="EXPCAL"      long_name="Export of Calcite"                       unit="molC/m2/s"  grid_ref="grid_T_3D" /> 
    503605       <field id="Cflx"        long_name="DIC flux"                                unit="mol/m2/s"                     /> 
    504606       <field id="Oflx"        long_name="Oxygen flux"                             unit="mol/m2/s"                     /> 
     
    510612       <field id="Ironsed"     long_name="Iron deposition from sediment"           unit="mol/m2/s"  grid_ref="grid_T_3D"/> 
    511613 
     614 
    512615       <!-- PISCES with Kriest parametisation : variables available with key_kriest --> 
    513        <field id="POCFlx"      long_name="Particulate organic C flux"              unit="mol/m2/s"   grid_ref="grid_T_3D" /> 
    514        <field id="NumFlx"      long_name="Particle number flux"                    unit="nbr/m2/s"   grid_ref="grid_T_3D" /> 
    515        <field id="SiFlx"       long_name="Biogenic Si flux"                        unit="mol/m2/s"   grid_ref="grid_T_3D" /> 
    516        <field id="CaCO3Flx"    long_name="CaCO3 flux"                              unit="mol/m2/s"   grid_ref="grid_T_3D" /> 
    517        <field id="xnum"        long_name="Number of particles in aggregats"        unit="-"          grid_ref="grid_T_3D" /> 
    518        <field id="W1"          long_name="sinking speed of mass flux"              unit="m2/s"       grid_ref="grid_T_3D" /> 
    519        <field id="W2"          long_name="sinking speed of number flux"            unit="m2/s"       grid_ref="grid_T_3D" /> 
     616       <field id="EPN100"      long_name="Particulate number flux at 100 m"        unit="mol/m2/s"                     /> 
     617       <field id="EXPN"        long_name="Particulate number flux"                 unit="mol/m2/s"   grid_ref="grid_T_3D" /> 
     618       <field id="XNUM"        long_name="Number of particles in aggregats"        unit="-"          grid_ref="grid_T_3D" /> 
     619       <field id="WSC"         long_name="sinking speed of mass flux"              unit="m2/s"       grid_ref="grid_T_3D" /> 
     620       <field id="WSN"         long_name="sinking speed of number flux"            unit="m2/s"       grid_ref="grid_T_3D" /> 
     621 
     622      <!-- dbio_T on T grid : variables available with key_diaar5 --> 
     623       <field id="TPP"         long_name="Total Primary production of phyto"                   unit="molC/m3/s"  grid_ref="grid_T_3D" /> 
     624       <field id="TPNEW"       long_name="New Primary production of phyto"                     unit="molC/m3/s"  grid_ref="grid_T_3D" /> 
     625       <field id="TPBFE"       long_name="Total biogenic iron production"                      unit="molFe/m3/s" grid_ref="grid_T_3D" /> 
     626       <field id="INTDIC"      long_name="DIC content"                                         unit="kg/m2"   /> 
     627       <field id="O2MIN"       long_name="Oxygen minimum concentration"                        unit="molC/L" /> 
     628       <field id="ZO2MIN"      long_name="Depth of oxygen minimum concentration"               unit="m"       /> 
     629       <field id="INTNFIX"     long_name="Nitrogen fixation rate : vert. integrated"           unit="molN/m2/s"  /> 
     630       <field id="INTPPPHY"    long_name="Vertically integrated primary production by nanophy" unit="molC/m2/s"  /> 
     631       <field id="INTPPPHY2"   long_name="Vertically integrated primary production by diatom"  unit="molC/m2/s"  /> 
     632       <field id="INTPP"       long_name="Vertically integrated primary production by phyto"   unit="molC/m2/s"  /> 
     633       <field id="INTPNEW"     long_name="Vertically integrated new primary production"        unit="molC/m2/s"  /> 
     634       <field id="INTPBFE"     long_name="Vertically integrated of biogenic iron production"   unit="molFe/m2/s" /> 
     635       <field id="INTPBSI"     long_name="Vertically integrated of biogenic Si production"     unit="molSi/m2/s" /> 
     636       <field id="INTPCAL"     long_name="Vertically integrated of calcite production"         unit="molC/m2/s" /> 
    520637 
    521638       <!-- PISCES light : variables available with key_pisces_reduced --> 
     
    549666       <field id="qintC14b"    long_name="Cumulative air-sea flux of Bomb C14"      unit="mol/m2"     /> 
    550667       <field id="fdecay"      long_name="Radiactive decay of Bomb C14"             unit="mol/m3"  grid_ref="grid_T_3D"  /> 
     668     </field_group> 
     669 
     670    <field_group id="PISCES_scalar"  domain_ref="1point" > 
     671       <field id="pno3tot"         long_name="global mean nitrate concentration"                  unit="mol/l"   /> 
     672       <field id="psiltot"         long_name="global mean silicate concentration"                 unit="mol/l"   /> 
     673       <field id="palktot"         long_name="global mean alkalinity concentration"               unit="mol/l"   /> 
     674       <field id="pfertot"         long_name="global mean iron concentration"                     unit="mol/l"   /> 
     675       <field id="tcflx"           long_name="total Flux of Carbon out of the ocean"              unit="molC/s" /> 
     676       <field id="tcflxcum"        long_name="cumulative total Flux of Carbon out of the ocean"   unit="molC/s" /> 
     677       <field id="tcexp"           long_name="total Carbon export at 100m"                        unit="molC/s" /> 
     678       <field id="tintpp"          long_name="global total integrated primary production"         unit="molC/s" /> 
     679       <field id="tnfix"           long_name="global total nitrogen fixation"                     unit="molN/s" /> 
     680       <field id="tdenit"          long_name="Total denitrification"                              unit="molN/s" /> 
    551681     </field_group> 
    552682 
     
    615745    <field_group id="groupU" > 
    616746      <field field_ref="uoce"         name="uo"      long_name="sea_water_x_velocity"      /> 
    617       <field field_ref="suoce"        name="uos"     long_name="sea_surface_x_velocity"    /> 
     747      <field field_ref="ssu"          name="uos"     long_name="sea_surface_x_velocity"    /> 
    618748      <field field_ref="utau"         name="tauuo"   long_name="surface_downward_x_stress" /> 
    619749    </field_group> 
     
    621751    <field_group id="groupV" > 
    622752      <field field_ref="voce"         name="vo"      long_name="sea_water_y_velocity"      /> 
    623       <field field_ref="svoce"        name="vos"     long_name="sea_surface_y_velocity"    /> 
     753      <field field_ref="ssv"          name="vos"     long_name="sea_surface_y_velocity"    /> 
    624754      <field field_ref="vtau"         name="tauvo"   long_name="surface_downward_y_stress" /> 
    625755    </field_group> 
     
    629759    </field_group> 
    630760     
     761    <!--  
     762============================================================================================================ 
     763                     Trend diagnostics : temperature, KE, PE, momentum 
     764============================================================================================================ 
     765    --> 
     766 
     767    <field_group id="trendT" grid_ref="grid_T_3D"> 
     768      <!-- variables available with ln_tra_trd --> 
     769      <field id="ttrd_xad"      long_name="temperature-trend: i-advection"                 unit="degC/s" /> 
     770      <field id="strd_xad"      long_name="salinity   -trend: i-advection"                 unit="psu/s"  /> 
     771      <field id="ttrd_yad"      long_name="temperature-trend: j-advection"                 unit="degC/s" /> 
     772      <field id="strd_yad"      long_name="salinity   -trend: j-advection"                 unit="psu/s"  /> 
     773      <field id="ttrd_zad"      long_name="temperature-trend: k-advection"                 unit="degC/s" /> 
     774      <field id="strd_zad"      long_name="salinity   -trend: k-advection"                 unit="psu/s"  /> 
     775      <field id="ttrd_sad"      long_name="temperature-trend: surface adv. (no-vvl)"       unit="degC/s" grid_ref="grid_T_2D" /> 
     776      <field id="strd_sad"      long_name="salinity   -trend: surface adv. (no-vvl)"       unit="psu/s"  grid_ref="grid_T_2D" /> 
     777      <field id="ttrd_ldf"      long_name="temperature-trend: lateral  diffusion"          unit="degC/s" /> 
     778      <field id="strd_ldf"      long_name="salinity   -trend: lateral  diffusion"          unit="psu/s"  /> 
     779      <field id="ttrd_zdf"      long_name="temperature-trend: vertical diffusion"          unit="degC/s" /> 
     780      <field id="strd_zdf"      long_name="salinity   -trend: vertical diffusion"          unit="psu/s"  /> 
     781      <!-- ln_traldf_iso=T only (iso-neutral diffusion) --> 
     782      <field id="ttrd_zdfp"     long_name="temperature-trend: pure vert. diffusion"        unit="degC/s" /> 
     783      <field id="strd_zdfp"     long_name="salinity   -trend: pure vert. diffusion"        unit="psu/s"  /> 
     784      <!-- --> 
     785      <field id="ttrd_dmp"      long_name="temperature-trend: interior restoring"          unit="degC/s" /> 
     786      <field id="strd_dmp"      long_name="salinity   -trend: interior restoring"          unit="psu/s"  /> 
     787      <field id="ttrd_bbl"      long_name="temperature-trend: bottom boundary layer"       unit="degC/s" /> 
     788      <field id="strd_bbl"      long_name="salinity   -trend: bottom boundary layer"       unit="psu/s"  /> 
     789      <field id="ttrd_npc"      long_name="temperature-trend: non-penetrative conv."       unit="degC/s" /> 
     790      <field id="strd_npc"      long_name="salinity   -trend: non-penetrative conv."       unit="psu/s"  /> 
     791      <field id="ttrd_qns"      long_name="temperature-trend: non-solar flux + runoff"     unit="degC/s" /> 
     792      <field id="strd_cdt"      long_name="salinity   -trend: C/D term       + runoff"     unit="degC/s" /> 
     793      <field id="ttrd_qsr"      long_name="temperature-trend: solar penetr. heating"       unit="degC/s" /> 
     794      <field id="ttrd_bbc"      long_name="temperature-trend: geothermal heating"          unit="degC/s" /> 
     795      <field id="ttrd_atf"      long_name="temperature-trend: asselin time filter"         unit="degC/s" /> 
     796      <field id="strd_atf"      long_name="salinity   -trend: asselin time filter"         unit="psu/s"  /> 
     797      <!-- variables available with ln_KE_trd --> 
     798      <field id="ketrd_hpg"     long_name="ke-trend: hydrostatic pressure gradient"        unit="W/s^3"  /> 
     799      <field id="ketrd_spg"     long_name="ke-trend: surface     pressure gradient"        unit="W/s^3"  /> 
     800      <field id="ketrd_spgexp"  long_name="ke-trend: surface pressure gradient (explicit)" unit="W/s^3"  /> 
     801      <field id="ketrd_spgflt"  long_name="ke-trend: surface pressure gradient (filter)"   unit="W/s^3"  /> 
     802      <field id="ssh_flt"       long_name="filtered contribution to ssh (dynspg_flt)"      unit="m"     grid_ref="grid_T_2D"   /> 
     803      <field id="w0"            long_name="surface vertical velocity"                      unit="m/s"   grid_ref="grid_T_2D"   /> 
     804      <field id="pw0_exp"       long_name="surface pressure flux due to ssh"               unit="W/s^2" grid_ref="grid_T_2D"   /> 
     805      <field id="pw0_flt"       long_name="surface pressure flux due to filtered ssh"      unit="W/s^2" grid_ref="grid_T_2D"   /> 
     806      <field id="ketrd_keg"     long_name="ke-trend: KE gradient         or hor. adv."     unit="W/s^3"  /> 
     807      <field id="ketrd_rvo"     long_name="ke-trend: relative  vorticity or metric term"   unit="W/s^3"  /> 
     808      <field id="ketrd_pvo"     long_name="ke-trend: planetary vorticity"                  unit="W/s^3"  /> 
     809      <field id="ketrd_zad"     long_name="ke-trend: vertical  advection"                  unit="W/s^3"  /> 
     810      <field id="ketrd_udx"     long_name="ke-trend: U.dx[U]"                              unit="W/s^3"  /> 
     811      <field id="ketrd_ldf"     long_name="ke-trend: lateral   diffusion"                  unit="W/s^3"  /> 
     812      <field id="ketrd_zdf"     long_name="ke-trend: vertical  diffusion"                  unit="W/s^3"  /> 
     813      <field id="ketrd_tau"     long_name="ke-trend: wind stress "                         unit="W/s^3" grid_ref="grid_T_2D"   /> 
     814      <field id="ketrd_bfr"     long_name="ke-trend: bottom friction (explicit)"           unit="W/s^3"  />    
     815      <field id="ketrd_bfri"    long_name="ke-trend: bottom friction (implicit)"           unit="W/s^3"  />    
     816      <field id="ketrd_atf"     long_name="ke-trend: asselin time filter trend"            unit="W/s^3"  />   
     817      <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)"    unit="W/s^3"  /> 
     818      <field id="KE"            long_name="kinetic energy: u(n)*u(n+1)/2"                  unit="W/s^2"  />    
     819      <!-- variables available with ln_PE_trd --> 
     820      <field id="petrd_xad"     long_name="pe-trend: i-advection"                          unit="W/m^3"  /> 
     821      <field id="petrd_yad"     long_name="pe-trend: j-advection"                          unit="W/m^3"  /> 
     822      <field id="petrd_zad"     long_name="pe-trend: k-advection"                          unit="W/m^3"  /> 
     823      <field id="petrd_sad"     long_name="pe-trend: surface adv. (no-vvl)"                unit="W/m^3"  grid_ref="grid_T_2D" /> 
     824      <field id="petrd_ldf"     long_name="pe-trend: lateral  diffusion"                   unit="W/m^3"  /> 
     825      <field id="petrd_zdf"     long_name="pe-trend: vertical diffusion"                   unit="W/m^3"  /> 
     826      <field id="petrd_zdfp"    long_name="pe-trend: pure vert. diffusion"                 unit="W/m^3"  /> 
     827      <field id="petrd_dmp"     long_name="pe-trend: interior restoring"                   unit="W/m^3"  /> 
     828      <field id="petrd_bbl"     long_name="pe-trend: bottom boundary layer"                unit="W/m^3"  /> 
     829      <field id="petrd_npc"     long_name="pe-trend: non-penetrative conv."                unit="W/m^3"  /> 
     830      <field id="petrd_nsr"     long_name="pe-trend: surface forcing + runoff"             unit="W/m^3"  /> 
     831      <field id="petrd_qsr"     long_name="pe-trend: solar penetr. heating"                unit="W/m^3"  /> 
     832      <field id="petrd_bbc"     long_name="pe-trend: geothermal heating"                   unit="W/m^3"  /> 
     833      <field id="petrd_atf"     long_name="pe-trend: asselin time filter"                  unit="W/m^3"  /> 
     834      <field id="PEanom"        long_name="potential energy anomaly"                       unit="SI"     />    
     835      <field id="alphaPE"       long_name="- partial deriv. of PEanom wrt T"               unit="/degC"  />    
     836      <field id="betaPE"        long_name="partial deriv. of PEanom wrt S"                 unit="/psu"   />    
     837    </field_group> 
     838 
     839    <field_group id="trendU" grid_ref="grid_U_3D"> 
     840     <!-- variables available with ln_dyn_trd --> 
     841     <field id="utrd_hpg"       long_name="i-trend: hydrostatic pressure gradient"         unit="m/s^2"                      /> 
     842     <field id="utrd_spg"       long_name="i-trend: surface     pressure gradient"         unit="m/s^2"                      /> 
     843     <field id="utrd_spgexp"    long_name="i-trend: surface pressure gradient (explicit)"  unit="m/s^2"                      /> 
     844     <field id="utrd_spgflt"    long_name="i-trend: surface pressure gradient (filtered)"  unit="m/s^2"                      /> 
     845     <field id="utrd_keg"       long_name="i-trend: KE gradient         or hor. adv."      unit="m/s^2"                      /> 
     846     <field id="utrd_rvo"       long_name="i-trend: relative  vorticity or metric term"    unit="m/s^2"                      /> 
     847     <field id="utrd_pvo"       long_name="i-trend: planetary vorticity"                   unit="m/s^2"                      /> 
     848     <field id="utrd_zad"       long_name="i-trend: vertical  advection"                   unit="m/s^2"                      /> 
     849     <field id="utrd_udx"       long_name="i-trend: U.dx[U]"                               unit="m/s^2"                      /> 
     850     <field id="utrd_ldf"       long_name="i-trend: lateral   diffusion"                   unit="m/s^2"                      /> 
     851     <field id="utrd_zdf"       long_name="i-trend: vertical  diffusion"                   unit="m/s^2"                      /> 
     852     <field id="utrd_tau"       long_name="i-trend: wind stress "                          unit="m/s^2" grid_ref="grid_U_2D" /> 
     853     <field id="utrd_bfr"       long_name="i-trend: bottom friction (explicit)"            unit="m/s^2"                      />    
     854     <field id="utrd_bfri"      long_name="i-trend: bottom friction (implicit)"            unit="m/s^2"                      />    
     855     <field id="utrd_tot"       long_name="i-trend: total momentum trend before atf"       unit="m/s^2"                      />    
     856     <field id="utrd_atf"       long_name="i-trend: asselin time filter trend"             unit="m/s^2"                      />    
     857    </field_group> 
     858 
     859    <field_group id="trendV" grid_ref="grid_V_3D"> 
     860     <!-- variables available with ln_dyn_trd --> 
     861     <field id="vtrd_hpg"       long_name="j-trend: hydrostatic pressure gradient"         unit="m/s^2"                      /> 
     862     <field id="vtrd_spg"       long_name="j-trend: surface     pressure gradient"         unit="m/s^2"                      /> 
     863     <field id="vtrd_spgexp"    long_name="j-trend: surface pressure gradient (explicit)"  unit="m/s^2"                      /> 
     864     <field id="vtrd_spgflt"    long_name="j-trend: surface pressure gradient (filtered)"  unit="m/s^2"                      /> 
     865     <field id="vtrd_keg"       long_name="j-trend: KE gradient         or hor. adv."      unit="m/s^2"                      /> 
     866     <field id="vtrd_rvo"       long_name="j-trend: relative  vorticity or metric term"    unit="m/s^2"                      /> 
     867     <field id="vtrd_pvo"       long_name="j-trend: planetary vorticity"                   unit="m/s^2"                      /> 
     868     <field id="vtrd_zad"       long_name="j-trend: vertical  advection"                   unit="m/s^2"                      /> 
     869     <field id="vtrd_vdy"       long_name="i-trend: V.dx[V]"                               unit="m/s^2"                      /> 
     870     <field id="vtrd_ldf"       long_name="j-trend: lateral   diffusion"                   unit="m/s^2"                      /> 
     871     <field id="vtrd_zdf"       long_name="j-trend: vertical  diffusion"                   unit="m/s^2"                      /> 
     872     <field id="vtrd_tau"       long_name="j-trend: wind stress "                          unit="m/s^2" grid_ref="grid_V_2D" /> 
     873     <field id="vtrd_bfr"       long_name="j-trend: bottom friction (explicit)"            unit="m/s^2"                      />    
     874     <field id="vtrd_bfri"      long_name="j-trend: bottom friction (implicit)"            unit="m/s^2"                      />    
     875     <field id="vtrd_tot"       long_name="j-trend: total momentum trend before atf"       unit="m/s^2"                      />    
     876     <field id="vtrd_atf"       long_name="j-trend: asselin time filter trend"            unit="m/s^2"                       />    
     877    </field_group> 
     878 
    631879    </field_definition> 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref

    r4205 r5034  
    1515!----------------------------------------------------------------------- 
    1616   cn_icerst_in  = "restart_ice"   !  suffix of ice restart name (input) 
    17    cn_icerst_out = "restart_ice"      !  suffix of ice restart name (output) 
    18    ln_limdyn   = .true.    !  ice dynamics (T) or thermodynamics only (F) 
    19    amax        = 0.999      !  maximum ice concentration 
    20    cai         =  1.40e-3  !  atmospheric drag over sea ice 
    21    cao         =  1.00e-3  !  atmospheric drag over ocean 
    22    ln_nicep    = .false.   !  Ice points output for debug (yes or no) 
    23    ln_limdiahsb  = .false.    !  check the heat and salt budgets (T) or not (F) 
    24    ln_limdiaout  = .false.    !  output the heat and salt budgets (T) or not (F) 
     17   cn_icerst_out = "restart_ice"   !  suffix of ice restart name (output) 
     18   ln_limdyn     = .true.          !  ice dynamics (T) or thermodynamics only (F) 
     19   amax          = 0.999           !  maximum ice concentration 
     20   cai           = 1.40e-3         !  atmospheric drag over sea ice (clio) 
     21   cao           = 1.00e-3         !  atmospheric drag over ocean   (clio) 
     22   ln_nicep      = .false.         !  Ice points output for debug (yes or no) 
     23   ln_limdiahsb  = .false.          !  check the heat and salt budgets (T) or not (F) 
     24   ln_limdiaout  = .true.          !  output the heat and salt budgets (T) or not (F) 
    2525/ 
    2626!----------------------------------------------------------------------- 
    2727&namiceini     !   ice initialisation 
    2828!----------------------------------------------------------------------- 
    29    ttest       =  2.0      !  threshold water temperature for initial sea ice 
    30    hninn       =  0.3      !  initial snow thickness in the north 
    31    hnins       =  0.1      !        "            "          south 
    32    hginn       =  3.50     !  initial undeformed ice thickness in the north 
    33    hgins       =  1.0      !        "            "              "     south 
    34    aginn       =  0.95     !  initial undeformed ice concentration in the north 
    35    agins       =  0.9      !        "            "              "         south 
    36    sinn        =  6.301    !  initial salinity in the north 
    37    sins        =  6.301    !        "            "    south 
     29   ln_limini   = .false.   !  activate ice initialization (T) or not (F) 
     30   thres_sst   =  0.0      !  threshold water temperature for initial sea ice 
     31   hts_ini_n   =  0.3      !  initial snow thickness in the north 
     32   hts_ini_s   =  0.3      !        "            "          south 
     33   hti_ini_n   =  1.0      !  initial ice thickness in the north 
     34   hti_ini_s   =  1.0      !        "            "         south 
     35   ati_ini_n   =  0.9      !  initial ice concentration in the north 
     36   ati_ini_s   =  0.9      !        "            "             south 
     37   smi_ini_n   =  6.301    !  initial ice salinity in the north 
     38   smi_ini_s   =  6.301    !        "            "    south 
     39   tmi_ini_n   =  270.     !  initial ice/snw temp in the north 
     40   tmi_ini_s   =  270.     !  initial ice/snw temp in the south 
    3841/ 
    3942!----------------------------------------------------------------------- 
     
    4144!----------------------------------------------------------------------- 
    4245   epsd        =   1.0e-20 !  tolerance parameter 
    43    alpha       =   0.5     !  coefficient for semi-implicit coriolis 
    44    dm          =   0.6e+03 !  diffusion constant for dynamics 
    45    nbiter      =   1       !  number of sub-time steps for relaxation 
    46    nbitdr      = 500       !  maximum number of iterations for relaxation 
    4746   om          =   0.5     !  relaxation constant  
    48    resl        =   5.0e-05 !  maximum value for the residual of relaxation 
    4947   cw          =   5.0e-03 !  drag coefficient for oceanic stress 
    50    angvg       =   0.0     !  turning angle for oceanic stress 
    5148   pstar       =   2.0e+04 !  1st bulk-rheology parameter 
    5249   c_rhg       =  20.0     !  2nd bulk-rhelogy parameter 
    53    etamn       =   0.0e+07 !  minimun value for viscosity 
    5450   creepl      =   1.0e-12 !  creep limit 
    5551   ecc         =   2.0     !  eccentricity of the elliptical yield curve 
    5652   ahi0        = 350.e0    !  horizontal eddy diffusivity coefficient for sea-ice [m2/s] 
    57    nevp        = 300       !  number of iterations for subcycling in EVP, SB, 400 
    58    telast      =9600.0     !  timescale for elastic waves, SB, 720.0 
     53   nevp        = 120       !  number of iterations for subcycling in EVP 
     54   relast      = 0.333     !  ratio of elastic timescale over ice time step (1/3 if nevp=120 ; 1/9 if nevp=300) 
    5955   alphaevp    =   1.0     !  coefficient for the solution of internal ice stresses 
    6056   hminrhg     =   0.001   !  ice volume (a*h in m) below which ice velocity equal ocean velocity 
    61 / 
    62 !----------------------------------------------------------------------- 
    63 &namicetrp     !   ice transport 
    64 !----------------------------------------------------------------------- 
    65    bound       =   0.      !  boundary conditions (=0.0 no-slip, =1.0 free-slip) 
    6657/ 
    6758!----------------------------------------------------------------------- 
     
    6960!----------------------------------------------------------------------- 
    7061   hmelt       = -0.15     !  maximum melting at the bottom 
    71    hiccrit     = 0.1 , 0.1 !  ice thickness for lateral accretion in the Northern (Southern) Hemisphere 
     62   hiccrit     = 0.1       !  ice thickness for lateral accretion  
    7263                           !         caution 1.0, 1.0 best value to be used!!! (gilles G.)  ???? 
    73    fraz_swi    = 0.0       !  use of frazil ice collection thickness in function of wind (1.0) or not (0.0) 
     64   fraz_swi    = 0         !  use of frazil ice collection thickness in function of wind (1.0) or not (0.0) 
    7465   maxfrazb    = 0.0       !  maximum portion of frazil ice collecting at the ice bottom 
    7566   vfrazb      = 0.4166667 !  thresold drift speed for frazil ice collecting at the ice bottom 
    7667   Cfrazb      = 5.0       !  squeezing coefficient for frazil ice collecting at the ice bottom 
    77    hicmin      = 0.2       !  ice thickness corr. to max. energy stored in brine pocket 
    7868   hiclim      = 0.10      !  minimum ice thickness 
    79    sbeta       = 1.        !  numerical caracteritic of the scheme for diffusion in ice 
    80                            !        Cranck-Nicholson (=0.5), implicit (=1), explicit (=0) 
    81    parlat      = 0.0       !  percentage of energy used for lateral ablation 
    82    hakspl      = 0.25      !  slope of distr. for Hakkinen-Mellor's lateral melting 
    83    hibspl      = 0.5       !  slope of distribution for Hibler's lateral melting 
    84    exld        = 2.0       !  exponent for leads-closure rate 
    85    hakdif      = 1.0       !  coefficient for diffusions of ice and snow 
    86    thth        = 0.2       !  threshold thickness for comp. of eq. thermal conductivity 
    8769   hnzst       = 0.1       !  thickness of the surf. layer in temp. computation 
    8870   parsub      = 1.0       !  switch for snow sublimation or not 
    89    alphs       = 1.0       !  coefficient for snow density when snow ice formation 
    9071   betas       = 0.6       !  exponent in lead-ice fractionation of snow precipitation 0.66 
    9172                           !        betas = 1 -> equipartition, betas < 1 -> more on leads 
     
    124105   astar            =   0.05 !  equivalent of gstar (0.05 for TH75 and 0.03 for weaker ice) 
    125106   Hstar            = 100.0  !  parameter determining the maximum thickness of ridged ice 
    126    raftswi          =   1    !  rafting or not 
     107   raft_swi         =   1    !  rafting or not 
    127108   hparmeter        =   0.75 !  threshold thickness for rafting or not 
    128109   Craft            =   5.0  !  coefficient used in the rafting function 
    129110   ridge_por        =   0.3  !  initial porosity of the ridged ice (typically 0.30) 
    130    sal_max_ridge    =  15.0  !  maximum ridged ice salinity 
    131111   partfun_swi      =   1    !  participation function linear, TH75 (0) or exponential Letal07 (1) 
    132    transfun_swi     =   0    !  transfer function uniform of H80 (0) or exponential Letal07 (1) 
    133112   brinstren_swi    =   0    !  (1) use brine volume to diminish ice strength 
    134113/ 
     
    144123!&namicehsb       !  Heat and salt budgets  
    145124!!----------------------------------------------------------------------- 
    146 ! 
    147125!/ 
    148 !----------------------------------------------------------------------- 
    149 &namiceout     !   parameters for outputs 
    150 !----------------------------------------------------------------------- 
    151    noumef      =   43      !  number of fields 
    152    add_diag_swi=    1      !  1 -> diagnose distribution in thickness space 
    153                            !  0 -> only simple diagnostics 
    154 ! 
    155 !           !         title of the field           !  name     !   units   !  save  ! multipl. ! additive ! 
    156 !           !                                      !           !           ! or not !  factor  !  factor  ! 
    157    field_1  = 'Ice concentration                  ', 'iiceconc', '%       ',    1   ,  1.0     ,    0.0 
    158    field_2  = 'Ice thickness                      ', 'iicethic', 'm       ',    1   ,  1.0     ,    0.0 
    159    field_3  = 'Snow thickness                     ', 'isnowthi', 'm       ',    1   ,  1.0     ,    0.0 
    160    field_4  = 'Daily bottom thermo ice production ', 'iicebopr', 'km3/day ',    1   , 1.0e-9   ,    0.0 
    161    field_5  = 'Daily dynamic ice production       ', 'iicedypr', 'km3/day ',    1   , 1.0e-9   ,    0.0 
    162    field_6  = 'Oceanic flux at the ice base       ', 'ioceflxb', 'w/m2    ',    1   ,  1.0     ,    0.0 
    163    field_7  = 'Ice velocity u                     ', 'iicevelu', 'm/s     ',    1   ,  1.0     ,    0.0 
    164    field_8  = 'Ice velocity v                     ', 'iicevelv', 'm/s     ',    1   ,  1.0     ,    0.0 
    165    field_9  = 'Sea surface temperature            ', 'isstempe', 'C       ',    1   ,  1.0     , -273.15 
    166    field_10 = 'Sea surface salinity               ', 'isssalin', 'PSU     ',    1   ,  1.0     ,    0.0 
    167    field_11 = 'Total flux at ocean surface        ', 'iocetflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
    168    field_12 = 'Solar flux at ocean surface        ', 'iocesflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
    169    field_13 = 'Non-solar flux at ocean surface    ', 'iocwnsfl', 'w/m2    ',    1   ,  1.0     ,    0.0 
    170    field_14 = 'Heat flux due to brine release     ', 'iocehebr', 'w/m2    ',    1   ,  1.0     ,    0.0 
    171    field_15 = 'Wind stress u                      ', 'iocestru', 'Pa      ',    1   ,  1.0     ,    0.0 
    172    field_16 = 'Wind stress v                      ', 'iocestrv', 'Pa      ',    1   ,  1.0     ,    0.0  
    173    field_17 = 'Solar flux at ice/ocean surface    ', 'iicesflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
    174    field_18 = 'Non-solar flux at ice/ocean surface', 'iicenflx', 'w/m2    ',    1   ,  1.0     ,    0.0 
    175    field_19 = 'Snow precipitation                 ', 'isnowpre', 'kg/m2/d ',    1   ,  1.0     ,    0.0 
    176    field_20 = 'Mean ice salinity                  ', 'iicesali', 'psu     ',    1   ,  1.0     ,    0.0 
    177    field_21 = 'Mean ice age                       ', 'iiceages', 'years   ',    1   ,  0.002739,    0.0 
    178    field_22 = 'Daily lateral thermo ice prod.     ', 'iicelapr', 'km3/day ',    1   ,1.0e-9    ,    0.0 
    179    field_23 = 'Daily snowice ice production       ', 'iicesipr', 'km3/day ',    1   ,1.0e-9    ,    0.0 
    180    field_24 = 'Mean ice temperature               ', 'iicetemp', 'C       ',    1   ,  1.0     , -273.15 
    181    field_25 = 'Ice total heat content             ', 'iiceheco', '10^9 J  ',    1   ,  1.0     ,    0.0 
    182    field_26 = 'Ice surface temperature            ', 'iicesurt', 'C       ',    1   ,  1.0     , -273.15 
    183    field_27 = 'Snow temperature                   ', 'isnotem2', 'C       ',    1   ,  1.0     , -273.15 
    184    field_28 = 'Fsbri - brine salt flux            ', 'iicefsbr', 'kg/m2/d ',    1   ,  1.0     ,    0.0 
    185    field_29 = 'Fseqv - equivalent FW salt flux    ', 'iicefseq', 'kg/m2/d ',    1   ,  1.0     ,    0.0 
    186    field_30 = 'Brine volume                       ', 'ibrinvol', '%       ',    1   ,  100.0   ,    0.0 
    187    field_31 = 'Frazil ice collection thickness    ', 'iicecolf', 'm       ',    1   ,  1.0     ,    0.0 
    188    field_32 = 'Ice strength                       ', 'iicestre', 'N/m     ',    1   ,  0.001   ,    0.0 
    189    field_33 = 'Ice velocity                       ', 'iicevelo', 'm/s     ',    1   ,  1.0     ,    0.0 
    190    field_34 = 'Surface melt                       ', 'iicesume', 'km3/day ',    1   ,1.0e-9    ,    0.0 
    191    field_35 = 'Bottom melt                        ', 'iicebome', 'km3/day ',    1   ,1.0e-9    ,    0.0 
    192    field_36 = 'Divergence                         ', 'iicedive', '10-8s-1 ',    1   ,  1.0e8   ,    0.0 
    193    field_37 = 'Shear                              ', 'iiceshea', '10-8s-1 ',    1   ,  1.0e8   ,    0.0 
    194    field_38 = 'Daily resultant ice prod/melt      ', 'iicerepr', 'km3/day ',    1   ,  1.0e-9  ,    0.0 
    195    field_39 = 'Ice volume                         ', 'iicevolu', 'km3     ',    1   ,  1.0e-9  ,    0.0 
    196    field_40 = 'Snow volume                        ', 'isnowvol', 'km3     ',    1   ,  1.0e-9  ,    0.0 
    197    field_41 = 'Fsrpo - salt flux from ridg/raft   ', 'iicefsrp', 'kg/m2/d ',    1   ,  1.0     ,    0.0 
    198    field_42 = 'Fsres - salt flux from limupdate   ', 'iicefsre', 'kg/m2/d ',    1   ,  1.0     ,    0.0 
    199    field_43 = 'Ice volume transport               ', 'iicevtrp', 'km3/day ',    1   ,1.0e-9    ,    0.0 
    200 /  
    201126 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/SHARED/namelist_ref

    r4384 r5034  
    11!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    2 !! NEMO/OPA  :  1 - run manager      (namrun, namcfg) 
    3 !! namelists    2 - Domain           (namzgr, namzgr_sco, namdom, namtsd) 
     2!! NEMO/OPA  :  1 - run manager      (namrun) 
     3!! namelists    2 - Domain           (namcfg, namzgr, namzgr_sco, namdom, namtsd) 
    44!!              3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core, namsbc_sas 
    55!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
     
    1818!!                   ***  Run management namelists  *** 
    1919!!====================================================================== 
    20 !!   namrun        parameters of the run 
     20!!   namrun       parameters of the run 
    2121!!====================================================================== 
    2222! 
     
    3131   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    3232   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
    33    nn_euler    =       1   !  = 0 : start with forward time step if ln_rstart=.true. 
    34    nn_rstctl   =       0   !  restart control => activated only if ln_rstart = T 
     33   nn_euler    =       1   !  = 0 : start with forward time step if ln_rstart=T 
     34   nn_rstctl   =       0   !  restart control ==> activated only if ln_rstart=T 
    3535                           !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
    3636                           !    = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart 
     
    4747/ 
    4848! 
    49 !----------------------------------------------------------------------- 
    50 &namcfg     !   default parameters of the configuration       
     49!!====================================================================== 
     50!!                      ***  Domain namelists  *** 
     51!!====================================================================== 
     52!!   namcfg       parameters of the configuration       
     53!!   namzgr       vertical coordinate 
     54!!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
     55!!   namdom       space and time domain (bathymetry, mesh, timestep) 
     56!!   namtsd       data: temperature & salinity 
     57!!====================================================================== 
     58! 
     59!----------------------------------------------------------------------- 
     60&namcfg     !   parameters of the configuration       
    5161!----------------------------------------------------------------------- 
    5262   cp_cfg      =  "default"            !  name of the configuration 
    53    cp_cfz      =         ''            !  name of the zoom of configuration 
     63   cp_cfz      =  "no zoom"            !  name of the zoom of configuration 
    5464   jp_cfg      =       0               !  resolution of the configuration 
    5565   jpidta      =      10               !  1st lateral dimension ( >= jpi ) 
     
    5767   jpkdta      =      31               !  number of levels      ( >= jpk ) 
    5868   jpiglo      =      10               !  1st dimension of global domain --> i =jpidta 
    59    jpjglo      =      12               !  2nd    -                  -    --> j  =jpjdta 
     69   jpjglo      =      12               !  2nd    -                  -    --> j =jpjdta 
    6070   jpizoom     =       1               !  left bottom (i,j) indices of the zoom 
    6171   jpjzoom     =       1               !  in data domain indices 
     
    6777                                       !  = 6 cyclic East-West AND North fold F-point pivot 
    6878/ 
    69 !!====================================================================== 
    70 !!                      ***  Domain namelists  *** 
    71 !!====================================================================== 
    72 !!   namzgr       vertical coordinate 
    73 !!   namzgr_sco   s-coordinate or hybrid z-s-coordinate 
    74 !!   namdom       space and time domain (bathymetry, mesh, timestep) 
    75 !!   namtsd       data: temperature & salinity 
    76 !!====================================================================== 
    77 ! 
    7879!----------------------------------------------------------------------- 
    7980&namzgr        !   vertical coordinate 
     
    8283   ln_zps      = .true.    !  z-coordinate - partial steps   (T/F) 
    8384   ln_sco      = .false.   !  s- or hybrid z-s-coordinate    (T/F) 
     85   ln_isfcav   = .false.   !  ice shelf cavity 
    8486/ 
    8587!----------------------------------------------------------------------- 
     
    183185/ 
    184186!----------------------------------------------------------------------- 
     187&namc1d        !   1D configuration options                             ("key_c1d") 
     188!----------------------------------------------------------------------- 
     189   rn_lat1d    =      50   !  Column latitude (default at PAPA station) 
     190   rn_lon1d    =    -145   !  Column longitude (default at PAPA station) 
     191   ln_c1d_locpt=  .true.   ! Localization of 1D config in a grid (T) or independant point (F) 
     192/ 
     193!----------------------------------------------------------------------- 
    185194&namtsd    !   data : Temperature  & Salinity 
    186195!----------------------------------------------------------------------- 
     
    208217!!   namtra_qsr      penetrative solar radiation 
    209218!!   namsbc_rnf      river runoffs 
     219!!   namsbc_isf      ice shelf melting/freezing 
    210220!!   namsbc_apr      Atmospheric Pressure 
    211221!!   namsbc_ssr      sea surface restoring term (for T and/or S) 
     
    223233   ln_blk_core = .true.    !  CORE bulk formulation                     (T => fill namsbc_core) 
    224234   ln_blk_mfs  = .false.   !  MFS bulk formulation                      (T => fill namsbc_mfs ) 
    225    ln_cpl      = .false.   !  Coupled formulation                       (T => fill namsbc_cpl ) 
    226235   ln_apr_dyn  = .false.   !  Patm gradient added in ocean & ice Eqs.   (T => fill namsbc_apr ) 
    227236   nn_ice      = 2         !  =0 no ice boundary condition   , 
    228237                           !  =1 use observed ice-cover      , 
    229                            !  =2 ice-model used                         ("key_lim3" or "key_lim2) 
     238                           !  =2 ice-model used                         ("key_lim3" or "key_lim2") 
    230239   nn_ice_embd = 1         !  =0 levitating ice (no mass exchange, concentration/dilution effect) 
    231240                           !  =1 levitating ice with mass and salt exchange but no presure effect 
    232241                           !  =2 embedded sea-ice (full salt and mass exchanges and pressure) 
    233242   ln_dm2dc    = .false.   !  daily mean to diurnal cycle on short wave 
    234    ln_rnf      = .true.    !  runoffs                                   (T => fill namsbc_rnf) 
     243   ln_rnf      = .true.    !  runoffs                                   (T   => fill namsbc_rnf) 
     244   nn_isf      = 0         !  ice shelf melting/freezing                (/=0 => fill namsbc_isf) 
     245                           !  0 =no isf                  1 = presence of ISF  
     246                           !  2 = bg03 parametrisation   3 = rnf file for isf    
     247                           !  4 = ISF fwf specified 
     248                           !  option 1 and 4 need ln_isfcav = .true. (domzgr) 
    235249   ln_ssr      = .true.    !  Sea Surface Restoring on T and/or S       (T => fill namsbc_ssr) 
    236250   nn_fwb      = 2         !  FreshWater Budget: =0 unchecked 
     
    240254   ln_cdgw = .false.       !  Neutral drag coefficient read from wave model (T => fill namsbc_wave) 
    241255   ln_sdw  = .false.       !  Computation of 3D stokes drift                (T => fill namsbc_wave) 
    242    cn_iceflx = 'linear'    !  redistribution of solar input into ice categories during coupling ice/atm. 
     256   nn_lsm  = 0             !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , 
     257                           !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) 
     258   nn_limflx = -1          !  LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used) 
     259                           !  =-1  Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled 
     260                           !  = 0  Average per-category fluxes (forced and coupled mode) 
     261                           !  = 1  Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled 
     262                           !  = 2  Redistribute a single flux over categories (coupled mode only) 
    243263/ 
    244264!----------------------------------------------------------------------- 
     
    296316 
    297317   cn_dir      = './'      !  root directory for the location of the bulk files 
    298    ln_2m       = .false.   !  air temperature and humidity referenced at 2m (T) instead 10m (F) 
    299318   ln_taudif   = .false.   !  HF tau contribution: use "mean of stress module - module of the mean stress" data 
    300    ln_bulk2z   = .false.   !  Air temperature/humidity and wind vectors are referenced at heights rn_zqt and rn_zu 
    301    rn_zqt      = 3.        !  Air temperature and humidity reference height (m) (ln_bulk2z) 
    302    rn_zu       = 4.        !  Wind vector reference height (m)                  (ln_bulk2z) 
     319   rn_zqt      = 10.        !  Air temperature and humidity reference height (m) 
     320   rn_zu       = 10.        !  Wind vector reference height (m)                  
    303321   rn_pfac     = 1.        !  multiplicative factor for precipitation (total & snow) 
    304322   rn_efac     = 1.        !  multiplicative factor for evaporation (0. or 1.) 
     
    327345!                    !                       ! categories !  reference  !    orientation       ! grids  ! 
    328346! send 
    329 sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
    330 sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
    331 sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   '' 
    332 sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
    333 sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
     347   sn_snd_temp   =       'weighted oce and ice' ,    'no'    ,     ''      ,         ''           ,   '' 
     348   sn_snd_alb    =       'weighted ice'         ,    'no'    ,     ''      ,         ''           ,   '' 
     349   sn_snd_thick  =       'none'                 ,    'no'   ,     ''      ,         ''           ,   '' 
     350   sn_snd_crt    =       'none'                 ,    'no'    , 'spherical' , 'eastward-northward' ,  'T' 
     351   sn_snd_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''           ,   '' 
    334352! receive 
    335 sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    336 sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    337 sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
    338 sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    339 sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
    340 sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
    341 sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
    342 sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    343 sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
    344 sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     353   sn_rcv_w10m   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
     354   sn_rcv_taumod =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     355   sn_rcv_tau    =       'oce only'             ,    'no'    , 'cartesian' , 'eastward-northward',  'U,V' 
     356   sn_rcv_dqnsdt =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     357   sn_rcv_qsr    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     358   sn_rcv_qns    =       'oce and ice'          ,    'no'    ,     ''      ,         ''          ,   '' 
     359   sn_rcv_emp    =       'conservative'         ,    'no'    ,     ''      ,         ''          ,   '' 
     360   sn_rcv_rnf    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     361   sn_rcv_cal    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     362   sn_rcv_co2    =       'coupled'              ,    'no'    ,     ''      ,         ''          ,   '' 
     363! 
     364   nn_cplmodel   =     1     !  Maximum number of models to/from which NEMO is potentialy sending/receiving data 
     365   ln_usecplmask = .false.   !  use a coupling mask file to merge data received from several models 
     366                             !   -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 
    345367/ 
    346368!----------------------------------------------------------------------- 
     
    396418   ln_rnf_tem   = .false.   !  read in temperature information for runoff 
    397419   ln_rnf_sal   = .false.   !  read in salinity information for runoff 
     420/ 
     421!----------------------------------------------------------------------- 
     422&namsbc_isf    !  Top boundary layer (ISF)  
     423!----------------------------------------------------------------------- 
     424!              ! file name ! frequency (hours) ! variable ! time interpol. !  clim   ! 'yearly'/ ! weights  ! rotation ! 
     425!              !           !  (if <0  months)  !   name   !    (logical)   !  (T/F)  ! 'monthly' ! filename ! pairing  ! 
     426! nn_isf == 4 
     427   sn_qisf      = 'rnfisf' ,         -12      ,'sohflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     428   sn_fwfisf    = 'rnfisf' ,         -12      ,'sowflisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     429! nn_isf == 3 
     430   sn_rnfisf    = 'runoffs' ,         -12      ,'sofwfisf',    .false.      , .true.  , 'yearly'  ,  ''      ,   '' 
     431! nn_isf == 2 and 3 
     432   sn_depmax_isf = 'runoffs' ,       -12        ,'sozisfmax' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
     433   sn_depmin_isf = 'runoffs' ,       -12        ,'sozisfmin' ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
     434! nn_isf == 2 
     435   sn_Leff_isf = 'rnfisf' ,       0          ,'Leff'         ,   .false.  , .true.  , 'yearly'  ,  ''      ,   '' 
     436! for all case 
     437   ln_divisf   = .true.  ! apply isf melting as a mass flux or in the salinity trend. (maybe I should remove this option as for runoff?) 
     438! only for nn_isf = 1 or 2 
     439   rn_gammat0  = 1.0e-4   ! gammat coefficient used in blk formula 
     440   rn_gammas0  = 1.0e-4   ! gammas coefficient used in blk formula 
     441! only for nn_isf = 1 
     442   nn_isfblk   =  1       ! 1 ISOMIP ; 2 conservative (3 equation formulation, Jenkins et al. 1991 ??) 
     443   rn_hisf_tbl =  30.      ! thickness of the top boundary layer           (Losh et al. 2008) 
     444                          ! 0 => thickness of the tbl = thickness of the first wet cell 
     445   ln_conserve = .true.   ! conservative case (take into account meltwater advection) 
     446   nn_gammablk = 1        ! 0 = cst Gammat (= gammat/s) 
     447                          ! 1 = velocity dependend Gamma (u* * gammat/s)  (Jenkins et al. 2010) 
     448                          !     if you want to keep the cd as in global config, adjust rn_gammat0 to compensate 
     449                          ! 2 = velocity and stability dependent Gamma    Holland et al. 1999 
    398450/ 
    399451!----------------------------------------------------------------------- 
     
    559611    nn_tra_dta    =  0                    !  = 0, bdy data are equal to the initial state 
    560612                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     613    cn_ice_lim      =  'none'             !   
     614    nn_ice_lim_dta  =  0                  !  = 0, bdy data are equal to the initial state 
     615                                          !  = 1, bdy data are read in 'bdydata   .nc' files 
     616    rn_ice_tem      = 270.                !  lim3 only: arbitrary temperature of incoming sea ice 
     617    rn_ice_sal      = 10.                 !  lim3 only:      --   salinity           -- 
     618    rn_ice_age      = 30.                 !  lim3 only:      --   age                -- 
     619 
    561620    ln_tra_dmp    =.false.                !  open boudaries conditions for tracers 
    562621    ln_dyn3d_dmp  =.false.                !  open boundary condition for baroclinic velocities 
     
    579638   bn_tem  =    'amm12_bdyT_tra' ,         24        , 'votemper' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    580639   bn_sal  =    'amm12_bdyT_tra' ,         24        , 'vosaline' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     640! for lim2 
     641!   bn_frld  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     642!   bn_hicif =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     643!   bn_hsnif =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     644! for lim3 
     645!   bn_a_i  =    'amm12_bdyT_ice' ,         24        , 'ileadfra' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     646!   bn_ht_i =    'amm12_bdyT_ice' ,         24        , 'iicethic' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
     647!   bn_ht_s =    'amm12_bdyT_ice' ,         24        , 'isnowthi' ,     .true.     , .false. ,  'daily'  ,    ''    ,   ''     , '' 
    581648   cn_dir  =    'bdydta/' 
    582649   ln_full_vel = .false. 
     
    607674   rn_bfeb2    =    2.5e-3 !  bottom turbulent kinetic energy background  (m2/s2) 
    608675   rn_bfrz0    =    3.e-3  !  bottom roughness [m] if ln_loglayer=T  
    609    ln_loglayer = .false.   !  logarithmic formulation (non linear case) 
    610676   ln_bfr2d    = .false.   !  horizontal variation of the bottom friction coef (read a 2D mask file ) 
    611677   rn_bfrien   =    50.    !  local multiplying factor of bfr (ln_bfr2d=T) 
     678   rn_tfri1    =    4.e-4  !  top drag coefficient (linear case) 
     679   rn_tfri2    =    2.5e-3 !  top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 
     680   rn_tfri2_max =   1.e-1  !  max. top drag coefficient (non linear case and ln_loglayer=T) 
     681   rn_tfeb2    =    0.0    !  top turbulent kinetic energy background  (m2/s2) 
     682   rn_tfrz0    =    3.e-3  !  top roughness [m] if ln_loglayer=T 
     683   ln_tfr2d    = .false.   !  horizontal variation of the top friction coef (read a 2D mask file ) 
     684   rn_tfrien   =    50.    !  local multiplying factor of tfr (ln_tfr2d=T) 
     685 
    612686   ln_bfrimp   = .true.    !  implicit bottom friction (requires ln_zdfexp = .false. if true) 
     687   ln_loglayer = .false.   !  logarithmic formulation (non linear case) 
    613688/ 
    614689!----------------------------------------------------------------------- 
     
    635710!!   nameos        equation of state 
    636711!!   namtra_adv    advection scheme 
     712!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.) 
    637713!!   namtra_ldf    lateral diffusion scheme 
    638714!!   namtra_dmp    T & S newtonian damping 
     
    642718&nameos        !   ocean physical parameters 
    643719!----------------------------------------------------------------------- 
    644    nn_eos      =   0       !  type of equation of state and Brunt-Vaisala frequency 
    645                            !     = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) 
    646                            !     = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T ) 
    647                            !     = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) 
    648    rn_alpha    =   2.0e-4  !  thermal expension coefficient (nn_eos= 1 or 2) 
    649    rn_beta     =   7.7e-4  !  saline  expension coefficient (nn_eos= 2) 
     720   nn_eos      =  -1     !  type of equation of state and Brunt-Vaisala frequency 
     721                                 !  =-1, TEOS-10  
     722                                 !  = 0, EOS-80  
     723                                 !  = 1, S-EOS   (simplified eos) 
     724   ln_useCT    = .true.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 
     725   !                             ! 
     726   !                     ! S-EOS coefficients : 
     727   !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
     728   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1) 
     729   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1) 
     730   rn_lambda1  =  5.9520e-2      !  cabbeling coeff in T^2  (=0 for linear eos) 
     731   rn_lambda2  =  7.4914e-4      !  cabbeling coeff in S^2  (=0 for linear eos) 
     732   rn_mu1      =  1.4970e-4      !  thermobaric coeff. in T (=0 for linear eos) 
     733   rn_mu2      =  1.1090e-5      !  thermobaric coeff. in S (=0 for linear eos) 
     734   rn_nu       =  2.4341e-3      !  cabbeling coeff in T*S  (=0 for linear eos) 
    650735/ 
    651736!----------------------------------------------------------------------- 
    652737&namtra_adv    !   advection scheme for tracer 
    653738!----------------------------------------------------------------------- 
    654    ln_traadv_cen2   =  .false.  !  2nd order centered scheme 
    655    ln_traadv_tvd    =  .true.   !  TVD scheme 
    656    ln_traadv_muscl  =  .false.  !  MUSCL scheme 
    657    ln_traadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
    658    ln_traadv_ubs    =  .false.  !  UBS scheme 
    659    ln_traadv_qck    =  .false.  !  QUICKEST scheme 
    660    ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl 
     739   ln_traadv_cen2   =  .false.   !  2nd order centered scheme 
     740   ln_traadv_tvd    =  .true.    !  TVD scheme 
     741   ln_traadv_muscl  =  .false.   !  MUSCL scheme 
     742   ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries 
     743   ln_traadv_ubs    =  .false.   !  UBS scheme 
     744   ln_traadv_qck    =  .false.   !  QUICKEST scheme 
     745   ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl 
     746   ln_traadv_tvd_zts=  .false.  !  TVD scheme with sub-timestepping of vertical tracer advection 
    661747/ 
    662748!----------------------------------------------------------------------- 
     
    733819   ln_dynadv_cen2= .false. !  flux form - 2nd order centered scheme 
    734820   ln_dynadv_ubs = .false. !  flux form - 3rd order UBS      scheme 
     821   ln_dynzad_zts = .false. !  Use (T) sub timestepping for vertical momentum advection 
    735822/ 
    736823!----------------------------------------------------------------------- 
     
    755842   ln_dynvor_mix = .false. !  mixed scheme 
    756843   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     844   ln_dynvor_een_old = .false.  !  energy & enstrophy scheme - original formulation 
    757845/ 
    758846!----------------------------------------------------------------------- 
     
    9141002!!                  ***  Miscellaneous namelists  *** 
    9151003!!====================================================================== 
     1004!!   namsol            elliptic solver / island / free surface 
    9161005!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi) 
    9171006!!   namctl            Control prints & Benchmark 
    918 !!   namsol            elliptic solver / island / free surface 
     1007!!   namc1d            1D configuration options                         ("key_c1d") 
     1008!!   namc1d_uvd        data: U & V currents                             ("key_c1d") 
     1009!!   namc1d_dyndmp     U & V newtonian damping                          ("key_c1d") 
    9191010!!====================================================================== 
    9201011! 
     
    9591050/ 
    9601051!----------------------------------------------------------------------- 
    961 &namc1d        !   1D configuration options                             ("key_c1d") 
    962 !----------------------------------------------------------------------- 
    963    rn_lat      =    50     !  Column latitude 
    964    rn_lon      =    -145   !  Column longitude 
    965 / 
    966 !----------------------------------------------------------------------- 
    9671052&namc1d_uvd    !   data: U & V currents                                 ("key_c1d") 
    9681053!----------------------------------------------------------------------- 
     
    9811066   ln_dyndmp   =  .false.  !  add a damping term (T) or not (F) 
    9821067/ 
     1068 
    9831069!!====================================================================== 
    9841070!!                  ***  Diagnostics namelists  *** 
    9851071!!====================================================================== 
    9861072!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4") 
    987 !!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld") 
     1073!!   namtrd       dynamics and/or tracer trends 
    9881074!!   namflo       float parameters                                      ("key_float") 
    9891075!!   namptr       Poleward Transport Diagnostics 
     
    10031089/ 
    10041090!----------------------------------------------------------------------- 
    1005 &namtrd        !   diagnostics on dynamics and/or tracer trends         ("key_trddyn" and/or "key_trdtra") 
    1006 !              !       or mixed-layer trends or barotropic vorticity    ("key_trdmld" or     "key_trdvor") 
    1007 !----------------------------------------------------------------------- 
    1008    nn_trd      = 365       !  time step frequency dynamics and tracers trends 
    1009    nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk) 
    1010    rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
    1011    cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input) 
    1012    cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output) 
    1013    ln_trdmld_restart = .false.         !  restart for ML diagnostics 
    1014    ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
    1015 / 
     1091&namtrd        !   diagnostics on dynamics and/or tracer trends 
     1092!              !       and/or mixed-layer trends and/or barotropic vorticity 
     1093!----------------------------------------------------------------------- 
     1094   ln_glo_trd  = .false.   ! (T) global domain averaged diag for T, T^2, KE, and PE 
     1095   ln_dyn_trd  = .false.   ! (T) 3D momentum trend output 
     1096   ln_dyn_mxl  = .FALSE.   ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 
     1097   ln_vor_trd  = .FALSE.   ! (T) 2D barotropic vorticity trends (not coded yet) 
     1098   ln_KE_trd   = .false.   ! (T) 3D Kinetic   Energy     trends 
     1099   ln_PE_trd   = .false.   ! (T) 3D Potential Energy     trends 
     1100   ln_tra_trd  = .FALSE.   ! (T) 3D tracer trend output 
     1101   ln_tra_mxl  = .false.   ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 
     1102   nn_trd      = 365       !  print frequency (ln_glo_trd=T) (unit=time step) 
     1103/ 
     1104!!gm   nn_ctls     =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk) 
     1105!!gm   rn_ucf      =   1.      !  unit conversion factor (=1 -> /seconds ; =86400. -> /day) 
     1106!!gm   cn_trdrst_in      = "restart_mld"   ! suffix of ocean restart name (input) 
     1107!!gm   cn_trdrst_out     = "restart_mld"   ! suffix of ocean restart name (output) 
     1108!!gm   ln_trdmld_restart = .false.         !  restart for ML diagnostics 
     1109!!gm   ln_trdmld_instant = .false.         !  flag to diagnose trends of instantantaneous or mean ML T/S 
     1110!!gm 
    10161111!----------------------------------------------------------------------- 
    10171112&namflo       !   float parameters                                      ("key_float") 
     
    10851180                           !     ln_ssh                  Logical switch for SSH observations 
    10861181 
    1087    ln_sst     = .true.     ! Logical switch for SST observations 
    1088    ln_reysst  = .true.     !     ln_reysst               Logical switch for Reynolds observations 
     1182   ln_sst     = .false.     ! Logical switch for SST observations 
     1183   ln_reysst  = .false.     !     ln_reysst               Logical switch for Reynolds observations 
    10891184   ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations       
    10901185 
    10911186   ln_sstfb   = .false.    ! Logical switch for feedback SST data 
    10921187                           !     ln_sss                  Logical switch for SSS observations 
    1093                            !     ln_seaice              Logical switch for Sea Ice observations 
     1188   ln_seaice  = .false.    ! Logical switch for Sea Ice observations 
    10941189                           !     ln_vel3d                Logical switch for velocity observations 
    10951190                           !     ln_velavcur             Logical switch for velocity daily av. cur. 
     
    11121207                           !     sstfiles                GHRSST input observation file name 
    11131208   !                       ! sstfbfiles: Feedback SST input observation file name 
    1114    sstfbfiles = 'sst_01.nc' 'sst_02.nc' 'sst_03.nc' 'sst_04.nc' 'sst_05.nc' 
    1115                            !     seaicefiles             Sea Ice input observation file name 
     1209   sstfbfiles = 'sst_01.nc' 
     1210                           !     seaicefiles             Sea Ice input observation file names 
     1211   seaicefiles = 'seaice_01.nc'   
    11161212                           !     velavcurfiles           Vel. cur. daily av. input file name 
    11171213                           !     velhvcurfiles           Vel. cur. high freq. input file name 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/cfg.txt

    r4529 r5034  
    11GYRE_PISCES OPA_SRC TOP_SRC 
    22ORCA2_LIM_CFC_C14b OPA_SRC LIM_SRC_2 NST_SRC TOP_SRC 
    3 GYRE OPA_SRC 
    43GYRE_XIOS OPA_SRC 
    5 ORCA2_OFF_PISCES OPA_SRC OFF_SRC TOP_SRC 
    64ORCA2_SAS_LIM OPA_SRC SAS_SRC LIM_SRC_2 NST_SRC 
    75C1D_PAPA OPA_SRC 
     6GYRE_BFM OPA_SRC TOP_SRC 
     7AMM12 OPA_SRC 
     8ORCA2_LIM_PISCES OPA_SRC LIM_SRC_2 NST_SRC TOP_SRC 
     9ORCA2_OFF_PISCES OPA_SRC OFF_SRC TOP_SRC 
     10ISOMIP OPA_SRC 
     11GYRE OPA_SRC 
    812ORCA2_LIM3 OPA_SRC LIM_SRC_3 NST_SRC 
    913ORCA2_LIM OPA_SRC LIM_SRC_2 NST_SRC 
    10 AMM12 OPA_SRC 
    11 GYRE_BFM OPA_SRC TOP_SRC 
    12 ORCA2_LIM_PISCES OPA_SRC LIM_SRC_2 NST_SRC TOP_SRC 
  • branches/2014/dev_r4650_UKMO14.12_STAND_ALONE_OBSOPER/NEMOGCM/CONFIG/makenemo

    r4148 r5034  
    3838# 
    3939# - NEW_CONF    : configuration to be created 
    40 # - REF_CONF    : reference configuration to build the new one 
     40# - REF_CONF    : reference configuration to build the new one from 
    4141# - CMP_NAM     : compiler name 
    4242# - NBR_PRC     : number of processes used to compile 
     43# - USP_CONF    : unsupported (external) configuration to build the new one from 
    4344# - NEM_SUBDIR  : NEMO subdirectory used (specified) 
    4445# 
     
    5152# - TOOLS_DIR   :   "    "    " 
    5253# - NEMO_DIR    :   "    "    " 
     54# - REMOTE_CTL  : URL link to a remote resource list for an external configuration  
     55#                 which is not part of the reference suite 
     56# - LOCAL_REF   : Nearest reference configuration to an external configuration  
     57#                 which is not part of the reference suite 
     58#                 (used to populate work directories if remote access is not available) 
    5359# 
    5460# EXAMPLES 
     
    8389x_n=""; 
    8490x_r=""; 
     91x_u=""; 
    8592x_m=""; 
    8693x_t=""; 
     
    106113export AGRIFUSE=10 
    107114declare -a TAB 
     115declare -a REMOTE_CTL 
     116declare -a LOCAL_REF 
    108117list_key=0 
    109118chk_key=1 
     
    114123#- 
    115124#- Choice of the options --- 
    116 while getopts :hd:n:r:m:j:e:s:v:t:k: V 
     125while getopts :hd:n:r:u:m:j:e:s:v:t:k: V 
    117126do 
    118127    case $V in 
    119128   (h) x_h=${OPTARG}; 
    120129        echo "Usage   : "${b_n} \ 
    121        " [-h] [-n name] [-m arch] [-d "dir1 dir2"] [-r conf] [-s Path] [-e Path] [-j No] [-v No] [-k 0/1]"; 
     130       " [-h] [-n name] [-m arch] [-d "dir1 dir2"] [-r conf] [-u conf] [-s Path] [-e Path] [-j No] [-v No] [-k 0/1]"; 
    122131   echo " -h           : help"; 
    123132   echo " -h institute : specific help for consortium members"; 
     
    126135   echo " -d dir       : choose NEMO sub-directories"; 
    127136   echo " -r conf      : choose reference configuration"; 
     137   echo " -u conf      : choose an unsupported (external) configuration"; 
    128138        echo " -s Path      : choose alternative location for NEMO main directory"; 
    129139        echo " -e Path      : choose alternative location for MY_SRC directory"; 
     
    139149   echo "Available configurations :"; cat ${CONFIG_DIR}/cfg.txt; 
    140150   echo ""; 
     151        echo "Available unsupported (external) configurations :"; cat ${CONFIG_DIR}/uspcfg.txt; 
     152   echo ""; 
    141153   echo "Example to remove bad configuration "; 
    142154   echo "./makenemo -n MY_CONFIG clean_config"; 
     
    161173   (n)  x_n=${OPTARG};; 
    162174   (r)  x_r=${OPTARG};; 
     175   (u)  x_u=${OPTARG};; 
    163176   (m)  x_m=${OPTARG};; 
    164177   (j)  x_j=${OPTARG};; 
     
    220233NEM_SUBDIR=${x_d} 
    221234REF_CONF=${x_r} 
     235USP_CONF=${x_u} 
    222236NEMO_TDIR=${x_t:-$NEMO_TDIR} 
    223237export NEMO_TDIR=${NEMO_TDIR:-$CONFIG_DIR} 
     
    228242    echo "Available configurations :" 
    229243    cat ${CONFIG_DIR}/cfg.txt 
     244    echo 
     245    echo "Available unsupported (external) configurations :"  
     246    cat ${CONFIG_DIR}/uspcfg.txt 
    230247    exit 
    231248fi 
     
    238255 
    239256if [ ${#NEW_CONF} -eq 0 ] ; then 
    240     if [ ${#NEM_SUBDIR} -eq 0 -a ${#REF_CONF} -eq 0 ]; then 
    241    echo "You are  installing a new configuration" 
     257    if [ ${#NEM_SUBDIR} -eq 0 ] && [ ${#REF_CONF} -eq 0 ] && [ ${#USP_CONF} -eq 0 ] ; then 
     258   echo "You are installing a new default (ORCA2_LIM) configuration" 
    242259   ind=0 
    243260   . ${COMPIL_DIR}/Fread_dir.sh OPA_SRC    YES 
     
    248265   . ${COMPIL_DIR}/Fread_dir.sh OFF_SRC    NO 
    249266   REF_CONF=ORCA2_LIM 
    250     elif [ ${#NEM_SUBDIR} -gt 0 ] && [ ${#REF_CONF} -eq 0 ]; then 
    251    echo "You are  installing a new configuration" 
     267    elif [ ${#NEM_SUBDIR} -gt 0 ] && [ ${#REF_CONF} -eq 0 ] && [ ${#USP_CONF} -eq 0 ] ; then 
     268   echo "You are  installing a new configuration based on ORCA2_LIM" 
    252269   TAB=( ${NEM_SUBDIR} ) 
    253270   REF_CONF=ORCA2_LIM 
     
    255272   echo "You are  installing a new configuration based on ${REF_CONF}" 
    256273   . ${COMPIL_DIR}/Fcopy_dir.sh ${REF_CONF} 
     274    elif [ ${#NEM_SUBDIR} -eq 0 ] && [ ${#USP_CONF} -gt 0 ]; then 
     275   echo "You are  installing a new configuration based on the unsupported (external) ${USP_CONF}" 
     276   . ${COMPIL_DIR}/Fcopy_extdir.sh ${USP_CONF}   
     277        #echo "TTT " ${TAB} 
     278        #echo "RRR " ${REMOTE_CTL} 
     279        #echo "LLL " ${LOCAL_REF} 
    257280    fi 
    258281    NEW_CONF=${x_n} 
    259     . ${COMPIL_DIR}/Fmake_config.sh ${NEW_CONF} ${REF_CONF} 
     282 
     283    if [ ${#USP_CONF} -gt 0 ]; then 
     284      . ${COMPIL_DIR}/Fmake_extconfig.sh ${NEW_CONF} ${LOCAL_REF} 
     285      . ${COMPIL_DIR}/Ffetch_extdir.sh ${NEW_CONF} ${REMOTE_CTL}   
     286    else 
     287      . ${COMPIL_DIR}/Fmake_config.sh ${NEW_CONF} ${REF_CONF} 
     288    fi 
    260289else 
    261290    sed -e "/${NEW_CONF} /d"  ${CONFIG_DIR}/cfg.txt >  ${COMPIL_DIR}/cfg.tmp 
Note: See TracChangeset for help on using the changeset viewer.