- Timestamp:
- 2017-12-19T13:00:24+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_GO6_conserv_Check/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r9128 r9129 281 281 !! as proxy of org matter from the ocean 282 282 !! -- needed for the coupling with atm 283 !! 07-12-2017 -- To make things cleaner, we want to store an 284 !! unscaled Chl field in the restart and only 285 !! scale it when reading it in. 286 283 287 IF( iom_varid( numrtr, 'N_CHL_srf', ldstop = .FALSE. ) > 0 ) THEN 284 IF(lwp) WRITE(numout,*) 'Chl surf concentration - reading in ...'288 IF(lwp) WRITE(numout,*) 'Chl cpl concentration - reading in ... - scale by ', scl_chl 285 289 CALL iom_get( numrtr, jpdom_autoglo, 'N_CHL_srf', zn_chl_srf(:,:) ) 286 290 ELSE 287 IF(lwp) WRITE(numout,*) ' Chl surf concentration - setting to zero ...'288 zn_chl_srf(:,:) = (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * 1.E-6291 IF(lwp) WRITE(numout,*) 'set Chl coupled concentration - scaled by ', scl_chl 292 zn_chl_srf(:,:) = MAX( 0.0, (trn(:,:,1,jpchn) + trn(:,:,1,jpchd)) * 1.E-6 ) 289 293 ENDIF 290 294 IF (lk_oasis) THEN 291 chloro_out_cpl(:,:) = zn_chl_srf(:,:) !! Coupling variable295 chloro_out_cpl(:,:) = zn_chl_srf(:,:) * scl_chl !! Coupling variable 292 296 END IF 293 297 !! … … 301 305 call trc_rst_dia_stat(zn_dms_srf(:,:), 'DMS surf') 302 306 call trc_rst_dia_stat(zn_co2_flx(:,:), 'CO2 flux') 303 call trc_rst_dia_stat(zn_chl_srf(:,:), 'CHL surf') 307 IF (lk_oasis) THEN 308 call trc_rst_dia_stat(chloro_out_cpl(:,:), 'CHL cpl') 309 END IF 304 310 !! 305 311 !! JPALM 14-06-2016 -- add Carbonate chenistry variables through the restart … … 461 467 CALL iom_rstput( kt, nitrst, numrtw, 'B_CO2_flx', zb_co2_flx(:,:) ) 462 468 CALL iom_rstput( kt, nitrst, numrtw, 'N_CO2_flx', zn_co2_flx(:,:) ) 469 !! JPALM 07-12-2017 -- To make things cleaner, we want to store an 470 !! unscaled Chl field in the restart and only 471 !! scale it when reading it in. 463 472 CALL iom_rstput( kt, nitrst, numrtw, 'N_CHL_srf', zn_chl_srf(:,:) ) 464 473 !! … … 472 481 call trc_rst_dia_stat(zn_dms_srf(:,:), 'DMS surf') 473 482 call trc_rst_dia_stat(zn_co2_flx(:,:), 'CO2 flux') 474 call trc_rst_dia_stat(zn_chl_srf(:,:), ' CHL surf')483 call trc_rst_dia_stat(zn_chl_srf(:,:), 'unscaled CHL cpl') 475 484 !! 476 485 IF(lwp) WRITE(numout,*) ' MEDUSA averaged prop. for dust and iron dep.'
Note: See TracChangeset
for help on using the changeset viewer.