Changeset 3875

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/ARCH/CNRS/arch-macport_osx.fcm

r3756	r3875
27	27	%NCDF_INC -I/opt/local/include
28	28	%NCDF_LIB -L/opt/local/lib -lnetcdf -lnetcdff
29		%XIOS_ROOT /Users/~~smasson~~/XIOS
	29	%XIOS_ROOT /Users/clairelevy/Documents/WORK/XIOS
30	30	%FC mpif90
31	31	%FCFLAGS -fdefault-real-8 -O3 -funroll-all-loops -fcray-pointer

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/AMM12/cpp_AMM12.fcm

r3833	r3875
1		bld::tool::fppkeys key_bdy key_tide key_vectopt_loop key_amm_12km key_dynspg_ts key_ldfslp key_zdfgls key_vvl key_diainstant key_mpp_mpi key_iomput
	1	bld::tool::fppkeys key_bdy key_tide key_vectopt_loop key_amm_12km key_dynspg_ts key_ldfslp key_zdfgls key_vvl key_diainstant key_mpp_mpi key_iomput key_nosignedzero

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/GYRE/EXP00/iodef.xml

-                      r3771
+                      r3875
     -->
     <file_definition type="one_file" sync_freq="1d" min_digits="4">
+    <file_definition type="multiple_file" sync_freq="1d" min_digits="4">
       <file_group id="1h" output_freq="1h"  output_level="10" enabled=".TRUE."/> <!-- 1h files -->
 …
      <variable id="buffer_server_factor_size" type="integer">2</variable>
      <variable id="info_level"                type="integer">0</variable>
      <variable id="using_server"              type="boolean">true</variable>
+     <variable id="using_server"              type="boolean">false</variable>
      <variable id="using_oasis"               type="boolean">false</variable>
      <variable id="oasis_codes_id"            type="string" >oceanx</variable>

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/GYRE/cpp_GYRE.fcm

r3695	r3875
1		bld::tool::fppkeys key_gyre key_dynspg_flt key_ldfslp key_zdftke key_iomput key_mpp_mpi
	1	bld::tool::fppkeys key_gyre key_dynspg_flt key_ldfslp key_zdftke key_iomput key_mpp_mpi key_nosignedzero

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/GYRE_BFM/EXP00/solver.stat

