Changeset 6201 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM
- Timestamp:
- 2016-01-04T14:06:11+01:00 (9 years ago)
- Location:
- branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/CFC/trcsms_cfc.F90
r6164 r6201 200 200 WRITE(NUMOUT,*) ' ' 201 201 END DO 202 # if defined key_debug_medusa 203 CALL flush(numout) 204 # endif 202 205 ENDIF 203 206 ! -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/IDTRA/trcsms_idtra.F90
r6164 r6201 82 82 !!---------------------------------------------------------------------- 83 83 IF (kt == nittrc000) THEN 84 IF(lwp) WRITE(numout,*) ' - JPALM - verif:'84 IF(lwp) WRITE(numout,*) ' trcsms_idtra :' 85 85 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~' 86 86 IF(lwp) WRITE(numout,*) ' - idtra decay factor : ', FDEC 87 # if defined key_debug_medusa 88 CALL flush(numout) 89 # endif 87 90 ENDIF 88 91 … … 160 163 IF( nn_timing == 1 ) CALL timing_stop('trc_sms_idtra') 161 164 ! 165 # if defined key_debug_medusa 166 IF(lwp) WRITE(numout,*) ' IDTRA DONE trc_sms_idtra -- ' 167 CALL flush(numout) 168 # endif 169 ! 162 170 END SUBROUTINE trc_sms_idtra 163 171 -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcnam.F90
r6164 r6201 270 270 ELSE ; IF(lwp) WRITE(numout,*) ' CFC not used' 271 271 ENDIF 272 272 ! 273 # if defined key_debug_medusa 274 CALL flush(numout) 275 IF (lwp) write (numout,*) '------------------------------' 276 IF (lwp) write (numout,*) 'Jpalm - debug' 277 IF (lwp) write (numout,*) 'CALL trc_nam_cfc -- OK' 278 IF (lwp) write (numout,*) 'in trc_nam - CALL trc_nam OK' 279 IF (lwp) write (numout,*) ' ' 280 # endif 281 ! 273 282 IF( lk_c14b ) THEN ; CALL trc_nam_c14b ! C14 bomb tracers 274 283 ELSE ; IF(lwp) WRITE(numout,*) ' C14 not used' -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r6174 r6201 86 86 IF( kt == nitrst - 2*nn_dttrc .OR. nstock == nn_dttrc .OR. ( kt == nitend - nn_dttrc .AND. .NOT. lrst_trc ) ) THEN 87 87 IF ( ln_rstdate ) THEN 88 !! JPALM -- 22-12-2015 -- modif to get the good date on restart trc 89 !file name 90 !! -- the condition to open the rst file is not 91 !the same than for the dynamic rst. 92 !! -- here it - for an obscure reason - is open 93 !2 time-step before the restart writing process 88 !! JPALM -- 22-12-2015 -- modif to get the good date on restart trc file name 89 !! -- the condition to open the rst file is not the same than for the dynamic rst. 90 !! -- here it - for an obscure reason - is open 2 time-step before the restart writing process 94 91 !! instead of 1. 95 !! -- i am not sure if someone forgot +1 in the 96 !if loop condition as 97 !! it is writen in all comments nitrst - 98 !2*nn_dttrc + 1 and the condition is nitrst - 2*nn_dttrc 99 !! -- nevertheless we didn't wanted to broke 100 !something already working 92 !! -- i am not sure if someone forgot +1 in the if loop condition as 93 !! it is writen in all comments nitrst -2*nn_dttrc + 1 and the condition is 94 !! nitrst - 2*nn_dttrc 95 !! -- nevertheless we didn't wanted to broke something already working 101 96 !! and just adapted the part we added. 102 !! -- So instead of calling ju2ymds( fjulday + 103 !(rdttra(1)) .... 104 !! we call ju2ymds( fjulday + (2*rdttra(1)) 105 !..... 97 !! -- So instead of calling ju2ymds( fjulday + (rdttra(1)) 98 !! we call ju2ymds( fjulday + (2*rdttra(1)) 106 99 !!-------------------------------------------------------------------- 107 100 CALL ju2ymds( fjulday + (2*rdttra(1)) / rday, iyear, imonth, iday, zsec ) … … 219 212 if (lwp) write (numout,'(a,3f15.5)') 'Sediment Ca ', & 220 213 & fq0, fq1, fq2 214 !! 215 !! AXY (07/07/15): read in temporally averaged fields for DMS 216 !! calculations 217 !! 218 IF( iom_varid( numrtr, 'B_DMS_CHN', ldstop = .FALSE. ) > 0 ) THEN 219 !! YES; in which case read them 220 !! 221 IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS present - reading in ...' 222 CALL iom_get( numrtr, jpdom_autoglo, 'B_DMS_CHN', zb_dms_chn(:,:) ) 223 CALL iom_get( numrtr, jpdom_autoglo, 'N_DMS_CHN', zn_dms_chn(:,:) ) 224 CALL iom_get( numrtr, jpdom_autoglo, 'B_DMS_CHD', zb_dms_chd(:,:) ) 225 CALL iom_get( numrtr, jpdom_autoglo, 'N_DMS_CHD', zn_dms_chd(:,:) ) 226 CALL iom_get( numrtr, jpdom_autoglo, 'B_DMS_MLD', zb_dms_mld(:,:) ) 227 CALL iom_get( numrtr, jpdom_autoglo, 'N_DMS_MLD', zn_dms_mld(:,:) ) 228 CALL iom_get( numrtr, jpdom_autoglo, 'B_DMS_QSR', zb_dms_qsr(:,:) ) 229 CALL iom_get( numrtr, jpdom_autoglo, 'N_DMS_QSR', zn_dms_qsr(:,:) ) 230 CALL iom_get( numrtr, jpdom_autoglo, 'B_DMS_DIN', zb_dms_din(:,:) ) 231 CALL iom_get( numrtr, jpdom_autoglo, 'N_DMS_DIN', zn_dms_din(:,:) ) 232 ELSE 233 !! NO; in which case set them to zero 234 !! 235 IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS absent - setting to zero ...' 236 zb_dms_chn(:,:) = 0.0 !! CHN 237 zn_dms_chn(:,:) = 0.0 238 zb_dms_chd(:,:) = 0.0 !! CHD 239 zn_dms_chd(:,:) = 0.0 240 zb_dms_mld(:,:) = 0.0 !! MLD 241 zn_dms_mld(:,:) = 0.0 242 zb_dms_qsr(:,:) = 0.0 !! QSR 243 zn_dms_qsr(:,:) = 0.0 244 zb_dms_din(:,:) = 0.0 !! DIN 245 zn_dms_din(:,:) = 0.0 246 ENDIF 247 !! 248 !! calculate stats on these fields 249 IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS stats (min, max, sum) ...' 250 fq0 = MINVAL(zn_dms_chn(:,:)) 251 fq1 = MAXVAL(zn_dms_chn(:,:)) 252 fq2 = SUM(zn_dms_chn(:,:)) 253 if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHN ', fq0, fq1, fq2 254 fq0 = MINVAL(zn_dms_chd(:,:)) 255 fq1 = MAXVAL(zn_dms_chd(:,:)) 256 fq2 = SUM(zn_dms_chd(:,:)) 257 if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHD ', fq0, fq1, fq2 258 fq0 = MINVAL(zn_dms_mld(:,:)) 259 fq1 = MAXVAL(zn_dms_mld(:,:)) 260 fq2 = SUM(zn_dms_mld(:,:)) 261 if (lwp) write (numout,'(a,3f15.5)') 'DMS, MLD ', fq0, fq1, fq2 262 fq0 = MINVAL(zn_dms_qsr(:,:)) 263 fq1 = MAXVAL(zn_dms_qsr(:,:)) 264 fq2 = SUM(zn_dms_qsr(:,:)) 265 if (lwp) write (numout,'(a,3f15.5)') 'DMS, QSR ', fq0, fq1, fq2 266 fq0 = MINVAL(zn_dms_din(:,:)) 267 fq1 = MAXVAL(zn_dms_din(:,:)) 268 fq2 = SUM(zn_dms_din(:,:)) 269 if (lwp) write (numout,'(a,3f15.5)') 'DMS, DIN ', fq0, fq1, fq2 270 221 271 #endif 222 272 … … 300 350 if (lwp) write (numout,'(a,3f15.5)') 'Sediment Ca ', & 301 351 & fq0, fq1, fq2 352 !! 353 !! AXY (07/07/15): write out temporally averaged fields for DMS 354 !! calculations 355 !! 356 IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS - writing out ...' 357 CALL iom_rstput( kt, nitrst, numrtw, 'B_DMS_CHN', zb_dms_chn(:,:) ) 358 CALL iom_rstput( kt, nitrst, numrtw, 'N_DMS_CHN', zn_dms_chn(:,:) ) 359 CALL iom_rstput( kt, nitrst, numrtw, 'B_DMS_CHD', zb_dms_chd(:,:) ) 360 CALL iom_rstput( kt, nitrst, numrtw, 'N_DMS_CHD', zn_dms_chd(:,:) ) 361 CALL iom_rstput( kt, nitrst, numrtw, 'B_DMS_MLD', zb_dms_mld(:,:) ) 362 CALL iom_rstput( kt, nitrst, numrtw, 'N_DMS_MLD', zn_dms_mld(:,:) ) 363 CALL iom_rstput( kt, nitrst, numrtw, 'B_DMS_QSR', zb_dms_qsr(:,:) ) 364 CALL iom_rstput( kt, nitrst, numrtw, 'N_DMS_QSR', zn_dms_qsr(:,:) ) 365 CALL iom_rstput( kt, nitrst, numrtw, 'B_DMS_DIN', zb_dms_din(:,:) ) 366 CALL iom_rstput( kt, nitrst, numrtw, 'N_DMS_DIN', zn_dms_din(:,:) ) 367 !! 368 !! calculate stats on these fields 369 IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS stats (min, max, sum) ...' 370 fq0 = MINVAL(zn_dms_chn(:,:)) 371 fq1 = MAXVAL(zn_dms_chn(:,:)) 372 fq2 = SUM(zn_dms_chn(:,:)) 373 if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHN ', fq0, fq1, fq2 374 fq0 = MINVAL(zn_dms_chd(:,:)) 375 fq1 = MAXVAL(zn_dms_chd(:,:)) 376 fq2 = SUM(zn_dms_chd(:,:)) 377 if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHD ', fq0, fq1, fq2 378 fq0 = MINVAL(zn_dms_mld(:,:)) 379 fq1 = MAXVAL(zn_dms_mld(:,:)) 380 fq2 = SUM(zn_dms_mld(:,:)) 381 if (lwp) write (numout,'(a,3f15.5)') 'DMS, MLD ', fq0, fq1, fq2 382 fq0 = MINVAL(zn_dms_qsr(:,:)) 383 fq1 = MAXVAL(zn_dms_qsr(:,:)) 384 fq2 = SUM(zn_dms_qsr(:,:)) 385 if (lwp) write (numout,'(a,3f15.5)') 'DMS, QSR ', fq0, fq1, fq2 386 fq0 = MINVAL(zn_dms_din(:,:)) 387 fq1 = MAXVAL(zn_dms_din(:,:)) 388 fq2 = SUM(zn_dms_din(:,:)) 389 if (lwp) write (numout,'(a,3f15.5)') 'DMS, DIN ', fq0, fq1, fq2 390 !! 302 391 #endif 303 392 -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcsms.F90
r6164 r6201 52 52 IF( lk_pisces ) CALL trc_sms_pisces ( kt ) ! main program of PISCES 53 53 IF( lk_medusa ) CALL trc_sms_medusa ( kt ) ! MEDUSA tracers 54 # if defined key_debug_medusa 55 IF(lwp) WRITE(numout,*) '--trcsms : MEDUSA OK -- next IDTRA -- ' 56 CALL flush(numout) 57 # endif 54 58 IF( lk_idtra ) CALL trc_sms_idtra ( kt ) ! radioactive decay of Id. tracer 59 # if defined key_debug_medusa 60 IF(lwp) WRITE(numout,*) '--trcsms : IDTRA OK -- next CFC -- ' 61 CALL flush(numout) 62 # endif 55 63 IF( lk_cfc ) CALL trc_sms_cfc ( kt ) ! surface fluxes of CFC 64 # if defined key_debug_medusa 65 IF(lwp) WRITE(numout,*) '--trcsms : CFC OK -- continue -- ' 66 CALL flush(numout) 67 # endif 56 68 IF( lk_c14b ) CALL trc_sms_c14b ( kt ) ! surface fluxes of C14 57 69 IF( lk_my_trc ) CALL trc_sms_my_trc ( kt ) ! MY_TRC tracers
Note: See TracChangeset
for help on using the changeset viewer.