Changeset 773 for branches/dev_001_GM/NEMO/TOP_SRC/trcdtr.F90
- Timestamp:
- 2007-12-17T17:55:32+01:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_001_GM/NEMO/TOP_SRC/trcdtr.F90
r772 r773 41 41 !! ** Method : 42 42 !!--------------------------------------------------------------------- 43 INTEGER :: ji,jj,jk,jn 44 #if defined key_pisces 43 INTEGER :: jn 44 # if defined key_pisces 45 INTEGER :: ji, jj, jk 45 46 REAL(wp) :: alka0, oxyg0, calc0, bioma0 46 47 REAL(wp) :: silic1, po4, no3, caralk, bicarb, co3 47 # endif48 # endif 48 49 !!--------------------------------------------------------------------- 49 50 … … 52 53 IF(lwp) WRITE(numout,*) '~~~~~~~' 53 54 54 #if defined key_cfc 55 ! CFC initialisation 56 trn(:,:,:,:) = 0.e0 55 IF( lk_cfc ) THEN 56 DO jn = jp_cfc0,jp_cfc1 57 trn(:,:,:,jn) = 0.e0 ! CFC initialisation to zera 58 END DO 59 ENDIF 57 60 58 # elif defined key_pisces61 # if defined key_pisces 59 62 ! PISCES initialisation 60 63 ! --------------------- … … 75 78 trn(:,:,:,jppo4) = po4 76 79 trn(:,:,:,jppoc) = bioma0 77 # if ! defined key_kriest80 # if ! defined key_kriest 78 81 trn(:,:,:,jpgoc) = bioma0 79 82 trn(:,:,:,jpbfe) = bioma0 * 5.e-6 80 # else83 # else 81 84 trn(:,:,:,jpnum) = bioma0 / ( 6. *xkr_massp ) 82 # endif85 # endif 83 86 trn(:,:,:,jpsil) = silic1 84 87 trn(:,:,:,jpbsi) = bioma0 * 0.15 … … 119 122 IF(lwp) WRITE(numout,*) ' ' 120 123 121 # elif defined key_lobster && ( defined key_eel_r6 || defined key_eel_r2 )124 # elif defined key_lobster && ( defined key_eel_r6 || defined key_eel_r2 ) 122 125 ! LOBSTER initialisation for EEL 123 126 ! ---------------------- … … 195 198 trn(:,:,30,jpno3) = 20.01 * tmask(:,:,30) 196 199 197 # elif defined key_lobster && defined key_gyre200 # elif defined key_lobster && defined key_gyre 198 201 ! LOBSTER initialisation for GYRE 199 202 ! ---------------------- … … 216 219 END DO 217 220 218 # else221 # else 219 222 ! Default case 220 223 ! ------------ … … 223 226 END DO 224 227 225 # endif226 227 # if defined key_dtatrc228 # endif 229 230 # if defined key_dtatrc 228 231 ! Initialization of tracer from a file that may also be used for damping 229 232 CALL dta_trc( nittrc000 ) … … 231 234 IF( lutini(jn) ) trn(:,:,:,jk) = trdta(:,:,:,jn) * tmask(:,:,:) ! initialisation from file if required 232 235 END DO 233 # endif236 # endif 234 237 235 238 ! before field : … … 237 240 trb(:,:,:,:) = trn(:,:,:,:) 238 241 239 # if defined key_lobster242 # if defined key_lobster 240 243 ! initialize the POC in sediments 241 244 sedpocb(:,:) = 0.e0 242 245 sedpocn(:,:) = 0.e0 243 246 sedpoca(:,:) = 0.e0 244 # endif247 # endif 245 248 ! 246 249 END SUBROUTINE trc_dtr
Note: See TracChangeset
for help on using the changeset viewer.