-                      r3399
+                      r3875
+it:       1 iter: 230 r: 0.7440246784E-05 b: 0.1840836680E+00
+it:       2 iter: 250 r: 0.9359445851E-05 b: 0.2352830875E+00
+it:       3 iter: 230 r: 0.7006458826E-05 b: 0.3162751066E+00
+it:       4 iter: 250 r: 0.7963375184E-05 b: 0.7600135454E+00
+it:       5 iter: 230 r: 0.6853491637E-05 b: 0.8308380404E+00
+it:       6 iter: 230 r: 0.7072073043E-05 b: 0.9601908004E+00
+it:       7 iter: 230 r: 0.7993938546E-05 b: 0.8968937663E+00
+it:       8 iter: 220 r: 0.9452523537E-05 b: 0.9286046333E+00
+it:       9 iter: 220 r: 0.7643866533E-05 b: 0.8678317378E+00
+it:      10 iter: 230 r: 0.6734816843E-05 b: 0.8320320445E+00
+it:      11 iter: 230 r: 0.6555766837E-05 b: 0.7763847661E+00
+it:      12 iter: 230 r: 0.7881189707E-05 b: 0.7649864886E+00
+it:      13 iter: 230 r: 0.8218451842E-05 b: 0.7173828982E+00
+it:      14 iter: 230 r: 0.6144547574E-05 b: 0.7251532602E+00
+it:      15 iter: 230 r: 0.8895658888E-05 b: 0.7578256449E+00
+it:      16 iter: 240 r: 0.6909441439E-05 b: 0.7022615898E+00
+it:      17 iter: 230 r: 0.8711982405E-05 b: 0.7820485476E+00
+it:      18 iter: 240 r: 0.5646937198E-05 b: 0.6278761617E+00
+it:      19 iter: 220 r: 0.7750888632E-05 b: 0.7304352525E+00
+it:      20 iter: 230 r: 0.7665403697E-05 b: 0.6090718109E+00
+it:      21 iter: 220 r: 0.9206169267E-05 b: 0.6046990692E+00
+it:      22 iter: 220 r: 0.7880544716E-05 b: 0.6708382256E+00
+it:      23 iter: 240 r: 0.6805896028E-05 b: 0.6918136672E+00
+it:      24 iter: 240 r: 0.5257969079E-05 b: 0.7100011136E+00
+it:      25 iter: 230 r: 0.7255731129E-05 b: 0.6299497651E+00
+it:      26 iter: 230 r: 0.7955511730E-05 b: 0.5981168369E+00
+it:      27 iter: 230 r: 0.8723250269E-05 b: 0.7065752333E+00
+it:      28 iter: 230 r: 0.9084870960E-05 b: 0.6811209792E+00
+it:      29 iter: 240 r: 0.6265972753E-05 b: 0.6438025850E+00
+it:      30 iter: 240 r: 0.6756411483E-05 b: 0.7380066457E+00
+it:      31 iter: 240 r: 0.8267499945E-05 b: 0.5718579407E+00
+it:      32 iter: 230 r: 0.7822761428E-05 b: 0.5468158515E+00
+it:      33 iter: 230 r: 0.8554659234E-05 b: 0.5944574379E+00
+it:      34 iter: 230 r: 0.9814944368E-05 b: 0.6066986782E+00
+it:      35 iter: 230 r: 0.8281265862E-05 b: 0.6257082797E+00
+it:      36 iter: 240 r: 0.6434775828E-05 b: 0.7161136630E+00
+it:      37 iter: 240 r: 0.8493419909E-05 b: 0.6469647840E+00
+it:      38 iter: 230 r: 0.9514775389E-05 b: 0.6252465293E+00
+it:      39 iter: 230 r: 0.9696091963E-05 b: 0.6405217223E+00
+it:      40 iter: 240 r: 0.7487372084E-05 b: 0.6986117635E+00
+it:      41 iter: 240 r: 0.7147871952E-05 b: 0.7762620324E+00
+it:      42 iter: 240 r: 0.7872673471E-05 b: 0.6105754302E+00
+it:      43 iter: 240 r: 0.6699980937E-05 b: 0.5967819085E+00
+it:      44 iter: 240 r: 0.6723050857E-05 b: 0.6588837375E+00
+it:      45 iter: 240 r: 0.7533196699E-05 b: 0.5549833564E+00
+it:      46 iter: 220 r: 0.6344734570E-05 b: 0.5860203472E+00
+it:      47 iter: 230 r: 0.9878074607E-05 b: 0.6014543803E+00
+it:      48 iter: 240 r: 0.7035058015E-05 b: 0.5783782087E+00
+it:      49 iter: 220 r: 0.6529639860E-05 b: 0.7361370508E+00
+it:      50 iter: 240 r: 0.7729340369E-05 b: 0.5842760800E+00
+it:      51 iter: 240 r: 0.6788462114E-05 b: 0.6143198588E+00
+it:      52 iter: 220 r: 0.9937168124E-05 b: 0.6243364139E+00
+it:      53 iter: 240 r: 0.6534750798E-05 b: 0.5595607447E+00
+it:      54 iter: 230 r: 0.8678988428E-05 b: 0.7792975280E+00
+it:      55 iter: 250 r: 0.5709524264E-05 b: 0.5643273114E+00
+it:      56 iter: 240 r: 0.6868924242E-05 b: 0.6748093667E+00
+it:      57 iter: 220 r: 0.8802456851E-05 b: 0.4833796887E+00
+it:      58 iter: 230 r: 0.7963865360E-05 b: 0.5811727569E+00
+it:      59 iter: 240 r: 0.7604935477E-05 b: 0.6757019066E+00
+it:      60 iter: 240 r: 0.8249974860E-05 b: 0.6491309310E+00
+it:      61 iter: 240 r: 0.6962681203E-05 b: 0.5942964677E+00
+it:      62 iter: 240 r: 0.7920516441E-05 b: 0.5333356047E+00
+it:      63 iter: 230 r: 0.9979996401E-05 b: 0.5654869154E+00
+it:      64 iter: 230 r: 0.9440808495E-05 b: 0.5511726149E+00
+it:      65 iter: 230 r: 0.9721137141E-05 b: 0.6463726156E+00
+it:      66 iter: 240 r: 0.8200359436E-05 b: 0.7128297479E+00
+it:      67 iter: 240 r: 0.8549759019E-05 b: 0.7316743242E+00
+it:      68 iter: 240 r: 0.9628430286E-05 b: 0.5951733959E+00
+it:      69 iter: 240 r: 0.7366094543E-05 b: 0.6633287622E+00
+it:      70 iter: 240 r: 0.8231034076E-05 b: 0.6088735400E+00
+it:      71 iter: 220 r: 0.8624631868E-05 b: 0.5750411699E+00
+it:      72 iter: 240 r: 0.6294860018E-05 b: 0.6310657497E+00
+it:      73 iter: 240 r: 0.9365148279E-05 b: 0.7392710383E+00
+it:      74 iter: 250 r: 0.5364158338E-05 b: 0.5537547290E+00
+it:      75 iter: 230 r: 0.9773008643E-05 b: 0.6682402833E+00
+it:      76 iter: 240 r: 0.8100467942E-05 b: 0.4151078256E+00
+it:      77 iter: 230 r: 0.7799429163E-05 b: 0.7316424154E+00
+it:      78 iter: 240 r: 0.7895404910E-05 b: 0.5546671882E+00
+it:      79 iter: 240 r: 0.8022559953E-05 b: 0.6427847760E+00
+it:      80 iter: 240 r: 0.7277566511E-05 b: 0.6194711553E+00
+it:      81 iter: 240 r: 0.7281194538E-05 b: 0.5302877545E+00
+it:      82 iter: 220 r: 0.9419193804E-05 b: 0.5830701993E+00
+it:      83 iter: 240 r: 0.6201472391E-05 b: 0.5753005254E+00
+it:      84 iter: 240 r: 0.8012278476E-05 b: 0.7371921321E+00
+it:      85 iter: 250 r: 0.5491651366E-05 b: 0.5417053834E+00
+it:      86 iter: 240 r: 0.5899673560E-05 b: 0.6119654171E+00
+it:      87 iter: 230 r: 0.9885748439E-05 b: 0.4452258232E+00
+it:      88 iter: 230 r: 0.8059891888E-05 b: 0.5150581792E+00
+it:      89 iter: 240 r: 0.6474085695E-05 b: 0.7299448326E+00
+it:      90 iter: 240 r: 0.9823416885E-05 b: 0.5900060073E+00
+it:      91 iter: 230 r: 0.9877402224E-05 b: 0.5782497516E+00
+it:      92 iter: 240 r: 0.7243880097E-05 b: 0.6371540604E+00
+it:      93 iter: 250 r: 0.6636279056E-05 b: 0.6822048243E+00
+it:      94 iter: 240 r: 0.7952782380E-05 b: 0.6007558897E+00
+it:      95 iter: 230 r: 0.9479950980E-05 b: 0.4812863791E+00
+it:      96 iter: 230 r: 0.9716970475E-05 b: 0.5460524242E+00
+it:      97 iter: 240 r: 0.7071442197E-05 b: 0.7167370757E+00
+it:      98 iter: 240 r: 0.9140383178E-05 b: 0.7007191239E+00
+it:      99 iter: 240 r: 0.9325496852E-05 b: 0.5653204254E+00
+it:     100 iter: 240 r: 0.7224897062E-05 b: 0.6429090388E+00
+it:     101 iter: 250 r: 0.6076157518E-05 b: 0.5903397131E+00
+it:     102 iter: 240 r: 0.7709458679E-05 b: 0.5631843117E+00
+it:     103 iter: 230 r: 0.9483458449E-05 b: 0.4630274284E+00
+it:     104 iter: 220 r: 0.9163801266E-05 b: 0.5297015960E+00
+it:     105 iter: 240 r: 0.7057520098E-05 b: 0.6752365938E+00
+it:     106 iter: 240 r: 0.8819884636E-05 b: 0.7966358545E+00
+it:     107 iter: 250 r: 0.6716847991E-05 b: 0.5455663238E+00
+it:     108 iter: 240 r: 0.6937130501E-05 b: 0.6438496359E+00
+it:     109 iter: 240 r: 0.8127329934E-05 b: 0.4770795941E+00
+it:     110 iter: 230 r: 0.9027417641E-05 b: 0.5278952315E+00
+it:     111 iter: 240 r: 0.6025113978E-05 b: 0.5847697136E+00
+it:     112 iter: 240 r: 0.8995445301E-05 b: 0.6485093354E+00
+it:     113 iter: 240 r: 0.8035579774E-05 b: 0.5015333593E+00
+it:     114 iter: 230 r: 0.9945631243E-05 b: 0.6028582927E+00
+it:     115 iter: 240 r: 0.9530143332E-05 b: 0.5302317353E+00
+it:     116 iter: 240 r: 0.7375366429E-05 b: 0.6999526779E+00
+it:     117 iter: 240 r: 0.8084203070E-05 b: 0.5968685018E+00
+it:     118 iter: 240 r: 0.9186455567E-05 b: 0.5347836484E+00
+it:     119 iter: 240 r: 0.6828632897E-05 b: 0.6006343873E+00
+it:     120 iter: 240 r: 0.8126692683E-05 b: 0.5644510961E+00
+it:     121 iter: 240 r: 0.8506885049E-05 b: 0.5134367555E+00
+it:     122 iter: 240 r: 0.6041055138E-05 b: 0.5791725527E+00
+it:     123 iter: 240 r: 0.8966067574E-05 b: 0.5183401294E+00
+it:     124 iter: 230 r: 0.9630371415E-05 b: 0.4865587842E+00
+it:     125 iter: 230 r: 0.9901262506E-05 b: 0.5654058301E+00
+it:     126 iter: 240 r: 0.9839887281E-05 b: 0.6406980159E+00
+it:     127 iter: 240 r: 0.8660444831E-05 b: 0.6779025257E+00
+it:     128 iter: 240 r: 0.7588719353E-05 b: 0.6064939748E+00
+it:     129 iter: 240 r: 0.9979402524E-05 b: 0.5718319709E+00
+it:     130 iter: 240 r: 0.7130939446E-05 b: 0.5808342373E+00
+it:     131 iter: 240 r: 0.7794682727E-05 b: 0.6034740220E+00
+it:     132 iter: 240 r: 0.9500862178E-05 b: 0.5277992927E+00
+it:     133 iter: 240 r: 0.6528002864E-05 b: 0.5651071124E+00
+it:     134 iter: 240 r: 0.9017936782E-05 b: 0.5644476204E+00
+it:     135 iter: 240 r: 0.8436198587E-05 b: 0.4824425524E+00
+it:     136 iter: 240 r: 0.5647234665E-05 b: 0.5667481086E+00
+it:     137 iter: 240 r: 0.9586446435E-05 b: 0.5090651378E+00
+it:     138 iter: 230 r: 0.9471783464E-05 b: 0.4579065362E+00
+it:     139 iter: 230 r: 0.9087123846E-05 b: 0.5621798442E+00
+it:     140 iter: 250 r: 0.6164278853E-05 b: 0.6400463775E+00
+it:     141 iter: 240 r: 0.8405453366E-05 b: 0.6408662750E+00
+it:     142 iter: 240 r: 0.7106874117E-05 b: 0.4394102239E+00
+it:     143 iter: 230 r: 0.9492975149E-05 b: 0.5603112032E+00
+it:     144 iter: 240 r: 0.7017532542E-05 b: 0.4947342810E+00
+it:     145 iter: 230 r: 0.9625966339E-05 b: 0.6755118834E+00
+it:     146 iter: 240 r: 0.9014760944E-05 b: 0.4974747733E+00
+it:     147 iter: 240 r: 0.7186872796E-05 b: 0.6999869410E+00
+it:     148 iter: 250 r: 0.6162130523E-05 b: 0.5810448777E+00
+it:     149 iter: 240 r: 0.8121650155E-05 b: 0.5307579148E+00
+it:     150 iter: 240 r: 0.5841702424E-05 b: 0.4475905272E+00
+it:     151 iter: 220 r: 0.9300228120E-05 b: 0.5222786160E+00
+it:     152 iter: 230 r: 0.9744052153E-05 b: 0.4525549979E+00
+it:     153 iter: 230 r: 0.9273651310E-05 b: 0.8101635553E+00
+it:     154 iter: 250 r: 0.7361411044E-05 b: 0.6602067514E+00
+it:     155 iter: 240 r: 0.9034822204E-05 b: 0.7450641854E+00
+it:     156 iter: 240 r: 0.9522222592E-05 b: 0.5221485717E+00
+it:     157 iter: 240 r: 0.8452617503E-05 b: 0.5829628861E+00
+it:     158 iter: 240 r: 0.7316529130E-05 b: 0.6265474727E+00
+it:     159 iter: 240 r: 0.8853389481E-05 b: 0.5616413830E+00
+it:     160 iter: 240 r: 0.7884514240E-05 b: 0.5455198798E+00
+it:     161 iter: 240 r: 0.6296424346E-05 b: 0.5814948170E+00
+it:     162 iter: 240 r: 0.9250432780E-05 b: 0.5206685079E+00
+it:     163 iter: 230 r: 0.9972290864E-05 b: 0.5190431048E+00
+it:     164 iter: 240 r: 0.6340621003E-05 b: 0.5620862805E+00
+it:     165 iter: 240 r: 0.9368773786E-05 b: 0.6772619087E+00
+it:     166 iter: 240 r: 0.8238385360E-05 b: 0.5066514750E+00
+it:     167 iter: 230 r: 0.9946777196E-05 b: 0.6054272709E+00
+it:     168 iter: 240 r: 0.9499294421E-05 b: 0.5244922683E+00
+it:     169 iter: 240 r: 0.7113592963E-05 b: 0.7106784558E+00
+it:     170 iter: 240 r: 0.9228364886E-05 b: 0.5905304516E+00
+it:     171 iter: 240 r: 0.8746740247E-05 b: 0.5373804561E+00
+it:     172 iter: 240 r: 0.6849933501E-05 b: 0.6218696094E+00
+it:     173 iter: 250 r: 0.6035027619E-05 b: 0.5378961005E+00
+it:     174 iter: 240 r: 0.7412329206E-05 b: 0.5297645900E+00
+it:     175 iter: 230 r: 0.9615917480E-05 b: 0.4269978634E+00
+it:     176 iter: 220 r: 0.9097078405E-05 b: 0.4851428437E+00
+it:     177 iter: 240 r: 0.6753351497E-05 b: 0.6551058178E+00
+it:     178 iter: 240 r: 0.8992624057E-05 b: 0.7759910549E+00
+it:     179 iter: 240 r: 0.9910899945E-05 b: 0.5123581884E+00
+it:     180 iter: 240 r: 0.6677946482E-05 b: 0.6326675073E+00
+it:     181 iter: 240 r: 0.9345725231E-05 b: 0.6187279633E+00
+it:     182 iter: 240 r: 0.9236700508E-05 b: 0.5015068284E+00
+it:     183 iter: 240 r: 0.5895912167E-05 b: 0.6520692342E+00
+it:     184 iter: 250 r: 0.6157986971E-05 b: 0.5409347688E+00
+it:     185 iter: 240 r: 0.7743859030E-05 b: 0.4909818976E+00
+it:     186 iter: 230 r: 0.8962161695E-05 b: 0.4552372908E+00
+it:     187 iter: 220 r: 0.9705527833E-05 b: 0.4732608271E+00
+it:     188 iter: 240 r: 0.6635165760E-05 b: 0.6111113062E+00
+it:     189 iter: 240 r: 0.8751641668E-05 b: 0.8354887705E+00
+it:     190 iter: 250 r: 0.7272106340E-05 b: 0.4686360906E+00
+it:     191 iter: 230 r: 0.9998260643E-05 b: 0.6417323699E+00
+it:     192 iter: 240 r: 0.8352764235E-05 b: 0.4649399091E+00
+it:     193 iter: 240 r: 0.7693050177E-05 b: 0.6751386667E+00
+it:     194 iter: 240 r: 0.7837567569E-05 b: 0.5896402728E+00
+it:     195 iter: 240 r: 0.8298605400E-05 b: 0.5325820463E+00
+it:     196 iter: 240 r: 0.7578099155E-05 b: 0.5453175304E+00
+it:     197 iter: 240 r: 0.6203395717E-05 b: 0.5563993510E+00
+it:     198 iter: 240 r: 0.8656046136E-05 b: 0.5160587698E+00
+it:     199 iter: 240 r: 0.6870280862E-05 b: 0.5042090074E+00
+it:     200 iter: 240 r: 0.6214822717E-05 b: 0.5353508505E+00
+it:     201 iter: 220 r: 0.9806381626E-05 b: 0.4642334713E+00
+it:     202 iter: 230 r: 0.9562177608E-05 b: 0.4913729447E+00
+it:     203 iter: 240 r: 0.8077155329E-05 b: 0.8520973000E+00
+it:     204 iter: 250 r: 0.7388277397E-05 b: 0.6191102859E+00
+it:     205 iter: 240 r: 0.7485371726E-05 b: 0.6372750347E+00
+it:     206 iter: 240 r: 0.8769830863E-05 b: 0.5000407239E+00
+it:     207 iter: 240 r: 0.7912570954E-05 b: 0.4915069139E+00
+it:     208 iter: 240 r: 0.5592815535E-05 b: 0.6046107420E+00
+it:     209 iter: 240 r: 0.9006566098E-05 b: 0.5153272632E+00
+it:     210 iter: 240 r: 0.7145306568E-05 b: 0.4604280214E+00
+it:     211 iter: 230 r: 0.9176495722E-05 b: 0.5733940681E+00
+it:     212 iter: 240 r: 0.9035218314E-05 b: 0.4540280089E+00
+it:     213 iter: 240 r: 0.6528802051E-05 b: 0.6632050333E+00
+it:     214 iter: 240 r: 0.8343957676E-05 b: 0.5438310607E+00
+it:     215 iter: 240 r: 0.8628323766E-05 b: 0.4650558834E+00
+it:     216 iter: 240 r: 0.6040918862E-05 b: 0.5784589612E+00
+it:     217 iter: 240 r: 0.8276901507E-05 b: 0.5032676750E+00
+it:     218 iter: 230 r: 0.9607985475E-05 b: 0.4584739098E+00
+it:     219 iter: 230 r: 0.9822106038E-05 b: 0.5610157955E+00
+it:     220 iter: 240 r: 0.9442636799E-05 b: 0.6225038684E+00
+it:     221 iter: 240 r: 0.8561098712E-05 b: 0.6635010554E+00
+it:     222 iter: 240 r: 0.7477402926E-05 b: 0.6009413263E+00
+it:     223 iter: 240 r: 0.9987679875E-05 b: 0.5492499510E+00
+it:     224 iter: 240 r: 0.7248148842E-05 b: 0.5634770560E+00
+it:     225 iter: 240 r: 0.7133084535E-05 b: 0.5988824784E+00
+it:     226 iter: 240 r: 0.9522823608E-05 b: 0.4996940676E+00
+it:     227 iter: 240 r: 0.6193545299E-05 b: 0.5514029879E+00
+it:     228 iter: 240 r: 0.8249961230E-05 b: 0.5585023455E+00
+it:     229 iter: 240 r: 0.8774167477E-05 b: 0.4502157386E+00
+it:     230 iter: 240 r: 0.5467874459E-05 b: 0.5559360104E+00
+it:     231 iter: 240 r: 0.9140465089E-05 b: 0.5049968419E+00
+it:     232 iter: 230 r: 0.9579902860E-05 b: 0.4279161802E+00
+it:     233 iter: 230 r: 0.9202819655E-05 b: 0.5545415983E+00
+it:     234 iter: 250 r: 0.6285138875E-05 b: 0.6304842881E+00
+it:     235 iter: 240 r: 0.8343570272E-05 b: 0.6154112359E+00
+it:     236 iter: 240 r: 0.6505122201E-05 b: 0.4400921324E+00
+it:     237 iter: 220 r: 0.9310317421E-05 b: 0.5431094757E+00
+it:     238 iter: 240 r: 0.6784936983E-05 b: 0.4669812808E+00
+it:     239 iter: 230 r: 0.9368649385E-05 b: 0.8479074480E+00
+it:     240 iter: 250 r: 0.7017894667E-05 b: 0.4686029721E+00
+it:     241 iter: 240 r: 0.6868518574E-05 b: 0.7473331886E+00
+it:     242 iter: 240 r: 0.9601335762E-05 b: 0.4644419933E+00
+it:     243 iter: 240 r: 0.7403614387E-05 b: 0.4990407123E+00
+it:     244 iter: 240 r: 0.6035952120E-05 b: 0.6348711003E+00
+it:     245 iter: 240 r: 0.9543698964E-05 b: 0.4667730597E+00
+it:     246 iter: 230 r: 0.9244068279E-05 b: 0.5106593689E+00
+it:     247 iter: 230 r: 0.9675379433E-05 b: 0.5823567227E+00
+it:     248 iter: 240 r: 0.9850007517E-05 b: 0.6032244784E+00
+it:     249 iter: 240 r: 0.8414484256E-05 b: 0.7379488667E+00
+it:     250 iter: 240 r: 0.9359864294E-05 b: 0.5992844742E+00
+it:     251 iter: 240 r: 0.9805924128E-05 b: 0.5437868254E+00
+it:     252 iter: 240 r: 0.7200118062E-05 b: 0.6431333476E+00
+it:     253 iter: 240 r: 0.9497951929E-05 b: 0.5727311966E+00
+it:     254 iter: 240 r: 0.8608618146E-05 b: 0.5247496678E+00
+it:     255 iter: 240 r: 0.6477773688E-05 b: 0.6198132050E+00
+it:     256 iter: 250 r: 0.6167985594E-05 b: 0.5147142200E+00
+it:     257 iter: 230 r: 0.9873526591E-05 b: 0.5107478304E+00
+it:     258 iter: 230 r: 0.9242400980E-05 b: 0.4175446613E+00
+it:     259 iter: 220 r: 0.9628045476E-05 b: 0.6609339294E+00
+it:     260 iter: 240 r: 0.9025707695E-05 b: 0.6319611640E+00
+it:     261 iter: 240 r: 0.8795250954E-05 b: 0.8830745300E+00
+it:     262 iter: 250 r: 0.7693624131E-05 b: 0.5665069240E+00
+it:     263 iter: 240 r: 0.7726791671E-05 b: 0.6297963525E+00
+it:     264 iter: 240 r: 0.8326853590E-05 b: 0.5189333421E+00
+it:     265 iter: 240 r: 0.8199421374E-05 b: 0.5259497410E+00
+it:     266 iter: 240 r: 0.6449113208E-05 b: 0.5886194174E+00
+it:     267 iter: 240 r: 0.8675422182E-05 b: 0.5338163661E+00
+it:     268 iter: 240 r: 0.7591974802E-05 b: 0.4720143198E+00
+it:     269 iter: 230 r: 0.9997893412E-05 b: 0.5626420956E+00
+it:     270 iter: 240 r: 0.8914510176E-05 b: 0.4817405806E+00
+it:     271 iter: 240 r: 0.7291755721E-05 b: 0.6766204399E+00
+it:     272 iter: 240 r: 0.8136831411E-05 b: 0.5450107949E+00
+it:     273 iter: 240 r: 0.8786746638E-05 b: 0.5019283584E+00
+it:     274 iter: 240 r: 0.6716926754E-05 b: 0.5685486814E+00
+it:     275 iter: 240 r: 0.7767047962E-05 b: 0.5160557965E+00
+it:     276 iter: 230 r: 0.9965265242E-05 b: 0.4854558213E+00
+it:     277 iter: 240 r: 0.5982356324E-05 b: 0.5422752979E+00
+it:     278 iter: 240 r: 0.9054533402E-05 b: 0.6510528018E+00
+it:     279 iter: 240 r: 0.8480006691E-05 b: 0.4571427144E+00
+it:     280 iter: 240 r: 0.5504541532E-05 b: 0.5922938569E+00
+it:     281 iter: 240 r: 0.9844504652E-05 b: 0.5175970175E+00
+it:     282 iter: 240 r: 0.7407536180E-05 b: 0.4525291246E+00
+it:     283 iter: 230 r: 0.8579166437E-05 b: 0.6006559903E+00
+it:     284 iter: 240 r: 0.9724668911E-05 b: 0.4354556592E+00
+it:     285 iter: 240 r: 0.6023890176E-05 b: 0.6591813277E+00
+it:     286 iter: 240 r: 0.8502411837E-05 b: 0.5651243431E+00
+it:     287 iter: 240 r: 0.9512739362E-05 b: 0.4421835563E+00
+it:     288 iter: 240 r: 0.5685678142E-05 b: 0.5911641339E+00
+it:     289 iter: 240 r: 0.9268166632E-05 b: 0.5152137872E+00
+it:     290 iter: 230 r: 0.9724484178E-05 b: 0.4423407239E+00
+it:     291 iter: 230 r: 0.9396760769E-05 b: 0.5813430382E+00
+it:     292 iter: 250 r: 0.6577516502E-05 b: 0.6288732492E+00
+it:     293 iter: 240 r: 0.8451461695E-05 b: 0.6518149542E+00
+it:     294 iter: 240 r: 0.6781848792E-05 b: 0.4501801403E+00
+it:     295 iter: 220 r: 0.9959637262E-05 b: 0.5363171884E+00
+it:     296 iter: 240 r: 0.6778822387E-05 b: 0.4825648186E+00
+it:     297 iter: 240 r: 0.7594642927E-05 b: 0.8625428827E+00
+it:     298 iter: 250 r: 0.7483566491E-05 b: 0.4695307145E+00
+it:     299 iter: 240 r: 0.6253892640E-05 b: 0.5902785545E+00
+it:     300 iter: 240 r: 0.8294760274E-05 b: 0.4798512636E+00
+it:     301 iter: 230 r: 0.9687031296E-05 b: 0.4219920052E+00
+it:     302 iter: 230 r: 0.9225592427E-05 b: 0.5648579003E+00
+it:     303 iter: 240 r: 0.9388302120E-05 b: 0.6700980624E+00
+it:     304 iter: 240 r: 0.9153470558E-05 b: 0.5937010051E+00
+it:     305 iter: 240 r: 0.6806149190E-05 b: 0.6181465279E+00
+it:     306 iter: 250 r: 0.6545420125E-05 b: 0.5756197731E+00
+it:     307 iter: 240 r: 0.7744927602E-05 b: 0.5092369391E+00
+it:     308 iter: 230 r: 0.8744916633E-05 b: 0.4573424476E+00
+it:     309 iter: 230 r: 0.9442396183E-05 b: 0.5066584846E+00
+it:     310 iter: 240 r: 0.7084839957E-05 b: 0.6300566890E+00
+it:     311 iter: 240 r: 0.8618669529E-05 b: 0.6884236710E+00
+it:     312 iter: 240 r: 0.9275503474E-05 b: 0.5036022021E+00
+it:     313 iter: 240 r: 0.6919289411E-05 b: 0.5728547455E+00
+it:     314 iter: 240 r: 0.8942676909E-05 b: 0.5844163973E+00
+it:     315 iter: 240 r: 0.8388772295E-05 b: 0.4846939204E+00
+it:     316 iter: 240 r: 0.5739100066E-05 b: 0.5831024730E+00
+it:     317 iter: 240 r: 0.9434372940E-05 b: 0.5178686374E+00
+it:     318 iter: 230 r: 0.9861691195E-05 b: 0.4656964406E+00
+it:     319 iter: 230 r: 0.8971064070E-05 b: 0.5733799575E+00
+it:     320 iter: 240 r: 0.9903829223E-05 b: 0.6593697397E+00
+it:     321 iter: 240 r: 0.9014966680E-05 b: 0.6684955387E+00
+it:     322 iter: 240 r: 0.8375983376E-05 b: 0.6251140670E+00
+it:     323 iter: 250 r: 0.6447050802E-05 b: 0.5607851626E+00
+it:     324 iter: 240 r: 0.7281354067E-05 b: 0.5898233933E+00
+it:     325 iter: 240 r: 0.7891438283E-05 b: 0.4257134507E+00
+it:     326 iter: 230 r: 0.9111914991E-05 b: 0.5140988199E+00
+it:     327 iter: 240 r: 0.6421472217E-05 b: 0.5071089201E+00
+it:     328 iter: 230 r: 0.9888393396E-05 b: 0.6476089032E+00
+it:     329 iter: 240 r: 0.8199863171E-05 b: 0.4681191792E+00
+it:     330 iter: 240 r: 0.6580117926E-05 b: 0.7135228391E+00
+it:     331 iter: 250 r: 0.6961575827E-05 b: 0.5216087612E+00
+it:     332 iter: 240 r: 0.6795216420E-05 b: 0.5342925318E+00
+it:     333 iter: 230 r: 0.9272336578E-05 b: 0.4464092554E+00
+it:     334 iter: 230 r: 0.9358211024E-05 b: 0.4543204929E+00
+it:     335 iter: 240 r: 0.6494875853E-05 b: 0.6897801097E+00
+it:     336 iter: 240 r: 0.9398593229E-05 b: 0.6483492695E+00
+it:     337 iter: 240 r: 0.8970751734E-05 b: 0.4861260074E+00
+it:     338 iter: 240 r: 0.6492204211E-05 b: 0.6283166015E+00
+it:     339 iter: 250 r: 0.6282521245E-05 b: 0.5356801975E+00
+it:     340 iter: 240 r: 0.7866565313E-05 b: 0.4963618785E+00
+it:     341 iter: 230 r: 0.9110858297E-05 b: 0.4354634931E+00
+it:     342 iter: 220 r: 0.9622046653E-05 b: 0.4655574988E+00
+it:     343 iter: 240 r: 0.6678022741E-05 b: 0.6184513419E+00
+it:     344 iter: 240 r: 0.8829474468E-05 b: 0.8104726331E+00
+it:     345 iter: 250 r: 0.7157722604E-05 b: 0.4727655304E+00
+it:     346 iter: 240 r: 0.6433293473E-05 b: 0.6331781769E+00
+it:     347 iter: 240 r: 0.8376333111E-05 b: 0.4483322391E+00
+it:     348 iter: 230 r: 0.9032626249E-05 b: 0.4573054905E+00
+it:     349 iter: 230 r: 0.9923279836E-05 b: 0.5803391117E+00
+it:     350 iter: 240 r: 0.9158589477E-05 b: 0.6182516123E+00
+it:     351 iter: 240 r: 0.8742408528E-05 b: 0.6654006406E+00
+it:     352 iter: 240 r: 0.7742233892E-05 b: 0.6065475129E+00
+it:     353 iter: 240 r: 0.9911840287E-05 b: 0.5503476040E+00
+it:     354 iter: 240 r: 0.7566892071E-05 b: 0.5592587298E+00
+it:     355 iter: 240 r: 0.7058886898E-05 b: 0.6070216169E+00
+it:     356 iter: 240 r: 0.9532244842E-05 b: 0.5021220588E+00
+it:     357 iter: 240 r: 0.6581905373E-05 b: 0.5428098796E+00
+it:     358 iter: 240 r: 0.7739991638E-05 b: 0.5710056940E+00
+it:     359 iter: 240 r: 0.8758601072E-05 b: 0.4513770944E+00
+it:     360 iter: 240 r: 0.5612180673E-05 b: 0.5456855029E+00
+it:     361 iter: 240 r: 0.8818081887E-05 b: 0.5219130511E+00
+it:     362 iter: 240 r: 0.7693135936E-05 b: 0.4230324205E+00
+it:     363 iter: 230 r: 0.9330402960E-05 b: 0.5466741815E+00
+it:     364 iter: 240 r: 0.9329723621E-05 b: 0.4571170068E+00
+it:     365 iter: 230 r: 0.9597716286E-05 b: 0.6113816086E+00
+it:     366 iter: 240 r: 0.6862418219E-05 b: 0.5346482659E+00
+it:     367 iter: 240 r: 0.9579590312E-05 b: 0.6636934893E+00
+it:     368 iter: 240 r: 0.8884686379E-05 b: 0.5200160853E+00
+it:     369 iter: 240 r: 0.6925065330E-05 b: 0.5806576814E+00
+it:     370 iter: 240 r: 0.9298083171E-05 b: 0.5320130285E+00
+it:     371 iter: 240 r: 0.6729450195E-05 b: 0.5107096618E+00
+it:     372 iter: 240 r: 0.8300216838E-05 b: 0.5725637345E+00
+it:     373 iter: 240 r: 0.8463547174E-05 b: 0.4607061379E+00
+it:     374 iter: 240 r: 0.5818753748E-05 b: 0.5183287506E+00
+it:     375 iter: 240 r: 0.9023857877E-05 b: 0.5269511380E+00
+it:     376 iter: 240 r: 0.7381909967E-05 b: 0.4197384345E+00
+it:     377 iter: 230 r: 0.9174465056E-05 b: 0.5253456364E+00
+it:     378 iter: 240 r: 0.9177801147E-05 b: 0.4661442161E+00
+it:     379 iter: 240 r: 0.6739015755E-05 b: 0.5961511088E+00
+it:     380 iter: 240 r: 0.7138913606E-05 b: 0.5186325957E+00
+it:     381 iter: 220 r: 0.9957458341E-05 b: 0.4708350683E+00
+it:     382 iter: 240 r: 0.6214118912E-05 b: 0.5084594626E+00
+it:     383 iter: 240 r: 0.8884779870E-05 b: 0.8302238398E+00
+it:     384 iter: 250 r: 0.6888254275E-05 b: 0.4363495398E+00
+it:     385 iter: 230 r: 0.9720238700E-05 b: 0.6293956176E+00
+it:     386 iter: 240 r: 0.9602391250E-05 b: 0.4306623454E+00
+it:     387 iter: 240 r: 0.7496942815E-05 b: 0.6284278745E+00
+it:     388 iter: 240 r: 0.7372637604E-05 b: 0.5847502653E+00
+it:     389 iter: 240 r: 0.9468369013E-05 b: 0.4885956905E+00
+it:     390 iter: 240 r: 0.7094294922E-05 b: 0.5274549895E+00
+it:     391 iter: 240 r: 0.6769266082E-05 b: 0.5442979043E+00
+it:     392 iter: 240 r: 0.8962624793E-05 b: 0.4623300794E+00
+it:     393 iter: 230 r: 0.9990717937E-05 b: 0.5163570594E+00
+it:     394 iter: 240 r: 0.7721915098E-05 b: 0.5096005808E+00
+it:     395 iter: 240 r: 0.8505263053E-05 b: 0.6318946869E+00
+it:     396 iter: 240 r: 0.8106982358E-05 b: 0.5139391077E+00
+it:     397 iter: 240 r: 0.8482016367E-05 b: 0.5357793469E+00
+it:     398 iter: 240 r: 0.7605386113E-05 b: 0.4951472641E+00
+it:     399 iter: 240 r: 0.6418019887E-05 b: 0.5151271467E+00
+it:     400 iter: 230 r: 0.9969223169E-05 b: 0.5131293456E+00
+it:     401 iter: 230 r: 0.9992956814E-05 b: 0.4492217122E+00
+it:     402 iter: 230 r: 0.9493837033E-05 b: 0.6823456869E+00
+it:     403 iter: 250 r: 0.7125462203E-05 b: 0.6722935400E+00
+it:     404 iter: 240 r: 0.9957566837E-05 b: 0.7253749806E+00
+it:     405 iter: 240 r: 0.8891937286E-05 b: 0.4636747592E+00
+it:     406 iter: 240 r: 0.8185325290E-05 b: 0.6021556031E+00
+it:     407 iter: 240 r: 0.8514384324E-05 b: 0.5698450922E+00
+it:     408 iter: 240 r: 0.8268173790E-05 b: 0.5223330648E+00
+it:     409 iter: 230 r: 0.9915960192E-05 b: 0.5589034893E+00
+it:     410 iter: 240 r: 0.7578172603E-05 b: 0.5145736742E+00
+it:     411 iter: 240 r: 0.9362076096E-05 b: 0.7068986819E+00
+it:     412 iter: 240 r: 0.8840618595E-05 b: 0.5057077629E+00
+it:     413 iter: 240 r: 0.6776175784E-05 b: 0.5686641264E+00
+it:     414 iter: 240 r: 0.9850784521E-05 b: 0.5733874559E+00
+it:     415 iter: 240 r: 0.7719911188E-05 b: 0.4785166967E+00
+it:     416 iter: 230 r: 0.9253049781E-05 b: 0.5869789134E+00
+it:     417 iter: 240 r: 0.9727638040E-05 b: 0.4936785828E+00
+it:     418 iter: 240 r: 0.7357747291E-05 b: 0.6731020200E+00
+it:     419 iter: 240 r: 0.8864878948E-05 b: 0.5607498834E+00
+it:     420 iter: 240 r: 0.9042773412E-05 b: 0.5081433387E+00
+it:     421 iter: 240 r: 0.7122642675E-05 b: 0.5837552421E+00
+it:     422 iter: 240 r: 0.9742336129E-05 b: 0.5209029834E+00
+it:     423 iter: 240 r: 0.7562697089E-05 b: 0.4980639854E+00
+it:     424 iter: 240 r: 0.6408566905E-05 b: 0.5642943617E+00
+it:     425 iter: 240 r: 0.9824615815E-05 b: 0.4688218575E+00
+it:     426 iter: 230 r: 0.9129894730E-05 b: 0.4883165764E+00
+it:     427 iter: 230 r: 0.9232905961E-05 b: 0.5424591442E+00
+it:     428 iter: 240 r: 0.9885656084E-05 b: 0.6080138351E+00
+it:     429 iter: 240 r: 0.8684206709E-05 b: 0.6998539087E+00
+it:     430 iter: 240 r: 0.9333254608E-05 b: 0.5672755834E+00
+it:     431 iter: 240 r: 0.9584027329E-05 b: 0.5469060250E+00
+it:     432 iter: 240 r: 0.7875190959E-05 b: 0.6072683806E+00
+it:     433 iter: 250 r: 0.5683773505E-05 b: 0.5484402581E+00
+it:     434 iter: 240 r: 0.8036333302E-05 b: 0.5256307673E+00
+it:     435 iter: 240 r: 0.6285860484E-05 b: 0.3885253755E+00
+it:     436 iter: 220 r: 0.9006940402E-05 b: 0.4995674972E+00
+it:     437 iter: 230 r: 0.9768021627E-05 b: 0.4121814020E+00
+it:     438 iter: 230 r: 0.8855905043E-05 b: 0.7424947274E+00
+it:     439 iter: 250 r: 0.6999985699E-05 b: 0.6453379767E+00
+it:     440 iter: 240 r: 0.9307017519E-05 b: 0.6903100806E+00
+it:     441 iter: 240 r: 0.8981222548E-05 b: 0.4821343709E+00
+it:     442 iter: 230 r: 0.9948547069E-05 b: 0.5717358533E+00
+it:     443 iter: 240 r: 0.7786240736E-05 b: 0.5668128653E+00
+it:     444 iter: 240 r: 0.8896388689E-05 b: 0.7276095222E+00
+it:     445 iter: 240 r: 0.9285725823E-05 b: 0.5178496122E+00
+it:     446 iter: 240 r: 0.7112559669E-05 b: 0.5962758277E+00
+it:     447 iter: 250 r: 0.6208478821E-05 b: 0.6037386886E+00
+it:     448 iter: 240 r: 0.8469099610E-05 b: 0.4830438756E+00
+it:     449 iter: 230 r: 0.9201679690E-05 b: 0.4470009097E+00
+it:     450 iter: 220 r: 0.9410464807E-05 b: 0.5274449367E+00
+it:     451 iter: 240 r: 0.7766112732E-05 b: 0.6044952597E+00
+it:     452 iter: 240 r: 0.8527821750E-05 b: 0.8504510693E+00
+it:     453 iter: 250 r: 0.7240964669E-05 b: 0.5143883743E+00
+it:     454 iter: 240 r: 0.7489859460E-05 b: 0.5922246735E+00
+it:     455 iter: 240 r: 0.8150684278E-05 b: 0.4857339643E+00
+it:     456 iter: 240 r: 0.7724457443E-05 b: 0.4876586048E+00
+it:     457 iter: 240 r: 0.6042020297E-05 b: 0.5549232420E+00
+it:     458 iter: 240 r: 0.8374412856E-05 b: 0.4953621744E+00
+it:     459 iter: 230 r: 0.9573243419E-05 b: 0.4447831730E+00
+it:     460 iter: 230 r: 0.9381329858E-05 b: 0.5309644787E+00
+it:     461 iter: 240 r: 0.9096529844E-05 b: 0.6416171934E+00
+it:     462 iter: 240 r: 0.9356532044E-05 b: 0.6327577836E+00
+it:     463 iter: 240 r: 0.7628339945E-05 b: 0.5861280263E+00
+it:     464 iter: 240 r: 0.9310552466E-05 b: 0.5580416643E+00
+it:     465 iter: 240 r: 0.7485669403E-05 b: 0.5403008993E+00
+it:     466 iter: 240 r: 0.7590805052E-05 b: 0.5849930918E+00
+it:     467 iter: 240 r: 0.8711915552E-05 b: 0.5042455237E+00
+it:     468 iter: 240 r: 0.6769849941E-05 b: 0.5298753819E+00
+it:     469 iter: 240 r: 0.8625549459E-05 b: 0.5461276100E+00
+it:     470 iter: 240 r: 0.7776347800E-05 b: 0.4534267118E+00
+it:     471 iter: 240 r: 0.5593618373E-05 b: 0.5323264584E+00
+it:     472 iter: 240 r: 0.8997664200E-05 b: 0.4923070869E+00
+it:     473 iter: 230 r: 0.9508383995E-05 b: 0.4317418978E+00
+it:     474 iter: 230 r: 0.8611632240E-05 b: 0.5294936990E+00
+it:     475 iter: 240 r: 0.9258400228E-05 b: 0.6313671325E+00
+it:     476 iter: 240 r: 0.8830978793E-05 b: 0.6126531403E+00
+it:     477 iter: 240 r: 0.7799803557E-05 b: 0.5809797260E+00
+it:     478 iter: 240 r: 0.9417847817E-05 b: 0.5402032508E+00
+it:     479 iter: 240 r: 0.7268906785E-05 b: 0.5385576541E+00
+it:     480 iter: 240 r: 0.8569757263E-05 b: 0.5725175500E+00
+it:     481 iter: 240 r: 0.8475435149E-05 b: 0.4905396017E+00
+it:     482 iter: 240 r: 0.6424251242E-05 b: 0.5364602675E+00
+it:     483 iter: 240 r: 0.9134335064E-05 b: 0.5240236375E+00
+it:     484 iter: 240 r: 0.7416399646E-05 b: 0.4484277866E+00
+it:     485 iter: 230 r: 0.9218489387E-05 b: 0.5391608057E+00
+it:     486 iter: 240 r: 0.9223369078E-05 b: 0.4644560614E+00
+it:     487 iter: 240 r: 0.7005374647E-05 b: 0.6378454473E+00
+it:     488 iter: 240 r: 0.7907276860E-05 b: 0.5260457717E+00
+it:     489 iter: 240 r: 0.9134300298E-05 b: 0.4802181676E+00
+it:     490 iter: 240 r: 0.6511731237E-05 b: 0.5403219967E+00
+it:     491 iter: 240 r: 0.8293826991E-05 b: 0.5000780339E+00
+it:     492 iter: 230 r: 0.9743260463E-05 b: 0.4640877349E+00
+it:     493 iter: 240 r: 0.6066313601E-05 b: 0.5200489203E+00
+it:     494 iter: 240 r: 0.9602688039E-05 b: 0.6348974306E+00
+it:     495 iter: 240 r: 0.8320959887E-05 b: 0.4386132859E+00
+it:     496 iter: 230 r: 0.9443219904E-05 b: 0.5712801852E+00
+it:     497 iter: 250 r: 0.6223852704E-05 b: 0.4987189619E+00
+it:     498 iter: 240 r: 0.7078328173E-05 b: 0.6403667777E+00
+it:     499 iter: 240 r: 0.6915966497E-05 b: 0.4042775883E+00
+it:     500 iter: 220 r: 0.8664600046E-05 b: 0.4916638907E+00
+it:     501 iter: 240 r: 0.6711099272E-05 b: 0.4742347645E+00
+it:     502 iter: 230 r: 0.9650599205E-05 b: 0.7509994968E+00
+it:     503 iter: 240 r: 0.8934242190E-05 b: 0.4543437950E+00
+it:     504 iter: 240 r: 0.7328752331E-05 b: 0.7420407966E+00
+it:     505 iter: 250 r: 0.6598827075E-05 b: 0.5884958251E+00
+it:     506 iter: 240 r: 0.8212692098E-05 b: 0.5053207308E+00
+it:     507 iter: 230 r: 0.9641290906E-05 b: 0.4774920508E+00
+it:     508 iter: 240 r: 0.8446481906E-05 b: 0.5142522888E+00
+it:     509 iter: 240 r: 0.7764252095E-05 b: 0.6518201897E+00
+it:     510 iter: 240 r: 0.7247433046E-05 b: 0.5365145975E+00
+it:     511 iter: 240 r: 0.8753242134E-05 b: 0.5110914352E+00
+it:     512 iter: 240 r: 0.7214231487E-05 b: 0.5192317487E+00
+it:     513 iter: 240 r: 0.6794496960E-05 b: 0.5265835226E+00
+it:     514 iter: 240 r: 0.8257783113E-05 b: 0.4778625755E+00
+it:     515 iter: 240 r: 0.6438038030E-05 b: 0.5007288407E+00
+it:     516 iter: 240 r: 0.7891620437E-05 b: 0.4941909213E+00
+it:     517 iter: 230 r: 0.9681402627E-05 b: 0.4341059740E+00
+it:     518 iter: 230 r: 0.9372592410E-05 b: 0.4966207185E+00
+it:     519 iter: 240 r: 0.8916130859E-05 b: 0.6315138857E+00
+it:     520 iter: 240 r: 0.9111939995E-05 b: 0.6113427916E+00
+it:     521 iter: 240 r: 0.7566852609E-05 b: 0.5552510727E+00
+it:     522 iter: 240 r: 0.9400566307E-05 b: 0.5644593327E+00
+it:     523 iter: 240 r: 0.7735477877E-05 b: 0.5041612002E+00
+it:     524 iter: 230 r: 0.9910217652E-05 b: 0.5639666282E+00
+it:     525 iter: 240 r: 0.9020769555E-05 b: 0.5112527467E+00
+it:     526 iter: 240 r: 0.7816155122E-05 b: 0.6997629593E+00
+it:     527 iter: 240 r: 0.9073213369E-05 b: 0.5421948127E+00
+it:     528 iter: 240 r: 0.8376702758E-05 b: 0.5392260510E+00

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/GYRE_BFM/cpp_GYRE_BFM.fcm

