diff trcini_pisces.F90 $PISCES
333,334c333
< IF(lwm) CALL ctl_opn( numonp , 'output.namelist.pis' &
< & , 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
---
> IF(lwm) CALL ctl_opn( numonp , 'output.namelist.pis' , 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
diff step3d_t.F $OCEAN
539,543d538
< ! modif SPOUS ASAP
< if (dRz .eq. 0) then
< ! write(*,*) 'step3d_t mynode', mynode
< cff = -1.E-14
< else
545,547d539
< endif
< ! cff = min( 1./dRz,-1.E-14 ) ! minimum stratification imposed
< ! modif SPOUS ASAP
cpv2 // cpv1 // Julien
wsbio = 2 // 3 // 2 :: POC sinking speed
ferat3 = 10e-6 // 3e-6 // 10e-6 :: Fe/C in zooplankton
kmort = 2e-7 // 1e-7 // 2e-7 :: half saturation constant for mortality
niter2max = 50 // 15 // 1
conc1, concnfer ... ??
concdnh4 = 3e-7 // 5e-7 // 3e-7 :: NH4 half saturation for phyto, diatoms
xksi1 = 2e-6 // 1.5e-6 // 2e-6 :: half saturation constant for Si uptake
xksi2 = 20e-6 // 4e-6 // 20e-6:: half saturation constant for Si/C
pislopen(pisloped) = 2 // (pislope2) 3 // 2 :: initial slope of P-I curve
excret = 0.05 // 0.1 // 0.05 :: delta' in iron in nano and diatoms
fecnm (fecdm) = 80e-6 // 15 (20) e-6 // 40e-6 :: iron quota in nano (diatoms)
grosip = 0.159 // 0.06 // 0.159 :: Optimal Si/C uptake ratio of diatoms
wchl = 0.01 // 0.001 // 0.01 :: nanophytoplankton aggregation term
wchld = 0.01 // 0.02 // 0.01 :: diatoms aggregation term
xlam1 = 0.15 // 0.005 // 0.005 :: Slope of the scavenging rate of iron
