Changeset 617 for trunk/NEMO/TOP_SRC/SMS/trclsm.pisces.h90

Timestamp:

2007-02-21T14:03:10+01:00 (17 years ago)

Author:

opalod

Message:

* empty log message *

File:

• trunk/NEMO/TOP_SRC/SMS/trclsm.pisces.h90 (modified) (7 diffs)

trunk/NEMO/TOP_SRC/SMS/trclsm.pisces.h90

@@ -34 +34 @@
       ! ------------------
       !
-      NAMELIST/natgas/ gasfac, igaswind, icice
       NAMELIST/natext/ atcco2
-      NAMELIST/natbio/caco3r,                                          &
+      NAMELIST/natbio/caco3r, kdca, nca, part,                       &
          &          dispo0,conc0,oxymin,grosip, sedlam,              &
          &          sedlostpoc,sedlostcal,sedlostsil,nrdttrc,        &
@@ -49 +48 @@
          &          chlcnm,chlcdm, sedfeinput
       NAMELIST/natsms/bdustfer, briver, bndepo, bsedinput
+#if defined key_trc_kriest
+      NAMELIST/natkriest/xkr_eta,xkr_zeta,xkr_sfact,xkr_mass_min,xkr_mass_max,&
+         &               xkr_dnano,xkr_ddiat,xkr_dmeso,xkr_daggr,xkr_stick
+#endif
 
       ! initialize the number of LOGICAL UNIT used
@@ -68 +71 @@
 
 
-      ! 1 Namelist natgas :
-      ! -------------------
-
-      READ(numnat,natgas)
-
-      IF(lwp) THEN
-         WRITE(numout,*) ' '
-         WRITE(numout,*) 'natgas'
-         write(numout,*) 'gasfac = ',gasfac
-         WRITE(numout,*) ' '
-         write(numout,*) 'igaswind = ',igaswind
-         WRITE(numout,*) ' '
-         write(numout,*) 'icice = ',icice
-         WRITE(numout,*) ' '
-      ENDIF
-
-      ! 2 Namelist natext :
+      ! 1 Namelist natext :
       ! -------------------
       READ(numnat,natext)
@@ -104 +91 @@
             &   ' mean rainratio                             =', caco3r
          WRITE(numout,*)      &
+            &   ' diss. rate constant calcite (per month)    =', kdca
+         WRITE(numout,*)      &
+            &   ' order of reaction for calcite dissolution  =', nca
+         WRITE(numout,*)      &
+            &   ' part of calcite not dissolved in guts      =', part
+         WRITE(numout,*)      &
             &   ' mean Si/C ratio                            =', grosip
          WRITE(numout,*)      &
@@ -112 +105 @@
             &   ' Sediment bioturbation factor               =', sedlam
          WRITE(numout,*)      &
-            &   ' Sediment burying ratio for POC         =', sedlostpoc
-         WRITE(numout,*)      &
-            &   ' Sediment burying ratio for CACO3       =', sedlostcal
-         WRITE(numout,*)      &
-            &   ' Sediment burying ratio for SI          =', sedlostsil
+            &   ' Sediment burying ratio for POC             =', sedlostpoc
+         WRITE(numout,*)      &
+            &   ' Sediment burying ratio for CACO3           =', sedlostcal
+         WRITE(numout,*)      &
+            &   ' Sediment burying ratio for SI              =', sedlostsil
          WRITE(numout,*)      &
             &   ' frequence pour la biologie                 =', nrdttrc
@@ -125 +118 @@
          WRITE(numout,*)      &
             &   ' POC sinking speed                          =', wsbio
+         WRITE(numout,*)      &
+            &   ' quadratic mortality of phytoplankton       =', wchl
+         WRITE(numout,*)      &
+            &   ' maximum quadratic mortality of diatoms     =', wchld
          WRITE(numout,*)      &
             &   ' exsudation rate of zooplankton             =', resrat
@@ -240 +237 @@
       ENDIF
 
+#if defined key_trc_kriest
+
+      READ(numnat,natkriest)
+      IF(lwp) THEN
+         WRITE(numout,*) ' '
+         WRITE(numout,*) 'natkriest'
+         WRITE(numout,*) ' '
+         WRITE(numout,*) 'Sinking  exponent                    xkr_eta       = ', xkr_eta 
+         WRITE(numout,*) 'N content exponent                   xkr_zeta      = ', xkr_zeta
+         WRITE(numout,*) 'Sinking factor                       xkr_sfact     = ', xkr_sfact
+         WRITE(numout,*) 'Stickiness                           xkr_stick     = ', xkr_stick
+         WRITE(numout,*) 'Minimum mass for Aggregates          xkr_mass_min  = ', xkr_mass_min
+         WRITE(numout,*) 'Maximum mass for Aggregates          xkr_mass_max  = ', xkr_mass_max
+         WRITE(numout,*) 'Size of particles in nano pool       xkr_dnano     = ', xkr_dnano
+         WRITE(numout,*) 'Size of particles in diatoms pool    xkr_ddiat     = ', xkr_ddiat
+         WRITE(numout,*) 'Size of particles in mesozoo pool    xkr_dmeso     = ', xkr_dmeso
+         WRITE(numout,*) 'Size of particles in aggregates pool xkr_daggr     = ', xkr_daggr
+     ENDIF
+
+
+     ! Computation of some variables
+     xkr_massp     = 5.7E-6 * 7.6 * xkr_mass_min**xkr_zeta
+     ! max and min vertical particle speed
+     xkr_wsbio_min = xkr_sfact * xkr_mass_min**xkr_eta
+     xkr_wsbio_max = xkr_sfact * xkr_mass_max**xkr_eta
+     WRITE(numout,*) ' max and min vertical particle speed ',xkr_wsbio_min,xkr_wsbio_max
+     !
+     !    effect of the sizes of the different living pools on particle numbers
+     !    nano = 2um-20um -> mean size=6.32 um -> ws=2.596 -> xnum=xnnano=2.337
+     !    diat and microzoo = 10um-200um -> 44.7 -> 8.732 -> xnum=xndiat=3.718
+     !    mesozoo = 200um-2mm -> 632.45 -> 45.14 -> xnum=xnmeso=7.147
+     !    aggregates = 200um-10mm -> 1414 -> 74.34 -> xnum=xnaggr=9.877
+     !    doc aggregates = 1um
+     ! ----------------------------------------------------------
+
+     xkr_nnano = 1. / ( xkr_massp * xkr_dnano )
+     xkr_ndiat = 1. / ( xkr_massp * xkr_ddiat )
+     xkr_nmeso = 1. / ( xkr_massp * xkr_dmeso )
+     xkr_naggr = 1. / ( xkr_massp * xkr_daggr )
+    
+
+#endif
+
    END SUBROUTINE trc_lsm