r3695	r3875
1		bld::tool::fppkeys key_gyre key_dynspg_flt key_ldfslp key_zdftke key_vectopt_loop key_top key_my_trc key_mpp_mpi key_iomput
	1	bld::tool::fppkeys key_gyre key_dynspg_flt key_ldfslp key_zdftke key_vectopt_loop key_top key_my_trc key_mpp_mpi key_iomput key_nosignedzero
2	2	inc $BFMDIR/src/nemo/bfm.fcm

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/GYRE_PISCES/cpp_GYRE_PISCES.fcm

r3695	r3875
1		bld::tool::fppkeys key_gyre key_dynspg_flt key_ldfslp key_zdftke key_top key_pisces_reduced key_iomput key_mpp_mpi
	1	bld::tool::fppkeys key_gyre key_dynspg_flt key_ldfslp key_zdftke key_top key_pisces_reduced key_iomput key_mpp_mpi key_nosignedzero

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/ORCA2_LIM/cpp_ORCA2_LIM.fcm

r3695	r3875
1		bld::tool::fppkeys key_trabbl key_orca_r2 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
	1	bld::tool::fppkeys key_trabbl key_orca_r2 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_nosignedzero

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/STRAIT.dat

-                      r3567
+                      r3875
  Net freshwater budget
    fwf    = -3.43323E-06 m3 =   0.00000 Sv
+   fwf    =  5.90515E-05 m3 =   0.00000 Sv
    zarea = 3.57005E+14
    zvol  = 1.32653E+18
  Mean sea level :
    at nit000 =  0.00000E+00 m3
    at nitend = -3.24688E+12 m3
    diff      = -3.24688E+12 m3 =  -1.25266 Sv
    mean sea level elevation    =  -0.00909 m
+   at nitend = -3.98438E-01 m3
+   diff      = -3.98438E-01 m3 =  -0.00000 Sv
+   mean sea level elevation    =  -0.00000 m
  Anomaly of salinity content :
    at nit000 =  5.65470E+14 psu.m3
    at nitend = -3.35735E+13 psu.m3
    diff      = -5.99044E+14 psu.m3
+   at nit000 = -2.20495E+15 psu.m3
+   at nitend = -2.76348E+15 psu.m3
+   diff      = -5.58530E+14 psu.m3
  Mean salinity :
    at nit000 =  34.72697 psu
    at nitend =  34.72651 psu
    diff      =  -0.00045 psu
    S-SLevitus=  -0.00003 psu
+   at nit000 =  34.72488 psu
+   at nitend =  34.72446 psu
+   diff      =  -0.00042 psu
+   S-SLevitus=  -0.00208 psu
  Gibraltar :
    Flux entrant (Sv) :   0.4037
    Flux sortant (Sv) :  -0.8196
    T entrant (deg)   :  13.1923
    T sortant (deg)   :  13.0994
    S entrant (psu)   :  31.3817
    S sortant (psu)   :  37.7129
+   Flux entrant (Sv) :   0.5713
+   Flux sortant (Sv) :  -0.5958
+   T entrant (deg)   :  16.1815
+   T sortant (deg)   :  12.6934
+   S entrant (psu)   :  36.0673
+   S sortant (psu)   :  37.8858
  Cadiz :
    Flux entrant (Sv) :   4.5633
    Flux sortant (Sv) :  -4.9890
    T entrant (deg)   :  13.2328
    T sortant (deg)   :  11.6841
    S entrant (psu)   :  35.8927
    S sortant (psu)   :  35.9800
+   Flux entrant (Sv) :   2.5583
+   Flux sortant (Sv) :  -2.5828
+   T entrant (deg)   :  14.3696
+   T sortant (deg)   :  11.9981
+   S entrant (psu)   :  36.0057
+   S sortant (psu)   :  36.0549
  Bab el Mandeb :
    Flux entrant (Sv) :   0.4526
    Flux sortant (Sv) :  -0.4375
    T entrant (deg)   :  10.1702
    T sortant (deg)   :  25.0599
    S entrant (psu)   :  37.7353
    S sortant (psu)   :  36.2338
+   Flux entrant (Sv) :   0.2526
+   Flux sortant (Sv) :  -0.2675
+   T entrant (deg)   :  15.8806
+   T sortant (deg)   :  23.5775
+   S entrant (psu)   :  36.9325
+   S sortant (psu)   :  36.1224
  Baltic :
    Flux entrant (Sv) :   0.0775
    Flux sortant (Sv) :  -0.0876
    T entrant (deg)   :   1.5469
    T sortant (deg)   :   2.6595
    S entrant (psu)   :   8.1471
    S sortant (psu)   :  14.5510
+   Flux entrant (Sv) :   0.0457
+   Flux sortant (Sv) :  -0.0617
+   T entrant (deg)   :   3.8966
+   T sortant (deg)   :   5.5943
+   S entrant (psu)   :   7.0356
+   S sortant (psu)   :   9.7869

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/cpp_ORCA2_LIM_CFC_C14b.fcm

r3695	r3875
1		bld::tool::fppkeys key_trabbl key_orca_r2 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_top key_cfc key_c14b key_iomput key_mpp_mpi
	1	bld::tool::fppkeys key_trabbl key_orca_r2 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_top key_cfc key_c14b key_iomput key_mpp_mpi key_nosignedzero

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/cpp_ORCA2_LIM_PISCES.fcm

r3695	r3875
1		bld::tool::fppkeys key_trabbl key_orca_r2 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_top key_pisces key_iomput key_mpp_mpi
	1	bld::tool::fppkeys key_trabbl key_orca_r2 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_top key_pisces key_iomput key_mpp_mpi key_nosignedzero

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/cpp_ORCA2_OFF_PISCES.fcm

r3695	r3875
1		bld::tool::fppkeys key_trabbl key_orca_r2 key_ldfslp key_traldf_c2d key_traldf_eiv key_top key_offline key_pisces key_iomput key_mpp_mpi
	1	bld::tool::fppkeys key_trabbl key_orca_r2 key_ldfslp key_traldf_c2d key_traldf_eiv key_top key_offline key_pisces key_iomput key_mpp_mpi key_nosignedzero

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/ORCA2_SAS_LIM/EXP00/opa

r3331	r3875
1		link /~~noc/msm/working/nemo/sga/NEMO/FORGE/SAS~~/NEMOGCM/CONFIG/ORCA2_SAS_LIM/BLD/bin/nemo.exe
	1	link /Users/clairelevy/Documents/WORK/CONFIG_SETTING2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/ORCA2_SAS_LIM/BLD/bin/nemo.exe

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/ORCA2_SAS_LIM/cpp_ORCA2_SAS_LIM.fcm

r3695	r3875
1		bld::tool::fppkeys key_trabbl key_orca_r2 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
	1	bld::tool::fppkeys key_trabbl key_orca_r2 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_nosignedzero

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/CONFIG/cfg.txt

-                      r3769
+                      r3875
+ORCA2_LIM3 OPA_SRC LIM_SRC_3
+ORCA2_LIM OPA_SRC LIM_SRC_2 NST_SRC
+ORCA2_LIM_PISCES OPA_SRC LIM_SRC_2 NST_SRC TOP_SRC
 AMM12 OPA_SRC
+GYRE_PISCES OPA_SRC TOP_SRC
 GYRE OPA_SRC
 GYRE_BFM OPA_SRC TOP_SRC
+GYRE_PISCES OPA_SRC TOP_SRC
+ORCA2_LIM3 OPA_SRC LIM_SRC_3
+ORCA2_LIM_PISCES OPA_SRC LIM_SRC_2 NST_SRC TOP_SRC
+O2L3 OPA_SRC LIM_SRC_3 TOP_SRC
 ORCA2_OFF_PISCES OPA_SRC OFF_SRC TOP_SRC
 ORCA2_SAS_LIM OPA_SRC SAS_SRC LIM_SRC_2 NST_SRC
 ORCA2_LIM_CFC_C14b OPA_SRC LIM_SRC_2 NST_SRC TOP_SRC
-ORCA2_LIM OPA_SRC LIM_SRC_2 NST_SRC

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_2/iceini_2.F90

-                      r3625
+                      r3875
          &     'use more ocean levels or less ice layers/categories.' )
