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 336 for trunk/NEMO/TOP_SRC/SMS/trclsm.pisces.h90 – NEMO

Ignore:
Timestamp:
2005-11-14T13:13:26+01:00 (18 years ago)
Author:
opalod
Message:

nemo_v1_update_024 : CE + RB + CT : new evolution of modules

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/SMS/trclsm.pisces.h90

    r274 r336  
    1 !! 
    2 !!                       trclsm.pisces.h 
    3 !!                       **************** 
    4 !! 
    5 !!  PURPOSE : 
    6 !!  --------- 
    7 !!     READs and PRINT options for PISCES namelist 
    8 !! 
    9 !!   MODIFICATIONS: 
    10 !!   -------------- 
    11 !!      original  : 99-10 (M.A. Foujols, M. Levy) passive tracer 
    12 !!      addition  : 00-01 (L. Bopp) hamocc3,p3zd 
    13 !!      
    141!!---------------------------------------------------------------------- 
     2!!                    ***  trclsm.pisces.h90 *** 
    153!!---------------------------------------------------------------------- 
    16 !! local declarations 
    17 !! ================== 
    18  
    19 #if defined key_passivetrc && defined key_trc_pisces 
     4CONTAINS 
     5 
     6   SUBROUTINE trc_lsm 
     7      !!---------------------------------------------------------------------- 
     8      !! 
     9      !!                       trclsm.pisces.h 
     10      !!                       **************** 
     11      !! 
     12      !!  PURPOSE : 
     13      !!  --------- 
     14      !!     READs and PRINT options for PISCES namelist 
     15      !! 
     16      !!   MODIFICATIONS: 
     17      !!   -------------- 
     18      !!      original  : 99-10 (M.A. Foujols, M. Levy) passive tracer 
     19      !!      addition  : 00-01 (L. Bopp) hamocc3,p3zd 
     20      !!      
     21      !!---------------------------------------------------------------------- 
     22      !!---------------------------------------------------------------------- 
     23      !! local declarations 
     24      !! ================== 
    2025      CHARACTER (len=32) clname 
    2126 
    22 !!--------------------------------------------------------------------- 
    23 !!  TOP 1.0,  LOCEAN-IPSL (2005) 
    24 !! $Header$ 
    25 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
    26 !!--------------------------------------------------------------------- 
    27  
    28 ! 0. initializations 
    29 ! ------------------ 
    30 ! 
    31        NAMELIST/natgas/ gasfac, igaswind, icice 
    32        NAMELIST/natext/ atcco2 
    33        NAMELIST/natbio/caco3r,                                          & 
    34        &               dispo0,conc0,oxymin,grosip, sedlam,              & 
    35        &               sedlostpoc,sedlostcal,sedlostsil,nrdttrc,        & 
    36        &               pislope, excret,wsbio,wchl,resrat,mprat,mzrat,   & 
    37        &               grazrat,xprefc,xprefp,unass,xkgraz,xkmort,xksi1, & 
    38        &               xksi2,sicmax,xremip,xremik,xsirem,xkdoc1,xkdoc2, & 
    39        &               excret2,resrat2,mprat2,mpratm,mzrat2,grazrat2,   & 
    40        &               xprefz,xprefpoc,unass2,wchl2,xkgraz2,xlam1,      & 
    41        &               ferat3,conc1,conc2,conc3,concnnh4,concdnh4,      & 
    42        &               nitrif,epsher,epsher2,pislope2,wsbio2,sigma1,    & 
    43        &               sigma2, zprefc, zprefp, zprefd,sedfeinput 
    44        NAMELIST/natsms/bdustfer, briver, bndepo, bsedinput 
    45  
    46 ! initialize the number of LOGICAL UNIT used 
    47 ! ------------------------------------------ 
    48      
    49       IF(lwp) THEN 
    50           WRITE(numout,*) ' ' 
    51           WRITE(numout,*) ' ROUTINE trclec' 
    52           WRITE(numout,*) ' **************' 
    53           WRITE(numout,*) ' ' 
    54           WRITE(numout,*) ' namelist for PISCES model' 
    55           WRITE(numout,*) ' ***********************' 
    56           WRITE(numout,*) ' ' 
     27      !!--------------------------------------------------------------------- 
     28      !!  OPA8, LODYC (15/11/96) 
     29      !!--------------------------------------------------------------------- 
     30 
     31      ! 0. initializations 
     32      ! ------------------ 
     33      ! 
     34      NAMELIST/natgas/ gasfac, igaswind, icice 
     35      NAMELIST/natext/ atcco2 
     36      NAMELIST/natbio/caco3r,                                          & 
     37         &          dispo0,conc0,oxymin,grosip, sedlam,              & 
     38         &          sedlostpoc,sedlostcal,sedlostsil,nrdttrc,        & 
     39         &          pislope, excret,wsbio,wchl,wchld,resrat,mprat,mzrat,   & 
     40         &          grazrat,xprefc,xprefp,unass,xkgraz,xkmort,xksi1, & 
     41         &          xksi2,xremip,xremik,xsirem,xkdoc1,xkdoc2, & 
     42         &          excret2,resrat2,mprat2,mpratm,mzrat2,grazrat2,   & 
     43         &          xprefz,xprefpoc,unass2,xkgraz2,xlam1,      & 
     44         &          ferat3,conc1,conc2,conc3,concnnh4,concdnh4,      & 
     45         &          nitrif,epsher,epsher2,pislope2,wsbio2,sigma1,    & 
     46         &          sigma2, zprefc, zprefp, zprefd,fecnm,fecdm,      & 
     47         &          chlcnm,chlcdm, sedfeinput 
     48      NAMELIST/natsms/bdustfer, briver, bndepo, bsedinput 
     49 
     50      ! initialize the number of LOGICAL UNIT used 
     51      ! ------------------------------------------ 
     52 
     53      IF(lwp) THEN 
     54         WRITE(numout,*) ' ' 
     55         WRITE(numout,*) ' ROUTINE trclec' 
     56         WRITE(numout,*) ' **************' 
     57         WRITE(numout,*) ' ' 
     58         WRITE(numout,*) ' namelist for PISCES model' 
     59         WRITE(numout,*) ' ***********************' 
     60         WRITE(numout,*) ' ' 
    5761      ENDIF 
    5862 
     
    6266 
    6367 
    64 ! 1 Namelist natgas : 
    65 ! ------------------- 
     68      ! 1 Namelist natgas : 
     69      ! ------------------- 
    6670 
    6771      READ(numnat,natgas) 
    6872 
    6973      IF(lwp) THEN 
    70           WRITE(numout,*) ' ' 
    71           WRITE(numout,*) 'natgas' 
    72           write(numout,*) 'gasfac = ',gasfac 
    73           WRITE(numout,*) ' ' 
    74           write(numout,*) 'igaswind = ',igaswind 
    75           WRITE(numout,*) ' ' 
    76           write(numout,*) 'icice = ',icice 
    77           WRITE(numout,*) ' ' 
    78       ENDIF 
    79  
    80 ! 2 Namelist natext : 
    81 ! ------------------- 
     74         WRITE(numout,*) ' ' 
     75         WRITE(numout,*) 'natgas' 
     76         write(numout,*) 'gasfac = ',gasfac 
     77         WRITE(numout,*) ' ' 
     78         write(numout,*) 'igaswind = ',igaswind 
     79         WRITE(numout,*) ' ' 
     80         write(numout,*) 'icice = ',icice 
     81         WRITE(numout,*) ' ' 
     82      ENDIF 
     83 
     84      ! 2 Namelist natext : 
     85      ! ------------------- 
    8286      READ(numnat,natext) 
    8387 
    8488      IF(lwp) THEN 
    85           WRITE(numout,*) ' ' 
    86           WRITE(numout,*) 'natext' 
    87           WRITE(numout,*) ' ' 
    88           WRITE(numout,*) 'atmospheric pCO2= ',atcco2 
    89           WRITE(numout,*) ' ' 
     89         WRITE(numout,*) ' ' 
     90         WRITE(numout,*) 'natext' 
     91         WRITE(numout,*) ' ' 
     92         WRITE(numout,*) 'atmospheric pCO2= ',atcco2 
     93         WRITE(numout,*) ' ' 
    9094      ENDIF 
    9195 
     
    9397      READ(numnat,natbio) 
    9498      IF(lwp) THEN 
    95           WRITE(numout,*) 'natbio' 
    96           WRITE(numout,*) ' ' 
    97           WRITE(numout,*)      & 
    98           &   ' mean rainratio                             =', caco3r 
    99           WRITE(numout,*)      & 
    100           &   ' mean Si/C ratio                            =', grosip 
    101           WRITE(numout,*)      & 
    102           &   ' Calcite dissolution half saturation        =', dispo0 
    103           WRITE(numout,*)      & 
    104           &   ' Phosphate half saturation                  =', conc0 
    105           WRITE(numout,*)      & 
    106           &   ' Sediment bioturbation factor               =', sedlam 
    107           WRITE(numout,*)      & 
    108           &   ' Sediment burying ratio for POC         =', sedlostpoc 
    109           WRITE(numout,*)      & 
    110           &   ' Sediment burying ratio for CACO3       =', sedlostcal 
    111           WRITE(numout,*)      & 
    112           &   ' Sediment burying ratio for SI          =', sedlostsil 
    113           WRITE(numout,*)      & 
    114           &   ' frequence pour la biologie                 =', nrdttrc 
    115           WRITE(numout,*)      & 
    116           &   ' P-I slope                                  =', pislope 
    117           WRITE(numout,*)      & 
    118           &   ' excretion ratio of phytoplankton           =', excret 
    119           WRITE(numout,*)      & 
    120           &   ' POC sinking speed                          =', wsbio 
    121           WRITE(numout,*)      & 
    122           &   ' exsudation rate of zooplankton             =', resrat 
    123           WRITE(numout,*)      & 
    124           &   ' phytoplankton mortality rate               =', mprat 
    125           WRITE(numout,*)      & 
    126           &   ' zooplankton mortality rate                 =', mzrat 
    127           WRITE(numout,*)      & 
    128           &   ' zoo preference for phyto                   =', xprefc 
    129           WRITE(numout,*)      & 
    130           &   ' zoo preference for POC                     =', xprefp 
    131           WRITE(numout,*)      & 
    132           &   ' maximal zoo grazing rate                   =', grazrat 
    133           WRITE(numout,*)      & 
    134           &   ' non assimilated fraction of phyto by zoo   =', unass 
    135           WRITE(numout,*)      & 
    136           &   ' half sturation constant for grazing        =', xkgraz 
    137           WRITE(numout,*)      & 
    138           &   ' half saturation constant for mortality     =', xkmort 
    139           WRITE(numout,*)      & 
    140           &   ' half saturation constant for Si uptake     =', xksi1 
    141           WRITE(numout,*)      & 
    142           &   ' half saturation constant for Si/C          =', xksi2 
    143           WRITE(numout,*)      & 
    144           &   ' maximum Si/C                               =', sicmax 
    145           WRITE(numout,*)      & 
    146           &   ' remineralisation rate of POC               =', xremip 
    147           WRITE(numout,*)      & 
    148           &   ' remineralization rate of DOC               =', xremik 
    149           WRITE(numout,*)      &  
    150           &   ' remineralization rate of Si                =', xsirem  
    151           WRITE(numout,*)      & 
    152           &   ' 1st half-sat. of DOC remineralization      =', xkdoc1 
    153           WRITE(numout,*)      & 
    154           &   ' 2nd half-sat. of DOC remineralization      =', xkdoc2 
    155           WRITE(numout,*)      & 
    156           &   ' excretion ratio of diatoms                 =', excret2 
    157           WRITE(numout,*)      & 
    158           &   ' exsudation rate of mesozooplankton         =', resrat2 
    159           WRITE(numout,*)      & 
    160           &   ' Diatoms mortality rate                     =', mprat2 
    161           WRITE(numout,*)      & 
    162           &   ' Phytoplankton minimum mortality rate       =', mpratm 
    163           WRITE(numout,*)      & 
    164           &   ' mesozooplankton mortality rate             =', mzrat2 
    165           WRITE(numout,*)      & 
    166           &   ' zoo preference for zoo                     =', xprefz 
    167           WRITE(numout,*)      & 
    168           &   ' zoo preference for poc                   =', xprefpoc 
    169           WRITE(numout,*)      & 
    170           &   ' maximal mesozoo grazing rate               =', grazrat2 
    171           WRITE(numout,*)      & 
    172           &   ' non assimilated fraction of P by mesozoo   =', unass2 
    173           WRITE(numout,*)      & 
    174           &   ' Efficicency of Mesozoo growth              =', epsher2  
    175           WRITE(numout,*)      & 
    176           &   ' Efficiency of microzoo growth              =', epsher 
    177           WRITE(numout,*)      & 
    178           &   ' half sturation constant for grazing 2      =', xkgraz2 
    179           WRITE(numout,*)      & 
    180           &   ' Aggregation rate for diatoms               =', wchl2 
    181           WRITE(numout,*)      & 
    182           &   ' scavenging rate of Iron                    =', xlam1 
    183           WRITE(numout,*)      & 
    184           &   ' Fe/C in zooplankton                        =', ferat3 
    185           WRITE(numout,*)      & 
    186           &   ' Phosphate half saturation for diatoms      =', conc1 
    187           WRITE(numout,*)      & 
    188           &   ' Iron half saturation for phyto             =', conc2 
    189           WRITE(numout,*)      & 
    190           &   ' Iron half saturation for diatoms           =', conc3 
    191           WRITE(numout,*)      & 
    192           &   ' NH4 half saturation for phyto              =', concnnh4 
    193           WRITE(numout,*)      & 
    194           &   ' NH4 half saturation for diatoms            =', concdnh4 
    195           WRITE(numout,*)      & 
    196           &   ' NH4 nitrification rate                     =', nitrif 
    197           WRITE(numout,*)      & 
    198           &   ' P-I slope  for diatoms                     =', pislope2 
    199           WRITE(numout,*)      & 
    200           &   ' Big particles sinking speed                =', wsbio2 
    201           WRITE(numout,*)      & 
    202           &   ' Fraction of microzoo excretion as DOM      =', sigma1 
    203           WRITE(numout,*)      & 
    204           &   ' Fraction of mesozoo excretion as DOM       =', sigma2 
    205           WRITE(numout,*)      & 
    206           &   ' Microzoo preference for POM                =', zprefc 
    207           WRITE(numout,*)      & 
    208           &   ' Microzoo preference for Nanophyto          =', zprefp 
    209           WRITE(numout,*)      & 
    210           &   ' Microzoo preference for Diatoms          =', zprefd 
    211           WRITE(numout,*)      & 
    212           &   ' Coastal release of Iron                 =', sedfeinput 
     99         WRITE(numout,*) 'natbio' 
     100         WRITE(numout,*) ' ' 
     101         WRITE(numout,*)      & 
     102            &   ' mean rainratio                             =', caco3r 
     103         WRITE(numout,*)      & 
     104            &   ' mean Si/C ratio                            =', grosip 
     105         WRITE(numout,*)      & 
     106            &   ' Calcite dissolution half saturation        =', dispo0 
     107         WRITE(numout,*)      & 
     108            &   ' Phosphate half saturation                  =', conc0 
     109         WRITE(numout,*)      & 
     110            &   ' Sediment bioturbation factor               =', sedlam 
     111         WRITE(numout,*)      & 
     112            &   ' Sediment burying ratio for POC         =', sedlostpoc 
     113         WRITE(numout,*)      & 
     114            &   ' Sediment burying ratio for CACO3       =', sedlostcal 
     115         WRITE(numout,*)      & 
     116            &   ' Sediment burying ratio for SI          =', sedlostsil 
     117         WRITE(numout,*)      & 
     118            &   ' frequence pour la biologie                 =', nrdttrc 
     119         WRITE(numout,*)      & 
     120            &   ' P-I slope                                  =', pislope 
     121         WRITE(numout,*)      & 
     122            &   ' excretion ratio of phytoplankton           =', excret 
     123         WRITE(numout,*)      & 
     124            &   ' POC sinking speed                          =', wsbio 
     125         WRITE(numout,*)      & 
     126            &   ' exsudation rate of zooplankton             =', resrat 
     127         WRITE(numout,*)      & 
     128            &   ' phytoplankton mortality rate               =', mprat 
     129         WRITE(numout,*)      & 
     130            &   ' zooplankton mortality rate                 =', mzrat 
     131         WRITE(numout,*)      & 
     132            &   ' zoo preference for phyto                   =', xprefc 
     133         WRITE(numout,*)      & 
     134            &   ' zoo preference for POC                     =', xprefp 
     135         WRITE(numout,*)      & 
     136            &   ' maximal zoo grazing rate                   =', grazrat 
     137         WRITE(numout,*)      & 
     138            &   ' non assimilated fraction of phyto by zoo   =', unass 
     139         WRITE(numout,*)      & 
     140            &   ' half sturation constant for grazing        =', xkgraz 
     141         WRITE(numout,*)      & 
     142            &   ' half saturation constant for mortality     =', xkmort 
     143         WRITE(numout,*)      & 
     144            &   ' half saturation constant for Si uptake     =', xksi1 
     145         WRITE(numout,*)      & 
     146            &   ' half saturation constant for Si/C          =', xksi2 
     147         WRITE(numout,*)      & 
     148            &   ' remineralisation rate of POC               =', xremip 
     149         WRITE(numout,*)      & 
     150            &   ' remineralization rate of DOC               =', xremik 
     151         WRITE(numout,*)      &  
     152            &   ' remineralization rate of Si                =', xsirem  
     153         WRITE(numout,*)      & 
     154            &   ' 1st half-sat. of DOC remineralization      =', xkdoc1 
     155         WRITE(numout,*)      & 
     156            &   ' 2nd half-sat. of DOC remineralization      =', xkdoc2 
     157         WRITE(numout,*)      & 
     158            &   ' excretion ratio of diatoms                 =', excret2 
     159         WRITE(numout,*)      & 
     160            &   ' exsudation rate of mesozooplankton         =', resrat2 
     161         WRITE(numout,*)      & 
     162            &   ' Diatoms mortality rate                     =', mprat2 
     163         WRITE(numout,*)      & 
     164            &   ' Phytoplankton minimum mortality rate       =', mpratm 
     165         WRITE(numout,*)      & 
     166            &   ' mesozooplankton mortality rate             =', mzrat2 
     167         WRITE(numout,*)      & 
     168            &   ' zoo preference for zoo                     =', xprefz 
     169         WRITE(numout,*)      & 
     170            &   ' zoo preference for poc                   =', xprefpoc 
     171         WRITE(numout,*)      & 
     172            &   ' maximal mesozoo grazing rate               =', grazrat2 
     173         WRITE(numout,*)      & 
     174            &   ' non assimilated fraction of P by mesozoo   =', unass2 
     175         WRITE(numout,*)      & 
     176            &   ' Efficicency of Mesozoo growth              =', epsher2  
     177         WRITE(numout,*)      & 
     178            &   ' Efficiency of microzoo growth              =', epsher 
     179         WRITE(numout,*)      & 
     180            &   ' half sturation constant for grazing 2      =', xkgraz2 
     181         WRITE(numout,*)      & 
     182            &   ' Maximum aggregation rate for diatoms       =', wchld 
     183         WRITE(numout,*)      & 
     184            &   ' scavenging rate of Iron                    =', xlam1 
     185         WRITE(numout,*)      & 
     186            &   ' Fe/C in zooplankton                        =', ferat3 
     187         WRITE(numout,*)      & 
     188            &   ' Phosphate half saturation for diatoms      =', conc1 
     189         WRITE(numout,*)      & 
     190            &   ' Iron half saturation for phyto             =', conc2 
     191         WRITE(numout,*)      & 
     192            &   ' Iron half saturation for diatoms           =', conc3 
     193         WRITE(numout,*)      & 
     194            &   ' NH4 half saturation for phyto              =', concnnh4 
     195         WRITE(numout,*)      & 
     196            &   ' NH4 half saturation for diatoms            =', concdnh4 
     197         WRITE(numout,*)      & 
     198            &   ' NH4 nitrification rate                     =', nitrif 
     199         WRITE(numout,*)      & 
     200            &   ' P-I slope  for diatoms                     =', pislope2 
     201         WRITE(numout,*)      & 
     202            &   ' Big particles sinking speed                =', wsbio2 
     203         WRITE(numout,*)      & 
     204            &   ' Fraction of microzoo excretion as DOM      =', sigma1 
     205         WRITE(numout,*)      & 
     206            &   ' Fraction of mesozoo excretion as DOM       =', sigma2 
     207         WRITE(numout,*)      & 
     208            &   ' Microzoo preference for POM                =', zprefc 
     209         WRITE(numout,*)      & 
     210            &   ' Microzoo preference for Nanophyto          =', zprefp 
     211         WRITE(numout,*)      & 
     212            &   ' Microzoo preference for Diatoms          =', zprefd 
     213         WRITE(numout,*)      & 
     214            &   ' Minimum Chl/C in nanophytoplankton         =', chlcnm 
     215         WRITE(numout,*)      & 
     216            &   ' Minimum Chl/C in diatoms                   =', chlcdm 
     217         WRITE(numout,*)      & 
     218            &   ' Maximum Fe/C in nanophytoplankton          =', fecnm 
     219         WRITE(numout,*)      & 
     220            &   ' Minimum Fe/C in diatoms                    =', fecdm 
     221         WRITE(numout,*)      & 
     222            &   ' Coastal release of Iron                 =', sedfeinput 
    213223      ENDIF 
    214224 
    215225      READ(numnat,natsms) 
    216226      IF(lwp) THEN 
    217           WRITE(numout,*) ' ' 
    218           WRITE(numout,*) 'natsms' 
    219           WRITE(numout,*) ' ' 
    220           WRITE(numout,*) 'Dust input from the atmosphere : ', bdustfer 
    221           WRITE(numout,*) ' ' 
    222           WRITE(numout,*) 'River input of nutrients : ', briver 
    223           WRITE(numout,*) ' ' 
    224           WRITE(numout,*) 'Atmospheric deposition of N : ', bndepo 
    225           WRITE(numout,*) ' ' 
    226           WRITE(numout,*) 'Fe input from sediments : ', bsedinput 
    227           WRITE(numout,*) ' ' 
    228       ENDIF 
    229  
    230  
    231 #else 
    232  
    233 ! no passive tracers 
    234  
    235 #endif 
     227         WRITE(numout,*) ' ' 
     228         WRITE(numout,*) 'natsms' 
     229         WRITE(numout,*) ' ' 
     230         WRITE(numout,*) 'Dust input from the atmosphere : ', bdustfer 
     231         WRITE(numout,*) ' ' 
     232         WRITE(numout,*) 'River input of nutrients : ', briver 
     233         WRITE(numout,*) ' ' 
     234         WRITE(numout,*) 'Atmospheric deposition of N : ', bndepo 
     235         WRITE(numout,*) ' ' 
     236         WRITE(numout,*) 'Fe input from sediments : ', bsedinput 
     237         WRITE(numout,*) ' ' 
     238      ENDIF 
     239 
     240   END SUBROUTINE trc_lsm 
Note: See TracChangeset for help on using the changeset viewer.