Changeset 336 for trunk/NEMO/TOP_SRC/SMS/trclsm.lobster1.h90
- Timestamp:
- 2005-11-14T13:13:26+01:00 (19 years ago)
- File:
-
- 1 edited
-
trunk/NEMO/TOP_SRC/SMS/trclsm.lobster1.h90 (modified) (9 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/SMS/trclsm.lobster1.h90
r274 r336 1 !!! 2 !!! trclsm.lobster1.h 3 !!! ********************** 4 !!! 5 !!! PURPOSE : 6 !!! --------- 7 !!! READS the specific NAMELIST for LOBSTER1 model 8 !!! 9 !! WORKSPACE : : no 10 !! --------- 11 !! 12 !! MODIFICATIONS: 13 !! -------------- 14 !! original : 99-10 (M.A. Foujols, M. Levy) passive tracer 15 !! additions : 00-12 (O. Aumont, E. Kestenare) add sediments 16 !! ---------------------------------------------------------------------- 17 !! local declarations 18 !! ================== 19 20 #if defined key_passivetrc && defined key_trc_lobster1 1 !!---------------------------------------------------------------------- 2 !! *** trclsm.lobster1.h90 *** 3 !!---------------------------------------------------------------------- 4 CONTAINS 5 6 SUBROUTINE trc_lsm 7 !!---------------------------------------------------------------------- 8 !! trclsm.lobster1.h 9 !! ********************** 10 !! 11 !! PURPOSE : 12 !! --------- 13 !! READS the specific NAMELIST for LOBSTER1 model 14 !! 15 !! WORKSPACE : : no 16 !! --------- 17 !! 18 !! MODIFICATIONS: 19 !! -------------- 20 !! original : 99-10 (M.A. Foujols, M. Levy) passive tracer 21 !! additions : 00-12 (O. Aumont, E. Kestenare) add sediments 22 !! ---------------------------------------------------------------------- 23 !! local declarations 24 !! ================== 21 25 INTEGER :: ji 22 26 CHARACTER (len=32) :: clname 23 27 24 !!--------------------------------------------------------------------- 25 !! TOP 1.0, LOCEAN-IPSL (2005) 26 !! $Header$ 27 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 28 !!--------------------------------------------------------------------- 29 30 ! 0. initializations 31 ! ------------------ 32 33 NAMELIST/natbio/apmin,azmin,anmin,admin, & 34 & redf,slopet,toptp,aknut,rcchl, & 28 !!--------------------------------------------------------------------- 29 !! OPA8, LODYC (15/11/96) 30 !!--------------------------------------------------------------------- 31 32 ! 0. initializations 33 ! ------------------ 34 NAMELIST/natbio/apmin,azmin,anmin,admin, & 35 & redf,reddom,slopet,toptp,psinut,akno3,aknh4,rcchl, & 35 36 & rgamma,toptgz,tmaxgz,rgz, & 36 37 & rppz,taus,aks,filmax,rpnaz,rdnaz,eggzoo,tauzn, & 37 38 & tmmaxp,tmminp,tmmaxz,tmminz,anumin,afdmin,taudn, & 38 & vsed,tmumax,aki,tmaxr,tminr, fdoml, taunn, taudomn,xhr, & 39 & sedlam, sedlostpoc 39 & vsed,tmumax,aki,tmaxr,tminr, taunn, taudomn,xhr, & 40 & sedlam, sedlostpoc, & 41 & fphylab,fzoolab,fdetlab,fdbod 40 42 NAMELIST/natopt/xkg0,xkr0,xkgp,xkrp,xlg,xlr,rpig 41 43 #if defined key_trc_diabio … … 43 45 #endif 44 46 45 46 IF(lwp) THEN 47 WRITE(numout,*) ' ' 48 WRITE(numout,*) ' ROUTINE trclsm' 49 WRITE(numout,*) ' **************' 50 WRITE(numout,*) ' ' 51 WRITE(numout,*) ' namelist for lobster1 model' 52 WRITE(numout,*) ' ***************************' 53 WRITE(numout,*) ' ' 47 IF(lwp) THEN 48 WRITE(numout,*) ' ' 49 WRITE(numout,*) ' ROUTINE trclsm' 50 WRITE(numout,*) ' **************' 51 WRITE(numout,*) ' ' 52 WRITE(numout,*) ' namelist for lobster1 model' 53 WRITE(numout,*) ' ***************************' 54 WRITE(numout,*) ' ' 54 55 ENDIF 55 56 … … 59 60 OPEN( numnat, FILE= clname, FORM='formatted', STATUS = 'old') 60 61 61 ! 1.4 namelist natbio : biological parameters 62 ! ------------------------------------------- 63 64 apmin = 0 65 azmin = 0 66 anmin = 0 67 admin = 0 68 redf = 0 69 slopet= 0 70 toptp = 0 71 aknut = 0 72 rcchl = 0 73 rgamma= 0 74 toptgz= 0 75 tmaxgz= 0 76 rgz = 0 77 rppz = 0 78 taus = 0 79 aks = 0 80 filmax= 0 81 rpnaz = 0 82 rdnaz = 0 83 eggzoo= 0 84 tauzn = 0 85 tmmaxp= 0 86 tmminp= 0 87 tmmaxz= 0 88 tmminz= 0 89 anumin= 0 90 afdmin= 0 91 taudn = 0 92 vsed = 0 93 tmumax= 0 94 aki = 0 62 ! 1.4 namelist natbio : biological parameters 63 ! ------------------------------------------- 64 65 apmin = 0. 66 azmin = 0. 67 anmin = 0. 68 admin = 0. 69 redf = 0. 70 reddom = 0. 71 slopet = 0. 72 toptp = 0. 73 psinut = 0. 74 akno3 = 0. 75 aknh4 = 0. 76 rcchl = 0. 77 rgamma= 0. 78 toptgz= 0. 79 tmaxgz= 0. 80 rgz = 0. 81 rppz = 0. 82 taus = 0. 83 aks = 0. 84 filmax= 0. 85 rpnaz = 0. 86 rdnaz = 0. 87 eggzoo= 0. 88 tauzn = 0. 89 tmmaxp= 0. 90 tmminp= 0. 91 tmmaxz= 0. 92 tmminz= 0. 93 anumin= 0. 94 afdmin= 0. 95 taudn = 0. 96 vsed = 0. 97 tmumax= 0. 98 aki = 0. 95 99 tmaxr = 1./( 4.*rday)*0. 96 100 tminr = 1./(24.*30.*rday)*0. 97 fdoml = 0 98 xhr=0 99 sedlam=0 100 sedlostpoc=0 101 taudomn = 0 102 taunn = 0 101 xhr=0. 102 sedlam=0. 103 sedlostpoc=0. 104 taudomn = 0. 105 taunn = 0. 106 fphylab = 0. 107 fzoolab = 0. 108 fdetlab = 0. 109 fdbod = 0. 103 110 104 111 READ(numnat,natbio) … … 122 129 & ' optimal photosynthesis temperature toptp =', toptp 123 130 WRITE(numout,*) & 124 & ' half-saturation nutrient aknut =', aknut 131 & ' inhibition of no3 uptake by nh4 psinut =', psinut 132 WRITE(numout,*) & 133 & ' half-saturation nutrient for no3 uptake akno3 =', akno3 134 WRITE(numout,*) & 135 & ' half-saturation nutrient for nh4 uptake aknh4 =', aknh4 125 136 WRITE(numout,*) & 126 137 & ' carbone/chlorophyl ratio rcchl =', rcchl … … 181 192 & ' damping-remineralisation rate tminr =', tminr 182 193 WRITE(numout,*) & 183 & ' damping-remineralisation rate fdoml =', fdoml184 WRITE(numout,*) &185 194 & ' nitrification rate taunn =', taunn 186 195 WRITE(numout,*) & … … 192 201 WRITE(numout,*) & 193 202 & ' Sediment geol loss for POC sedlostpoc =', sedlostpoc 203 WRITE(numout,*) & 204 & ' NH4 fraction of phytoplankton exsudation fphylab =', fphylab 205 WRITE(numout,*) & 206 & ' NH4 fraction of zooplankton excretion fzoolab =', fzoolab 207 WRITE(numout,*) & 208 & ' NH4 fraction of detritus dissolution fdetlab =', fdetlab 209 WRITE(numout,*) & 210 & ' Zooplankton mortality fraction that goes to detritus fdbod =', fdbod 194 211 ENDIF 195 212 196 ! 1.5 namelist natopt : parameters for optic197 ! ------------------------------------------213 ! 1.5 namelist natopt : parameters for optic 214 ! ------------------------------------------ 198 215 199 216 xkg0 = 0. … … 208 225 209 226 IF(lwp) THEN 210 211 212 213 214 215 216 217 218 219 227 WRITE(numout,*) 'natopt' 228 WRITE(numout,*) ' ' 229 WRITE(numout,*) ' green water absorption coeff xkg0 = ',xkg0 230 WRITE(numout,*) ' red water absorption coeff xkr0 = ',xkr0 231 WRITE(numout,*) ' pigment red absorption coeff xkrp = ',xkrp 232 WRITE(numout,*) ' pigment green absorption coeff xkgp = ',xkgp 233 WRITE(numout,*) ' green chl exposant xlg = ',xlg 234 WRITE(numout,*) ' red chl exposant xlr = ',xlr 235 WRITE(numout,*) ' chla/chla+phea ratio rpig = ',rpig 236 WRITE(numout,*) ' ' 220 237 221 238 ENDIF … … 223 240 #if defined key_trc_diabio 224 241 225 ! NAMELIST : natdbi226 227 ! default name for biological trends : short and long name, units242 ! NAMELIST : natdbi 243 244 ! default name for biological trends : short and long name, units 228 245 229 246 DO ji=1,jpdiabio 230 IF (ji < 10) THEN247 IF (ji < 10) THEN 231 248 WRITE (ctrbio(ji),'("BIO_",I1)') ji 232 ELSE IF (ji < 100) THEN249 ELSE IF (ji < 100) THEN 233 250 WRITE (ctrbio(ji),'("BIO_",I2)') ji 234 ELSE251 ELSE 235 252 WRITE (ctrbio(ji),'("BIO_",I3)') ji 236 ENDIF237 WRITE (ctrbil(ji),'("BIOLOGICAL TREND NUMBER ",I2)') ji238 ctrbiu(ji)='mmoleN/m3/s '239 END DO 253 ENDIF 254 WRITE (ctrbil(ji),'("BIOLOGICAL TREND NUMBER ",I2)') ji 255 ctrbiu(ji)='mmoleN/m3/s ' 256 END DO 240 257 241 258 nwritebio = 10 … … 244 261 245 262 IF(lwp) THEN 246 247 248 249 & ' frequency of outputs for biological outputs = ' &250 & ,nwritebio251 252 263 WRITE(numout,*) 'natdbi' 264 WRITE(numout,*) ' ' 265 WRITE(numout,*) & 266 & ' frequency of outputs for biological outputs = ' & 267 & ,nwritebio 268 WRITE(numout,*) ' ' 269 DO ji=1,jpdiabio 253 270 WRITE(numout,*) & 254 & 'name of biological trend number :',ji,' : ',ctrbio(ji)271 & 'name of biological trend number :',ji,' : ',ctrbio(ji) 255 272 WRITE(numout,*) ctrbil(ji) 256 273 WRITE(numout,*) ' in unit = ',ctrbiu(ji) 257 END DO258 END IF 274 END DO 275 END IF 259 276 #endif 260 277 261 #else 262 263 ! no passive tracers 264 265 #endif 278 END SUBROUTINE trc_lsm
Note: See TracChangeset
for help on using the changeset viewer.