+      !                                ! Open the namelist file
+      CALL ctl_opn( numnam_ice, 'namelist_ice', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
+      !                                ! Open the reference and configuration namelist files and namelist output file
+      CALL ctl_opn( numnam_ice_ref, 'namelist_ice_ref',    'OLD',     'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
+      CALL ctl_opn( numnam_ice_cfg, 'namelist_ice_cfg',    'OLD',     'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
+      CALL ctl_opn( numoni,         'output.namelist.ice', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
+      !
       CALL ice_run_2                   ! read in namelist some run parameters
 …
       !!-------------------------------------------------------------------
       NAMELIST/namicerun/ cn_icerst_in, cn_icerst_out, ln_limdyn, ln_limdmp, acrit, hsndif, hicdif
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       !!-------------------------------------------------------------------
+      !
+      REWIND( numnam_ice )                      ! Read Namelist namicerun
+      READ  ( numnam_ice , namicerun )
+      REWIND( numnam_ice_ref )              ! Namelist namicerun in reference namelist : Parameters for ice
+      READ  ( numnam_ice_ref, namicerun, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicerun in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namicerun in configuration namelist : Parameters for ice
+      READ  ( numnam_ice_cfg, namicerun, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicerun in configuration namelist', lwp )
+      WRITE ( numoni, namicerun )
+      !
       IF(lwp) THEN                              ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_2/limdia_2.F90

-                      r3625
+                      r3875
        INTEGER  ::   jv   ! dummy loop indice
        INTEGER  ::   ntot , ndeb, nv, ierr   ! local integer
+       INTEGER  ::   ios                     ! Local integer output status for namelist read
        REAL(wp) ::   zxx0, zxx1              ! local scalars
        NAMELIST/namicedia/fmtinf, nfrinf, ninfo, ntmoy
        !!-------------------------------------------------------------------
+       REWIND( numnam_ice )                     ! Read Namelist namicedia
+       READ  ( numnam_ice  , namicedia )
+       REWIND( numnam_ice_ref )              ! Namelist namicedia in reference namelist : Ice diagnostics in ice_evolu
+       READ  ( numnam_ice_ref, namicedia, IOSTAT = ios, ERR = 901)
+    IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicedia in reference namelist', lwp )
+       REWIND( numnam_ice_cfg )              ! Namelist namicedia in configuration namelist : Ice diagnostics in ice_evolu
+       READ  ( numnam_ice_cfg, namicedia, IOSTAT = ios, ERR = 902 )
+    IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicedia in configuration namelist', lwp )
+       WRITE ( numoni, namicedia )
        IF(lwp) THEN                             ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_2/limdmp_2.F90

-                      r3635
+                      r3875
       INTEGER  :: ji, jj, jk       ! dummy loop indices
       INTEGER  :: irelax, ierror   ! error flag for allocation
+      INTEGER  ::   ios            ! Local integer output status for namelist read
+      !
       REAL(wp) :: zdmpmax, zdmpmin, zfactor, zreltim ! temporary scalar
 …
       sn_hicif = FLD_N( 'ice_damping ', -1       , 'hicif'  ,  .true.    , .true.  ,  'yearly'   ,  ''      ,  ''      )
       sn_frld  = FLD_N( 'ice_damping ', -1       , 'frld'   ,  .true.    , .true.  ,  'yearly'   ,  ''      ,  ''      )
+      REWIND( numnam_ice_ref )              ! Namelist namice_dmp in reference namelist : Ice restoring
+      READ  ( numnam_ice_ref, namice_dmp, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_dmp in reference namelist', lwp )
+      REWIND( numnam_ice )                !* read in namelist_ice namicedmp
+      READ  ( numnam_ice, namice_dmp )
+      REWIND( numnam_ice_cfg )              ! Namelist  namice_dmp in configuration namelist : Ice restoring
+      READ  ( numnam_ice_cfg, namice_dmp, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_dmp in configuration namelist', lwp )
+      WRITE ( numoni, namice_dmp )
+      !
       IF ( lwp ) THEN                     !* control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_2/limdyn_2.F90

-                      r3625
+                      r3875
       !! ** input   :   Namelist namicedyn
       !!-------------------------------------------------------------------
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       NAMELIST/namicedyn/ epsd, alpha,     &
          &                dm, nbiter, nbitdr, om, resl, cw, angvg, pstar,   &
 …
          &                nevp, telast,alphaevp
       !!-------------------------------------------------------------------
+      REWIND ( numnam_ice )                       ! Read Namelist namicedyn
+      READ   ( numnam_ice  , namicedyn )
+      REWIND( numnam_ice_ref )              ! Namelist namicedyn in reference namelist : Ice dynamics
+      READ  ( numnam_ice_ref, namicedyn, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicedyn in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namicedyn in configuration namelist : Ice dynamics
+      READ  ( numnam_ice_cfg, namicedyn, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicedyn in configuration namelist', lwp )
+      WRITE ( numoni, namicedyn )
       IF(lwp) THEN                                ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_2/limistate_2.F90

-                      r3625
+                      r3875
       INTEGER :: inum_ice
       INTEGER :: ji,jj
+      INTEGER :: ios                 ! Local integer output status for namelist read
       NAMELIST/namiceini/ ln_limini, ttest, hninn, hginn, alinn, &
          &                hnins, hgins, alins
       !!-------------------------------------------------------------------
+      !
+      REWIND ( numnam_ice )               ! Read Namelist namiceini
+      READ   ( numnam_ice , namiceini )
+      REWIND( numnam_ice_ref )              ! Namelist namiceini in reference namelist : Ice initial state
+      READ  ( numnam_ice_ref, namiceini, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceini in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namiceini in configuration namelist : Ice initial state
+      READ  ( numnam_ice_cfg, namiceini, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceini in configuration namelist', lwp )
+      WRITE ( numoni, namiceini )
+      !
       IF(lwp) THEN

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_2/limthd_2.F90

-                      r3625
+                      r3875
       !! ** input   :   Namelist namicether
       !!-------------------------------------------------------------------
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       NAMELIST/namicethd/ hmelt , hiccrit, hicmin, hiclim, amax  ,        &
          &                swiqst, sbeta  , parlat, hakspl, hibspl, exld,  &
          &                hakdif, hnzst  , thth  , parsub, alphs
       !!-------------------------------------------------------------------
+      !
+      REWIND( numnam_ice )                  ! read namelist
+      READ  ( numnam_ice , namicethd )
+      REWIND( numnam_ice_ref )              ! Namelist namicethd in reference namelist : Ice thermodynamics
+      READ  ( numnam_ice_ref, namicethd, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicethd in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namicethd in configuration namelist : Ice thermodynamics
+      READ  ( numnam_ice_cfg, namicethd, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicethd in configuration namelist', lwp )
+      WRITE ( numoni, namicethd )
       IF( lk_cpl .AND. parsub /= 0.0 )   CALL ctl_stop( 'In coupled mode, use parsub = 0. or send dqla' )
+      !

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_2/limtrp_2.F90

-                      r3764
+                      r3875
       !! ** input   :   Namelist namicetrp
       !!-------------------------------------------------------------------
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       NAMELIST/namicetrp/ bound
       !!-------------------------------------------------------------------
+      !
+      REWIND ( numnam_ice )      ! Read Namelist namicetrp
+      READ   ( numnam_ice  , namicetrp )
+      REWIND( numnam_ice_ref )              ! Namelist namicetrp in reference namelist : Ice advection
+      READ  ( numnam_ice_ref, namicetrp, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicetrp in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namicetrp in configuration namelist : Ice advection
+      READ  ( numnam_ice_cfg, namicetrp, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicetrp in configuration namelist', lwp )
+      WRITE ( numoni, namicetrp )
       IF(lwp) THEN
          WRITE(numout,*)

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_2/limwri_2.F90

-                      r3764
+                      r3875
       !!-------------------------------------------------------------------
       INTEGER ::   nf      ! ???
+      INTEGER ::   ios     ! Local integer output status for namelist read
       TYPE FIELD
          CHARACTER(len = 35) :: ztitle
 …
          CALL ctl_stop( 'STOP', 'lim_wri_init_2 : unable to allocate standard arrays' )   ;   RETURN
       ENDIF
+      REWIND ( numnam_ice )                ! Read Namelist namicewri
+      READ   ( numnam_ice  , namiceout )
+      REWIND( numnam_ice_ref )              ! Namelist namiceout in reference namelist : Ice outputs
+      READ  ( numnam_ice_ref, namiceout, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceout in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namiceout in configuration namelist : Ice outputs
+      READ  ( numnam_ice_cfg, namiceout, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceout in configuration namelist', lwp )
+      WRITE ( numoni, namiceout )
       zfield( 1) = field_1

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_3/iceini.F90

-                      r3625
+                      r3875
          &     'use more ocean levels or less ice/snow layers/categories.' )
+      !                                ! Open the namelist file
+      CALL ctl_opn( numnam_ice, 'namelist_ice', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
+                                       ! Open the reference and configuration namelist files and namelist output file
+      CALL ctl_opn( numnam_ice_ref, 'namelist_ice_ref',    'OLD',     'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
+      CALL ctl_opn( numnam_ice_cfg, 'namelist_ice_cfg',    'OLD',     'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
+      CALL ctl_opn( numoni,         'output.namelist.ice', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp )
+      !
       CALL ice_run                     ! set some ice run parameters
 …
       !!-------------------------------------------------------------------
       NAMELIST/namicerun/ cn_icerst_in, cn_icerst_out, ln_limdyn, acrit, hsndif, hicdif, cai, cao, ln_nicep
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       !!-------------------------------------------------------------------
+      !
+      REWIND( numnam_ice )                ! Read Namelist namicerun
+      READ  ( numnam_ice , namicerun )
+      REWIND( numnam_ice_ref )              ! Namelist namicerun in reference namelist : Parameters for ice
+      READ  ( numnam_ice_ref, namicerun, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicerun in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namicerun in configuration namelist : Parameters for ice
+      READ  ( numnam_ice_cfg, namicerun, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicerun in configuration namelist', lwp )
+      WRITE ( numoni, namicerun )
+      !
       IF( lk_mpp .AND. ln_nicep ) THEN

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_3/limdia.F90

-                      r3625
+                      r3875
       INTEGER  ::   jv    ! dummy loop indice
       INTEGER  ::   ierr, ntot , ndeb , irecl   ! local integers
+      INTEGER  ::   ios                         ! Local integer output status for namelist read
       REAL(wp) ::   zxx0, zxx1    ! local scalars
       CHARACTER(len=jpchinf) ::   titinf
 …
       !!-------------------------------------------------------------------
+      !
+      REWIND( numnam_ice )             ! read namicedia namelist
+      READ  ( numnam_ice, namicedia )
+      REWIND( numnam_ice_ref )              ! Namelist namicedia in reference namelist : Ice diagnostics in ice_evolu
+      READ  ( numnam_ice_ref, namicedia, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicedia in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namicedia in configuration namelist : Ice diagnostics in ice_evolu
+      READ  ( numnam_ice_cfg, namicedia, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicedia in configuration namelist', lwp )
+      WRITE ( numoni, namicedia )
+      !
       IF(lwp) THEN                     ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_3/limdyn.F90

-                      r3625
+                      r3875
       !! ** input   :   Namelist namicedyn
       !!-------------------------------------------------------------------
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       NAMELIST/namicedyn/ epsd, alpha,     &
          &                dm, nbiter, nbitdr, om, resl, cw, angvg, pstar,   &
 …
       !!-------------------------------------------------------------------
+      REWIND( numnam_ice )                ! Read Namelist namicedyn
+      READ  ( numnam_ice  , namicedyn )
+      REWIND( numnam_ice_ref )              ! Namelist namicedyn in reference namelist : Ice dynamics
+      READ  ( numnam_ice_ref, namicedyn, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicedyn in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namicedyn in configuration namelist : Ice dynamics
+      READ  ( numnam_ice_cfg, namicedyn, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicedyn in configuration namelist', lwp )
+      WRITE ( numoni, namicedyn )
       IF(lwp) THEN                        ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_3/limistate.F90

-                      r3625
+                      r3875
       !! ** input  :   namelist namiceini
       !!-----------------------------------------------------------------------------
+      INTEGER :: ios                 ! Local integer output status for namelist read
       NAMELIST/namiceini/ ttest, hninn, hginn_u, aginn_u, hginn_d, aginn_d, hnins,   &
          &                hgins_u, agins_u, hgins_d, agins_d, sinn, sins
       !!-----------------------------------------------------------------------------
+      !
+      REWIND ( numnam_ice )               ! Read Namelist namiceini
+      READ   ( numnam_ice , namiceini )
+      REWIND( numnam_ice_ref )              ! Namelist namiceini in reference namelist : Ice initial state
+      READ  ( numnam_ice_ref, namiceini, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceini in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namiceini in configuration namelist : Ice initial state
+      READ  ( numnam_ice_cfg, namiceini, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceini in configuration namelist', lwp )
+      WRITE ( numoni, namiceini )
+      !
       IF(lwp) THEN                        ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90

-                      r3764
+                      r3875
       !! ** input   :   Namelist namiceitdme
       !!-------------------------------------------------------------------
+      INTEGER :: ios                 ! Local integer output status for namelist read
       NAMELIST/namiceitdme/ ridge_scheme_swi, Cs, Cf, fsnowrdg, fsnowrft,&
          Gstar, astar,                                &
 …
       !!-------------------------------------------------------------------
+      !
+      REWIND( numnam_ice )                   ! read namiceitdme namelist
+      READ  ( numnam_ice , namiceitdme)
+      REWIND( numnam_ice_ref )              ! Namelist namicetdme in reference namelist : Ice mechanical ice redistribution
+      READ  ( numnam_ice_ref, namiceitdme, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceitdme in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namiceitdme in configuration namelist : Ice mechanical ice redistribution
+      READ  ( numnam_ice_cfg, namiceitdme, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceitdme in configuration namelist', lwp )
+      WRITE ( numoni, namiceitdme )
+      !
       IF (lwp) THEN                          ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90

-                      r3625
+                      r3875
       !! ** input   :   Namelist namicether
       !!-------------------------------------------------------------------
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       NAMELIST/namicethd/ hmelt , hiccrit, fraz_swi, maxfrazb, vfrazb, Cfrazb,   &
          &                hicmin, hiclim, amax  ,                                &
 …
       ENDIF
+      !
+      REWIND( numnam_ice )                  ! read Namelist numnam_ice
+      READ  ( numnam_ice , namicethd )
+      REWIND( numnam_ice_ref )              ! Namelist namicethd in reference namelist : Ice thermodynamics
+      READ  ( numnam_ice_ref, namicethd, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicethd in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namicethd in configuration namelist : Ice thermodynamics
+      READ  ( numnam_ice_cfg, namicethd, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicethd in configuration namelist', lwp )
+      WRITE ( numoni, namicethd )
+      !
       IF(lwp) THEN                          ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_3/limthd_sal.F90

-                      r3625
+                      r3875
       !! ** input   :   Namelist namicesal
       !!-------------------------------------------------------------------
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       NAMELIST/namicesal/ num_sal, bulk_sal, sal_G, time_G, sal_F, time_F,   &
          &                s_i_max, s_i_min, s_i_0, s_i_1
       !!-------------------------------------------------------------------
+      !
+      REWIND( numnam_ice )                   ! Read Namelist namicesal
+      READ  ( numnam_ice  , namicesal )
+      REWIND( numnam_ice_ref )              ! Namelist namicesal in reference namelist : Ice salinity
+      READ  ( numnam_ice_ref, namicesal, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicesal in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namicesal in configuration namelist : Ice salinity
+      READ  ( numnam_ice_cfg, namicesal, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namicesal in configuration namelist', lwp )
+      WRITE ( numoni, namicesal )
+      !
       IF(lwp) THEN                           ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/LIM_SRC_3/limwri.F90

-                      r3853
+                      r3875
       !!-------------------------------------------------------------------
       INTEGER ::   nf      ! ???
+      INTEGER ::   ios     ! Local integer output status for namelist read
       TYPE FIELD
 …
          field_37, add_diag_swi
       !!-------------------------------------------------------------------
+      REWIND( numnam_ice )                ! Read Namelist namicewri
+      READ  ( numnam_ice  , namiceout )
+      REWIND( numnam_ice_ref )              ! Namelist namiceout in reference namelist : Ice outputs
+      READ  ( numnam_ice_ref, namiceout, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceout in reference namelist', lwp )
+      REWIND( numnam_ice_cfg )              ! Namelist namiceout in configuration namelist : Ice outputs
+      READ  ( numnam_ice_cfg, namiceout, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceout in configuration namelist', lwp )
+      WRITE ( numoni, namiceout )
       zfield(1)  = field_1

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/NST_SRC/agrif_user.F90

-                      r3680
+                      r3875
    IMPLICIT NONE
+   !
+   INTEGER  ::   ios                 ! Local integer output status for namelist read
    NAMELIST/namagrif/ nn_cln_update, rn_sponge_tra, rn_sponge_dyn, ln_spc_dyn
    !!----------------------------------------------------------------------
+   !
+   REWIND( numnam )                ! Read namagrif namelist
+   READ  ( numnam, namagrif )
+      REWIND( numnam_ref )              ! Namelist namagrif in reference namelist : AGRIF zoom
+      READ  ( numnam_ref, namagrif, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namagrif in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namagrif in configuration namelist : AGRIF zoom
+      READ  ( numnam_cfg, namagrif, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namagrif in configuration namelist', lwp )
+      WRITE ( numond, namagrif )
+   !
    IF(lwp) THEN                    ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OFF_SRC/domain.F90

-                      r3632
+                      r3875
       !!----------------------------------------------------------------------
       USE ioipsl
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       NAMELIST/namrun/ nn_no   , cn_exp    , cn_ocerst_in, cn_ocerst_out, ln_rstart , nn_rstctl,   &
          &             nn_it000, nn_itend  , nn_date0    , nn_leapy     , nn_istate , nn_stock ,   &
 …
       !!----------------------------------------------------------------------
+      REWIND( numnam )              ! Namelist namrun : parameters of the run
+      READ  ( numnam, namrun )
+      REWIND( numnam_ref )              ! Namelist namrun in reference namelist : Parameters of the run
+      READ  ( numnam_ref, namrun, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namrun in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namrun in configuration namelist : Parameters of the run
+      READ  ( numnam_cfg, namrun, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namrun in configuration namelist', lwp )
+      WRITE ( numond, namrun )
+      !
       IF(lwp) THEN                  ! control print
 …
 #endif
+      REWIND( numnam )             ! Domain
+      READ  ( numnam, namdom )
+      REWIND( numnam_ref )              ! Namelist namdom in reference namelist : space & time domain (bathymetry, mesh, timestep)
+      READ  ( numnam_ref, namdom, IOSTAT = ios, ERR = 903)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdom in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namdom in configuration namelist : space & time domain (bathymetry, mesh, timestep)
+      READ  ( numnam_cfg, namdom, IOSTAT = ios, ERR = 904 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdom in configuration namelist', lwp )
+      WRITE ( numond, namdom )
       IF(lwp) THEN
 …
       rdth      = rn_rdth
+      REWIND( numnam )             ! Namelist cross land advection
+      READ  ( numnam, namcla )
+      REWIND( numnam_ref )              ! Namelist namcla in reference namelist : Cross land advection
+      READ  ( numnam_ref, namcla, IOSTAT = ios, ERR = 905)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namcla in configuration namelist : Cross land advection
+      READ  ( numnam_cfg, namcla, IOSTAT = ios, ERR = 906 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in configuration namelist', lwp )
+      WRITE( numond, namcla )
       IF(lwp) THEN
          WRITE(numout,*)
 …
 #if defined key_netcdf4
       !                             ! NetCDF 4 case   ("key_netcdf4" defined)
+      REWIND( numnam )                    ! Namelist namnc4 : netcdf4 chunking parameters
+      READ  ( numnam, namnc4 )
+      REWIND( numnam_ref )              ! Namelist namnc4 in reference namelist : NETCDF
+      READ  ( numnam_ref, namnc4, IOSTAT = ios, ERR = 907)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namnc4 in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namnc4 in configuration namelist : NETCDF
+      READ  ( numnam_cfg, namnc4, IOSTAT = ios, ERR = 908 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namnc4 in configuration namelist', lwp )
+      WRITE( numond, namnc4 )
       IF(lwp) THEN                        ! control print
          WRITE(numout,*)
 …
       !!----------------------------------------------------------------------
       INTEGER ::   ioptio = 0   ! temporary integer
+      INTEGER ::   ios
       !!
       NAMELIST/namzgr/ ln_zco, ln_zps, ln_sco
       !!----------------------------------------------------------------------
+      REWIND ( numnam )                ! Read Namelist namzgr : vertical coordinate'
+      READ   ( numnam, namzgr )
+      REWIND( numnam_ref )              ! Namelist namzgr in reference namelist : Vertical coordinate
+      READ  ( numnam_ref, namzgr, IOSTAT = ios, ERR = 901 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzgr in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namzgr in configuration namelist : Vertical coordinate
+      READ  ( numnam_cfg, namzgr, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzgr in configuration namelist', lwp )
+      WRITE ( numond, namzgr )
       IF(lwp) THEN                     ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90

-                      r3827
+                      r3875
       INTEGER  :: jfld                               ! dummy loop arguments
       INTEGER  :: inum, idv, idimv                   ! local integer
+      INTEGER  :: ios                                ! Local integer output status for namelist read
       !!
       CHARACTER(len=100)            ::  cn_dir   !   Root directory for location of core files
 …
       sn_eiw  = FLD_N( 'dyna_grid_W' ,    120    , 'voveaeiw' ,  .true.    , .true. ,   'yearly'  , ''       , ''         )
+      !
+      REWIND( numnam )                          ! read in namlist namdta_dyn
+      READ  ( numnam, namdta_dyn )
+      REWIND( numnam_ref )              ! Namelist namdta_dyn in reference namelist : Offline: init. of dynamical data
+      READ  ( numnam_ref, namdta_dyn, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdta_dyn in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namdta_dyn in configuration namelist : Offline: init. of dynamical data
+      READ  ( numnam_cfg, namdta_dyn, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdta_dyn in configuration namelist', lwp )
+      WRITE ( numond, namdta_dyn )
       !                                         ! store namelist information in an array
       !                                         ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OFF_SRC/nemogcm.F90

-                      r3827
+                      r3875
       INTEGER ::   ji            ! dummy loop indices
       INTEGER ::   ilocal_comm   ! local integer
+      INTEGER ::   ios
       CHARACTER(len=80), DIMENSION(16) ::   cltxt
       !!
 …
       cltxt = ''
+      !
+      !                             ! open Namelist file
+      CALL ctl_opn( numnam, 'namelist', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. )
+      !
+      READ( numnam, namctl )        ! Namelist namctl : Control prints & Benchmark
+      !                             ! Open reference namelist and configuration namelist files
+      CALL ctl_opn( numnam_ref, 'namelist_ref', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. )
+      CALL ctl_opn( numnam_cfg, 'namelist_cfg', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. )
+      CALL ctl_opn( numond, 'output.namelist.dyn', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea )
+      !
+      !
+      REWIND( numnam_ref )              ! Namelist namctl in reference namelist : Control prints & Benchmark
+      READ  ( numnam_ref, namctl, IOSTAT = ios, ERR = 901 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namctl in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namctl in confguration namelist : Control prints & Benchmark
+      READ  ( numnam_cfg, namctl, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namctl in configuration namelist', lwp )
+      WRITE( numond, namctl )
+      !
       !                             !--------------------------------------------!
 …
 #if defined key_iomput
          CALL  xios_initialize( "nemo",return_comm=ilocal_comm )
       narea = mynode( cltxt, numnam, nstop, ilocal_comm )   ! Nodes selection
+      narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop, ilocal_comm )   ! Nodes selection
 #else
       ilocal_comm = 0
       narea = mynode( cltxt, numnam, nstop )                 ! Nodes selection (control print return in cltxt)
+      narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop )                ! Nodes selection (control print return in cltxt)
 #endif
 …
+      !
       IF( numstp     /= -1 )   CLOSE( numstp     )   ! time-step file
+      IF( numnam     /= -1 )   CLOSE( numnam     )   ! oce namelist
+      IF( numnam_ref /= -1 )   CLOSE( numnam_ref )   ! oce reference namelist
+      IF( numnam_cfg /= -1 )   CLOSE( numnam_cfg )   ! oce configuration namelist
       IF( numout     /=  6 )   CLOSE( numout     )   ! standard model output file
       numout = 6                                     ! redefine numout in case it is used after this point...

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

-                      r3785
+                      r3875
       INTEGER :: iitiaustr_date  ! Date YYYYMMDD of IAU interval start time step
       INTEGER :: iitiaufin_date  ! Date YYYYMMDD of IAU interval final time step
+      INTEGER :: ios             ! Local integer output status for namelist read
       REAL(wp) :: znorm        ! Normalization factor for IAU weights
 …
       niaufn    = 0
+      REWIND ( numnam )
+      READ   ( numnam, nam_asminc )
+      REWIND( numnam_ref )              ! Namelist nam_asminc in reference namelist : Assimilation increment
+      READ  ( numnam_ref, nam_asminc, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_asminc in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist nam_asminc in configuration namelist : Assimilation increment
+      READ  ( numnam_cfg, nam_asminc, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_asminc in configuration namelist', lwp )
+      WRITE ( numond, nam_asminc )
       ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90

-                      r3851
+                      r3875
       !!
       INTEGER     ::  ib_bdy, jfld, jstart, jend, ierror  ! local indices
+      INTEGER      ::   ios                               ! Local integer output status for namelist read
       !!
       CHARACTER(len=100)                     ::   cn_dir        ! Root directory for location of data files
 …
       ! Read namelists
       ! --------------
       REWIND(numnam)
+      REWIND(numnam_cfg)
       jfld = 0
       DO ib_bdy = 1, nb_bdy
 …
 #endif
+            ! Read configuration namelist only to avoid unsuccessful overwrite
             ! Important NOT to rewind here.
+            READ( numnam, nambdy_dta )
+!!          REWIND( numnam_ref )              ! Namelist nambdy_dta in reference namelist :
+!!          READ  ( numnam_ref, nambdy_dta, IOSTAT = ios, ERR = 901)
+!!901       IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_dta in reference namelist', lwp )
+!!          REWIND( numnam_cfg )              ! Namelist nambdy_dta in configuration namelist :
+            READ  ( numnam_cfg, nambdy_dta, IOSTAT = ios, ERR = 902 )
+         IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_dta in configuration namelist', lwp )
+            WRITE ( numond, nambdy_dta )
             cn_dir_array(ib_bdy) = cn_dir

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

-                      r3703
+                      r3875
       !!
       NAMELIST/nambdy_index/ ctypebdy, nbdyind, nbdybeg, nbdyend
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       !!----------------------------------------------------------------------
 …
       nn_rimwidth(:)    = -1  ! uninitialised flag
+      REWIND( numnam )
+      READ  ( numnam, nambdy )
+      REWIND( numnam_ref )              ! Namelist nambdy in reference namelist :Unstructured open boundaries
+      READ  ( numnam_ref, nambdy, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist nambdy in configuration namelist :Unstructured open boundaries
+      READ  ( numnam_cfg, nambdy, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy in configuration namelist', lwp )
+      WRITE ( numond, nambdy )
       ! -----------------------------------------
 …
       ! Work out global dimensions of boundary data
       ! ---------------------------------------------
+      REWIND( numnam )
+      REWIND( numnam_cfg )
+      !!----------------------------------------------------------------------
       nblendta(:,:) = 0
 …
             icount = icount + 1
             ! No REWIND here because may need to read more than one nambdy_index namelist.
+            READ  ( numnam, nambdy_index )
+            ! Read only namelist_cfg to avoid unseccessfull overwrite
+!!          REWIND( numnam_ref )              ! Namelist nambdy_index in reference namelist : Open boundaries indexes
+!!          READ  ( numnam_ref, namrun, IOSTAT = ios, ERR = 903)
+!!903       IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_index in reference namelist', lwp )
+!!          REWIND( numnam_cfg )              ! Namelist nambdy_index in configuration namelist : Open boundaries indexes
+            READ  ( numnam_cfg, nambdy_index, IOSTAT = ios, ERR = 904 )
+         IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_index in configuration namelist', lwp )
+            WRITE ( numond, nambdy_index )
             SELECT CASE ( TRIM(ctypebdy) )

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/BDY/bdytides.F90

-                      r3651
+                      r3875
       INTEGER, DIMENSION(3)                     ::   ilen0       !: length of boundary data (from OBC arrays)
       INTEGER, POINTER, DIMENSION(:)            ::   nblen, nblenrim     ! short cuts
+      INTEGER                                   ::   ios                 ! Local integer output status for namelist read
       CHARACTER(len=80)                         ::   clfile              !: full file name for tidal input file
       REAL(wp),ALLOCATABLE, DIMENSION(:,:,:)    ::   dta_read            !: work space to read in tidal harmonics data
 …
       ln_bdytide_conj  = .FALSE.
+      REWIND(numnam)
+      REWIND(numnam_cfg)
       DO ib_bdy = 1, nb_bdy
          IF( nn_dyn2d_dta(ib_bdy) .ge. 2 ) THEN
 …
             ! Don't REWIND here - may need to read more than one of these namelists.
+            READ  ( numnam, nambdy_tide )
+            ! Read only configuration namelist to avoid unsecessful overwrite
+!!          REWIND( numnam_ref )              ! Namelist nambdy_tide in reference namelist :
+!!          READ  ( numnam_ref, nambdy_tide, IOSTAT = ios, ERR = 901)
+!!    901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_tide in reference namelist', lwp )
+!!          REWIND( numnam_cfg )              ! Namelist nambdy_tide in configuration namelist :
+            READ  ( numnam_cfg, nambdy_tide, IOSTAT = ios, ERR = 902 )
+         IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy_tide in configuration namelist', lwp )
+            WRITE ( numond, nambdy_tide )
             !                                               ! Parameter control and print
             IF(lwp) WRITE(numout,*) '  '

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90

-                      r3680
+                      r3875
      !!---------------------------------------------------------------------
      NAMELIST/namdct/nn_dct,nn_dctwri,nn_secdebug
+     INTEGER  ::   ios                 ! Local integer output status for namelist read
      IF( nn_timing == 1 )   CALL timing_start('dia_dct_init')
+     !read namelist
+     REWIND( numnam )
+     READ  ( numnam, namdct )
+     REWIND( numnam_ref )              ! Namelist namdct in reference namelist : Diagnostic: transport through sections
+     READ  ( numnam_ref, namdct, IOSTAT = ios, ERR = 901)
+  IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdct in reference namelist', lwp )
+     REWIND( numnam_cfg )              ! Namelist namdct in configuration namelist : Diagnostic: transport through sections
+     READ  ( numnam_cfg, namdct, IOSTAT = ios, ERR = 902 )
+  IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdct in configuration namelist', lwp )
+     WRITE ( numond, namdct )
      IF( lwp ) THEN

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DIA/diaharm.F90

-                      r3294
+                      r3875
       !! * Local declarations
       INTEGER :: jh, nhan, jk, ji
+      INTEGER ::   ios                 ! Local integer output status for namelist read
       NAMELIST/nam_diaharm/ nit000_han, nitend_han, nstep_han, tname
 …
       tname(:)=''
+      !
+      ! Read Namelist nam_diaharm
+      REWIND ( numnam )
+      READ   ( numnam, nam_diaharm )
+      REWIND( numnam_ref )              ! Namelist nam_diaharm in reference namelist : Tidal harmonic analysis
+      READ  ( numnam_ref, nam_diaharm, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_diaharm in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist nam_diaharm in configuration namelist : Tidal harmonic analysis
+      READ  ( numnam_cfg, nam_diaharm, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_diaharm in configuration namelist', lwp )
+      WRITE ( numond, nam_diaharm )
+      !
       IF(lwp) THEN

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90

-                      r3625
+                      r3875
       INTEGER            ::   jk       ! dummy loop indice
       INTEGER            ::   ierror   ! local integer
+      INTEGER            ::   ios      ! Local integer output status for namelist read
       !!
       NAMELIST/namhsb/ ln_diahsb
       !!----------------------------------------------------------------------
+      !
+      REWIND ( numnam )              ! Read Namelist namhsb
+      READ   ( numnam, namhsb )
+      REWIND( numnam_ref )              ! Namelist namhsb in reference namelist : Heat & salt budget
+      READ  ( numnam_ref, namhsb, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namhsb in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namhsb in configuration namelist : Heat & salt budget
+      READ  ( numnam_cfg, namhsb, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namhsb in configuration namelist', lwp )
+      WRITE ( numond, namhsb )
+      !
       IF(lwp) THEN                   ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DIA/diaptr.F90

-                      r3764
+                      r3875
       INTEGER ::   jn           ! dummy loop indices
       INTEGER ::   inum, ierr   ! local integers
+      INTEGER ::   ios          ! Local integer output status for namelist read
 #if defined key_mpp_mpi
       INTEGER, DIMENSION(1) :: iglo, iloc, iabsf, iabsl, ihals, ihale, idid
 …
       !!----------------------------------------------------------------------
+      REWIND( numnam )                 ! Read Namelist namptr : poleward transport parameters
+      READ  ( numnam, namptr )
+      REWIND( numnam_ref )              ! Namelist namptr in reference namelist : Poleward transport
+      READ  ( numnam_ref, namptr, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namptr in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namptr in configuration namelist : Poleward transport
+      READ  ( numnam_cfg, namptr, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namptr in configuration namelist', lwp )
+      WRITE ( numond, namptr )
       IF(lwp) THEN                     ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

-                      r3764
+                      r3875
       NAMELIST/namnc4/ nn_nchunks_i, nn_nchunks_j, nn_nchunks_k, ln_nc4zip
 #endif
+      !!----------------------------------------------------------------------
+      REWIND( numnam )              ! Namelist namrun : parameters of the run
+      READ  ( numnam, namrun )
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
+      !!----------------------------------------------------------------------
+      REWIND( numnam_ref )              ! Namelist namrun in reference namelist : Parameters of the run
+      READ  ( numnam_ref, namrun, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namrun in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namrun in configuration namelist : Parameters of the run
+      READ  ( numnam_cfg, namrun, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namrun in configuration namelist', lwp )
+      WRITE ( numond, namrun )
+      !
       IF(lwp) THEN                  ! control print
 …
 #endif
+      REWIND( numnam )              ! Namelist namdom : space & time domain (bathymetry, mesh, timestep)
+      READ  ( numnam, namdom )
+      REWIND( numnam_ref )              ! Namelist namdom in reference namelist : space & time domain (bathymetry, mesh, timestep)
+      READ  ( numnam_ref, namdom, IOSTAT = ios, ERR = 903)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdom in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namdom in configuration namelist : space & time domain (bathymetry, mesh, timestep)
+      READ  ( numnam_cfg, namdom, IOSTAT = ios, ERR = 904 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdom in configuration namelist', lwp )
+      WRITE ( numond, namdom )
       IF(lwp) THEN
 …
       rdth      = rn_rdth
+      REWIND( numnam )              ! Namelist cross land advection
+      READ  ( numnam, namcla )
+      REWIND( numnam_ref )              ! Namelist namcla in reference namelist : Cross land advection
+      READ  ( numnam_ref, namcla, IOSTAT = ios, ERR = 905)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namcla in configuration namelist : Cross land advection
+      READ  ( numnam_cfg, namcla, IOSTAT = ios, ERR = 906 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in configuration namelist', lwp )
+      WRITE( numond, namcla )
       IF(lwp) THEN
          WRITE(numout,*)
 …
 #if defined key_netcdf4
       !                             ! NetCDF 4 case   ("key_netcdf4" defined)
+      REWIND( numnam )                    ! Namelist namnc4 : netcdf4 chunking parameters
+      READ  ( numnam, namnc4 )
+      REWIND( numnam_ref )              ! Namelist namnc4 in reference namelist : NETCDF
+      READ  ( numnam_ref, namnc4, IOSTAT = ios, ERR = 907)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namnc4 in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namnc4 in configuration namelist : NETCDF
+      READ  ( numnam_cfg, namnc4, IOSTAT = ios, ERR = 908 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namnc4 in configuration namelist', lwp )
+      WRITE( numond, namnc4 )
       IF(lwp) THEN                        ! control print
          WRITE(numout,*)

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90

-                      r3294
+                      r3875
       INTEGER  ::   iif, iil, ii0, ii1, ii   ! local integers
       INTEGER  ::   ijf, ijl, ij0, ij1       !   -       -
+      INTEGER  ::   ios
       INTEGER , POINTER, DIMENSION(:,:) ::  imsk
       REAL(wp), POINTER, DIMENSION(:,:) ::  zwf
 …
       CALL wrk_alloc( jpi, jpj, zwf  )
+      !
+      REWIND( numnam )              ! Namelist namlbc : lateral momentum boundary condition
+      READ  ( numnam, namlbc )
+      REWIND( numnam_ref )              ! Namelist namlbc in reference namelist : Lateral momentum boundary condition
+      READ  ( numnam_ref, namlbc, IOSTAT = ios, ERR = 901 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlbc in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namlbc in configuration namelist : Lateral momentum boundary condition
+      READ  ( numnam_cfg, namlbc, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlbc in configuration namelist', lwp )
+      WRITE ( numond, namlbc )
       IF(lwp) THEN                  ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

-                      r3764
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER ::   ioptio, ibat   ! local integer
+      INTEGER ::   ios
+      !
       NAMELIST/namzgr/ ln_zco, ln_zps, ln_sco
 …
       IF( nn_timing == 1 )   CALL timing_start('dom_zgr')
+      !
+      REWIND( numnam )                 ! Read Namelist namzgr : vertical coordinate'
+      READ  ( numnam, namzgr )
+      REWIND( numnam_ref )              ! Namelist namzgr in reference namelist : Vertical coordinate
+      READ  ( numnam_ref, namzgr, IOSTAT = ios, ERR = 901 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzgr in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namzgr in configuration namelist : Vertical coordinate
+      READ  ( numnam_cfg, namzgr, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzgr in configuration namelist', lwp )
+      WRITE ( numond, namzgr )
       IF(lwp) THEN                     ! Control print
 …
       INTEGER  ::   ji, jj, jk, jl           ! dummy loop argument
       INTEGER  ::   iip1, ijp1, iim1, ijm1   ! temporary integers
+      INTEGER  ::   ios                      ! Local integer output status for namelist read
       REAL(wp) ::   zrmax, ztaper   ! temporary scalars
+      !
 …
       CALL wrk_alloc( jpi, jpj,      zenv, ztmp, zmsk, zri, zrj, zhbat                           )
+      !
+      REWIND( numnam )                       ! Read Namelist namzgr_sco : sigma-stretching parameters
+      READ  ( numnam, namzgr_sco )
+      REWIND( numnam_ref )              ! Namelist namzgr_sco in reference namelist : Sigma-stretching parameters
+      READ  ( numnam_ref, namzgr_sco, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzgr_sco in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namzgr_sco in configuration namelist : Sigma-stretching parameters
+      READ  ( numnam_cfg, namzgr_sco, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzgr_sco in configuration namelist', lwp )
+      WRITE ( numond, namzgr_sco )
       IF(lwp) THEN                           ! control print
          WRITE(numout,*)
          WRITE(numout,*) 'dom:zgr_sco : s-coordinate or hybrid z-s-coordinate'
+         WRITE(numout,*) 'domzgr_sco : s-coordinate or hybrid z-s-coordinate'
          WRITE(numout,*) '~~~~~~~~~~~'
          WRITE(numout,*) '   Namelist namzgr_sco'

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DOM/dtatsd.F90

-                      r3294
+                      r3875
       !!
       NAMELIST/namtsd/   ln_tsd_init, ln_tsd_tradmp, cn_dir, sn_tem, sn_sal
+      INTEGER  ::   ios
       !!----------------------------------------------------------------------
+      !
 …
       sn_tem = FLD_N( 'temperature',  -1.  , 'votemper',  .false.   , .true.  ,  'monthly'  , ''       , ''       )
       sn_sal = FLD_N( 'salinity'   ,  -1.  , 'vosaline',  .false.   , .true.  ,  'monthly'  , ''       , ''       )
+      REWIND( numnam )              ! read in namlist namdta_tsd
+      READ  ( numnam, namtsd )
+      REWIND( numnam_ref )              ! Namelist namtsd in reference namelist :
+      READ  ( numnam_ref, namtsd, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtsd in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namtsd in configuration namelist : Parameters of the run
+      READ  ( numnam_cfg, namtsd, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtsd in configuration namelist', lwp )
+      WRITE ( numond, namtsd )
       IF( PRESENT( ld_tradmp ) )   ln_tsd_tradmp = .TRUE.     ! forces the initialization when tradmp is used

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DYN/dynadv.F90

-                      r3294
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER ::   ioptio
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       !!
       NAMELIST/namdyn_adv/ ln_dynadv_vec, ln_dynadv_cen2 , ln_dynadv_ubs
       !!----------------------------------------------------------------------
+      REWIND ( numnam )               ! Read Namelist namdyn_adv : momentum advection scheme
+      READ   ( numnam, namdyn_adv )
+      REWIND( numnam_ref )              ! Namelist namdyn_adv in reference namelist : Momentum advection scheme
+      READ  ( numnam_ref, namdyn_adv, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_adv in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namdyn_adv in configuration namelist : Momentum advection scheme
+      READ  ( numnam_cfg, namdyn_adv, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_adv in configuration namelist', lwp )
+      WRITE ( numond, namdyn_adv )
       IF(lwp) THEN                    ! Namelist print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DYN/dynhpg.F90

-                      r3764
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER ::   ioptio = 0      ! temporary integer
+      INTEGER ::   ios             ! Local integer output status for namelist read
       !!
       NAMELIST/namdyn_hpg/ ln_hpg_zco, ln_hpg_zps, ln_hpg_sco,     &
 …
       !!----------------------------------------------------------------------
+      !
+      REWIND( numnam )               ! Read Namelist namdyn_hpg
+      READ  ( numnam, namdyn_hpg )
+      REWIND( numnam_ref )              ! Namelist namdyn_hpg in reference namelist : Hydrostatic pressure gradient
+      READ  ( numnam_ref, namdyn_hpg, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_hpg in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namdyn_hpg in configuration namelist : Hydrostatic pressure gradient
+      READ  ( numnam_cfg, namdyn_hpg, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_hpg in configuration namelist', lwp )
+      WRITE ( numond, namdyn_hpg )
+      !
       IF(lwp) THEN                   ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DYN/dynnept.F90

-                      r3723
+                      r3875
       NAMELIST/namdyn_nept/ ln_neptsimp, ln_smooth_neptvel, rn_tslse, rn_tslsp,      &
                             ln_neptramp, rn_htrmin, rn_htrmax
+      INTEGER  ::   ios
       !!----------------------------------------------------------------------
       ! Define the (simplified) Neptune parameters
       ! ==========================================
+      REWIND( numnam )                  ! Read Namelist namdyn_nept:  Simplified Neptune
+      READ  ( numnam, namdyn_nept )
+      REWIND( numnam_ref )              ! Namelist namdyn_nept in reference namelist : Simplified Neptune
+      READ  ( numnam_ref, namdyn_nept, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_nept in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namdyn_nept in reference namelist : Simplified Neptune
+      READ  ( numnam_cfg, namdyn_nept, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_nept in configuration namelist', lwp )
+      WRITE ( numond, namdyn_nept )
       IF(lwp) THEN                      ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/DYN/dynvor.F90

-                      r3802
+                      r3875
       INTEGER ::   ioptio          ! local integer
       INTEGER ::   ji, jj, jk      ! dummy loop indices
+      INTEGER ::   ios             ! Local integer output status for namelist read
       !!
       NAMELIST/namdyn_vor/ ln_dynvor_ens, ln_dynvor_ene, ln_dynvor_mix, ln_dynvor_een
       !!----------------------------------------------------------------------
+      REWIND ( numnam )               ! Read Namelist namdyn_vor : Vorticity scheme options
+      READ   ( numnam, namdyn_vor )
+      REWIND( numnam_ref )              ! Namelist namdyn_vor in reference namelist : Vorticity scheme options
+      READ  ( numnam_ref, namdyn_vor, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_vor in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namdyn_vor in configuration namelist : Vorticity scheme options
+      READ  ( numnam_cfg, namdyn_vor, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_vor in configuration namelist', lwp )
+      WRITE ( numond, namdyn_vor )
       IF(lwp) THEN                    ! Namelist print
 …
          WRITE(numout,*) 'dyn_vor_init : vorticity term : read namelist and control the consistency'
          WRITE(numout,*) '~~~~~~~~~~~~'
          WRITE(numout,*) '        Namelist namdyn_vor : oice of the vorticity term scheme'
+         WRITE(numout,*) '        Namelist namdyn_vor : choice of the vorticity term scheme'
          WRITE(numout,*) '           energy    conserving scheme                ln_dynvor_ene = ', ln_dynvor_ene
          WRITE(numout,*) '           enstrophy conserving scheme                ln_dynvor_ens = ', ln_dynvor_ens

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/FLO/floats.F90

-                      r3294
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER :: jfl
+      INTEGER :: ios                 ! Local integer output status for namelist read
+      !
       NAMELIST/namflo/ jpnfl, jpnnewflo, ln_rstflo, nn_writefl, nn_stockfl, ln_argo, ln_flork4, ln_ariane, ln_flo_ascii
 …
       IF(lwp) WRITE(numout,*) '~~~~~~~'
+      REWIND( numnam )              ! Namelist namflo : floats
+      READ  ( numnam, namflo )
+      REWIND( numnam_ref )              ! Namelist namflo in reference namelist : Floats
+      READ  ( numnam_ref, namflo, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namflo in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namflo in configuration namelist : Floats
+      READ  ( numnam_cfg, namflo, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namflo in configuration namelist', lwp )
+      WRITE ( numond, namflo )
+      !
       IF(lwp) THEN                  ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/ICB/icbini.F90

-                      r3785
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER  ::   jn      ! dummy loop indices
+      INTEGER  ::   ios     ! Local integer output status for namelist read
       REAL(wp) ::   zfact   ! local scalar
+      !
 …
       sn_icb = FLD_N( 'calving' ,    -1     , 'calving'  ,  .TRUE.    , .TRUE. ,   'yearly'  , ''       , ''         )
+      REWIND( numnam )              ! Namelist namrun : iceberg parameters
+      READ  ( numnam, namberg )
+      REWIND( numnam_ref )              ! Namelist namberg in reference namelist : Iceberg parameters
+      READ  ( numnam_ref, namberg, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namberg in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namberg in configuration namelist : Iceberg parameters
+      READ  ( numnam_cfg, namberg, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namberg in configuration namelist', lwp )
+      WRITE ( numond, namberg )
       IF( .NOT. ln_icebergs ) THEN   ! no icebergs
          IF(lwp) THEN
             WRITE(numout,*)
             WRITE(numout,*) 'icb_nam :   ln_icebergs = F , NO icebergs used'
+            WRITE(numout,*) 'icbini :   Namelist namberg ln_icebergs = F , NO icebergs used'
             WRITE(numout,*) '~~~~~~~~ '
          ENDIF

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/IOM/in_out_manager.F90

-                      r3680
+                      r3875
    !!                        logical units
    !!----------------------------------------------------------------------
+   INTEGER ::   numstp      =   -1      !: logical unit for time step
+   INTEGER ::   numtime     =   -1      !: logical unit for timing
+   INTEGER ::   numout      =    6      !: logical unit for output print
+   INTEGER ::   numnam      =   -1      !: logical unit for namelist
+   INTEGER ::   numnam_ice  =   -1      !: logical unit for ice namelist
+   INTEGER ::   numevo_ice  =   -1      !: logical unit for ice variables (temp. evolution)
+   INTEGER ::   numsol      =   -1      !: logical unit for solver statistics
+   INTEGER ::   numdct_in   =   -1      !: logical unit for transports computing
+   INTEGER ::   numdct_vol  =   -1      !: logical unit for voulume transports output
+   INTEGER ::   numdct_heat =   -1      !: logical unit for heat    transports output
+   INTEGER ::   numdct_salt =   -1      !: logical unit for salt    transports output
+   INTEGER ::   numfl      =   -1      !: logical unit for floats ascii output
+   INTEGER ::   numflo     =   -1      !: logical unit for floats ascii output
+   INTEGER ::   numstp          =   -1      !: logical unit for time step
+   INTEGER ::   numtime         =   -1      !: logical unit for timing
+   INTEGER ::   numout          =    6      !: logical unit for output print
+   INTEGER ::   numnam_ref      =   -1      !: logical unit for reference namelist
+   INTEGER ::   numnam_cfg      =   -1      !: logical unit for configuration specific namelist
+   INTEGER ::   numond          =    7      !: logical unit for Output Namelist Dynamics
+   INTEGER ::   numnam_ice_ref  =   -1      !: logical unit for ice reference namelist
+   INTEGER ::   numnam_ice_cfg  =   -1      !: logical unit for ice reference namelist
+   INTEGER ::   numoni          =    8      !: logical unit for Output Namelist Ice
+   INTEGER ::   numevo_ice      =   -1      !: logical unit for ice variables (temp. evolution)
+   INTEGER ::   numsol          =   -1      !: logical unit for solver statistics
+   INTEGER ::   numdct_in       =   -1      !: logical unit for transports computing
+   INTEGER ::   numdct_vol      =   -1      !: logical unit for voulume transports output
+   INTEGER ::   numdct_heat     =   -1      !: logical unit for heat    transports output
+   INTEGER ::   numdct_salt     =   -1      !: logical unit for salt    transports output
+   INTEGER ::   numfl           =   -1      !: logical unit for floats ascii output
+   INTEGER ::   numflo          =   -1      !: logical unit for floats ascii output
    !!----------------------------------------------------------------------

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90

r3680	r3875
156	156	WRITE(numout,*)
157	157	SELECT CASE ( jprstlib )
158		CASE ( jpnf90 ) ; WRITE(numout,*) 'rst_read : read oce NetCDF restart file'
	158	CASE ( jpnf90 ) ; WRITE(numout,*) 'rst_read : read oce NetCDF restart file ',TRIM(cn_ocerst_in)//'.nc'
159	159	CASE ( jprstdimg ) ; WRITE(numout,*) 'rst_read : read oce binary restart file'
160	160	END SELECT

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/LBC/lib_mpp.F90

-                      r3799
+                      r3875
    !!   ctl_warn   : initialization, namelist read, and parameters control
    !!   ctl_opn    : Open file and check if required file is available.
+   !!   get_unit    : give the index of an unused logical unit
+   !!   ctl_nam    : Prints informations when an error occurs while reading a namelist
+   !!   get_unit   : give the index of an unused logical unit
    !!----------------------------------------------------------------------
 #if   defined key_mpp_mpi
 …
    IMPLICIT NONE
    PRIVATE
    PUBLIC   ctl_stop, ctl_warn, get_unit, ctl_opn
+   PUBLIC   ctl_stop, ctl_warn, get_unit, ctl_opn, ctl_nam
    PUBLIC   mynode, mppstop, mppsync, mpp_comm_free
    PUBLIC   mpp_ini_north, mpp_lbc_north, mpp_lbc_north_e
 …
    FUNCTION mynode( ldtxt, kumnam, kstop, localComm )
+   FUNCTION mynode( ldtxt, kumnam_ref , kumnam_cfg , kumond , kstop, localComm )
       !!----------------------------------------------------------------------
       !!                  ***  routine mynode  ***
 …
       !!----------------------------------------------------------------------
       CHARACTER(len=*),DIMENSION(:), INTENT(  out) ::   ldtxt
+      INTEGER                      , INTENT(in   ) ::   kumnam       ! namelist logical unit
+      INTEGER                      , INTENT(inout) ::   kstop        ! stop indicator
+      INTEGER                      , INTENT(in   ) ::   kumnam_ref     ! logical unit for reference namelist
+      INTEGER                      , INTENT(in   ) ::   kumnam_cfg     ! logical unit for configuration namelist
+      INTEGER                      , INTENT(in   ) ::   kumond         ! logical unit for namelist output
+      INTEGER                      , INTENT(inout) ::   kstop          ! stop indicator
       INTEGER, OPTIONAL            , INTENT(in   ) ::   localComm
+      !
       INTEGER ::   mynode, ierr, code, ji, ii
+      INTEGER ::   mynode, ierr, code, ji, ii, ios
       LOGICAL ::   mpi_was_called
+      !
 …
+      !
       jpni = -1; jpnj = -1; jpnij = -1
+      REWIND( kumnam )               ! Namelist namrun : parameters of the run
+      READ  ( kumnam, nammpp )
+      REWIND( kumnam_ref )              ! Namelist nammpp in reference namelist: mpi variables
+      READ  ( kumnam_ref, nammpp, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nammpp in reference namelist', lwp )
+      REWIND( kumnam_cfg )              ! Namelist nammpp in configuration namelist: mpi variables
+      READ  ( kumnam_cfg, nammpp, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nammpp in configuration namelist', lwp )
+      WRITE(kumond, nammpp)
       !                              ! control print
       WRITE(ldtxt(ii),*) '   Namelist nammpp'                                                     ;   ii = ii + 1
 …
    END FUNCTION lib_mpp_alloc
    FUNCTION mynode( ldtxt, kumnam, kstop, localComm ) RESULT (function_value)
+   FUNCTION mynode( ldtxt, kumnam_ref, knumnam_cfg, kstop, localComm ) RESULT (function_value)
       INTEGER, OPTIONAL            , INTENT(in   ) ::   localComm
       CHARACTER(len=*),DIMENSION(:) ::   ldtxt
       INTEGER ::   kumnam, kstop
+      INTEGER ::   kumnam_ref, knumnam_cfg , kstop
       IF( PRESENT( localComm ) .OR. .NOT.PRESENT( localComm ) )   function_value = 0
       IF( .FALSE. )   ldtxt(:) = 'never done'
 …
    !!----------------------------------------------------------------------
    !!   All cases:         ctl_stop, ctl_warn, get_unit, ctl_opn   routines
+   !!   All cases:         ctl_stop, ctl_warn, get_unit, ctl_opn, ctl_nam   routines
    !!----------------------------------------------------------------------
 …
    END SUBROUTINE ctl_opn
+   SUBROUTINE ctl_nam ( kios, cdnam, ldwp )
+      !!----------------------------------------------------------------------
+      !!                  ***  ROUTINE ctl_nam  ***
+      !!
+      !! ** Purpose :   Informations when error while reading a namelist
+      !!
+      !! ** Method  :   Fortan open
+      !!----------------------------------------------------------------------
+      INTEGER          , INTENT(inout) ::   kios      ! IO status after reading the namelist
+      CHARACTER(len=*) , INTENT(in   ) ::   cdnam     ! group name of namelist for which error occurs
+      CHARACTER(len=4)                 ::   clios     ! string to convert iostat in character for print
+      LOGICAL          , INTENT(in   ) ::   ldwp      ! boolean term for print
+      !!----------------------------------------------------------------------
+      !
+      ! ----------------
+      WRITE (clios, '(I4.0)') kios
+      IF( kios < 0 ) THEN
+         CALL ctl_warn( 'W A R N I N G:  end of record or file while reading namelist ' &
+ &           // TRIM(cdnam) // ' iostat = ' // TRIM(clios) )
+      ENDIF
+      IF( kios > 0 ) THEN
+         CALL ctl_stop( 'E R R O R :   misspelled variable in namelist ' &
+ &           // TRIM(cdnam) // ' iostat = ' // TRIM(clios) )
+      ENDIF
+      kios = 0
+      RETURN
+   END SUBROUTINE ctl_nam
    INTEGER FUNCTION get_unit()

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/LBC/mppini_2.h90

-                      r3818
+                      r3875
       INTEGER ::  inum                        ! temporary logical unit
       INTEGER ::  idir                        ! temporary integers
+      INTEGER ::   ios                        ! Local integer output status for namelist read
       INTEGER ::   &
          ii, ij, ifreq, il1, il2,          &  ! temporary integers
 …
       !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
       !!----------------------------------------------------------------------
+      REWIND ( numnam )              ! Read Namelist namzgr : vertical coordinate'
+      READ   ( numnam, namzgr )
+      REWIND( numnam_ref )              ! Namelist namzgr in reference namelist : Vertical coordinate
+      READ  ( numnam_ref, namzgr, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzgr in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namzgr in configuration namelist : Vertical coordinate
+      READ  ( numnam_cfg, namzgr, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzgr in configuration namelist', lwp )
+      WRITE ( numond, namzgr )
       IF(lwp)WRITE(numout,*)

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90

-                      r3634
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER ::   ioptio         ! ???
+      INTEGER ::   ios            ! Local : output status for namelist read
       LOGICAL ::   ll_print = .FALSE.    ! Logical flag for printing viscosity coef.
       !!
 …
    !!----------------------------------------------------------------------
+      REWIND( numnam )                  ! Read Namelist namdyn_ldf : Lateral physics
+      READ  ( numnam, namdyn_ldf )
+      REWIND( numnam_ref )              ! Namelist namdyn_ldf in reference namelist : Lateral physics
+      READ  ( numnam_ref, namdyn_ldf, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_ldf in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namdyn_ldf in configuration namelist : Lateral physics
+      READ  ( numnam_cfg, namdyn_ldf, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdyn_ldf in configuration namelist', lwp )
+      WRITE ( numond, namdyn_ldf )
       IF(lwp) THEN                      ! Parameter print
 …
          WRITE(numout,*) 'ldf_dyn : lateral momentum physics'
          WRITE(numout,*) '~~~~~~~'
          WRITE(numout,*) '   Namelist nam_dynldf : set lateral mixing parameters'
+         WRITE(numout,*) '   Namelist namdyn_ldf : set lateral mixing parameters'
          WRITE(numout,*) '      laplacian operator                      ln_dynldf_lap   = ', ln_dynldf_lap
          WRITE(numout,*) '      bilaplacian operator                    ln_dynldf_bilap = ', ln_dynldf_bilap

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90

-                      r3634
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER ::   ioptio               ! temporary integer
+      INTEGER ::   ios                  ! temporary integer
       LOGICAL ::   ll_print = .FALSE.   ! =T print eddy coef. in numout
       !!
 …
       ! =============================================
+      REWIND( numnam )                  ! Read Namelist namtra_ldf : Lateral physics on tracers
+      READ  ( numnam, namtra_ldf )
+      REWIND( numnam_ref )              ! Namelist namtra_ldf in reference namelist : Lateral physics on tracers
+      READ  ( numnam_ref, namtra_ldf, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtra_ldf in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namtra_ldf in configuration namelist : Lateral physics on tracers
+      READ  ( numnam_cfg, namtra_ldf, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtra_ldf in configuration namelist', lwp )
+      WRITE ( numond, namtra_ldf )
       IF(lwp) THEN                      ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/OBC/obcini.F90

-                      r2715
+                      r3875
       !!
       INTEGER  ::   ji, jj, istop , inumfbc
+      INTEGER  ::   ios                     ! Local integer output status for namelist read
       INTEGER, DIMENSION(4) ::   icorner
       REAL(wp), DIMENSION(2) ::   ztestmask
 …
       !!----------------------------------------------------------------------
+      REWIND( numnam )              ! Namelist namobc : open boundaries
+      READ  ( numnam, namobc )
+      REWIND( numnam_ref )              ! Namelist namobc in reference namelist : Open boundaries
+      READ  ( numnam_ref, namobc, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namobc in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namobc in configuration namelist : Open boundaries
+      READ  ( numnam_cfg, namobc, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namobc in configuration namelist', lwp )
+      WRITE ( numond, namobc )
       ! convert DOCTOR namelist name into the OLD names

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/OBS/diaobs.F90

-                      r3651
+                      r3875
       INTEGER :: ji
       INTEGER :: jset
+      INTEGER :: ios                 ! Local integer output status for namelist read
       LOGICAL :: lmask(MaxNumFiles), ll_u3d, ll_v3d
 …
       ! Read Namelist namobs : control observation diagnostics
+      REWIND( numnam )
+      READ  ( numnam, namobs )
+      REWIND( numnam_ref )              ! Namelist namobs in reference namelist : Diagnostic: control observation
+      READ  ( numnam_ref, namobs, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namobs in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namobs in configuration namelist : Diagnostic: control observation
+      READ  ( numnam_cfg, namobs, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namobs in configuration namelist', lwp )
+      WRITE ( numond, namobs )
       ! Count number of files for each type

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/albedo.F90

-                      r3625
+                      r3875
       !! ** Method  :   Read the namelist namsbc_alb
       !!----------------------------------------------------------------------
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       NAMELIST/namsbc_alb/ rn_cloud, rn_albice, rn_alphd, rn_alphdi, rn_alphc
       !!----------------------------------------------------------------------
 …
       albd_init = 1                     ! indicate that the initialization has been done
+      !
+      REWIND( numnam )                  ! Read Namelist namsbc_alb : albedo parameters
+      READ  ( numnam, namsbc_alb )
+      REWIND( numnam_ref )              ! Namelist namsbc_alb in reference namelist : Albedo parameters
+      READ  ( numnam_ref, namsbc_alb, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_alb in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namsbc_alb in configuration namelist : Albedo parameters
+      READ  ( numnam_cfg, namsbc_alb, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_alb in configuration namelist', lwp )
+      WRITE ( numond, namsbc_alb )
+      !
       IF(lwp) THEN                      ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/cyclone.F90

-                      r3680
+                      r3875
          sn_tc = FLD_N( 'tc_track',     6     ,  'tc'     ,  .true.    , .false. ,   'yearly'  , ''       , ''         )
+         !
+         !         REWIND( numnam )                    ! ... read in namlist namsbc_core
+         !         READ  ( numnam, namsbc_tc )
+         !
+         !  Namelist is read in namsbc_core
          ! set sf structure
          ALLOCATE( sf(1), STAT=ierror )

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/sbcana.F90

-                      r3625
+                      r3875
       INTEGER, INTENT(in) ::   kt       ! ocean time step
+      !
+      INTEGER  ::   ios                   ! Local integer output status for namelist read
       REAL(wp) ::   zfacto                ! local scalar
       REAL(wp) ::   zrhoa  = 1.22_wp      ! Air density kg/m3
 …
       IF( kt == nit000 ) THEN
+         !
+         REWIND( numnam )                    ! Read Namelist namsbc : surface fluxes
+         READ  ( numnam, namsbc_ana )
+         REWIND( numnam_ref )              ! Namelist namsbc_ana in reference namelist : Analytical surface fluxes
+         READ  ( numnam_ref, namsbc_ana, IOSTAT = ios, ERR = 901)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_ana in reference namelist', lwp )
+         REWIND( numnam_cfg )              ! Namelist namsbc_ana in configuration namelist : Analytical surface fluxes
+         READ  ( numnam_cfg, namsbc_ana, IOSTAT = ios, ERR = 902 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_ana in configuration namelist', lwp )
+         WRITE ( numond, namsbc_ana )
+         !
          IF(lwp) WRITE(numout,*)' '

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/sbcapr.F90

-                      r3795
+                      r3875
       !!
       INTEGER            ::   ierror  ! local integer
+      INTEGER            ::   ios     ! Local integer output status for namelist read
       !!
       CHARACTER(len=100) ::  cn_dir   ! Root directory for location of ssr files
 …
          cn_dir  = './'          ! directory in which the Patm data are
+         REWIND( numnam )                             !* read in namlist namsbc_apr
+         READ  ( numnam, namsbc_apr )
+         REWIND( numnam_ref )              ! Namelist namsbc_apr in reference namelist : File for atmospheric pressure forcing
+         READ  ( numnam_ref, namsbc_apr, IOSTAT = ios, ERR = 901)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_apr in reference namelist', lwp )
+         REWIND( numnam_cfg )              ! Namelist namsbc_apr in configuration namelist : File for atmospheric pressure forcing
+         READ  ( numnam_cfg, namsbc_apr, IOSTAT = ios, ERR = 902 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_apr in configuration namelist', lwp )
+         WRITE ( numond, namsbc_apr )
+         !
          ALLOCATE( sf_apr(1), STAT=ierror )           !* allocate and fill sf_sst (forcing structure) with sn_sst

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_clio.F90

-                      r3625
+                      r3875
       INTEGER, INTENT( in  ) ::   kt   ! ocean time step
       !!
       INTEGER  ::   ifpr, jfpr   ! dummy indices
+      INTEGER  ::   ifpr, jfpr                   ! dummy indices
       INTEGER  ::   ierr0, ierr1, ierr2, ierr3   ! return error code
+      INTEGER  ::   ios                          ! Local integer output status for namelist read
       !!
       CHARACTER(len=100) ::  cn_dir                            !   Root directory for location of CLIO files
 …
          sn_prec = FLD_N( 'precip' ,    -1     , 'precip' ,  .true.    , .false. ,   'yearly'  , ''       , ''       )
+         REWIND( numnam )                    ! ... read in namlist namsbc_clio
+         READ  ( numnam, namsbc_clio )
+         REWIND( numnam_ref )              ! Namelist namsbc_clio in reference namelist : CLIO files
+         READ  ( numnam_ref, namsbc_clio, IOSTAT = ios, ERR = 901)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_clio in reference namelist', lwp )
+         REWIND( numnam_cfg )              ! Namelist namsbc_clio in configuration namelist : CLIO files
+         READ  ( numnam_cfg, namsbc_clio, IOSTAT = ios, ERR = 902 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_clio in configuration namelist', lwp )
+         WRITE ( numond, namsbc_clio )
          ! store namelist information in an array

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

-                      r3772
+                      r3875
       INTEGER  ::   ifpr     ! dummy loop indice
       INTEGER  ::   jfld     ! dummy loop arguments
+      INTEGER  ::   ios      ! Local integer output status for namelist read
       !!
       CHARACTER(len=100) ::  cn_dir   !   Root directory for location of core files
 …
          sn_tdif = FLD_N( 'taudif' ,    24     , 'taudif' ,  .true.    , .false. ,   'yearly'  , ''       , ''       )
+         !
+         REWIND( numnam )                          ! read in namlist namsbc_core
+         READ  ( numnam, namsbc_core )
+         REWIND( numnam_ref )              ! Namelist namsbc_core in reference namelist : CORE bulk parameters
+         READ  ( numnam_ref, namsbc_core, IOSTAT = ios, ERR = 901)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_core in reference namelist', lwp )
+         REWIND( numnam_cfg )              ! Namelist namsbc_core in configuration namelist : CORE bulk parameters
+         READ  ( numnam_cfg, namsbc_core, IOSTAT = ios, ERR = 902 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_core in configuration namelist', lwp )
+         WRITE ( numond, namsbc_core )
          !                                         ! check: do we plan to use ln_dm2dc with non-daily forcing?
          IF( ln_dm2dc .AND. sn_qsr%nfreqh /= 24 )   &

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90

-                      r3625
+                      r3875
       INTEGER  :: ifpr     ! dummy loop indice
       INTEGER  :: jj,ji    ! dummy loop arguments
+      INTEGER  ::   ios    ! Local integer output status for namelist read
       REAL(wp) :: act_hour
       !!--------------------------------------------------------------------
 …
             sn_prec = FLD_N( 'precip_cmap' ,  -1   ,  'precip'  ,  .true.    ,  .true. ,   'yearly'  , ''       , ''         )
+            !
+            REWIND( numnam )                    ! ... read in namlist namsbc_mfs
+            READ  ( numnam, namsbc_mfs )
+            REWIND( numnam_ref )              ! Namelist namsbc_msf in reference namelist : MFS files
+            READ  ( numnam_ref, namsbc_mfs, IOSTAT = ios, ERR = 901)
+         IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_mfs in reference namelist', lwp )
+            REWIND( numnam_cfg )              ! Namelist namsbc_msf in configuration namelist : MFS files
+            READ  ( numnam_cfg, namsbc_mfs, IOSTAT = ios, ERR = 902 )
+         IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_mfs in configuration namelist', lwp )
+            WRITE ( numond, namsbc_mfs )
+            !
             ! store namelist information in an array

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

-                      r3680
+                      r3875
       !!
       INTEGER ::   jn   ! dummy loop index
+      INTEGER ::   ios  ! Local integer output status for namelist read
       REAL(wp), POINTER, DIMENSION(:,:) ::   zacs, zaos
       !!
 …
       sn_rcv_co2    = FLD_C( 'none'                ,    'no'    ,     ''      ,         ''          ,   ''    )
+      REWIND( numnam )                    ! ... read namlist namsbc_cpl
+      READ  ( numnam, namsbc_cpl )
+      REWIND( numnam_ref )              ! Namelist namsbc_cpl in reference namelist : Variables for OASIS coupling
+      READ  ( numnam_ref, namsbc_cpl, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_cpl in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namsbc_cpl in configuration namelist : Variables for OASIS coupling
+      READ  ( numnam_cfg, namsbc_cpl, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_cpl in configuration namelist', lwp )
+      WRITE ( numond, namsbc_cpl )
       IF(lwp) THEN                        ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/sbcflx.F90

-                      r3625
+                      r3875
       INTEGER  ::   ji, jj, jf            ! dummy indices
       INTEGER  ::   ierror                ! return error code
+      INTEGER  ::   ios                   ! Local integer output status for namelist read
       REAL(wp) ::   zfact                 ! temporary scalar
       REAL(wp) ::   zrhoa  = 1.22         ! Air density kg/m3
 …
          sn_emp  = FLD_N(  'emp'  ,    24     ,  'emp'     ,  .false.   , .false. ,   'yearly'  , ''       , ''        )
+         !
+         REWIND ( numnam )                         ! read in namlist namflx
+         READ   ( numnam, namsbc_flx )
+         REWIND( numnam_ref )              ! Namelist namsbc_flx in reference namelist : Files for fluxes
+         READ  ( numnam_ref, namsbc_flx, IOSTAT = ios, ERR = 901)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_flx in reference namelist', lwp )
+         REWIND( numnam_cfg )              ! Namelist namsbc_flx in configuration namelist : Files for fluxes
+         READ  ( numnam_cfg, namsbc_flx, IOSTAT = ios, ERR = 902 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_flx in configuration namelist', lwp )
+         WRITE ( numond, namsbc_flx )
+         !
          !                                         ! check: do we plan to use ln_dm2dc with non-daily forcing?

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

-                      r3625
+                      r3875
          sn_bot5 = FLD_N( 'botmeltn5_1m' ,    -1.    ,  'botmeltn5' ,  .true.    , .true.  ,  ' yearly'  , ''       , ''         )
+!         REWIND ( numnam )               ! ... at some point might read in from NEMO namelist?
+!         READ   ( numnam, namsbc_cice )
+! ... at some point might read in from NEMO namelist?
+!!$      REWIND( numnam_ref )              ! Namelist namsbc_cice in reference namelist :
+!!$      READ  ( numnam_ref, namsbc_cice, IOSTAT = ios, ERR = 901)
+!!$901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_cice in reference namelist', lwp )
+!!$
+!!$      REWIND( numnam_cfg )              ! Namelist namsbc_cice in configuration namelist : Parameters of the run
+!!$      READ  ( numnam_cfg, namsbc_cice, IOSTAT = ios, ERR = 902 )
+!!$902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_cice in configuration namelist', lwp )
+!!$      WRITE ( numond, namsbc_cice )
          ! store namelist information in an array

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_if.F90

-                      r3625
+                      r3875
       INTEGER  ::   ji, jj     ! dummy loop indices
       INTEGER  ::   ierror     ! return error code
+      INTEGER  ::   ios        ! Local integer output status for namelist read
       REAL(wp) ::   ztrp, zsice, zt_fzp, zfr_obs
       REAL(wp) ::   zqri, zqrj, zqrp, zqi
 …
          sn_ice = FLD_N('ice_cover',    -1    ,  'ice_cov' ,  .true.    , .true. ,   'yearly'  , ''       , ''         )
+         REWIND ( numnam )               ! ... read in namlist namiif
+         READ   ( numnam, namsbc_iif )
+         REWIND( numnam_ref )              ! Namelist namsbc_iif in reference namelist : Ice if file
+         READ  ( numnam_ref, namsbc_iif, IOSTAT = ios, ERR = 901)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_iif in reference namelist', lwp )
+         REWIND( numnam_cfg )              ! Namelist Namelist namsbc_iif in configuration namelist : Ice if file
+         READ  ( numnam_cfg, namsbc_iif, IOSTAT = ios, ERR = 902 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_iif in configuration namelist', lwp )
+         WRITE ( numond, namsbc_iif )
          ALLOCATE( sf_ice(1), STAT=ierror )

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

-                      r3764
+                      r3875
          &             ln_blk_mfs, ln_apr_dyn, nn_ice,  nn_ice_embd, ln_dm2dc   , ln_rnf,   &
          &             ln_ssr    , nn_fwb    , ln_cdgw , ln_wave , ln_sdw
+      INTEGER  ::   ios
       !!----------------------------------------------------------------------
 …
       ENDIF
+      REWIND( numnam )           ! Read Namelist namsbc
+      READ  ( numnam, namsbc )
+      REWIND( numnam_ref )              ! Namelist namsbc in reference namelist : Surface boundary
+      READ  ( numnam_ref, namsbc, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namsbc in configuration namelist : Parameters of the run
+      READ  ( numnam_cfg, namsbc, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc in configuration namelist', lwp )
+      WRITE ( numond, namsbc )
       !                          ! overwrite namelist parameter using CPP key information

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/sbcrnf.F90

-                      r3832
+                      r3875
       INTEGER           ::   ji, jj, jk    ! dummy loop indices
       INTEGER           ::   ierror, inum  ! temporary integer
+      INTEGER           ::   ios           ! Local integer output status for namelist read
+      !
       NAMELIST/namsbc_rnf/ cn_dir, ln_rnf_emp, ln_rnf_depth, ln_rnf_tem, ln_rnf_sal,   &
 …
       sn_dep_rnf = FLD_N( 'runoffs',   0.  , 'rodepth'  ,  .FALSE.   , .true. ,   'yearly'  , ''    , ''  )
+      !
+      REWIND ( numnam )                         ! Read Namelist namsbc_rnf
+      READ   ( numnam, namsbc_rnf )
+      REWIND( numnam_ref )              ! Namelist namsbc_rnf in reference namelist : Runoffs
+      READ  ( numnam_ref, namsbc_rnf, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_rnf in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namsbc_rnf in configuration namelist : Runoffs
+      READ  ( numnam_cfg, namsbc_rnf, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_rnf in configuration namelist', lwp )
+      WRITE ( numond, namsbc_rnf )
+      !
       !                                         ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssr.F90

-                      r3764
+                      r3875
       CHARACTER(len=100) ::  cn_dir          ! Root directory for location of ssr files
       TYPE(FLD_N) ::   sn_sst, sn_sss        ! informations about the fields to be read
-      NAMELIST/namsbc_ssr/ cn_dir, nn_sstr, nn_sssr, rn_dqdt, rn_deds, sn_sst, sn_sss, ln_sssr_bnd, rn_sssr_bnd
       !!----------------------------------------------------------------------
+      !
 …
       TYPE(FLD_N) ::   sn_sst, sn_sss        ! informations about the fields to be read
       NAMELIST/namsbc_ssr/ cn_dir, nn_sstr, nn_sssr, rn_dqdt, rn_deds, sn_sst, sn_sss, ln_sssr_bnd, rn_sssr_bnd
+      INTEGER     ::  ios
       !!----------------------------------------------------------------------
+      !
 …
       sn_sst = FLD_N( 'sst'    ,    24     ,  'sst'     ,  .false.   , .false. ,   'yearly'  , ''       , ''         )
       sn_sss = FLD_N( 'sss'    ,    -1     ,  'sss'     ,  .true.    , .false. ,   'yearly'  , ''       , ''         )
+      REWIND( numnam )             !* read in namlist namflx
+      READ  ( numnam, namsbc_ssr )
+      REWIND( numnam_ref )              ! Namelist namsbc_ssr in reference namelist :
+      READ  ( numnam_ref, namsbc_ssr, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_ssr in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namsbc_ssr in configuration namelist :
+      READ  ( numnam_cfg, namsbc_ssr, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_ssr in configuration namelist', lwp )
+      WRITE ( numond, namsbc_ssr )
       IF(lwp) THEN                 !* control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/sbcwave.F90

-                      r3680
+                      r3875
       INTEGER                ::  ierror   ! return error code
       INTEGER                ::  ifpr, jj,ji,jk
+      INTEGER                ::   ios     ! Local integer output status for namelist read
       REAL(wp),DIMENSION(:,:,:),POINTER             ::  udummy,vdummy,hdivdummy,rotdummy
       REAL                                          ::  z2dt,z1_2dt
 …
          cn_dir = './'          ! directory in which the wave data are
+         REWIND( numnam )                             !* read in namlist namsbc_wave
+         READ  ( numnam, namsbc_wave )
+         REWIND( numnam_ref )              ! Namelist namsbc_wave in reference namelist : File for drag coeff. from wave model
+         READ  ( numnam_ref, namsbc_wave, IOSTAT = ios, ERR = 901)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_wave in reference namelist', lwp )
+         REWIND( numnam_cfg )              ! Namelist namsbc_wave in configuration namelist : File for drag coeff. from wave model
+         READ  ( numnam_cfg, namsbc_wave, IOSTAT = ios, ERR = 902 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_wave in configuration namelist', lwp )
+         WRITE ( numond, namsbc_wave )
+         !

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SBC/tideini.F90

-                      r3651
+                      r3875
     INTEGER, INTENT( in ) ::   kt     ! ocean time-step
     CHARACTER(LEN=4), DIMENSION(jpmax_harmo) :: clname
+    INTEGER  ::   ios                 ! Local integer output status for namelist read
+    !
     NAMELIST/nam_tide/ln_tide_pot, ln_tide_ramp, rdttideramp, clname
 …
+       !
        ! Read Namelist nam_tide
+       REWIND ( numnam )
+       READ   ( numnam, nam_tide )
+       REWIND( numnam_ref )              ! Namelist nam_tide in reference namelist : Tides
+       READ  ( numnam_ref, nam_tide, IOSTAT = ios, ERR = 901)
+    IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_tide in reference namelist', lwp )
+       REWIND( numnam_cfg )              ! Namelist nam_tide in configuration namelist : Tides
+       READ  ( numnam_cfg, nam_tide, IOSTAT = ios, ERR = 902 )
+    IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_tide in configuration namelist', lwp )
+       WRITE ( numond, nam_tide )
+       !
        nb_harmo=0

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/SOL/solver.F90

-                      r3764
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER, INTENT(in) :: kt
+      INTEGER             ::   ios       ! Local integer output status for namelist read
       !!
       NAMELIST/namsol/ nn_solv, nn_sol_arp, nn_nmin, nn_nmax, nn_nmod, rn_eps, rn_resmax, rn_sor
 …
       ENDIF
+      REWIND( numnam )              !* Namelist namsol : elliptic solver / free surface
+      READ  ( numnam, namsol )
+      REWIND( numnam_ref )              ! Namelist namsol in reference namelist : Elliptic solver / free surface
+      READ  ( numnam_ref, namsol, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsol in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namsol in configuration namelist : Elliptic solver / free surface
+      READ  ( numnam_cfg, namsol, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsol in configuration namelist', lwp )
+      WRITE ( numond, namsol )
       IF(lwp) THEN                  !* Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90

-                      r3625
+                      r3875
       NAMELIST/nameos/ nn_eos, rn_alpha, rn_beta
       !!----------------------------------------------------------------------
+      !
+      REWIND( numnam )            ! Read Namelist nameos : equation of state
+      READ  ( numnam, nameos )
+      INTEGER  ::   ios
+      !
+      REWIND( numnam_ref )              ! Namelist nameos in reference namelist : equation of state
+      READ  ( numnam_ref, nameos, IOSTAT = ios, ERR = 901 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nameos in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist nameos in configuration namelist : equation of state
+      READ  ( numnam_cfg, nameos, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nameos in configuration namelist', lwp )
+      WRITE( numond, nameos )
+      !
       IF(lwp) THEN                ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90

-                      r3718
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER ::   ioptio
+      INTEGER ::   ios                 ! Local integer output status for namelist read
       !!
       NAMELIST/namtra_adv/ ln_traadv_cen2 , ln_traadv_tvd,     &
 …
       !!----------------------------------------------------------------------
+      REWIND( numnam )                ! Read Namelist namtra_adv : tracer advection scheme
+      READ  ( numnam, namtra_adv )
+      REWIND( numnam_ref )              ! Namelist namtra_adv in reference namelist : Tracer advection scheme
+      READ  ( numnam_ref, namtra_adv, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtra_adv in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namtra_adv in configuration namelist : Tracer advection scheme
+      READ  ( numnam_cfg, namtra_adv, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtra_adv in configuration namelist', lwp )
+      WRITE ( numond, namtra_adv )
       IF(lwp) THEN                    ! Namelist print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/TRA/trabbc.F90

-                      r3625
+                      r3875
       INTEGER  ::   ji, jj              ! dummy loop indices
       INTEGER  ::   inum                ! temporary logical unit
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       !!
       NAMELIST/nambbc/ln_trabbc, nn_geoflx, rn_geoflx_cst
       !!----------------------------------------------------------------------
+      REWIND( numnam )                 ! Read Namelist nambbc : bottom momentum boundary condition
+      READ  ( numnam, nambbc )
+      REWIND( numnam_ref )              ! Namelist nambbc in reference namelist : Bottom momentum boundary condition
+      READ  ( numnam_ref, nambbc, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambbc in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist nambbc in configuration namelist : Bottom momentum boundary condition
+      READ  ( numnam_cfg, nambbc, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambbc in configuration namelist', lwp )
+      WRITE ( numond, nambbc )
       IF(lwp) THEN                     ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/TRA/trabbl.F90

-                      r3764
+                      r3875
       INTEGER ::   ji, jj               ! dummy loop indices
       INTEGER ::   ii0, ii1, ij0, ij1   ! temporary integer
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       REAL(wp), POINTER, DIMENSION(:,:) :: zmbk
       !!
 …
+      !
+      REWIND ( numnam )              !* Read Namelist nambbl : bottom boundary layer scheme
+      READ   ( numnam, nambbl )
+      REWIND( numnam_ref )              ! Namelist nambbl in reference namelist : Bottom boundary layer scheme
+      READ  ( numnam_ref, nambbl, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambbl in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist nambbl in configuration namelist : Bottom boundary layer scheme
+      READ  ( numnam_cfg, nambbl, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambbl in configuration namelist', lwp )
+      WRITE ( numond, nambbl )
+      !
       l_bbl = .TRUE.                 !* flag to compute bbl coef and transport

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/TRA/tradmp.F90

-                      r3294
+                      r3875
       !!----------------------------------------------------------------------
       NAMELIST/namtra_dmp/ ln_tradmp, nn_hdmp, nn_zdmp, rn_surf, rn_bot, rn_dep, nn_file
+      !!----------------------------------------------------------------------
+      REWIND ( numnam )                  ! Read Namelist namtra_dmp : temperature and salinity damping term
+      READ   ( numnam, namtra_dmp )
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
+      !!----------------------------------------------------------------------
+      REWIND( numnam_ref )              ! Namelist namtra_dmp in reference namelist : Temperature and salinity damping term
+      READ  ( numnam_ref, namtra_dmp, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtra_dmp in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namtra_dmp in configuration namelist : Temperature and salinity damping term
+      READ  ( numnam_cfg, namtra_dmp, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtra_dmp in configuration namelist', lwp )
+      WRITE ( numond, namtra_dmp )
       IF( lzoom )   nn_zdmp = 0          ! restoring to climatology at closed north or south boundaries

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/TRA/traqsr.F90

-                      r3680
+                      r3875
       !!----------------------------------------------------------------------
+      !
       INTEGER  ::   ji, jj, jk     ! dummy loop indices
+      INTEGER  ::   ji, jj, jk                   ! dummy loop indices
       INTEGER  ::   irgb, ierror, ioptio, nqsr   ! local integer
+      INTEGER  ::   ios                          ! Local integer output status for namelist read
       REAL(wp) ::   zz0, zc0  , zc1, zcoef       ! local scalars
       REAL(wp) ::   zz1, zc2  , zc3, zchl        !   -      -
 …
       sn_chl = FLD_N( 'chlorophyll' ,    -1     ,  'CHLA'    ,  .true.     , .true.  ,   'yearly'  , ''       , ''         )
+      !
+      REWIND( numnam )            ! Read Namelist namtra_qsr : ratio and length of penetration
+      READ  ( numnam, namtra_qsr )
+      REWIND( numnam_ref )              ! Namelist namtra_qsr in reference namelist : Ratio and length of penetration
+      READ  ( numnam_ref, namtra_qsr, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtra_qsr in reference namelist', lwp )
+      REWIND( numnam_cfg )              !  Namelist namtra_qsr in configuration namelist : Ratio and length of penetration
+      READ  ( numnam_cfg, namtra_qsr, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtra_qsr in configuration namelist', lwp )
+      WRITE ( numond, namtra_qsr )
+      !
       IF(lwp) THEN                ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/TRD/trdmod.F90

-                      r3294
+                      r3875
       !!----------------------------------------------------------------------
       USE in_out_manager          ! I/O manager
+      USE lib_mpp                 ! MPP library
       !!
       NAMELIST/namtrd/ nn_trd, nn_ctls, cn_trdrst_in, cn_trdrst_out, ln_trdmld_restart, rn_ucf, ln_trdmld_instant
+      INTEGER  ::   ios           ! Local integer output status for namelist read
       !!----------------------------------------------------------------------
       IF( l_trdtra .OR. l_trddyn )   THEN
+         REWIND( numnam )
+         READ  ( numnam, namtrd )      ! namelist namtrd : trends diagnostic
+         REWIND( numnam_ref )              ! Namelist namtrd in reference namelist : Diagnostics: trends
+         READ  ( numnam_ref, namtrd, IOSTAT = ios, ERR = 901)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrd in reference namelist', lwp )
+         REWIND( numnam_cfg )              ! Namelist namtrd in configuration namelist : Diagnostics: trends
+         READ  ( numnam_cfg, namtrd, IOSTAT = ios, ERR = 902 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrd in configuration namelist', lwp )
+         WRITE ( numond, namtrd )
          IF(lwp) THEN

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfbfr.F90

-                      r3820
+                      r3875
       INTEGER   ::   ikbu, ikbv   ! temporary integers
       INTEGER   ::   ictu, ictv   !    -          -
+      INTEGER   ::   ios
       REAL(wp)  ::   zminbfr, zmaxbfr   ! temporary scalars
       REAL(wp)  ::   zfru, zfrv         !    -         -
 …
       IF( nn_timing == 1 )  CALL timing_start('zdf_bfr_init')
+      !
+      REWIND ( numnam )               !* Read Namelist nam_bfr : bottom momentum boundary condition
+      READ   ( numnam, nambfr )
+      REWIND( numnam_ref )              ! Namelist nambfr in reference namelist : Bottom momentum boundary condition
+      READ  ( numnam_ref, nambfr, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambfr in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist nambfr in configuration namelist : Bottom momentum boundary condition
+      READ  ( numnam_cfg, nambfr, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambfr in configuration namelist', lwp )
+      WRITE ( numond, nambfr )
       !                               !* Parameter control and print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfddm.F90

-                      r3610
+                      r3875
       !!----------------------------------------------------------------------
       NAMELIST/namzdf_ddm/ rn_avts, rn_hsbfr
+      !!----------------------------------------------------------------------
+      !
+      REWIND( numnam )                ! Read Namelist namzdf_ddm : double diffusion mixing scheme
+      READ  ( numnam, namzdf_ddm )
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
+      !!----------------------------------------------------------------------
+      !
+      REWIND( numnam_ref )              ! Namelist namzdf_ddm in reference namelist : Double diffusion mixing scheme
+      READ  ( numnam_ref, namzdf_ddm, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf_ddm in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namzdf_ddm in configuration namelist : Double diffusion mixing scheme
+      READ  ( numnam_cfg, namzdf_ddm, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf_ddm in configuration namelist', lwp )
+      WRITE ( numond, namzdf_ddm )
+      !
       IF(lwp) THEN                    ! Parameter print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90

-                      r3798
+                      r3875
+      !
       INTEGER ::   jk    ! dummy loop indices
+      INTEGER ::   ios   ! Local integer output status for namelist read
       REAL(wp)::   zcr   ! local scalar
       !!
 …
       IF( nn_timing == 1 )  CALL timing_start('zdf_gls_init')
+      !
+      REWIND( numnam )                 !* Read Namelist namzdf_gls
+      READ  ( numnam, namzdf_gls )
+      REWIND( numnam_ref )              ! Namelist namzdf_gls in reference namelist : Vertical eddy diffivity and viscosity using gls turbulent closure scheme
+      READ  ( numnam_ref, namzdf_gls, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf_gls in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namzdf_gls in configuration namelist : Vertical eddy diffivity and viscosity using gls turbulent closure scheme
+      READ  ( numnam_cfg, namzdf_gls, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf_gls in configuration namelist', lwp )
+      WRITE ( numond, namzdf_gls )
       IF(lwp) THEN                     !* Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfini.F90

-                      r3680
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER ::   ioptio       ! temporary scalar
+      INTEGER ::   ios
       !!
       NAMELIST/namzdf/ rn_avm0, rn_avt0, nn_avb, nn_havtb, ln_zdfexp, nn_zdfexp,   &
 …
       !!----------------------------------------------------------------------
+      REWIND( numnam )           !* Read namzdf namelist : vertical mixing parameters
+      READ  ( numnam, namzdf )
+      REWIND( numnam_ref )              ! Namelist namzdf in reference namelist : Vertical mixing parameters
+      READ  ( numnam_ref, namzdf, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namzdf in reference namelist : Vertical mixing parameters
+      READ  ( numnam_cfg, namzdf, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf in configuration namelist', lwp )
+      WRITE ( numond, namzdf )
       IF(lwp) THEN               !* Parameter print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfkpp.F90

-                      r3792
+                      r3875
       REAL(wp) ::   zustar, zucube, zustvk, zeta, zehat   ! tempory scalars
 #endif
+      INTEGER  ::   ios            ! Local integer output status for namelist read
       REAL(wp) ::   zhbf           ! tempory scalars
       LOGICAL  ::   ll_kppcustom   ! 1st ocean level taken as surface layer
 …
       IF( nn_timing == 1 )  CALL timing_start('zdf_kpp_init')
+      !
+      REWIND ( numnam )               ! Read Namelist namkpp : K-Profile Parameterisation
+      READ   ( numnam, namzdf_kpp )
+      REWIND( numnam_ref )              ! Namelist namzdf_kpp in reference namelist : Vertical eddy diffivity and viscosity using kpp turbulent closure scheme
+      READ  ( numnam_ref, namzdf_kpp, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf_kpp in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namzdf_kpp in configuration namelist : Vertical eddy diffivity and viscosity using kpp turbulent closure scheme
+      READ  ( numnam_cfg, namzdf_kpp, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf_kpp in configuration namelist', lwp )
+      WRITE ( numond, namzdf_kpp )
       IF(lwp) THEN                    ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90

-                      r3625
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER :: ji, jj, jk   ! dummy loop indices
+      INTEGER ::   ios        ! Local integer output status for namelist read
       !!
       NAMELIST/namzdf_ric/ rn_avmri, rn_alp   , nn_ric  , rn_ekmfc,  &
 …
       !!----------------------------------------------------------------------
+      !
+      REWIND( numnam )               ! Read Namelist namzdf_ric : richardson number dependent Kz
+      READ  ( numnam, namzdf_ric )
+      REWIND( numnam_ref )              ! Namelist namzdf_ric in reference namelist : Vertical diffusion Kz depends on Richardson number
+      READ  ( numnam_ref, namzdf_ric, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf_ric in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namzdf_ric in configuration namelist : Vertical diffusion Kz depends on Richardson number
+      READ  ( numnam_cfg, namzdf_ric, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf_ric in configuration namelist', lwp )
+      WRITE ( numond, namzdf_ric )
+      !
       IF(lwp) THEN                   ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90

-                      r3680
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER ::   ji, jj, jk   ! dummy loop indices
+      INTEGER ::   ios
       !!
       NAMELIST/namzdf_tke/ rn_ediff, rn_ediss , rn_ebb , rn_emin  ,   &
 …
       !!----------------------------------------------------------------------
+      !
+      REWIND ( numnam )               !* Read Namelist namzdf_tke : Turbulente Kinetic Energy
+      READ   ( numnam, namzdf_tke )
+      REWIND( numnam_ref )              ! Namelist namzdf_tke in reference namelist : Turbulent Kinetic Energy
+      READ  ( numnam_ref, namzdf_tke, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf_tke in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namzdf_tke in configuration namelist : Turbulent Kinetic Energy
+      READ  ( numnam_cfg, namzdf_tke, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf_tke in configuration namelist', lwp )
+      WRITE ( numond, namzdf_tke )
+      !
       ri_cri   = 2._wp    / ( 2._wp + rn_ediss / rn_ediff )   ! resulting critical Richardson number

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftmx.F90

-                      r3625
+                      r3875
       INTEGER  ::   ji, jj, jk   ! dummy loop indices
       INTEGER  ::   inum         ! local integer
+      INTEGER  ::   ios
       REAL(wp) ::   ztpc, ze_z   ! local scalars
       REAL(wp), DIMENSION(:,:)  , POINTER ::  zem2, zek1   ! read M2 and K1 tidal energy
 …
       CALL wrk_alloc( jpi,jpj,jpk, zpc )
+      REWIND( numnam )               ! Read Namelist namtmx : Tidal Mixing
+      READ  ( numnam, namzdf_tmx )
+      REWIND( numnam_ref )              ! Namelist namzdf_tmx in reference namelist : Tidal Mixing
+      READ  ( numnam_ref, namzdf_tmx, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf_tmx in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namzdf_tmx in configuration namelist : Tidal Mixing
+      READ  ( numnam_cfg, namzdf_tmx, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namzdf_tmx in configuration namelist', lwp )
+      WRITE ( numond, namzdf_tmx )
       IF(lwp) THEN                   ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/module_example

-                      r2737
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER ::   ji, jj, jk, jit   ! dummy loop indices
+      INTEGER  ::   ios              ! Local integer output status for namelist read
       !!
       NAMELIST/namexa/ exa_v1, exa_v2, nexa_0, sn_ex
 …
       sn_ex%nfreqh = 2
+      !
+      REWIND( numnam )                          ! Read Namelist namexa : example parameters
+      READ  ( numnam, namexa )
+      REWIND( numnam_ref )              ! Namelist namexa in reference namelist : Example
+      READ  ( numnam_ref, namexa, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namexa in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namexa in configuration namelist : Example
+      READ  ( numnam_cfg, namexa, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namexa in configuration namelist', lwp )
+   ! Output namelist for control
+      WRITE ( numond, namexa )
+      !
       IF(lwp) THEN                              ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90

-                      r3769
+                      r3875
       INTEGER ::   ji            ! dummy loop indices
       INTEGER ::   ilocal_comm   ! local integer
+      INTEGER ::   ios
       CHARACTER(len=80), DIMENSION(16) ::   cltxt
       !!
 …
       cltxt = ''
+      !
+      !                             ! open Namelist file
+      CALL ctl_opn( numnam, 'namelist', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. )
+      !
+      READ( numnam, namctl )        ! Namelist namctl : Control prints & Benchmark
+      !                             ! Open reference namelist and configuration namelist files
+      CALL ctl_opn( numnam_ref, 'namelist_ref', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. )
+      CALL ctl_opn( numnam_cfg, 'namelist_cfg', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. )
+      CALL ctl_opn( numond, 'output.namelist.dyn', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea )
+      !
+      REWIND( numnam_ref )              ! Namelist namctl in reference namelist : Control prints & Benchmark
+      READ  ( numnam_ref, namctl, IOSTAT = ios, ERR = 901 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namctl in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namctl in confguration namelist : Control prints & Benchmark
+      READ  ( numnam_cfg, namctl, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namctl in configuration namelist', lwp )
+      WRITE( numond, namctl )
+      !
       !                             !--------------------------------------------!
 …
 # endif
       ENDIF
       narea = mynode( cltxt, numnam, nstop, ilocal_comm )   ! Nodes selection
+      narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop, ilocal_comm )   ! Nodes selection
 #else
 # if defined key_oasis3 || defined key_oasis4
 …
          CALL cpl_prism_init( ilocal_comm )                 ! nemo local communicator given by oasis
       ENDIF
       narea = mynode( cltxt, numnam, nstop, ilocal_comm )   ! Nodes selection (control print return in cltxt)
+      narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop, ilocal_comm )   ! Nodes selection (control print return in cltxt)
 # else
       ilocal_comm = 0
       narea = mynode( cltxt, numnam, nstop )                 ! Nodes selection (control print return in cltxt)
+      narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop )                 ! Nodes selection (control print return in cltxt)
 # endif
 #endif
 …
       CALL iom_close                                 ! close all input/output files managed by iom_*
+      !
+      IF( numstp      /= -1 )   CLOSE( numstp      )   ! time-step file
+      IF( numsol      /= -1 )   CLOSE( numsol      )   ! solver file
+      IF( numnam      /= -1 )   CLOSE( numnam      )   ! oce namelist
+      IF( numnam_ice  /= -1 )   CLOSE( numnam_ice  )   ! ice namelist
+      IF( numevo_ice  /= -1 )   CLOSE( numevo_ice  )   ! ice variables (temp. evolution)
+      IF( numout      /=  6 )   CLOSE( numout      )   ! standard model output file
+      IF( numdct_vol  /= -1 )   CLOSE( numdct_vol  )   ! volume transports
+      IF( numdct_heat /= -1 )   CLOSE( numdct_heat )   ! heat transports
+      IF( numdct_salt /= -1 )   CLOSE( numdct_salt )   ! salt transports
+      IF( numstp          /= -1 )   CLOSE( numstp          )   ! time-step file
+      IF( numsol          /= -1 )   CLOSE( numsol          )   ! solver file
+      IF( numnam_ref      /= -1 )   CLOSE( numnam_ref      )   ! oce reference namelist
+      IF( numnam_cfg      /= -1 )   CLOSE( numnam_cfg      )   ! oce configuration namelist
+      IF( numnam_ice_ref  /= -1 )   CLOSE( numnam_ice_ref  )   ! ice reference namelist
+      IF( numnam_ice_cfg  /= -1 )   CLOSE( numnam_ice_cfg  )   ! ice configuration namelist
+      IF( numevo_ice      /= -1 )   CLOSE( numevo_ice      )   ! ice variables (temp. evolution)
+      IF( numout          /=  6 )   CLOSE( numout          )   ! standard model output file
+      IF( numdct_vol      /= -1 )   CLOSE( numdct_vol      )   ! volume transports
+      IF( numdct_heat     /= -1 )   CLOSE( numdct_heat     )   ! heat transports
+      IF( numdct_salt     /= -1 )   CLOSE( numdct_salt     )   ! salt transports
+      !

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/SAS_SRC/nemogcm.F90

-                      r3769
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER ::   ji            ! dummy loop indices
+      INTEGER ::   ios
       INTEGER ::   ilocal_comm   ! local integer
       CHARACTER(len=80), DIMENSION(16) ::   cltxt
 …
       cltxt = ''
+      !
+      !                             ! open Namelist file
+      CALL ctl_opn( numnam, 'namelist', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. )
+      !
+      READ( numnam, namctl )        ! Namelist namctl : Control prints & Benchmark
+      !                             ! Open reference namelist and configuration namelist files
+      CALL ctl_opn( numnam_ref, 'namelist_ref', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. )
+      CALL ctl_opn( numnam_cfg, 'namelist_cfg', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE. )
+      CALL ctl_opn( numond, 'output.namelist.dyn', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, 6, .FALSE., narea )
+      !
+      REWIND( numnam_ref )              ! Namelist namctl in reference namelist : Control prints & Benchmark
+      READ  ( numnam_ref, namctl, IOSTAT = ios, ERR = 901 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namctl in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namctl in confguration namelist : Control prints & Benchmark
+      READ  ( numnam_cfg, namctl, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namctl in configuration namelist', lwp )
+      !
       !                             !--------------------------------------------!
 …
          CALL  xios_initialize( "nemo",return_comm=ilocal_comm )
       ENDIF
       narea = mynode( cltxt, numnam, nstop, ilocal_comm )   ! Nodes selection
+      narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop, ilocal_comm )   ! Nodes selection
 #else
       ilocal_comm = 0
       narea = mynode( cltxt, numnam, nstop )                 ! Nodes selection (control print return in cltxt)
+      narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop, ilocal_comm )   ! Nodes selection (control print return in cltxt)
 #endif
       narea = narea + 1                                     ! mynode return the rank of proc (0 --> jpnij -1 )
 …
       CALL iom_close                                 ! close all input/output files managed by iom_*
+      !
+      IF( numstp      /= -1 )   CLOSE( numstp      )   ! time-step file
+      IF( numnam      /= -1 )   CLOSE( numnam      )   ! oce namelist
+      IF( numnam_ice  /= -1 )   CLOSE( numnam_ice  )   ! ice namelist
+      IF( numevo_ice  /= -1 )   CLOSE( numevo_ice  )   ! ice variables (temp. evolution)
+      IF( numout      /=  6 )   CLOSE( numout      )   ! standard model output file
+      IF( numstp          /= -1 )   CLOSE( numstp         )   ! time-step file
+      IF( numnam_ref      /= -1 )   CLOSE( numnam_ref     )   ! oce reference namelist
+      IF( numnam_cfg      /= -1 )   CLOSE( numnam_cfg     )   ! oce configuration namelist
+      IF( numnam_ice_ref  /= -1 )   CLOSE( numnam_ice_ref )   ! ice reference namelist
+      IF( numnam_ice_cfg  /= -1 )   CLOSE( numnam_ice_cfg )   ! ice configuration namelist
+      IF( numevo_ice      /= -1 )   CLOSE( numevo_ice     )   ! ice variables (temp. evolution)
+      IF( numout          /=  6 )   CLOSE( numout         )   ! standard model output file
+      !
       numout = 6                                     ! redefine numout in case it is used after this point...

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/SAS_SRC/sbcssm.F90

-                      r3364
+                      r3875
       INTEGER  :: ifpr                               ! dummy loop indice
       INTEGER  :: inum, idv, idimv, jpm              ! local integer
+      INTEGER  :: ios                                ! Local integer output status for namelist read
       !!
       CHARACTER(len=100)                     ::  cn_dir       ! Root directory for location of core files
 …
       sn_ssh  = FLD_N( 'ssm_grid_T' ,    120    , 'sossheig' ,  .true.    , .true. ,   'yearly'  , ''       , ''         )
+      !
+      REWIND( numnam )                          ! read in namlist namsbc_ssm
+      READ  ( numnam, namsbc_ssm )
+      REWIND( numnam_ref )              ! Namelist namsbc_ssm in reference namelist : SAS files for dynamical data
+      READ  ( numnam_ref, namsbc_ssm, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_ssm in reference namelist', lwp )
+      REWIND( numnam_cfg )              ! Namelist namsbc_ssm in configuration namelist : SAS files for dynamical data
+      READ  ( numnam_cfg, namsbc_ssm, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_ssm in configuration namelist', lwp )
+      WRITE ( numond, namsbc_ssm )
       !                                         ! store namelist information in an array
       !                                         ! Control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/C14b/trcnam_c14b.F90

-                      r3294
+                      r3875
       !! ** input   :   Namelist namelist_c14b
       !!----------------------------------------------------------------------
+      INTEGER ::   numnatb
+      INTEGER ::  numnatb_ref = -1   ! Logical unit for reference c14b namelist
+      INTEGER ::  numnatb_cfg = -1   ! Logical unit for configuration c14b namelist
+      INTEGER ::  numonb      = -1   ! Logical unit for output namelist
+      INTEGER :: ios                 ! Local integer output status for namelist read
       ! definition of additional diagnostic as a structure
 …
       !                             ! Open namelist file
+      CALL ctl_opn( numnatb, 'namelist_c14b', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      READ( numnatb , namc14date )     ! read namelist
+      CALL ctl_opn( numnatb_ref, 'namelist_c14b_ref'  ,     'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      CALL ctl_opn( numnatb_cfg, 'namelist_c14b_cfg'  ,     'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      CALL ctl_opn( numonb     , 'output.namelist.c14', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      REWIND( numnatb_ref )              ! Namelist namc14date in reference namelist : c14b parameters
+      READ  ( numnatb_ref, namc14date, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namc14date in reference namelist', lwp )
+      REWIND( numnatb_cfg )              ! Namelist namc14date in configuration namelist : c14b parameters
+      READ  ( numnatb_cfg, namc14date, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namc14date in configuration namelist', lwp )
+      WRITE ( numonb, namc14date )
       IF(lwp) THEN                  ! control print
          WRITE(numout,*)
 …
          END DO
+         REWIND( numnatb )               !
+         READ  ( numnatb, namc14dia )
+         REWIND( numnatb_ref )              ! Namelist namc14dia in reference namelist : c14b diagnostics
+         READ  ( numnatb_ref, namc14dia, IOSTAT = ios, ERR = 903)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namc14dia in reference namelist', lwp )
+         REWIND( numnatb_cfg )              ! Namelist namc14dia in configuration namelist : c14b diagnostics
+         READ  ( numnatb_cfg, namc14dia, IOSTAT = ios, ERR = 904 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namc14dia in configuration namelist', lwp )
+         WRITE ( numonb, namc14dia )
          DO jl = 1, jp_c14b_2d

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/CFC/trcnam_cfc.F90

-                      r3294
+                      r3875
       !! ** input   :   Namelist namcfc
       !!----------------------------------------------------------------------
+      INTEGER ::  numnatc
+      INTEGER ::  numnatc_ref = -1   ! Logical unit for reference CFC namelist
+      INTEGER ::  numnatc_cfg = -1   ! Logical unit for configuration CFC namelist
+      INTEGER ::  numonc      = -1   ! Logical unit for output namelist
+      INTEGER :: ios                 ! Local integer output status for namelist read
       INTEGER :: jl, jn
       TYPE(DIAG), DIMENSION(jp_cfc_2d) :: cfcdia2d
 …
       nyear_res = 1950
+      !                             ! Open namelist file
+      CALL ctl_opn( numnatc, 'namelist_cfc', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      READ( numnatc , namcfcdate )     ! read namelist
+      !                             ! Open namelist files
+      CALL ctl_opn( numnatc_ref, 'namelist_cfc_ref'   ,     'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      CALL ctl_opn( numnatc_cfg, 'namelist_cfc_cfg'   ,     'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      CALL ctl_opn( numonc     , 'output.namelist.cfc', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      REWIND( numnatc_ref )              ! Namelist namcfcdate in reference namelist : CFC parameters
+      READ  ( numnatc_ref, namcfcdate, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcfcdate in reference namelist', lwp )
+      REWIND( numnatc_cfg )              ! Namelist namcfcdate in configuration namelist : CFC parameters
+      READ  ( numnatc_cfg, namcfcdate, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcfcdate in configuration namelist', lwp )
+      WRITE ( numonc, namcfcdate )
       IF(lwp) THEN                  ! control print
 …
          END DO
+         REWIND( numnatc )               ! read natrtd
+         READ  ( numnatc, namcfcdia )
+         REWIND( numnatc_ref )              ! Namelist namcfcdia in reference namelist : CFC diagnostics
+         READ  ( numnatc_ref, namcfcdia, IOSTAT = ios, ERR = 903)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcfcdia in reference namelist', lwp )
+         REWIND( numnatc_cfg )              ! Namelist namcfcdia in configuration namelist : CFC diagnostics
+         READ  ( numnatc_cfg, namcfcdia, IOSTAT = ios, ERR = 904 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcfcdia in configuration namelist', lwp )
+         WRITE ( numonc, namcfcdia )
          DO jl = 1, jp_cfc_2d

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zbio.F90

-                      r3446
+                      r3875
       !!                  ***  ROUTINE p2z_bio_init  ***
       !!
       !! ** Purpose :  bilogical parameters
+      !! ** Purpose :  biological parameters
       !!
       !! ** Method  :   Read namelist and check the parameters
 …
       NAMELIST/namlobdet/  taudn, fdetlab
       NAMELIST/namlobdom/ taudomn
+      INTEGER :: ios                 ! Local integer output status for namelist read
       !!----------------------------------------------------------------------
+      REWIND( numnatp )
+      READ  ( numnatp, namlobphy )
+      REWIND( numnatp_ref )              ! Namelist namlobphy in reference namelist : Lobster biological parameters
+      READ  ( numnatp_ref, namlobphy, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobphy in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist namlobphy in configuration namelist : Lobster biological parameters
+      READ  ( numnatp_cfg, namlobphy, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobphy in configuration namelist', lwp )
+      WRITE ( numonp, namlobphy )
       IF(lwp) THEN
 …
       ENDIF
+      REWIND( numnatp )
+      READ  ( numnatp, namlobnut )
+      REWIND( numnatp_ref )              ! Namelist namlobnut in reference namelist : Lobster nutriments parameters
+      READ  ( numnatp_ref, namlobnut, IOSTAT = ios, ERR = 903)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobnut in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist namlobnut in configuration namelist : Lobster nutriments parameters
+      READ  ( numnatp_cfg, namlobnut, IOSTAT = ios, ERR = 904 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobnut in configuration namelist', lwp )
+      WRITE ( numonp, namlobnut )
       IF(lwp) THEN
           WRITE(numout,*) ' Namelist namlobnut'
 …
       ENDIF
+      REWIND( numnatp )
+      READ  ( numnatp, namlobzoo )
+      REWIND( numnatp_ref )              ! Namelist namlobzoo in reference namelist : Lobster zooplankton parameters
+      READ  ( numnatp_ref, namlobzoo, IOSTAT = ios, ERR = 905)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobzoo in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist namlobzoo in configuration namelist : Lobster zooplankton parameters
+      READ  ( numnatp_cfg, namlobzoo, IOSTAT = ios, ERR = 906 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobzoo in configuration namelist', lwp )
+      WRITE ( numonp, namlobzoo )
       IF(lwp) THEN
 …
       ENDIF
+      REWIND( numnatp )
+      READ  ( numnatp, namlobdet )
+      REWIND( numnatp_ref )              ! Namelist namlobdet in reference namelist : Lobster detritus parameters
+      READ  ( numnatp_ref, namlobdet, IOSTAT = ios, ERR = 907)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobdet in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist namlobdet in configuration namelist : Lobster detritus parameters
+      READ  ( numnatp_cfg, namlobdet, IOSTAT = ios, ERR = 908 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobdet in configuration namelist', lwp )
+      WRITE ( numonp, namlobdet )
       IF(lwp) THEN
 …
       ENDIF
+      REWIND( numnatp )
+      READ  ( numnatp, namlobdom )
+      REWIND( numnatp_ref )              ! Namelist namlobdom in reference namelist : Lobster DOM breakdown rate
+      READ  ( numnatp_ref, namlobdom, IOSTAT = ios, ERR = 909)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobdom in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist namlobdom in configuration namelist : Lobster DOM breakdown rate
+      READ  ( numnatp_cfg, namlobdom, IOSTAT = ios, ERR = 910 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobdom in configuration namelist', lwp )
+      WRITE ( numonp, namlobdom )
       IF(lwp) THEN

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zopt.F90

-                      r3443
+                      r3875
       NAMELIST/namlobopt/ xkg0, xkr0, xkgp, xkrp, xlg, xlr, rpig
       NAMELIST/namlobrat/ rcchl, redf, reddom
+      INTEGER :: ios                 ! Local integer output status for namelist read
       !!----------------------------------------------------------------------
+      REWIND( numnatp )
+      READ  ( numnatp, namlobopt )
+      REWIND( numnatp_ref )              ! Namelist namlobopt in reference namelist : Lobster options
+      READ  ( numnatp_ref, namlobopt, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobopt in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist namlobopt in configuration namelist : Lobster options
+      READ  ( numnatp_cfg, namlobopt, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobopt in configuration namelist', lwp )
+      WRITE ( numonp, namlobopt )
       IF(lwp) THEN
 …
       ENDIF
+      !
+      REWIND( numnatp )
+      READ  ( numnatp, namlobrat )
+      REWIND( numnatp_ref )              ! Namelist namlobrat in reference namelist : Lobster ratios
+      READ  ( numnatp_ref, namlobrat, IOSTAT = ios, ERR = 903)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobrat in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist namlobrat in configuration namelist : Lobster ratios
+      READ  ( numnatp_cfg, namlobrat, IOSTAT = ios, ERR = 904 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobrat in configuration namelist', lwp )
+      WRITE ( numonp, namlobrat )
       IF(lwp) THEN

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsed.F90

-                      r3443
+                      r3875
       !!----------------------------------------------------------------------
       NAMELIST/namlobsed/ sedlam, sedlostpoc, vsed, xhr
+      INTEGER :: ios                 ! Local integer output status for namelist read
+      REWIND( numnatp )
+      READ  ( numnatp, namlobsed )
+      REWIND( numnatp_ref )              ! Namelist namlobsed in reference namelist : Lobster sediments
+      READ  ( numnatp_ref, namlobsed, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlosed in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist namlobsed in configuration namelist : Lobster sediments
+      READ  ( numnatp_cfg, namlobsed, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobsed in configuration namelist', lwp )
+      WRITE ( numonp, namlobsed )
       IF(lwp) THEN

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zfechem.F90

-                      r3780
+                      r3875
       !!----------------------------------------------------------------------
       NAMELIST/nampisfer/ ln_fechem, ln_ligvar, xlam1, xlamdust, ligand
+      REWIND( numnatp )                     ! read numnatp
+      READ  ( numnatp, nampisfer )
+      INTEGER :: ios                 ! Local integer output status for namelist read
+      REWIND( numnatp_ref )              ! Namelist nampisfer in reference namelist : Pisces iron chemistry
+      READ  ( numnatp_ref, nampisfer, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisfer in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampisfer in configuration namelist : Pisces iron chemistry
+      READ  ( numnatp_cfg, nampisfer, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisfer in configuration namelist', lwp )
+      WRITE ( numonp, nampisfer )
       IF(lwp) THEN                         ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90

-                      r3481
+                      r3875
       NAMELIST/nampisext/ln_co2int, atcco2, clname, nn_offset
       INTEGER :: jm
+      !!----------------------------------------------------------------------
+      !
+      REWIND( numnatp )                     ! read numnatp
+      READ  ( numnatp, nampisext )
+      INTEGER :: ios                 ! Local integer output status for namelist read
+      !!----------------------------------------------------------------------
+      !
+      REWIND( numnatp_ref )              ! Namelist nampisext in reference namelist : Pisces atm. conditions
+      READ  ( numnatp_ref, nampisext, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisext in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampisext in configuration namelist : Pisces atm. conditions
+      READ  ( numnatp_cfg, nampisext, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisext in configuration namelist', lwp )
+      WRITE ( numonp, nampisext )
+      !
       IF(lwp) THEN                         ! control print
 …
+      !
       INTEGER            ::  ierr
+      INTEGER            ::  ios      ! Local integer output status for namelist read
       CHARACTER(len=100) ::  cn_dir   ! Root directory for location of ssr files
       TYPE(FLD_N)        ::  sn_patm  ! informations about the fields to be read
 …
          cn_dir  = './'          ! directory in which the Patm data are
+         REWIND( numnatp )                             !* read in namlist nampisatm
+         READ  ( numnatp, nampisatm )
+         REWIND( numnatp_ref )              ! Namelist nampisatm in reference namelist : Pisces atm. sea level pressure file
+         READ  ( numnatp_ref, nampisatm, IOSTAT = ios, ERR = 901)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisatm in reference namelist', lwp )
+         REWIND( numnatp_cfg )              ! Namelist nampisatm in configuration namelist : Pisces atm. sea level pressure file
+         READ  ( numnatp_cfg, nampisatm, IOSTAT = ios, ERR = 902 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisatm in configuration namelist', lwp )
+         WRITE ( numonp, nampisatm )
+         !
+         !

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlim.F90

-                      r3446
+                      r3875
          &                concbno3, concbnh4, xsizedia, xsizephy, xsizern, xsizerd,          &
          &                xksi1, xksi2, xkdoc, qnfelim, qdfelim, caco3r
+      REWIND( numnatp )                     ! read numnat
+      READ  ( numnatp, nampislim )
+      INTEGER :: ios                 ! Local integer output status for namelist read
+      REWIND( numnatp_ref )              ! Namelist nampislim in reference namelist : Pisces nutrient limitation parameters
+      READ  ( numnatp_ref, nampislim, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampislim in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampislim in configuration namelist : Pisces nutrient limitation parameters
+      READ  ( numnatp_cfg, nampislim, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampislim in configuration namelist', lwp )
+      WRITE ( numonp, nampislim )
       IF(lwp) THEN                         ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlys.F90

-                      r3443
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER  ::  ji, jj, jk
+      INTEGER  ::  ios                 ! Local integer output status for namelist read
       REAL(wp) ::  zcaralk, zbicarb, zco3
       REAL(wp) ::  ztmas, ztmas1
 …
       !!----------------------------------------------------------------------
+      REWIND( numnatp )                     ! read numnatp
+      READ  ( numnatp, nampiscal )
+      REWIND( numnatp_ref )              ! Namelist nampiscal in reference namelist : Pisces CaCO3 dissolution
+      READ  ( numnatp_ref, nampiscal, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampiscal in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampiscal in configuration namelist : Pisces CaCO3 dissolution
+      READ  ( numnatp_cfg, nampiscal, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampiscal in configuration namelist', lwp )
+      WRITE ( numonp, nampiscal )
       IF(lwp) THEN                         ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90

-                      r3829
+                      r3875
          &                xprefpoc, xthresh2dia, xthresh2phy, xthresh2zoo, xthresh2poc, &
          &                xthresh2, xkgraz2, epsher2, sigma2, unass2, grazflux
+      REWIND( numnatp )                     ! read numnatp
+      READ  ( numnatp, nampismes )
+      INTEGER :: ios                 ! Local integer output status for namelist read
+      REWIND( numnatp_ref )              ! Namelist nampismes in reference namelist : Pisces mesozooplankton
+      READ  ( numnatp_ref, nampismes, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampismes in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampismes in configuration namelist : Pisces mesozooplankton
+      READ  ( numnatp_cfg, nampismes, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampismes in configuration namelist', lwp )
+      WRITE ( numonp, nampismes )

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90

-                      r3829
+                      r3875
          &                xpref2d,  xthreshdia,  xthreshphy,  xthreshpoc, &
          &                xthresh, xkgraz, epsher, sigma1, unass
+      REWIND( numnatp )                     ! read numnatp
+      READ  ( numnatp, nampiszoo )
+      INTEGER :: ios                 ! Local integer output status for namelist read
+      REWIND( numnatp_ref )              ! Namelist nampiszoo in reference namelist : Pisces microzooplankton
+      READ  ( numnatp_ref, nampiszoo, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampiszoo in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampiszoo in configuration namelist : Pisces microzooplankton
+      READ  ( numnatp_cfg, nampiszoo, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampiszoo in configuration namelist', lwp )
+      WRITE ( numonp, nampiszoo )
       IF(lwp) THEN                         ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmort.F90

-                      r3475
+                      r3875
       NAMELIST/nampismort/ wchl, wchld, wchldm, mprat, mprat2
+      REWIND( numnatp )                     ! read numnatp
+      READ  ( numnatp, nampismort )
+      INTEGER :: ios                 ! Local integer output status for namelist read
+      REWIND( numnatp_ref )              ! Namelist nampismort in reference namelist : Pisces phytoplankton
+      READ  ( numnatp_ref, nampismort, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampismort in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampismort in configuration namelist : Pisces phytoplankton
+      READ  ( numnatp_cfg, nampismort, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampismort in configuration namelist', lwp )
+      WRITE ( numonp, nampismort )
       IF(lwp) THEN                         ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zopt.F90

-                      r3475
+                      r3875
       INTEGER :: numpar
       INTEGER :: ierr
+      INTEGER :: ios                 ! Local integer output status for namelist read
       REAL(wp), DIMENSION(nbtimes) :: zsteps                 ! times records
+      !
 …
       sn_par     = FLD_N( 'par_fraction',    24     ,  'fr_par'     ,  .true.    , .true.  ,   'yearly'  , ''       , ''         )
+      REWIND( numnatp )                     ! read numnatp
+      READ  ( numnatp, nampisopt )
+      REWIND( numnatp_ref )              ! Namelist nampisopt in reference namelist : Pisces attenuation coef. and PAR
+      READ  ( numnatp_ref, nampisopt, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisopt in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampisopt in configuration namelist : Pisces attenuation coef. and PAR
+      READ  ( numnatp_cfg, nampisopt, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisopt in configuration namelist', lwp )
+      WRITE ( numonp, nampisopt )
       IF(lwp) THEN

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

-                      r3686
+                      r3875
       NAMELIST/nampisprod/ pislope, pislope2, ln_newprod, bresp, excret, excret2,  &
          &                 chlcnm, chlcdm, chlcmin, fecnm, fecdm, grosip
+      INTEGER :: ios                 ! Local integer output status for namelist read
       !!----------------------------------------------------------------------
+      REWIND( numnatp )                     ! read numnatp
+      READ  ( numnatp, nampisprod )
+      REWIND( numnatp_ref )              ! Namelist nampisprod in reference namelist : Pisces phytoplankton production
+      READ  ( numnatp_ref, nampisprod, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisprod in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampisprod in configuration namelist : Pisces phytoplankton production
+      READ  ( numnatp_cfg, nampisprod, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisprod in configuration namelist', lwp )
+      WRITE ( numonp, nampisprod )
       IF(lwp) THEN                         ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90

-                      r3496
+                      r3875
       NAMELIST/nampisrem/ xremik, xremip, nitrif, xsirem, xsiremlab, xsilab,   &
       &                   oxymin
+      REWIND( numnatp )                     ! read numnatp
+      READ  ( numnatp, nampisrem )
+      INTEGER :: ios                 ! Local integer output status for namelist read
+      REWIND( numnatp_ref )              ! Namelist nampisrem in reference namelist : Pisces remineralization
+      READ  ( numnatp_ref, nampisrem, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisrem in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampisrem in configuration namelist : Pisces remineralization
+      READ  ( numnatp_cfg, nampisrem, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisrem in configuration namelist', lwp )
+      WRITE ( numonp, nampisrem )
       IF(lwp) THEN                         ! control print

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90

-                      r3557
+                      r3875
       INTEGER  :: numdust, numsolub, numriv, numiron, numdepo, numhydro
       INTEGER  :: ierr, ierr1, ierr2, ierr3
+      INTEGER  :: ios                 ! Local integer output status for namelist read
       REAL(wp) :: zexpide, zdenitide, zmaskt
       REAL(wp) :: ztimes_dust, ztimes_riv, ztimes_ndep
 …
       sn_hydrofe  = FLD_N( 'hydrofe'    ,   -12     ,  'hydro'    ,  .false.   , .true.  ,   'yearly'  , ''       , ''         )
+      REWIND( numnatp )                     ! read numnatp
+      READ  ( numnatp, nampissbc )
+      REWIND( numnatp_ref )              ! Namelist nampissbc in reference namelist : Pisces external sources of nutrients
+      READ  ( numnatp_ref, nampissbc, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampissbc in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampissbc in configuration namelist : Pisces external sources of nutrients
+      READ  ( numnatp_cfg, nampissbc, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampissbc in configuration namelist', lwp )
+      WRITE ( numonp, nampissbc )
       IF(lwp) THEN

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsink.F90

-                      r3829
+                      r3875
       !!----------------------------------------------------------------------
       INTEGER  ::   jk, jn, kiter
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
       REAL(wp) ::   znum, zdiv
       REAL(wp) ::   zws, zwr, zwl,wmax, znummax
 …
       IF( nn_timing == 1 )  CALL timing_start('p4z_sink_init')
+      !
+      REWIND( numnatp )                     ! read nampiskrs
+      READ  ( numnatp, nampiskrs )
+      REWIND( numnatp_ref )              ! Namelist nampiskrs in reference namelist : Pisces sinking Kriest
+      READ  ( numnatp_ref, nampiskrs, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampiskrs in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampiskrs in configuration namelist : Pisces sinking Kriest
+      READ  ( numnatp_cfg, nampiskrs, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampiskrs in configuration namelist', lwp )
+      WRITE ( numonp, nampiskrs )
       IF(lwp) THEN

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90

-                      r3780
+                      r3875
       NAMELIST/nampisdmp/ ln_pisdmp, nn_pisdmp, ln_pisclo
       NAMELIST/nampismass/ ln_check_mass
+      !!----------------------------------------------------------------------
+      REWIND( numnatp )
+      READ  ( numnatp, nampisbio )
+      INTEGER :: ios                 ! Local integer output status for namelist read
+      !!----------------------------------------------------------------------
+      REWIND( numnatp_ref )              ! Namelist nampisbio in reference namelist : Pisces variables
+      READ  ( numnatp_ref, nampisbio, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisbio in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampisbio in configuration namelist : Pisces variables
+      READ  ( numnatp_cfg, nampisbio, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisbio in configuration namelist', lwp )
+      WRITE ( numonp, nampisbio )
       IF(lwp) THEN                         ! control print
 …
       xkr_mass_max = 1.
+      REWIND( numnatp )                     ! read natkriest
+      READ  ( numnatp, nampiskrp )
+      REWIND( numnatp_ref )              ! Namelist nampiskrp in reference namelist : Pisces Kriest
+      READ  ( numnatp_ref, nampiskrp, IOSTAT = ios, ERR = 903)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampiskrp in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampiskrp in configuration namelist : Pisces Kriest
+      READ  ( numnatp_cfg, nampiskrp, IOSTAT = ios, ERR = 904 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampiskrp in configuration namelist', lwp )
+      WRITE ( numonp, nampiskrp )
       IF(lwp) THEN
 …
       ln_pisclo = .false.
+      REWIND( numnatp )
+      READ  ( numnatp, nampisdmp )
+      REWIND( numnatp_ref )              ! Namelist nampisdmp in reference namelist : Pisces damping
+      READ  ( numnatp_ref, nampisdmp, IOSTAT = ios, ERR = 905)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisdmp in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampisdmp in configuration namelist : Pisces damping
+      READ  ( numnatp_cfg, nampisdmp, IOSTAT = ios, ERR = 906 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisdmp in configuration namelist', lwp )
+      WRITE ( numonp, nampisdmp )
       IF(lwp) THEN                         ! control print
 …
       ln_check_mass = .false.
+      REWIND( numnatp )
+      READ  ( numnatp, nampismass )
+      REWIND( numnatp_ref )              ! Namelist nampismass in reference namelist : Pisces mass conservation check
+      READ  ( numnatp_ref, nampismass, IOSTAT = ios, ERR = 907)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampismass in reference namelist', lwp )
+      REWIND( numnatp_cfg )              ! Namelist nampismass in configuration namelist : Pisces mass conservation check
+      READ  ( numnatp_cfg, nampismass, IOSTAT = ios, ERR = 908 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampismass in configuration namelist', lwp )
+      WRITE ( numonp, nampismass )
       IF(lwp) THEN                         ! control print
          WRITE(numout,*) ' '

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/sms_pisces.F90

-                      r3780
+                      r3875
    PUBLIC
+   INTEGER ::   numnatp
+   INTEGER ::   numnatp_ref = -1           !! Logical units for namelist pisces
+   INTEGER ::   numnatp_cfg = -1           !! Logical units for namelist pisces
+   INTEGER ::   numonp      = -1           !! Logical unit for namelist pisces output
    !!*  Biological fluxes for light : variables shared by pisces & lobster

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/trcnam_pisces.F90

-                      r3680
+                      r3875
       !!
       INTEGER :: jl, jn
+      INTEGER :: ios                 ! Local integer output status for namelist read
       TYPE(DIAG), DIMENSION(jp_pisces_2d)  :: pisdia2d
       TYPE(DIAG), DIMENSION(jp_pisces_3d)  :: pisdia3d
 …
 #endif
       IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~'
+      CALL ctl_opn( numnatp, TRIM( clname ), 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      CALL ctl_opn( numnatp_ref, TRIM( clname )//'_ref', 'OLD'    , 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      CALL ctl_opn( numnatp_cfg, TRIM( clname )//'_cfg', 'OLD'    , 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      CALL ctl_opn( numonp     , 'output.namelist.pis' , 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
+      !
       IF( .NOT.lk_iomput .AND. ln_diatrc ) THEN
 …
          END DO
+         REWIND( numnatp )               !
+         READ  ( numnatp, nampisdia )
+         REWIND( numnatp_ref )              ! Namelist nampisdia in reference namelist : Pisces diagnostics
+         READ  ( numnatp_ref, nampisdia, IOSTAT = ios, ERR = 901)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisdia in reference namelist', lwp )
+         REWIND( numnatp_cfg )              ! Namelist nampisdia in configuration namelist : Pisces diagnostics
+         READ  ( numnatp_cfg, nampisdia, IOSTAT = ios, ERR = 902 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisdia in configuration namelist', lwp )
+         WRITE ( numonp, nampisdia )
          DO jl = 1, jp_pisces_2d
 …
          END DO
+         REWIND( numnatp )
+         READ  ( numnatp, nampisdbi )
+         REWIND( numnatp_ref )              ! Namelist nampisdbi in reference namelist : Pisces add. diagnostics
+         READ  ( numnatp_ref, nampisdbi, IOSTAT = ios, ERR = 903)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisdbi in reference namelist', lwp )
+         REWIND( numnatp_cfg )              ! Namelist nampisdbi in configuration namelist : Pisces add. diagnostics
+         READ  ( numnatp_cfg, nampisdbi, IOSTAT = ios, ERR = 904 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisdbi in configuration namelist', lwp )
+         WRITE ( numonp, nampisdbi )
          DO jl = 1, jp_pisces_trd

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/TRP/trcnam_trp.F90

-                      r3718
+                      r3875
    USE trc                 ! ocean passive tracers variables
    USE in_out_manager      ! ocean dynamics and active tracers variables
+   USE lib_mpp           ! distributed memory computing library
    IMPLICIT NONE
 …
       !! ** Purpose :   Read Namelist for tracer transport option
       !!----------------------------------------------------------------------
+      INTEGER ::  ios                 ! Local integer output status for namelist read
       NAMELIST/namtrc_adv/ ln_trcadv_cen2 , ln_trcadv_tvd   ,    &
          &                 ln_trcadv_muscl, ln_trcadv_muscl2,    &
 …
       IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~'
+      REWIND ( numnat )               ! Read Namelist namtrc_adv : tracer advection scheme
+      READ   ( numnat, namtrc_adv )
+      REWIND( numnat_ref )              ! Namelist namtrc_adv in reference namelist : Tracer advection scheme
+      READ  ( numnat_ref, namtrc_adv, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_adv in reference namelist', lwp )
+      REWIND( numnat_cfg )              ! Namelist namtrc_adv in configuration namelist : Tracer advection scheme
+      READ  ( numnat_cfg, namtrc_adv, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_adv in configuration namelist', lwp )
+      WRITE ( numont, namtrc_adv )
       IF(lwp) THEN                    ! Namelist print
 …
       ENDIF
+      !
+      REWIND( numnat )                ! Namelist namtrc_ldf
+      READ  ( numnat, namtrc_ldf )
+      REWIND( numnat_ref )              ! Namelist namtrc_ldf in reference namelist : Tracer lateral diffusive operator
+      READ  ( numnat_ref, namtrc_ldf, IOSTAT = ios, ERR = 903)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_ldf in reference namelist', lwp )
+      REWIND( numnat_cfg )              ! Namelist namtrc_ldf in configuration namelist : Tracer lateral diffusive operator
+      READ  ( numnat_cfg, namtrc_ldf, IOSTAT = ios, ERR = 904 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_ldf in configuration namelist', lwp )
+      WRITE ( numont, namtrc_ldf )
       IF(lwp) THEN                    ! Namelist print
 …
       !                                ! Vertical mixing
+      REWIND( numnat )                 !   ! Read namtopzdf namelist
+      READ  ( numnat, namtrc_zdf )
+      REWIND( numnat_ref )              ! Namelist namtrc_zdf in reference namelist : Tracer vertical mixing
+      READ  ( numnat_ref, namtrc_zdf, IOSTAT = ios, ERR = 905)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_zdf in reference namelist', lwp )
+      REWIND( numnat_cfg )              ! Namelist namtrc_zdf in configuration namelist : Tracer vertical mixing
+      READ  ( numnat_cfg, namtrc_zdf, IOSTAT = ios, ERR = 906 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_zdf in configuration namelist', lwp )
+      WRITE ( numont, namtrc_zdf )
       IF(lwp) THEN                     !   ! Control print
 …
+      !
+      REWIND( numnat )                 !   Read Namelist namtoprad
+      READ  ( numnat, namtrc_rad )
+      REWIND( numnat_ref )              ! Namelist namtrc_rad in reference namelist : Tracer negative concentrations
+      READ  ( numnat_ref, namtrc_rad, IOSTAT = ios, ERR = 907)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_rad in reference namelist', lwp )
+      REWIND( numnat_cfg )              ! Namelist namtrc_rad in configuration namelist : Tracer negative concentrations
+      READ  ( numnat_cfg, namtrc_rad, IOSTAT = ios, ERR = 908 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_rad in configuration namelist', lwp )
+      WRITE ( numont, namtrc_rad )
       IF(lwp) THEN                     !   ! Control print
 …
 # if defined key_trcdmp
+      REWIND ( numnat )                  ! Read Namelist namtra_dmp : temperature and salinity damping term
+      READ   ( numnat, namtrc_dmp )
+      REWIND( numnat_ref )              ! Namelist namtrc_dmp in reference namelist : Passive tracers newtonian damping
+      READ  ( numnat_ref, namtrc_dmp, IOSTAT = ios, ERR = 909)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_dmp in reference namelist', lwp )
+      REWIND( numnat_cfg )              ! Namelist namtrc_dmp in configuration namelist : Passive tracers newtonian damping
+      READ  ( numnat_cfg, namtrc_dmp, IOSTAT = ios, ERR = 910)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_dmp in configuration namelist', lwp )
+      WRITE ( numont, namtrc_dmp )
       IF( lzoom )   nn_zdmp_trc = 0           ! restoring to climatology at closed north or south boundaries

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/trc.F90

-                      r3770
+                      r3875
    !! parameters for the control of passive tracers
+   !! --------------------------------------------------
+   INTEGER, PUBLIC                                                 ::   numnat        !: logicla unit for the passive tracer NAMELIST
+   !! ---------------------------------------------
+   INTEGER, PUBLIC                                                 ::   numnat_ref = -1   !: logical unit for the reference passive tracer namelist_top_ref
+   INTEGER, PUBLIC                                                 ::   numnat_cfg = -1   !: logical unit for the reference passive tracer namelist_top_cfg
+   INTEGER, PUBLIC                                                 ::   numont     = -1   !: logical unit for the reference passive tracer namelist output output.namelist.top
    INTEGER, PUBLIC                                                 ::   numstr        !: logical unit for tracer statistics
    INTEGER, PUBLIC                                                 ::   numrtr        !: logical unit for trc restart (read )

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/trcdta.F90

-                      r3827
+                      r3875
       !!----------------------------------------------------------------------
+      !
       INTEGER            :: jl, jn                   ! dummy loop indicies
+      INTEGER            :: jl, jn                           ! dummy loop indicies
       INTEGER            :: ierr0, ierr1, ierr2, ierr3       ! temporary integers
+      INTEGER            :: ios                              ! Local integer output status for namelist read
       CHARACTER(len=100) :: clndta, clntrc
       REAL(wp)           :: zfact
 …
       END DO
+      !
+      REWIND( numnat )               ! read nattrc
+      READ  ( numnat, namtrc_dta )
+      REWIND( numnat_ref )              ! Namelist namtrc_dta in reference namelist : Passive tracer input data
+      READ  ( numnat_ref, namtrc_dta, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_dta in reference namelist', lwp )
+      REWIND( numnat_cfg )              ! Namelist namtrc_dta in configuration namelist : Passive tracer input data
+      READ  ( numnat_cfg, namtrc_dta, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_dta in configuration namelist', lwp )
+      WRITE ( numont, namtrc_dta )
       IF( lwp ) THEN

branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/trcnam.F90

-                      r3749
+                      r3875
       !!---------------------------------------------------------------------
       INTEGER ::  jn, ierr
+      INTEGER ::  ios                 ! Local integer output status for namelist read
       ! Definition of a tracer as a structure
       TYPE(PTRACER), DIMENSION(jptra) :: sn_tracer  ! type of tracer for saving if not key_iomput
 …
       IF(lwp) WRITE(numout,*) '~~~~~~~'
+      CALL ctl_opn( numnat, 'namelist_top', 'OLD', 'FORMATTED', 'SEQUENTIAL', 1, numout, .FALSE. )
+      CALL ctl_opn( numnat_ref, 'namelist_top_ref'   , 'OLD'    , 'FORMATTED', 'SEQUENTIAL', 1, numout, .FALSE. )
+      CALL ctl_opn( numnat_cfg, 'namelist_top_cfg'   , 'OLD'    , 'FORMATTED', 'SEQUENTIAL', 1, numout, .FALSE. )
+      CALL ctl_opn( numont    , 'output.namelist.top', 'REPLACE', 'FORMATTED', 'SEQUENTIAL',-1, numout, .FALSE. )
       ! Namelist nattrc (files)
 …
       ln_trcdmp = .FALSE.
+      REWIND( numnat )               ! read nattrc
+      READ  ( numnat, namtrc )
+      REWIND( numnat_ref )              ! Namelist namtrc in reference namelist : Passive tracer variables
+      READ  ( numnat_ref, namtrc, IOSTAT = ios, ERR = 901)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc in reference namelist', lwp )
+      REWIND( numnat_cfg )              ! Namelist namtrc in configuration namelist : Passive tracer variables
+      READ  ( numnat_cfg, namtrc, IOSTAT = ios, ERR = 902 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc in configuration namelist', lwp )
+      WRITE ( numont, namtrc )
       DO jn = 1, jptra
 …
       nn_writebio = 10
+      REWIND( numnat )               !  namelist namtoptrd : passive tracer trends diagnostic
+      READ  ( numnat, namtrc_dia )
+      REWIND( numnat_ref )              ! Namelist namtrc_dia in reference namelist : Passive tracer diagnostics
+      READ  ( numnat_ref, namtrc_dia, IOSTAT = ios, ERR = 903)
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_dia in reference namelist', lwp )
+      REWIND( numnat_cfg )              ! Namelist namtrc_dia in configuration namelist : Passive tracer diagnostics
+      READ  ( numnat_cfg, namtrc_dia, IOSTAT = ios, ERR = 904 )
+   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_dia in configuration namelist', lwp )
+      WRITE ( numont, namtrc_dia )
       IF(lwp) THEN
 …
          ln_trdtrc(:) = .FALSE.
+         REWIND( numnat )               !  namelist namtoptrd : passive tracer trends diagnostic
+         READ  ( numnat, namtrc_trd )
+         REWIND( numnat_ref )              ! Namelist namtrc_trd in reference namelist : Passive tracer trends
+         READ  ( numnat_ref, namtrc_trd, IOSTAT = ios, ERR = 905)
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_trd in reference namelist', lwp )
+         REWIND( numnat_cfg )              ! Namelist namtrc_trd in configuration namelist : Passive tracer trends
+         READ  ( numnat_cfg, namtrc_trd, IOSTAT = ios, ERR = 906 )
+      IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_trd in configuration namelist', lwp )
+         WRITE ( numont, namtrc_trd )
          IF(lwp) THEN

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

Context Navigation

Legend:

Download in other formats: