- Timestamp:
- 2018-06-21T11:58:42+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_nemo2cice_prints/NEMOGCM/NEMO/TOP_SRC/trcnam.F90
r9816 r9817 11 11 !! - ! 2001-01 (E Kestenare) suppress ndttrc=1 for CEN2 and TVD schemes 12 12 !! 1.0 ! 2005-03 (O. Aumont, A. El Moussaoui) F90 13 !! - ! 2014-06 (A. Yool, J. Palmieri) adding MEDUSA-2 13 14 !!---------------------------------------------------------------------- 14 15 #if defined key_top … … 25 26 USE trcnam_c14b ! C14 SMS namelist 26 27 USE trcnam_my_trc ! MY_TRC SMS namelist 28 USE trcnam_medusa ! MEDUSA namelist 29 USE trcnam_idtra ! Idealise tracer namelist 30 USE trcnam_age ! AGE SMS namelist 27 31 USE trd_oce 28 32 USE trdtrc_oce … … 54 58 !! ** Method : - read passive tracer namelist 55 59 !! - read namelist of each defined SMS model 56 !! ( (PISCES, CFC, MY_TRC )57 !!--------------------------------------------------------------------- 58 INTEGER :: jn ! dummy loop indice60 !! ( (PISCES, CFC, MY_TRC, MEDUSA, IDTRA, Age ) 61 !!--------------------------------------------------------------------- 62 INTEGER :: jn, jk ! dummy loop indice 59 63 ! ! Parameters of the run 60 64 IF( .NOT. lk_offline ) CALL trc_nam_run 61 65 62 66 ! ! passive tracer informations 67 # if defined key_debug_medusa 68 CALL flush(numout) 69 IF (lwp) write (numout,*) '------------------------------' 70 IF (lwp) write (numout,*) 'Jpalm - debug' 71 IF (lwp) write (numout,*) 'in trc_nam - just before CALL trc_nam_trc' 72 IF (lwp) write (numout,*) ' ' 73 # endif 74 ! 63 75 CALL trc_nam_trc 64 76 65 77 ! ! Parameters of additional diagnostics 78 # if defined key_debug_medusa 79 CALL flush(numout) 80 IF (lwp) write (numout,*) '------------------------------' 81 IF (lwp) write (numout,*) 'Jpalm - debug' 82 IF (lwp) write (numout,*) 'CALL trc_nam_trc -- OK' 83 IF (lwp) write (numout,*) 'in trc_nam - just before CALL trc_nam_dia' 84 IF (lwp) write (numout,*) ' ' 85 # endif 86 ! 87 66 88 CALL trc_nam_dia 67 89 68 90 ! ! namelist of transport 91 # if defined key_debug_medusa 92 CALL flush(numout) 93 IF (lwp) write (numout,*) '------------------------------' 94 IF (lwp) write (numout,*) 'Jpalm - debug' 95 IF (lwp) write (numout,*) 'CALL trc_nam_dia -- OK' 96 IF (lwp) write (numout,*) 'in trc_nam - just before CALL trc_nam_trp' 97 IF (lwp) write (numout,*) ' ' 98 # endif 99 ! 69 100 CALL trc_nam_trp 101 ! 102 # if defined key_debug_medusa 103 CALL flush(numout) 104 IF (lwp) write (numout,*) '------------------------------' 105 IF (lwp) write (numout,*) 'Jpalm - debug' 106 IF (lwp) write (numout,*) 'CALL trc_nam_trp -- OK' 107 IF (lwp) write (numout,*) 'continue trc_nam ' 108 IF (lwp) write (numout,*) ' ' 109 CALL flush(numout) 110 # endif 111 ! 70 112 71 113 … … 89 131 END DO 90 132 WRITE(numout,*) ' ' 133 # if defined key_debug_medusa 134 CALL flush(numout) 135 # endif 91 136 ENDIF 92 137 … … 107 152 WRITE(numout,*) 108 153 ENDIF 109 ENDIF 110 154 # if defined key_debug_medusa 155 CALL flush(numout) 156 # endif 157 ENDIF 158 159 # if defined key_debug_medusa 160 DO jk = 1, jpk 161 WRITE(numout,*) ' level number: ', jk, 'rdttrc: ',rdttrc(jk),'rdttra: ', rdttra(jk),'nn_dttrc: ', nn_dttrc 162 END DO 163 CALL flush(numout) 164 # endif 111 165 112 166 rdttrc(:) = rdttra(:) * FLOAT( nn_dttrc ) ! vertical profile of passive tracer time-step … … 116 170 WRITE(numout,*) ' Passive Tracer time step rdttrc = ', rdttrc(1) 117 171 WRITE(numout,*) 172 # if defined key_debug_medusa 173 CALL flush(numout) 174 # endif 118 175 ENDIF 119 176 … … 143 200 IF( ln_trdtrc(jn) ) WRITE(numout,*) ' compute ML trends for tracer number :', jn 144 201 END DO 202 WRITE(numout,*) ' ' 203 CALL flush(numout) 145 204 ENDIF 146 205 #endif 147 206 207 # if defined key_debug_medusa 208 CALL flush(numout) 209 IF (lwp) write (numout,*) '------------------------------' 210 IF (lwp) write (numout,*) 'Jpalm - debug' 211 IF (lwp) write (numout,*) 'just before ice module for tracers call : ' 212 IF (lwp) write (numout,*) ' ' 213 # endif 214 ! 148 215 149 216 ! Call the ice module for tracers 150 217 ! ------------------------------- 151 218 CALL trc_nam_ice 219 220 # if defined key_debug_medusa 221 CALL flush(numout) 222 IF (lwp) write (numout,*) '------------------------------' 223 IF (lwp) write (numout,*) 'Jpalm - debug' 224 IF (lwp) write (numout,*) 'Will now read SMS namelists : ' 225 IF (lwp) write (numout,*) ' ' 226 # endif 227 ! 152 228 153 229 ! namelist of SMS … … 156 232 ELSE ; IF(lwp) WRITE(numout,*) ' PISCES not used' 157 233 ENDIF 158 234 ! 235 # if defined key_debug_medusa 236 CALL flush(numout) 237 IF (lwp) write (numout,*) '------------------------------' 238 IF (lwp) write (numout,*) 'Jpalm - debug' 239 IF (lwp) write (numout,*) 'CALL trc_nam_pisces -- OK' 240 IF (lwp) write (numout,*) 'in trc_nam - just before CALL trc_nam_medusa' 241 IF (lwp) write (numout,*) ' ' 242 # endif 243 ! 244 IF( lk_medusa ) THEN ; CALL trc_nam_medusa ! MEDUSA tracers 245 ELSE ; IF(lwp) WRITE(numout,*) ' MEDUSA not used' 246 ENDIF 247 ! 248 # if defined key_debug_medusa 249 CALL flush(numout) 250 IF (lwp) write (numout,*) '------------------------------' 251 IF (lwp) write (numout,*) 'Jpalm - debug' 252 IF (lwp) write (numout,*) 'CALL trc_nam_medusa -- OK' 253 IF (lwp) write (numout,*) 'in trc_nam - just before CALL trc_nam_idtra' 254 IF (lwp) write (numout,*) ' ' 255 # endif 256 ! 257 IF( lk_idtra ) THEN ; CALL trc_nam_idtra ! Idealize tracers 258 ELSE ; IF(lwp) WRITE(numout,*) ' Idealize tracers not used' 259 ENDIF 260 ! 261 # if defined key_debug_medusa 262 CALL flush(numout) 263 IF (lwp) write (numout,*) '------------------------------' 264 IF (lwp) write (numout,*) 'Jpalm - debug' 265 IF (lwp) write (numout,*) 'CALL trc_nam_idtra -- OK' 266 IF (lwp) write (numout,*) 'in trc_nam - just before CALL trc_nam_cfc' 267 IF (lwp) write (numout,*) ' ' 268 # endif 269 ! 159 270 IF( lk_cfc ) THEN ; CALL trc_nam_cfc ! CFC tracers 160 271 ELSE ; IF(lwp) WRITE(numout,*) ' CFC not used' 161 272 ENDIF 162 273 ! 274 # if defined key_debug_medusa 275 CALL flush(numout) 276 IF (lwp) write (numout,*) '------------------------------' 277 IF (lwp) write (numout,*) 'Jpalm - debug' 278 IF (lwp) write (numout,*) 'CALL trc_nam_cfc -- OK' 279 IF (lwp) write (numout,*) 'in trc_nam - just before CALL trc_nam_c14' 280 IF (lwp) write (numout,*) ' ' 281 # endif 282 ! 163 283 IF( lk_c14b ) THEN ; CALL trc_nam_c14b ! C14 bomb tracers 164 284 ELSE ; IF(lwp) WRITE(numout,*) ' C14 not used' 165 285 ENDIF 166 286 ! 287 # if defined key_debug_medusa 288 CALL flush(numout) 289 IF (lwp) write (numout,*) '------------------------------' 290 IF (lwp) write (numout,*) 'Jpalm - debug' 291 IF (lwp) write (numout,*) 'CALL trc_nam_c14 -- OK' 292 IF (lwp) write (numout,*) 'in trc_nam - just before CALL trc_nam_age' 293 IF (lwp) write (numout,*) ' ' 294 # endif 295 ! 296 IF( lk_age ) THEN ; CALL trc_nam_age ! AGE tracer 297 ELSE ; IF(lwp) WRITE(numout,*) ' AGE not used' 298 ENDIF 299 ! 300 # if defined key_debug_medusa 301 CALL flush(numout) 302 IF (lwp) write (numout,*) '------------------------------' 303 IF (lwp) write (numout,*) 'Jpalm - debug' 304 IF (lwp) write (numout,*) 'CALL trc_nam_age -- OK' 305 IF (lwp) write (numout,*) 'in trc_nam - CALL trc_nam -- OK' 306 IF (lwp) write (numout,*) ' ' 307 # endif 308 ! 167 309 IF( lk_my_trc ) THEN ; CALL trc_nam_my_trc ! MY_TRC tracers 168 310 ELSE ; IF(lwp) WRITE(numout,*) ' MY_TRC not used' 169 311 ENDIF 170 ! 312 313 IF(lwp) CALL flush(numout) 171 314 END SUBROUTINE trc_nam 172 315 … … 216 359 WRITE(numout,*) ' Use euler integration for TRC (y/n) ln_top_euler = ', ln_top_euler 217 360 WRITE(numout,*) ' ' 361 CALL flush(numout) 218 362 ENDIF 219 363 ! … … 306 450 ln_trc_wri(jn) = sn_tracer(jn)%llsave 307 451 END DO 308 452 IF(lwp) CALL flush(numout) 453 309 454 END SUBROUTINE trc_nam_trc 310 455 … … 357 502 WRITE(numout,*) ' frequency of outputs for biological trends nn_writebio = ', nn_writebio 358 503 WRITE(numout,*) ' ' 359 ENDIF 360 361 IF( ln_diatrc .AND. .NOT. lk_iomput ) THEN 504 CALL flush(numout) 505 ENDIF 506 !! 507 !! JPALM -- 17-07-2015 -- 508 !! MEDUSA is not yet up-to-date with the iom server. 509 !! we use it for the main tracer, but not fully with diagnostics. 510 !! will have to adapt it properly when visiting Christian Ethee 511 !! for now, we change 512 !! IF( ln_diatrc .AND. .NOT. lk_iomput ) THEN 513 !! to : 514 !! 515 IF( ( ln_diatrc .AND. .NOT. lk_iomput ) .OR. ( ln_diatrc .AND. lk_medusa ) ) THEN 362 516 ALLOCATE( trc2d(jpi,jpj,jpdia2d), trc3d(jpi,jpj,jpk,jpdia3d), & 363 517 & ctrc2d(jpdia2d), ctrc2l(jpdia2d), ctrc2u(jpdia2d) , & … … 368 522 trc3d(:,:,:,:) = 0._wp ; ctrc3d(:) = ' ' ; ctrc3l(:) = ' ' ; ctrc3u(:) = ' ' 369 523 ! 524 !! ELSE IF ( lk_iomput .AND. lk_medusa .AND. .NOT. ln_diatrc) THEN 525 !! CALL trc_nam_iom_medusa 370 526 ENDIF 371 527 … … 397 553 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 398 554 !!====================================================================== 399 END MODULE 555 END MODULE trcnam
Note: See TracChangeset
for help on using the changeset viewer.