- Timestamp:
- 2005-12-08T11:42:00+01:00 (19 years ago)
- Location:
- trunk/NEMO/TOP_SRC
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/SMS/trcctl.cfc.h90
r341 r350 1 !! 2 !! Modifications: 3 !! -------------- 4 !! 05-10 (C. Ethe): 5 !! assign a parameter to name individual tracers 6 !! 7 !!---------------------------------------------------------------------- 8 !! TOP 1.0 , LOCEAN-IPSL (2005) 9 !! $Header$ 10 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 11 !!---------------------------------------------------------------------- 1 ! $Id$ 2 !!! 3 !!! Modifications: 4 !!! -------------- 5 !!! 05-10 (C. Ethe): 6 !!! assign a parameter to name individual tracers 7 !!! 12 8 13 9 IF(lwp) THEN … … 16 12 ENDIF 17 13 18 19 14 ! Check number of tracers 15 ! ----------------------- 20 16 21 17 IF (jptra .GT. 2) THEN … … 29 25 END IF 30 26 31 ! Check tracer names 32 ! ------------------ 33 27 ! Check tracer names 28 ! ------------------ 34 29 IF ( jptra == 1 ) THEN 35 30 IF ( jp11 == 1 ) THEN 36 IF ( ctrcnm(jp11) /= 'CFC -11') THEN37 ctrcnm(jp11)='CFC -11'31 IF ( ctrcnm(jp11) /= 'CFC11') THEN 32 ctrcnm(jp11)='CFC11' 38 33 ctrcnl(jp11)='Carbonate 11 concentration' 39 34 ENDIF 40 35 ENDIF 41 36 IF ( jp12 == 1 ) THEN 42 IF ( ctrcnm(jp12) /= 'CFC -12') THEN43 ctrcnm(jp12)='CFC -12'37 IF ( ctrcnm(jp12) /= 'CFC12') THEN 38 ctrcnm(jp12)='CFC12' 44 39 ctrcnl(jp12)='Carbonate 12 concentration' 45 40 ENDIF … … 48 43 49 44 IF ( jptra == 2 ) THEN 50 IF ( ctrcnm(jp11) /= 'CFC -11' .OR. ctrcnm(jp12) /= 'CFC-12' ) THEN51 ctrcnm(jp11)='CFC -11'45 IF ( ctrcnm(jp11) /= 'CFC11' .OR. ctrcnm(jp12) /= 'CFC12' ) THEN 46 ctrcnm(jp11)='CFC11' 52 47 ctrcnl(jp11)='Carbonate 11 concentration' 53 ctrcnm(jp12)='CFC -12'48 ctrcnm(jp12)='CFC12' 54 49 ctrcnl(jp12)='Carbonate 12 concentration' 55 50 ENDIF … … 67 62 68 63 69 ! Check tracer units 70 ! ------------------ 64 ! Check tracer units 71 65 72 66 DO jn=1,jptra -
trunk/NEMO/TOP_SRC/trcrst.F90
r340 r350 30 30 CHARACTER (len=48) :: & 31 31 trestart = 'initial.trc.nc' ! restart file name 32 33 !! * Substitutions 34 # include "passivetrc_substitute.h90" 32 35 33 36 CONTAINS … … 85 88 86 89 REAL(wp), DIMENSION(3) :: zinfo 90 91 #if defined key_trc_pisces && ( defined key_orca_r4 || defined key_orca_r2 || defined key_orca_r05 || defined key_orca_r025 ) 92 REAL(wp) , DIMENSION(jpi,jpj,jpk) :: zvolk 93 REAL(wp) :: zareatot, zpo4tot 94 #endif 87 95 88 96 !!--------------------------------------------------------------------- … … 238 246 239 247 #if defined key_trc_pisces 248 249 #if defined key_orca_r4 || defined key_orca_r2 || defined key_orca_r05 || defined key_orca_r025 250 251 zvolk(:,:,:) = 0. 252 zareatot = 0. 253 DO jk = 1, jpkm1 254 DO jj = 2, jpjm1 255 DO ji = 2, jpim1 256 zvolk(ji,jj,jk) = e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) * tmask(ji,jj,jk) 257 zareatot = zareatot + zvolk(ji,jj,jk) 258 ENDDO 259 ENDDO 260 ENDDO 261 262 263 zpo4tot = 0. 264 DO jk = 1, jpkm1 265 DO jj = 2, jpjm1 266 DO ji = 2, jpim1 267 zpo4tot = zpo4tot + trn(ji,jj,jk,jptal) * zvolk(ji,jj,jk) 268 END DO 269 END DO 270 END DO 271 272 WRITE(0,*) 'TALK moyen ', zpo4tot/zareatot*1E6 273 zpo4tot = zpo4tot/zareatot*1E6 274 trn(:,:,:,jptal) = trn(:,:,:,jptal)*2391./zpo4tot 275 276 277 zpo4tot = 0. 278 DO jk = 1, jpkm1 279 DO jj = 2, jpjm1 280 DO ji = 2, jpim1 281 zpo4tot = zpo4tot + trn(ji,jj,jk,jppo4) * zvolk(ji,jj,jk) 282 END DO 283 END DO 284 END DO 285 286 WRITE(0,*) 'PO4 moyen ', zpo4tot/zareatot*1E6/122. 287 zpo4tot = zpo4tot/zareatot*1E6/122. 288 trn(:,:,:,jppo4) = trn(:,:,:,jppo4)*2.165/zpo4tot 289 290 291 zpo4tot = 0. 292 DO jk = 1, jpkm1 293 DO jj = 2, jpjm1 294 DO ji = 2, jpim1 295 zpo4tot = zpo4tot + trn(ji,jj,jk,jpno3) * zvolk(ji,jj,jk) 296 END DO 297 END DO 298 END DO 299 300 WRITE(0,*) 'NO3 moyen ', zpo4tot/zareatot*1E6/7.6 301 zpo4tot = zpo4tot/zareatot*1E6/7.6 302 trn(:,:,:,jpno3) = trn(:,:,:,jpno3)*30.9/zpo4tot 303 304 zpo4tot = 0. 305 DO jk = 1, jpkm1 306 DO jj = 2, jpjm1 307 DO ji = 2, jpim1 308 zpo4tot = zpo4tot + trn(ji,jj,jk,jpsil) * zvolk(ji,jj,jk) 309 END DO 310 END DO 311 END DO 312 313 WRITE(0,*) 'SiO3 moyen ', zpo4tot/zareatot*1E6 314 zpo4tot = zpo4tot/zareatot*1E6 315 trn(:,:,:,jpsil) = MIN( 400E-6,trn(:,:,:,jpsil)*91.51/zpo4tot) 316 317 #endif 240 318 !! Initialization of chemical variables of the carbon cycle 241 319 !! -------------------------------------------------------- 242 243 320 DO jk = 1,jpk 244 321 DO jj = 1,jpj -
trunk/NEMO/TOP_SRC/trcrst_fdir.h90
r340 r350 7 7 !!------------------------------------------------------------------------ 8 8 !! TOP 1.0, LOCEAN-IPSL (2005) 9 !! $Header$10 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt9 !! $Header$ 10 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 11 11 !!---------------------------------------------------------------------- 12 12
Note: See TracChangeset
for help on using the changeset viewer.