- Timestamp:
- 2006-09-12T13:15:31+02:00 (18 years ago)
- Location:
- trunk/NEMO/TOP_SRC/SMS
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/SMS/p4zflx.F
r341 r502 120 120 kgco2(ji,jj) = kgco2(ji,jj)/(100.*3600.) 121 121 & *(1-freeze(ji,jj))*tmask(ji,jj,1) 122 # if defined key_off_degrad 123 & *facvol(ji,jj,1) 124 # endif 122 125 C 123 126 END DO … … 137 140 kgo2(ji,jj) = (0.3*ws*ws + 2.5*(0.5246+ttc*(0.016256+ 138 141 & ttc*0.00049946)))*sqrt(660./schmitto2) 142 139 143 C 140 144 C CONVERT TO m/s AND APPLY SEA ICE COVER … … 142 146 C 143 147 kgo2(ji,jj) = kgo2(ji,jj)/(100.*3600.) 144 $ *(1-freeze(ji,jj))*tmask(ji,jj,1) 148 & *(1-freeze(ji,jj))*tmask(ji,jj,1) 149 # if defined key_off_degrad 150 & *facvol(ji,jj,1) 151 # endif 145 152 C 146 153 ENDDO -
trunk/NEMO/TOP_SRC/SMS/p4zsed.F
r341 r502 53 53 REAL xlim,xconctmp2,zstep,zfact 54 54 REAL irondep(jpi,jpj,jpk),sidep(jpi,jpj) 55 REAL zvol 55 56 CC 56 57 C … … 112 113 END DO 113 114 C 115 116 C 114 117 C Add the external input of iron which is 3D distributed 115 118 C (dust, river and sediment mobilization) … … 234 237 xlim=(1.-xnanono3(ji,jj,jk)-xnanonh4(ji,jj,jk)) 235 238 if (xlim.le.0.2) xlim=0.01 236 nitrpot(ji,jj,jk)=max(0.,(prmax(ji,jj,jk)-2.15/rjjss)) 239 nitrpot(ji,jj,jk)=max(0.,(0.6*tgfunc(ji,jj,jk)-2.15)/rjjss) 240 #if defined key_off_degrad 241 & *facvol(ji,jj,jk) 242 #endif 237 243 & *xlim*rfact2*trn(ji,jj,jk,jpfer)/(conc3 238 244 & +trn(ji,jj,jk,jpfer))*(1.-exp(-etot(ji,jj,jk)/50.)) … … 241 247 END DO 242 248 C 243 249 nitrpottot=0. 244 250 DO jk=1,jpkm1 245 251 DO jj=2,jpj-1 … … 259 265 DO jj=1,jpj 260 266 DO ji=1,jpi 267 #if ! defined key_cfg_1d && ( defined key_orca_r4 || defined key_orca_r2 || defined key_orca_r05 || defined key_orca_r025 ) 268 zfact=nitrpot(ji,jj,jk)*denitot/nitrpottot 269 #else 261 270 zfact=nitrpot(ji,jj,jk)*1.E-7 271 #endif 262 272 trn(ji,jj,jk,jpnh4)=trn(ji,jj,jk,jpnh4)+zfact 263 273 trn(ji,jj,jk,jpoxy)=trn(ji,jj,jk,jpoxy)+zfact*o2nit
Note: See TracChangeset
for help on using the changeset viewer.