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

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

Ignore:
Timestamp:
2008-06-20T17:17:41+02:00 (16 years ago)
Author:
cetlod
Message:

style of all top namelist has been modified ; update modules to take it into account, see ticket:196

Location:
trunk/NEMO/TOP_SRC/PISCES
Files:
18 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/TOP_SRC/PISCES/p4zbio.F90

    r1073 r1119  
    1515   !!---------------------------------------------------------------------- 
    1616   USE oce_trc         ! 
    17    USE trp_trc         !  
     17   USE trc         !  
    1818   USE sms_pisces      !  
    1919   USE p4zsink         !  
    … …  
    131131 
    132132# if defined key_trc_dia3d && defined key_kriest 
    133 !!gm potential bug  hard coded index on trc3d 
    134       trc3d(:,:,:,11) = tra(:,:,:,jpcal) * xnegtr(:,:,:) * 1.e3 * rfact2r 
     133      trc3d(:,:,:,jp_pcs0_3d + 10) = tra(:,:,:,jpcal) * xnegtr(:,:,:) * 1.e3 * rfact2r 
    135134# endif 
    136135      ! 

  • trunk/NEMO/TOP_SRC/PISCES/p4zche.F90

    r1073 r1119  
    1818   !!---------------------------------------------------------------------- 
    1919   USE oce_trc         ! 
    20    USE trp_trc         ! 
     20   USE trc         ! 
    2121   USE sms_pisces      !  
    2222 

  • trunk/NEMO/TOP_SRC/PISCES/p4zflx.F90

    r1073 r1119  
    1919   USE trc 
    2020   USE oce_trc         ! 
    21    USE trp_trc 
     21   USE trc 
    2222   USE sms_pisces 
    2323   USE prtctl_trc 
    … …  
    167167# if defined key_trc_diaadd 
    168168            ! Save diagnostics 
    169             trc2d(ji,jj,1) = ( zfld - zflu ) * 1000. 
    170             trc2d(ji,jj,2) = ( zfld16 - zflu16 ) * 1000. 
    171             trc2d(ji,jj,3) = zkgco2(ji,jj) 
    172             trc2d(ji,jj,4) = atcco2 - zh2co3(ji,jj) / ( chemc(ji,jj,1) + rtrn ) 
     169            trc2d(ji,jj,jp_pcs0_2d    ) = ( zfld - zflu ) * 1000. 
     170            trc2d(ji,jj,jp_pcs0_2d + 1) = ( zfld16 - zflu16 ) * 1000. 
     171            trc2d(ji,jj,jp_pcs0_2d + 2) = zkgco2(ji,jj) 
     172            trc2d(ji,jj,jp_pcs0_2d + 3) = atcco2 - zh2co3(ji,jj) / ( chemc(ji,jj,1) + rtrn ) 
    173173# endif 
    174174         END DO 
    … …  
    209209      !! ** Purpose :   Initialization of atmospheric conditions 
    210210      !! 
    211       !! ** Method  :   Read the natext namelist and check the parameters 
     211      !! ** Method  :   Read the nampisext namelist and check the parameters 
    212212      !!      called at the first timestep (nittrc000) 
    213       !! ** input   :   Namelist natext 
     213      !! ** input   :   Namelist nampisext 
    214214      !! 
    215215      !!---------------------------------------------------------------------- 
    216216 
    217       NAMELIST/natext/ atcco2 
     217      NAMELIST/nampisext/ atcco2 
    218218 
    219219      REWIND( numnat )                     ! read numnat 
    220       READ  ( numnat, natext ) 
     220      READ  ( numnat, nampisext ) 
    221221 
    222222      IF(lwp) THEN                         ! control print 
    223223         WRITE(numout,*) ' ' 
    224          WRITE(numout,*) ' Namelist parameters for air-sea exchange, natext' 
     224         WRITE(numout,*) ' Namelist parameters for air-sea exchange, nampisext' 
    225225         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    226226         WRITE(numout,*) '    Atmospheric pCO2      atcco2      =', atcco2 

  • trunk/NEMO/TOP_SRC/PISCES/p4zint.F90

    r1073 r1119  
    1414   !!---------------------------------------------------------------------- 
    1515   USE oce_trc         ! 
    16    USE trp_trc 
     16   USE trc 
    1717   USE sms_pisces 
    1818 

  • trunk/NEMO/TOP_SRC/PISCES/p4zlim.F90

    r1073 r1119  
    1616   USE trc 
    1717   USE oce_trc         ! 
    18    USE trp_trc         !  
     18   USE trc         !  
    1919   USE sms_pisces      !  
    2020 
    … …  
    180180      !! ** Purpose :   Initialization of nutrient limitation parameters 
    181181      !! 
    182       !! ** Method  :   Read the natlim namelist and check the parameters 
     182      !! ** Method  :   Read the nampislim namelist and check the parameters 
    183183      !!      called at the first timestep (nittrc000) 
    184184      !! 
    185       !! ** input   :   Namelist natlim 
     185      !! ** input   :   Namelist nampislim 
    186186      !! 
    187187      !!---------------------------------------------------------------------- 
    188188 
    189       NAMELIST/natlim/ conc0, conc1, conc2, conc2m, conc3, conc3m,   & 
     189      NAMELIST/nampislim/ conc0, conc1, conc2, conc2m, conc3, conc3m,   & 
    190190         &             concnnh4, concdnh4, xksi1, xksi2, xkdoc, caco3r 
    191191 
    192192      REWIND( numnat )                     ! read numnat 
    193       READ  ( numnat, natlim ) 
     193      READ  ( numnat, nampislim ) 
    194194 
    195195      IF(lwp) THEN                         ! control print 
    196196         WRITE(numout,*) ' ' 
    197          WRITE(numout,*) ' Namelist parameters for nutrient limitations, natlim' 
     197         WRITE(numout,*) ' Namelist parameters for nutrient limitations, nampislim' 
    198198         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    199199         WRITE(numout,*) '    mean rainratio                            caco3r    =', caco3r 

  • trunk/NEMO/TOP_SRC/PISCES/p4zlys.F90

    r1073 r1119  
    1919   USE trc 
    2020   USE oce_trc         ! 
    21    USE trp_trc 
     21   USE trc 
    2222   USE sms_pisces 
    2323   USE prtctl_trc 
    … …  
    151151 
    152152# if defined key_trc_dia3d 
    153          trc3d(:,:,:,1) = hi(:,:,:) 
    154          trc3d(:,:,:,2) = zco3(:,:,:) 
    155          trc3d(:,:,:,3) = aksp(:,:,:) / calcon 
     153         trc3d(:,:,:,jp_pcs0_3d    ) = hi(:,:,:) 
     154         trc3d(:,:,:,jp_pcs0_3d + 1) = zco3(:,:,:) 
     155         trc3d(:,:,:,jp_pcs0_3d + 2) = aksp(:,:,:) / calcon 
    156156# endif 
    157157      ! 
    … …  
    171171      !! ** Purpose :   Initialization of CaCO3 dissolution parameters 
    172172      !! 
    173       !! ** Method  :   Read the natcal namelist and check the parameters 
     173      !! ** Method  :   Read the nampiscal namelist and check the parameters 
    174174      !!      called at the first timestep (nittrc000) 
    175175      !! 
    176       !! ** input   :   Namelist natcal 
     176      !! ** input   :   Namelist nampiscal 
    177177      !! 
    178178      !!---------------------------------------------------------------------- 
    179179 
    180       NAMELIST/natcal/ kdca, nca 
     180      NAMELIST/nampiscal/ kdca, nca 
    181181 
    182182      REWIND( numnat )                     ! read numnat 
    183       READ  ( numnat, natcal ) 
     183      READ  ( numnat, nampiscal ) 
    184184 
    185185      IF(lwp) THEN                         ! control print 
    186186         WRITE(numout,*) ' ' 
    187          WRITE(numout,*) ' Namelist parameters for CaCO3 dissolution, natcal' 
     187         WRITE(numout,*) ' Namelist parameters for CaCO3 dissolution, nampiscal' 
    188188         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    189189         WRITE(numout,*) '    diss. rate constant calcite (per month)   kdca      =', kdca 

  • trunk/NEMO/TOP_SRC/PISCES/p4zmeso.F90

    r1073 r1119  
    1616   USE trc 
    1717   USE oce_trc         ! 
    18    USE trp_trc         !  
     18   USE trc         !  
    1919   USE sms_pisces      !  
    2020   USE prtctl_trc 
    … …  
    299299      !! ** Purpose :   Initialization of mesozooplankton parameters 
    300300      !! 
    301       !! ** Method  :   Read the natmes namelist and check the parameters 
     301      !! ** Method  :   Read the nampismes namelist and check the parameters 
    302302      !!      called at the first timestep (nittrc000) 
    303303      !! 
    304       !! ** input   :   Namelist natmes 
     304      !! ** input   :   Namelist nampismes 
    305305      !! 
    306306      !!---------------------------------------------------------------------- 
    307307 
    308       NAMELIST/natmes/ grazrat2,resrat2,mzrat2,xprefc, xprefp, & 
     308      NAMELIST/nampismes/ grazrat2,resrat2,mzrat2,xprefc, xprefp, & 
    309309         &             xprefz, xprefpoc, xkgraz2, epsher2, sigma2, unass2, grazflux 
    310310 
    311311      REWIND( numnat )                     ! read numnat 
    312       READ  ( numnat, natmes ) 
     312      READ  ( numnat, nampismes ) 
    313313 
    314314 
    315315      IF(lwp) THEN                         ! control print 
    316316         WRITE(numout,*) ' '  
    317          WRITE(numout,*) ' Namelist parameters for mesozooplankton, natmes' 
     317         WRITE(numout,*) ' Namelist parameters for mesozooplankton, nampismes' 
    318318         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    319319         WRITE(numout,*) '    zoo preference for phyto                  xprefc    =', xprefc 

  • trunk/NEMO/TOP_SRC/PISCES/p4zmicro.F90

    r1073 r1119  
    1616   USE trc 
    1717   USE oce_trc         ! 
    18    USE trp_trc         !  
     18   USE trc         !  
    1919   USE sms_pisces      !  
    2020   USE prtctl_trc 
    … …  
    229229      !! ** Purpose :   Initialization of microzooplankton parameters 
    230230      !! 
    231       !! ** Method  :   Read the natzoo namelist and check the parameters 
     231      !! ** Method  :   Read the nampiszoo namelist and check the parameters 
    232232      !!      called at the first timestep (nittrc000) 
    233233      !! 
    234       !! ** input   :   Namelist natzoo 
     234      !! ** input   :   Namelist nampiszoo 
    235235      !! 
    236236      !!---------------------------------------------------------------------- 
    237237 
    238       NAMELIST/natzoo/ grazrat,resrat,mzrat,xpref2c, xpref2p, & 
     238      NAMELIST/nampiszoo/ grazrat,resrat,mzrat,xpref2c, xpref2p, & 
    239239         &             xpref2d, xkgraz, epsher, sigma1, unass 
    240240 
    241241      REWIND( numnat )                     ! read numnat 
    242       READ  ( numnat, natzoo ) 
     242      READ  ( numnat, nampiszoo ) 
    243243 
    244244      IF(lwp) THEN                         ! control print 
    245245         WRITE(numout,*) ' ' 
    246          WRITE(numout,*) ' Namelist parameters for microzooplankton, natzoo' 
     246         WRITE(numout,*) ' Namelist parameters for microzooplankton, nampiszoo' 
    247247         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    248248         WRITE(numout,*) '    zoo preference for POC                    xpref2c    =', xpref2c 

  • trunk/NEMO/TOP_SRC/PISCES/p4zmort.F90

    r1073 r1119  
    1616   USE trc 
    1717   USE oce_trc         ! 
    18    USE trp_trc         !  
     18   USE trc         !  
    1919   USE sms_pisces      !  
    2020   USE p4zsink 
    … …  
    241241      !! ** Purpose :   Initialization of phytoplankton parameters 
    242242      !! 
    243       !! ** Method  :   Read the natmort namelist and check the parameters 
     243      !! ** Method  :   Read the nampismort namelist and check the parameters 
    244244      !!      called at the first timestep (nittrc000) 
    245245      !! 
    246       !! ** input   :   Namelist natmort 
     246      !! ** input   :   Namelist nampismort 
    247247      !! 
    248248      !!---------------------------------------------------------------------- 
    249249 
    250       NAMELIST/natmort/ wchl, wchld, mprat, mprat2, mpratm 
     250      NAMELIST/nampismort/ wchl, wchld, mprat, mprat2, mpratm 
    251251 
    252252      REWIND( numnat )                     ! read numnat 
    253       READ  ( numnat, natmort ) 
     253      READ  ( numnat, nampismort ) 
    254254 
    255255      IF(lwp) THEN                         ! control print 
    256256         WRITE(numout,*) ' ' 
    257          WRITE(numout,*) ' Namelist parameters for phytoplankton mortality, natmort' 
     257         WRITE(numout,*) ' Namelist parameters for phytoplankton mortality, nampismort' 
    258258         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    259259         WRITE(numout,*) '    quadratic mortality of phytoplankton      wchl      =', wchl 

  • trunk/NEMO/TOP_SRC/PISCES/p4zopt.F90

    r1073 r1119  
    1515   USE trc 
    1616   USE oce_trc         ! 
    17    USE trp_trc 
     17   USE trc 
    1818   USE sms_pisces 
    1919 
    … …  
    255255 
    256256# if defined key_trc_diaadd 
    257       trc2d(:,:,11) = heup(:,:) 
     257      trc2d(:,:,jp_pcs0_2d + 10) = heup(:,:) 
    258258# endif 
    259259      ! 

  • trunk/NEMO/TOP_SRC/PISCES/p4zprod.F90

    r1073 r1119  
    1515   USE trc 
    1616   USE oce_trc         ! 
    17    USE trp_trc         !  
    1817   USE sms_pisces      !  
    1918   USE prtctl_trc 
    … …  
    376375!   Supplementary diagnostics 
    377376!   ------------------------- 
    378       trc3d(:,:,:,4)  = etot(:,:,:) 
    379       trc3d(:,:,:,5)  = zprorca(:,:,:)  * zrfact2 
    380       trc3d(:,:,:,6)  = zprorcad(:,:,:) * zrfact2 
    381       trc3d(:,:,:,7)  = zpronew(:,:,:)  * zrfact2 
    382       trc3d(:,:,:,8)  = zpronewd(:,:,:) * zrfact2 
    383       trc3d(:,:,:,9)  = zprorcad(:,:,:) * zysopt(:,:,:) * zrfact2 
    384       trc3d(:,:,:,10) = zprofed(:,:,:) * zrfact2 
     377      trc3d(:,:,:,jp_pcs0_3d + 3)  = etot(:,:,:) 
     378      trc3d(:,:,:,jp_pcs0_3d + 4)  = zprorca(:,:,:)  * zrfact2 
     379      trc3d(:,:,:,jp_pcs0_3d + 5)  = zprorcad(:,:,:) * zrfact2 
     380      trc3d(:,:,:,jp_pcs0_3d + 6)  = zpronew(:,:,:)  * zrfact2 
     381      trc3d(:,:,:,jp_pcs0_3d + 7)  = zpronewd(:,:,:) * zrfact2 
     382      trc3d(:,:,:,jp_pcs0_3d + 8)  = zprorcad(:,:,:) * zysopt(:,:,:) * zrfact2 
     383      trc3d(:,:,:,jp_pcs0_3d + 9) = zprofed(:,:,:) * zrfact2 
    385384#if ! defined key_kriest 
    386       trc3d(:,:,:,11) = zprofen(:,:,:) * zrfact2 
     385      trc3d(:,:,:,jp_pcs0_3d + 10) = zprofen(:,:,:) * zrfact2 
    387386#endif 
    388387#endif 
    … …  
    403402      !! ** Purpose :   Initialization of phytoplankton production parameters 
    404403      !! 
    405       !! ** Method  :   Read the natprod namelist and check the parameters 
     404      !! ** Method  :   Read the nampisprod namelist and check the parameters 
    406405      !!      called at the first timestep (nittrc000) 
    407406      !! 
    408       !! ** input   :   Namelist natprod 
     407      !! ** input   :   Namelist nampisprod 
    409408      !! 
    410409      !!---------------------------------------------------------------------- 
    411410 
    412       NAMELIST/natprod/ pislope, pislope2, excret, excret2, chlcnm, chlcdm,   & 
     411      NAMELIST/nampisprod/ pislope, pislope2, excret, excret2, chlcnm, chlcdm,   & 
    413412         &              fecnm, fecdm, grosip 
    414413 
    415414      REWIND( numnat )                     ! read numnat 
    416       READ  ( numnat, natprod ) 
     415      READ  ( numnat, nampisprod ) 
    417416 
    418417      IF(lwp) THEN                         ! control print 
    419418         WRITE(numout,*) ' ' 
    420          WRITE(numout,*) ' Namelist parameters for phytoplankton growth, natprod' 
     419         WRITE(numout,*) ' Namelist parameters for phytoplankton growth, nampisprod' 
    421420         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    422421         WRITE(numout,*) '    mean Si/C ratio                           grosip    =', grosip 

  • trunk/NEMO/TOP_SRC/PISCES/p4zrem.F90

    r1073 r1119  
    1616   USE trc 
    1717   USE oce_trc         ! 
    18    USE trp_trc         !  
    1918   USE sms_pisces      !  
    2019   USE prtctl_trc 
    … …  
    409408      !! ** Purpose :   Initialization of remineralization parameters 
    410409      !! 
    411       !! ** Method  :   Read the natrem namelist and check the parameters 
     410      !! ** Method  :   Read the nampisrem namelist and check the parameters 
    412411      !!      called at the first timestep (nittrc000) 
    413412      !! 
    414       !! ** input   :   Namelist natrem 
     413      !! ** input   :   Namelist nampisrem 
    415414      !! 
    416415      !!---------------------------------------------------------------------- 
    417416 
    418       NAMELIST/natrem/ xremik, xremip, nitrif, xsirem, xlam1, oxymin 
     417      NAMELIST/nampisrem/ xremik, xremip, nitrif, xsirem, xlam1, oxymin 
    419418 
    420419      REWIND( numnat )                     ! read numnat 
    421       READ  ( numnat, natrem ) 
     420      READ  ( numnat, nampisrem ) 
    422421 
    423422      IF(lwp) THEN                         ! control print 
    424423         WRITE(numout,*) ' ' 
    425          WRITE(numout,*) ' Namelist parameters for remineralization, natrem' 
     424         WRITE(numout,*) ' Namelist parameters for remineralization, nampisrem' 
    426425         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    427426         WRITE(numout,*) '    remineralisation rate of POC              xremip    =', xremip 

  • trunk/NEMO/TOP_SRC/PISCES/p4zsed.F90

    r1073 r1119  
    1717   USE trc 
    1818   USE oce_trc         ! 
    19    USE trp_trc 
    2019   USE sms_pisces 
    2120   USE lib_mpp 
    … …  
    316315      DO jj = 1,jpj 
    317316         DO ji = 1,jpi 
    318             trc2d(ji,jj,13) = znitrpot(ji,jj,1) * 1.e-7 * fse3t(ji,jj,1) * 1.e+3 / rfact2 
    319             trc2d(ji,jj,12) = zirondep(ji,jj,1) * 1.e+3 * rfact2r * fse3t(ji,jj,1) 
     317            trc2d(ji,jj,jp_pcs0_2d + 11) = zirondep(ji,jj,1) * 1.e+3 * rfact2r * fse3t(ji,jj,1) 
     318            trc2d(ji,jj,jp_pcs0_2d + 12) = znitrpot(ji,jj,1) * 1.e-7 * fse3t(ji,jj,1) * 1.e+3 / rfact2 
    320319         END DO 
    321320      END DO 
    … …  
    454453      REAL(wp), DIMENSION(jpi,jpj,12)    ::   zdustmo 
    455454 
    456       NAMELIST/natsms/ bdustfer, briver, bndepo, bsedinput, sedfeinput, dustsolub 
     455      NAMELIST/nampissed/ bdustfer, briver, bndepo, bsedinput, sedfeinput, dustsolub 
    457456 
    458457 
    459458      REWIND( numnat )                     ! read numnat 
    460       READ  ( numnat, natsms ) 
     459      READ  ( numnat, nampissed ) 
    461460 
    462461      IF(lwp) THEN 
    463462         WRITE(numout,*) ' ' 
    464          WRITE(numout,*) ' Namelist : natsms ' 
     463         WRITE(numout,*) ' Namelist : nampissed ' 
    465464         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~ ' 
    466465         WRITE(numout,*) '    Dust input from the atmosphere           bdustfer  = ', bdustfer 

  • trunk/NEMO/TOP_SRC/PISCES/p4zsink.F90

    r1073 r1119  
    1212   USE trc 
    1313   USE oce_trc         ! 
    14    USE trp_trc 
    1514   USE sms_pisces 
    1615   USE prtctl_trc 
    … …  
    279278#    if defined key_trc_dia3d 
    280279      zrfact2 = 1.e3 * rfact2r 
    281       trc2d(:,:, 5)   = sinking (:,:,iksed+1) * zrfact2 
    282       trc2d(:,:, 6)   = sinking2(:,:,iksed+1) * zrfact2 
    283       trc2d(:,:, 7)   = sinkfer (:,:,iksed+1) * zrfact2 
    284       trc2d(:,:, 9)   = sinksil (:,:,iksed+1) * zrfact2 
    285       trc2d(:,:,10)   = sinkcal (:,:,iksed+1) * zrfact2 
    286       trc3d(:,:,:,12) = sinking (:,:,:)       * zrfact2 
    287       trc3d(:,:,:,13) = sinking2(:,:,:)       * zrfact2 
    288       trc3d(:,:,:,14) = sinksil (:,:,:)       * zrfact2 
    289       trc3d(:,:,:,15) = sinkcal (:,:,:)       * zrfact2 
    290       trc3d(:,:,:,16) = znum3d  (:,:,:) 
    291       trc3d(:,:,:,17) = wsbio3  (:,:,:) 
    292       trc3d(:,:,:,18) = wsbio4  (:,:,:) 
     280      trc2d(:,:,  jp_pcs0_2d + 4)   = sinking (:,:,iksed+1) * zrfact2 
     281      trc2d(:,:,  jp_pcs0_2d + 5)   = sinking2(:,:,iksed+1) * zrfact2 
     282      trc2d(:,:,  jp_pcs0_2d + 6)   = sinkfer (:,:,iksed+1) * zrfact2 
     283      trc2d(:,:,  jp_pcs0_2d + 8)   = sinksil (:,:,iksed+1) * zrfact2 
     284      trc2d(:,:,  jp_pcs0_2d + 9)   = sinkcal (:,:,iksed+1) * zrfact2 
     285      trc3d(:,:,:,jp_pcs0_3d + 11) = sinking (:,:,:)       * zrfact2 
     286      trc3d(:,:,:,jp_pcs0_3d + 12) = sinking2(:,:,:)       * zrfact2 
     287      trc3d(:,:,:,jp_pcs0_3d + 13) = sinksil (:,:,:)       * zrfact2 
     288      trc3d(:,:,:,jp_pcs0_3d + 14) = sinkcal (:,:,:)       * zrfact2 
     289      trc3d(:,:,:,jp_pcs0_3d + 15) = znum3d  (:,:,:) 
     290      trc3d(:,:,:,jp_pcs0_3d + 16) = wsbio3  (:,:,:) 
     291      trc3d(:,:,:,jp_pcs0_3d + 17) = wsbio4  (:,:,:) 
    293292#    endif 
    294293      ! 
    … …  
    308307      !!                Kriest parameterization only 
    309308      !! 
    310       !! ** Method  :   Read the natkriest namelist and check the parameters 
     309      !! ** Method  :   Read the nampiskrs namelist and check the parameters 
    311310      !!      called at the first timestep (nittrc000) 
    312311      !! 
    313       !! ** input   :   Namelist natkriest 
     312      !! ** input   :   Namelist nampiskrs 
    314313      !! 
    315314      !!---------------------------------------------------------------------- 
    … …  
    319318      REAL(wp) ::   zmin, zmax, zl, zr, xacc 
    320319 
    321       NAMELIST/natkrsize/ xkr_sfact, xkr_stick ,  & 
     320      NAMELIST/nampiskrs/ xkr_sfact, xkr_stick ,  & 
    322321         &                xkr_nnano, xkr_ndiat, xkr_nmeso, xkr_naggr 
    323322 
    324323      !!---------------------------------------------------------------------- 
    325       !                               ! natkriest : kriest parameters 
    326       !                               ! ----------------------------- 
    327       REWIND( numnat )                     ! read natkriest 
    328       READ  ( numnat, natkrsize ) 
     324      REWIND( numnat )                     ! read nampiskrs 
     325      READ  ( numnat, nampiskrs ) 
    329326 
    330327      IF(lwp) THEN 
    331328         WRITE(numout,*) 
    332          WRITE(numout,*) ' Namelist : natkrsize' 
     329         WRITE(numout,*) ' Namelist : nampiskrs' 
    333330         WRITE(numout,*) '    Sinking factor                           xkr_sfact    = ', xkr_sfact 
    334331         WRITE(numout,*) '    Stickiness                               xkr_stick    = ', xkr_stick 
    … …  
    594591      END DO 
    595592 
    596 # if defined key_trc_dia3d 
     593# if defined key_trc_diaadd 
    597594      zrfact2 = 1.e3 * rfact2r 
    598       trc2d(:,:, 5) = sinking (:,:,iksed+1) * zrfact2 
    599       trc2d(:,:, 6) = sinking2(:,:,iksed+1) * zrfact2 
    600       trc2d(:,:, 7) = sinkfer (:,:,iksed+1) * zrfact2 
    601       trc2d(:,:, 8) = sinkfer2(:,:,iksed+1) * zrfact2 
    602       trc2d(:,:, 9) = sinksil (:,:,iksed+1) * zrfact2 
    603       trc2d(:,:,10) = sinkcal (:,:,iksed+1) * zrfact2 
     595      trc2d(:,:, jp_pcs0_2d + 4) = sinking (:,:,iksed+1) * zrfact2 
     596      trc2d(:,:, jp_pcs0_2d + 5) = sinking2(:,:,iksed+1) * zrfact2 
     597      trc2d(:,:, jp_pcs0_2d + 6) = sinkfer (:,:,iksed+1) * zrfact2 
     598      trc2d(:,:, jp_pcs0_2d + 7) = sinkfer2(:,:,iksed+1) * zrfact2 
     599      trc2d(:,:, jp_pcs0_2d + 8) = sinksil (:,:,iksed+1) * zrfact2 
     600      trc2d(:,:, jp_pcs0_2d + 9) = sinkcal (:,:,iksed+1) * zrfact2 
    604601# endif 
    605602      ! 

  • trunk/NEMO/TOP_SRC/PISCES/par_pisces.F90

    r935 r1119  
    1010   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
    1111   !!---------------------------------------------------------------------- 
    12    USE par_lobster, ONLY : jp_lobster   !: number of tracers in LOBSTER 
     12   USE par_lobster, ONLY : jp_lobster      !: number of tracers in LOBSTER 
     13   USE par_lobster, ONLY : jp_lobster_2d   !: number of 2D diag in LOBSTER 
     14   USE par_lobster, ONLY : jp_lobster_3d   !: number of 3D diag in LOBSTER 
     15   USE par_lobster, ONLY : jp_lobster_trd  !: number of biological diag in LOBSTER 
    1316 
    1417   IMPLICIT NONE 
    1518   PUBLIC 
    1619 
    17    INTEGER, PUBLIC, PARAMETER ::   jp_l  = jp_lobster      !: cumulative number of already defined TRC 
     20   INTEGER, PUBLIC, PARAMETER ::   jp_l      = jp_lobster      !: cumulative number of already defined TRC 
     21   INTEGER, PUBLIC, PARAMETER ::   jp_l_2d   = jp_lobster_2d   !: 
     22   INTEGER, PUBLIC, PARAMETER ::   jp_l_3d   = jp_lobster_3d   !: 
     23   INTEGER, PUBLIC, PARAMETER ::   jp_l_trd  = jp_lobster_trd  !: 
    1824 
    1925#if defined key_pisces  &&  defined key_kriest 
    … …  
    107113 
    108114   ! Starting/ending PISCES do-loop indices (N.B. no PISCES : jpl_pcs < jpf_pcs the do-loop are never done) 
    109    INTEGER, PUBLIC, PARAMETER ::   jp_pcs0 = jp_l + 1           !: First index of CFC passive tracers 
    110    INTEGER, PUBLIC, PARAMETER ::   jp_pcs1 = jp_l + jp_pisces   !: Last  index of CFC passive tracers 
     115   INTEGER, PUBLIC, PARAMETER ::   jp_pcs0     = jp_l + 1                  !: First index of PISCES tracers 
     116   INTEGER, PUBLIC, PARAMETER ::   jp_pcs1     = jp_l + jp_pisces          !: Last  index of PISCES tracers 
     117   INTEGER, PUBLIC, PARAMETER ::   jp_pcs0_2d  = jp_l_2d + 1               !: First index of 2D diag 
     118   INTEGER, PUBLIC, PARAMETER ::   jp_pcs1_2d  = jp_l_2d + jp_pisces_2d    !: Last  index of 2D diag 
     119   INTEGER, PUBLIC, PARAMETER ::   jp_pcs0_3d  = jp_l_3d + 1               !: First index of 3D diag 
     120   INTEGER, PUBLIC, PARAMETER ::   jp_pcs1_3d  = jp_l_3d + jp_pisces_3d    !: Last  index of 3d diag 
     121   INTEGER, PUBLIC, PARAMETER ::   jp_pcs0_trd = jp_l_trd + 1              !: First index of bio diag 
     122   INTEGER, PUBLIC, PARAMETER ::   jp_pcs1_trd = jp_l_trd + jp_pisces_trd  !: Last  index of bio diag 
     123 
    111124 
    112125   !!====================================================================== 

  • trunk/NEMO/TOP_SRC/PISCES/trcini_pisces.F90

    r1073 r1119  
    2020   USE trc 
    2121   USE oce_trc         ! ocean variables 
    22    USE trp_trc         ! 
    2322   USE p4zche  
    2423   USE lib_mpp 

  • trunk/NEMO/TOP_SRC/PISCES/trclsm_pisces.F90

    r1089 r1119  
    4545      !!                       natkriest ("key_kriest") 
    4646      !!---------------------------------------------------------------------- 
    47       CHARACTER (len=32) ::   clname 
     47      CHARACTER (len=32) ::  clname 
     48      INTEGER            ::  jn 
    4849      !! 
    49       NAMELIST/natbio/ part, nrdttrc, wsbio, xkmort, ferat3, wsbio2 
     50#if defined key_trc_diaadd 
     51      ! definition of additional diagnostic as a structure 
     52      TYPE DIAG 
     53         CHARACTER(len = 20)  :: snamedia   !: short name 
     54         CHARACTER(len = 80 ) :: lnamedia   !: long name 
     55         CHARACTER(len = 20 ) :: unitdia    !: unit 
     56      END TYPE DIAG 
    5057 
     58      TYPE(DIAG) , DIMENSION(jp_pisces_2d) :: pisdia2d 
     59      TYPE(DIAG) , DIMENSION(jp_pisces_3d) :: pisdia3d 
     60#endif 
     61 
     62      NAMELIST/nampisbio/ part, nrdttrc, wsbio, xkmort, ferat3, wsbio2 
    5163#if defined key_kriest 
    52       NAMELIST/natkrpar/xkr_eta, xkr_zeta, xkr_mass_min, xkr_mass_max 
     64      NAMELIST/nampiskrp/xkr_eta, xkr_zeta, xkr_mass_min, xkr_mass_max 
    5365#endif 
     66#if defined key_trc_diaadd 
     67      NAMELIST/nampisdia/nwritedia, pisdia3d, pisdia2d     ! additional diagnostics 
     68#endif 
     69 
    5470      !!---------------------------------------------------------------------- 
    5571 
    … …  
    6278      !                               ! ---------------------- 
    6379      clname ='namelist_pisces' 
    64       CALL ctlopn( numnat, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL',   & 
    65          &           1, numout, .FALSE., 1 ) 
     80      CALL ctlopn( numnat, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', 1, numout, .FALSE., 1 ) 
    6681 
    67       !                               ! natbio : biological parameters 
    68       !                               ! ------------------------------ 
    69       REWIND( numnat )                     ! read natbio 
    70       READ  ( numnat, natbio ) 
     82      REWIND( numnat )                     
     83      READ  ( numnat, nampisbio ) 
    7184 
    7285      IF(lwp) THEN                         ! control print 
    73          WRITE(numout,*) ' Namelist : natbio' 
     86         WRITE(numout,*) ' Namelist : nampisbio' 
    7487         WRITE(numout,*) '    part of calcite not dissolved in guts     part      =', part 
    7588         WRITE(numout,*) '    frequence pour la biologie                nrdttrc   =', nrdttrc 
    … …  
    8295#if defined key_kriest 
    8396 
    84       !                               ! natkrpar : kriest parameters 
     97      !                               ! nampiskrp : kriest parameters 
    8598      !                               ! ----------------------------- 
    8699      xkr_eta      = 0.62         
    … …  
    90103 
    91104      REWIND( numnat )                     ! read natkriest 
    92       READ  ( numnat, natkrpar ) 
     105      READ  ( numnat, nampiskrp ) 
    93106 
    94107      IF(lwp) THEN 
    95108         WRITE(numout,*) 
    96          WRITE(numout,*) ' Namelist : natkrpar' 
     109         WRITE(numout,*) ' Namelist : nampiskrp' 
    97110         WRITE(numout,*) '    Sinking  exponent                        xkr_eta      = ', xkr_eta 
    98111         WRITE(numout,*) '    N content exponent                       xkr_zeta     = ', xkr_zeta 
    … …  
    108121#endif 
    109122      ! 
     123#if defined key_trc_diaadd 
     124 
     125      ! Namelist namlobdia 
     126      ! ------------------- 
     127      nwritedia = 10                   ! default values 
     128 
     129      DO jn = jp_pcs0_2d, jp_pcs1_2d 
     130         WRITE(ctrc2d(jn),'("2D_",I1)') jn                      ! short name 
     131         WRITE(ctrc2l(jn),'("2D DIAGNOSTIC NUMBER ",I2)') jn    ! long name 
     132         ctrc2u(jn) = ' '                                       ! units 
     133      END DO 
     134      !                                 ! 3D output arrays 
     135      DO jn = jp_pcs0_3d, jp_pcs1_3d 
     136         WRITE(ctrc3d(jn),'("3D_",I1)') jn                      ! short name 
     137         WRITE(ctrc3l(jn),'("3D DIAGNOSTIC NUMBER ",I2)') jn    ! long name 
     138         ctrc3u(jn) = ' '                                       ! units 
     139      END DO 
     140 
     141      REWIND( numnat )               ! read natrtd 
     142      READ  ( numnat, nampisdia ) 
     143 
     144      DO jn = jp_pcs0_2d, jp_pcs1_2d 
     145         ctrc2d(jn) = pisdia2d(jn)%snamedia 
     146         ctrc2l(jn) = pisdia2d(jn)%lnamedia 
     147         ctrc2u(jn) = pisdia2d(jn)%unitdia 
     148      END DO 
     149 
     150      DO jn = jp_pcs0_3d, jp_pcs1_3d 
     151         ctrc3d(jn) = pisdia3d(jn)%snamedia 
     152         ctrc3l(jn) = pisdia3d(jn)%lnamedia 
     153         ctrc3u(jn) = pisdia3d(jn)%unitdia 
     154      END DO 
     155 
     156      IF(lwp) THEN                   ! control print 
     157         WRITE(numout,*) 
     158         WRITE(numout,*) ' Namelist : natadd' 
     159         WRITE(numout,*) '    frequency of outputs for additional arrays nwritedia = ', nwritedia 
     160         DO jn = jp_pcs0_3d, jp_pcs1_3d 
     161            WRITE(numout,*) '   3d output field No : ',jn 
     162            WRITE(numout,*) '   short name         : ', TRIM(ctrc3d(jn)) 
     163            WRITE(numout,*) '   long name          : ', TRIM(ctrc3l(jn)) 
     164            WRITE(numout,*) '   unit               : ', TRIM(ctrc3u(jn)) 
     165            WRITE(numout,*) ' ' 
     166         END DO 
     167 
     168         DO jn = jp_pcs0_2d, jp_pcs1_2d 
     169            WRITE(numout,*) '   2d output field No : ',jn 
     170            WRITE(numout,*) '   short name         : ', TRIM(ctrc2d(jn)) 
     171            WRITE(numout,*) '   long name          : ', TRIM(ctrc2l(jn)) 
     172            WRITE(numout,*) '   unit               : ', TRIM(ctrc2u(jn)) 
     173            WRITE(numout,*) ' ' 
     174         END DO 
     175      ENDIF 
     176#endif 
     177 
    110178   END SUBROUTINE trc_lsm_pisces 
    111179 

  • trunk/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90

    r1073 r1119  
    1414   !!---------------------------------------------------------------------- 
    1515   USE oce_trc         ! 
    16    USE trp_trc 
     16   USE trc 
    1717   USE sms_pisces 
    1818   USE lbclnk 
Note: See TracChangeset for help on using the changeset viewer.